Statsmodels.api GLM中的偏移量论证

时间:2017-02-12 23:09:12

标签: python statsmodels

在建模泊松生存时,我可以简单地用这样的偏移关键字包装曝光预测器:

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

但是,返回代码会产生拦截和曝光的系数度量,这对我正在解决的示例问题是不准确的。

0 个答案:

没有答案