【问题标题】:GLM Residual in Python statsmodelPython statsmodel 中的 GLM 残差
【发布时间】:2021-06-04 04:49:46
【问题描述】:

如何在 Python 中为所有 303 个观察结果生成残差:

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

res.resid()

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

res$resid

【问题讨论】:

    标签: python r regression linear-regression


    【解决方案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 方法返回的结果实例的属性。

    【讨论】:

    • 谢谢约瑟夫。如果,Result=smf.glm('y~x', data=x)Result=Result.fit()。如何使用 resid_response?
    • 更清晰的是区分模型和结果实例:model=smf.glm('y~x', data=x)results=model.fit()resid=results.resid_response
    • 代替data=x,即使我放data=y或其他一些数据,我仍然得到相同的结果。那么,这里的数据是什么?文档不是很清楚。
    • Data 是一个 pandas DataFrame,其列名在公式中使用,例如您的示例中的“y”和“x”。 (它也可以是字典或其他通过名称访问的结构。)这种浮动有很多例子。 statsmodels.org/dev/examples/notebooks/generated/formulas.html)
    • 嗨 Josef,有没有办法获得 RLM 的响应类型残差? results.resid 不与来自 R 的 RLM 残差绑定
    猜你喜欢
    • 2012-09-27
    • 2019-05-02
    • 1970-01-01
    • 2020-10-22
    • 1970-01-01
    • 2015-05-23
    • 2013-11-15
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多