python - 如何使用类似 R 的公式在 Python 中运行 GLM Gamma 回归

Question

我正在使用以下代码使用 statsmodels 在 Python 中运行 GLM 回归。我特别想实现一个日志链接功能。我能够使用 Statsmodels 编写类似 R 的公式。

以下代码成功运行但抛出弃用警告。任何人都可以建议如何摆脱这个警告。感谢您的帮助。

代码:

mod = smf.glm(formula='y ~  C(x1) + C(x2) + C(x3) + x4 + x5', data=data,family=sm.families.Gamma(link=sm.families.links.log))
reg = mod.fit()
print(reg.summary())

警告:弃用警告:不推荐使用链接类作为参数调用 Family(..)。请改用链接类的实例。

Answer 1

我认为以上任何一个答案都不正确。按照警告的指示，正确的方法是使用括号内链接函数的相应实例:

mod = smf.glm(formula='y ~  C(x1) + C(x2) + C(x3) + x4 + x5', data=data,family=sm.families.Gamma(link=sm.families.links.log()))
reg = mod.fit()
print(reg.summary())