Python 中的 GLM 残差 statsmodel

Question

如何为 Python 中的所有 303 个观测值生成残差：

from statsmodels.stats.outliers_influence import OLSInfluence
OLSInfluence(resid)

或者

res.resid()

我正在尝试生成类似于我们在 R 中生成的残差：

res$resid

Answer 1

statsmodels 没有GLM的默认resid ，但它具有以下内容

resid_anscombe Anscombe 残差。

resid_anscombe_scaled 缩放的 Anscombe 残差。

resid_anscombe_unscaled 未缩放的 Anscombe 残差。

resid_deviance 偏差残差。

resid_pearson 皮尔逊残差。

resid_response 响应残差。

resid_working 工作残差。

https://www.statsmodels.org/stable/generated/statsmodels.genmod.generalized_linear_model.GLMResults.html

残差y - E(y|x)是响应残差resid_response

这些残差可用作fit方法返回的结果实例的属性。