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

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：广义线性模型 线性模型 模型预测 R语言 预测值

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

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

结果还是一样的

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

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

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")

1. newdata <- data.frame(sex.code = as.factor(1), class.code = as.factor(2))
   
2. predict.glm(model2, newdata=newdata, type = "response")

复制代码

输出结果
1: 52.6700000000001