如何从R中的ols对象中提取斜率的p值

Question

如果我这样做

data(mtcars)
m1 <- lm(mpg ~ cyl, data= mtcars, x= TRUE, y= TRUE)

然后我可以使用summary(m1)$coefficients[2, 4] 提取斜率的 p 值。

如果我这样做了

library(rms)
data(mtcars)
m2 <- ols(mpg ~ cyl, data= mtcars, x= TRUE, y= TRUE)

我需要做什么来提取斜率的 p 值？

Answer 1

可以使用对应的提取器功能，但需要调用summary.lm：

> coef(summary.lm(m2))
          Estimate Std. Error   t value     Pr(>|t|)
Intercept 37.88458  2.0738436 18.267808 8.369155e-18
cyl       -2.87579  0.3224089 -8.919699 6.112687e-10