2012-02-14 2 views
3

Если яКак извлечь значение р для наклона от МНК объекта в R

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

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

Но если я

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

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

ответ

4

Вы можете использовать соответствующую функцию экстрактор, но вам нужно позвонить 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