Как извлечь значение p для наклона из объекта ols в R

Question

Если я сделаю

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

, тогда я могу извлечь p-значение для наклона, используя summary(m1)$coefficients[2, 4].

Но если я сделаю

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

Что мне нужно сделать, чтобы извлечь значение р для наклона?

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