如何在具有R之类公式的Python中运行GLM伽玛回归

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

警告：DeprecationWarning：不建议使用链接类作为参数调用Family（..）。改用链接类的实例。

Answer 1

我正在处理相同的问题。我注意到，如果删除（）内的所有内容，结果将完全相同，而不会出现错误。模型摘要具有相同的值，并且AIC值相同。

带有警告的代码：

glm_poisson_log = sm.GLM(endog, exog, family = sm.families.Poisson(sm.families.links.log))
results = glm_poisson_log.fit()

没有警告的代码：

glm_poisson_log = sm.GLM(endog, exog, family = sm.families.Poisson **()**)
results = glm_poisson_log.fit()

尝试删除链接信息。

此站点也特别有用：https://www.statsmodels.org/stable/glm.html