我知道这可能是一个基本问题...但是我似乎找不到答案。
我正在为Poisson系列安装GLM,然后尝试查看预测,但是似乎确实考虑了偏移量:
model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003),
offset=(log(population)), data=data, subset=28:36, family=poisson())
predict (model_glm, type="response")
我得到的情况不是费率...
我也尝试过
model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+
offset(log(population)), data=data, subset=28:36, family=poisson())
结果相同。但是,当我使用mgcv从GAM进行预测时,这些预测会考虑偏移量(我得到比率)。
我想念什么吗?
1
请不要在此处和r-help列表上交叉发布...如果您打算在stackoverflow / stackexchange论坛上发布,我认为这样会更好(这是一个技术性R问题,而不是统计问题) ...)
—
Ben Bolker