GLM 的 R 摘要未显示残差偏差

Question

当我在 R 中的 GLM 上运行摘要时，我希望看到偏差残差的摘要。我相信默认情况下应该显示它，但它不会出现在我身上。我知道我可以使用类似的东西

summary(resid(model,type="deviance")) 

要访问此信息，但我希望残差偏差与所有其他信息一起显示在摘要中。

我正在遵循教程，这就是示例：

#fit logistic regression model
model <- glm(am ~ disp + hp, data=mtcars, family=binomial)

#view model summary
summary(model)

我从摘要中得到的输出是：

Call:
glm(formula = am ~ disp + hp, family = binomial, data = mtcars)

Coefficients:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)  1.40342    1.36757   1.026   0.3048  
disp        -0.09518    0.04800  -1.983   0.0474 *
hp           0.12170    0.06777   1.796   0.0725 .
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 43.230  on 31  degrees of freedom
Residual deviance: 16.713  on 29  degrees of freedom
AIC: 22.713

Number of Fisher Scoring iterations: 8

这不包含偏差残差。我相信输出应该是：

Call:
glm(formula = am ~ disp + hp, family = binomial, data = mtcars)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.9665  -0.3090  -0.0017   0.3934   1.3682  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)  1.40342    1.36757   1.026   0.3048  
disp        -0.09518    0.04800  -1.983   0.0474 *
hp           0.12170    0.06777   1.796   0.0725 .
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 43.230  on 31  degrees of freedom
Residual deviance: 16.713  on 29  degrees of freedom
AIC: 22.713

Number of Fisher Scoring iterations: 8

Answer 1

这是由于 R 版本 4.3.0 中实现的更改所致。

相关新闻项目说

“summary.glm”类的 print() 方法默认不再显示偏差残差的汇总统计信息。如果需要，可以使用其可选参数 show.residuals 来显示它们。

因此，您可以通过显式

print

调用查看偏差残差的摘要：

print(summary(model),show.residuals=TRUE)