[R] How to request AIC information from "lm" object?
ronggui
ronggui.huang at gmail.com
Fri Jun 9 08:41:18 CEST 2006
You can use AIC to get what you want.
#example from lm help page
> ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
> trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
> group <- gl(2,10,20, labels=c("Ctl","Trt"))
> weight <- c(ctl, trt)
> anova(lm.D9 <- lm(weight ~ group))
Analysis of Variance Table
Response: weight
Df Sum Sq Mean Sq F value Pr(>F)
group 1 0.6882 0.6882 1.4191 0.249
Residuals 18 8.7293 0.4850
#get the AIC of lm model.
> AIC(lm.D9)
[1] 46.17648
2006/6/9, Michael <comtech.usa at gmail.com>:
> Can "lm" return AIC information?
>
> Thanks a lot!
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
--
»ÆÈÙ¹ó
Deparment of Sociology
Fudan University
More information about the R-help
mailing list