df - data.frame( matrix( rnorm(1000), nc=10 ) )
colnames(df) - c(y, paste(x, 1:9, sep=))
ifit - glm( y ~ ., data=df ) # initial fit
a - stepAIC( ifit, keep=extractAIC )
a$keep
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 10.000 9. 8. 7. 6. 5. 4.
[2,] 319.356 317.3819 315.4327 314.3526 313.2192 312.3311 311.1450
[,8] [,9][,10]
[1,] 3. 2. 1.
[2,] 310.2517 309.1266 308.1171
On Tue, 2005-11-29 at 19:01 +0100, [EMAIL PROTECTED] wrote:
Hi all,
I'm fitting GLM's using the step or stepAIC procedures and I would
like to save the AIC of the intermediate models. I would appreciate
very much information about how todo this.
Best wishes
Germán López
__
R-help@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
__
R-help@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