2012-02-14 40 views
3

如果我如何提取p值用于从OLS斜率对象中的R

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值?

回答

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