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

Question

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

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,6),loconst = 1)
pred1 <- predict(m2, newdata = NewData1, type = "response", interval = "confidence")
pred2 <- predict(m2, newdata = NewData2, type = "response", 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, lwd = 1.5, 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