如何解决在泊松回归中绘制预测的置信区间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 举报,一经查实,本站将立刻删除。