我试图使用statsmodels进行负二项式回归。
我的代码是:
import statsmodels.api as sm
model = sm.NegativeBinomial
mod_nbin = model(Y, X)
res = mod_nbin.fit(disp=False, maxiter = 200)
res_nbin2 = mod_nbin.fit_regularized(start_params=res.params, alpha=10)
然而,当我看结果时:
print res_nbin2.summary()
我看到以下内容:
....
x45 -0.1876 0.024 -7.708 0.000 -0.235 -0.140
alpha 0.4742 0.011 41.510 0.000 0.452 0.497
所以不知何故,statsmodels正在拟合alpha值而不是使用我指定的值。关于如何解决这个问题的任何建议?
谢谢!