R语言广义线性模型预测值

meijinrong

2294

收藏 2020-05-08

已知模型是poisson线性模型，自变量为class 和sex，给定一对class和sex的值，如何用predict函数预测因变量的值？原来的数据里有10个值。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

bwchen

2020-5-9 04:08:57

你返回的值是in sample predictions
用下面的代码应该可以得到你想要的：
predict.glm(model2, newdata = newdata, type = 'response')

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

meijinrong

2020-5-9 13:00:56

bwchen 发表于 2020-5-9 04:08
你返回的值是in sample predictions
用下面的代码应该可以得到你想要的：
predict.glm(model2, newdata = ...

结果还是一样的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

bwchen

2020-5-9 16:51:28

meijinrong 发表于 2020-5-9 13:00
结果还是一样的

把你的数据和整个有关代码发上来

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

meijinrong

2020-5-9 19:07:23

bwchen 发表于 2020-5-9 16:51
把你的数据和整个有关代码发上来

class.code <- as.factor(class.code)
sex.code <- as.factor(sex.code)
model0<-glm(n.policies ~ 1, family= poisson(),data = policies_data)
model1<-glm(n.policies ~ class.code, family= poisson(),data = policies_data)
anova(model0, model1, test = "Chisq")
summary(model1)
confint(model1, level = 0.95)
model2<-glm(n.policies ~ class.code+sex.code, family= poisson(),data = policies_data)
summary(model2)
confint(model2, level = 0.95)
anova(model1, model2, test = "Chisq")
newdata <- data.frame(c(sex.code=1,class.code = 2))
predict.glm(model2, newdata=newdata, type = "response")

附件列表

微信图片_20200509190651.png

原图尺寸 25.42 KB