summary() 与 glance() 之间的 p 值差异？

Question

我注意到在线性模型对象上调用 summary() 与调用 broom::glance() 时的 p 值不同。我认为调用 summary() 时的浮点精度仅限于 2.2e-16 而 glance 可以达到 1e-100 以上。我的怀疑是正确的，还是这些值本质上不同？

x <- c(1, 2, 3, 4, 5, 6, 7, 8, 9)
y <- c(10, 20, 30, 40, 50, 60, 70, 80, 90)

mod <- lm(y~x)

summary(mod) # p < 2.2e-16

broom::glance(mod) # p = 4.66e-112

Answer 1

他们是一样的。请在我们打印系数后查看 summary 中的 P 值。

summary(mod) 

s$coefficients
#                 Estimate   Std. Error      t value      Pr(>|t|)
# (Intercept) 9.473903e-15 3.161050e-15 2.997075e+00  2.002483e-02
# x           1.000000e+01 5.617334e-16 1.780204e+16 4.661081e-112

Answer 2

这只是一个可视化问题，

summary(mod)表示p值小于2.2e-16（<2.2e-16），不等于。

例如，如果您运行

summary(mod)$coefficient 你得到 4.66e-112

最佳

summary() 与 glance() 之间的 p 值差异？

Difference in p-values between summary() vs glance()?

r

lm

summary

broom