我注意到在线性模型对象上调用 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
【问题讨论】：
它们是一样的。请在我们打印系数后查看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
【讨论】：
x 系数的 p 值，但不是整个模型的 p 值？那么这是两个不同的统计数据？
这只是一个可视化的问题，
summary(mod) 表示 p 值小于 2.2e-16 (
例如，如果你运行
summary(mod)$coefficient 你得到 4.66e-112
最好的
【讨论】：