在建模泊松生存时,我可以简单地用这样的偏移关键字包装曝光预测器:
glm(deaths ~ offset(exposure), family=poisson, data=data)
摘要只会产生截距估计。
In,Statsmodels.api,在曝光预测器上应用偏移的等效代码是什么? 我尝试了以下内容:
smf.glm('dead~exposure',\
data=datgrp,\
offset=np.log(datgrp['exposure']), \
family=sm.families.Poisson(\
link=sm.families.links.log)).fit().summary()
但是,返回代码会产生拦截和曝光的系数度量,这对我正在解决的示例问题是不准确的。