在线性模型对象上调用summary()
与调用broom::glance()
时,我注意到p值有所不同。我认为,调用summary()
时的浮点精度仅限于2.2e-16
,而扫视范围可以超过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
答案 0 :(得分:2)
它们是相同的。打印系数后,请查看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
答案 1 :(得分:2)
这只是可视化问题,
summary(mod)
说p值小于2.2e-16(<2.2e-16),不相等。
例如,如果您运行
summary(mod)$coefficient
得到4.66e-112
最佳