我想用置信区间绘制模型回归。我使用了ggplot,但是它仅显示我的两个类别之一的置信区间。可能是一类相对于另一类的标准误范围太小,这就是为什么ggplot只显示一类的置信区间的原因?
当我分别绘制两个类别时,我得到了置信区间,但是当我将它们合并到ggplot中的一个图中时,则只有其中一个显示有置信区间。
我想绘制出食物消耗率随温度(temp)升高如何变化的线性回归,并比较这两种类型的食物。
这是我的模特
m1 <- glmmPQL(totbeakfuls ~ scale(temp)*day.type + offset(log(duration/60)),
random = ~1|birdid,
family = quasipoisson(link="log"),
data = focaldata)
然后,我考虑到它是日类型和温度之间的交互模型,因此使用ggeffects来生成模型的预测。
dat<- ggeffects::ggpredict(m1, terms = c("temp", "day.type"), type = "re")
然后我使用ggplot进行这些预测
ggplot(data = focaldata, aes(x = temp, y = consumption)) +
geom_point(aes(color = day.type)) +
geom_smooth(data= dat, method = "lm",
aes(x = x, y= predicted, fill= group, color=day.type),
se = TRUE)
即使我指定了“ se = TRUE”,它也只会绘制一天类型类别之一的置信区间。