在泊松回归中绘制预测的置信区间r

问题描述

所以我有一个表格形式的模型

m2<- glm(totalcha ~ envp + loconst,data=capitalcontrols,family=poisson(link="log"))

我想在更改某些变量时绘制其预测值和置信区间

NewData1 <- data.frame(envp = seq(1,6),loconst = 0)
NewData2 <- data.frame(envp = seq(1,loconst = 1)
pred1 <- predict(m2,newdata = NewData1,type = "response",interval = "confidence")
pred2 <- predict(m2,newdata = NewData2,interval = "confidence")

`

但是，尽管我能够绘制线条，但无法使用置信区间进行绘制：

plot(t_wait,pred1,lwd = 1.5,col = "steelblue",type="l",ylim = c(0,10),main = " ",xlab = "s",ylab = "y")
 lines(t_wait,pred2,lty=2,col = "red")
 lines(t_wait,pred1[,2],lty = 'dashed',col = 'blue')

，最后一行给我错误“尺寸不正确”（即使它通常适用于线性模型）。 关于如何进行的任何想法？如果有人知道如何使用ggplot进行设置，那也很好

编辑：数据具有以下格式

head(capitalcontrols)

govcode   country year totalcha envp loconst envp_loconst
104 Australia 1951        0    2       0            0
105 Australia 1954        0    2       0            0
106 Australia 1956        0    2       0            0
107 Australia 1958        1    2       0            0
108 Australia 1961        0    2       1            2
109 Australia 1963        0    2       1            2

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）