多元glm Poisson回归中的曲线置信区间(IC)
多元glm Poisson回归中的曲线置信区间(IC)
提问于 2019-03-17 06:07:56
我喜欢在多重glm泊松回归中绘制置信区间(IC)，但没有成功。如果我做了：
#Artificial data set
Consumption <- c(501, 502, 503, 504, 26, 27, 55, 56, 68, 69, 72, 93)
Gender <- gl(n = 2, k = 6, length = 2*6, labels = c("Male", "Female"), ordered = FALSE)
Income <- c(5010, 5020, 5030, 5040, 260, 270, 550, 560, 680, 690, 720, 930)
df3 <- data.frame(Consumption, Gender, Income)
df3创建glm泊松多元回归模型：
fm1 <- glm(Consumption~Gender+Income, data=df3, family=poisson)
summary(fm1)请参阅变量的重要性：
# ANOVA
anova(fm1,test="Chi")预测值和计算IC：
df3 = cbind(df3, pred = predict(fm1, type = "response"))#Estimate values
df3 = cbind(df3, se = predict(fm1, type="link",se.fit = TRUE)) ## Confidence interval of estimated
df3 = cbind(df3, ucl=exp(df3$pred + 1.96*df3$se.fit))
df3 = cbind(df3, lcl=exp(df3$pred - 1.96*df3$se.fit))现在，如果我尝试绘制这个：
#Plot
ggplot(data=df3, mapping=aes(x=Income, y=Consumption, color=Gender)) +
geom_point() +
geom_line(mapping=aes(y=pred)) +
geom_smooth(data=df3, aes(ymin = lcl, ymax = ucl), stat="identity") 不起作用：
Error in if ((w[1] * sm + w[2] * cm + w[3] * dm + w[4]) < best$score) break :
missing value where TRUE/FALSE needed有谁能帮我吗？
提前感谢！
回答 1
Stack Overflow用户
回答已采纳
发布于 2019-03-17 06:38:04
我建议使用以下代码计算预测值的置信区间：
pred <- predict(fm1, type="link", se.fit = TRUE)
df3 = cbind(df3, pred = pred$fit)
df3 = cbind(df3, se = pred$se.fit)
df3 = cbind(df3, ucl=exp(df3$pred + 1.96*df3$se))
df3 = cbind(df3, lcl=exp(df3$pred - 1.96*df3$se))或
pred <- predict(fm1, type="response", se.fit = TRUE)
df3 = cbind(df3, pred = pred$fit)
df3 = cbind(df3, se = pred$se.fit)
df3 = cbind(df3, ucl=df3$pred + 1.96*df3$se)
df3 = cbind(df3, lcl=df3$pred - 1.96*df3$se)
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/55201976
复制相关文章
相似问题
腾讯云开发者
