Hi List
I'm running a series of mixed models using lme, and I wonder if there
is a way to sort them by AIC prior to testing using anova
(lme1,lme2,lme3,....lme7) other than by hand.
My current output looks like this.
anova
(lme.T97NULL.ml,lme.T97FULL.ml,lme.T97NOINT.ml,lme.T972way.ml,lme.T97fc.
ml, lme.T97ns.ml, lme.T97min.ml)
Model df AIC BIC logLik Test
L.Ratio p-value
lme.T97NULL.ml 1 3 624.3053 635.0671 -309.1527
lme.T97FULL.ml 2 14 508.7594 558.9809 -240.3797 1 vs 2
137.54588 <.0001
lme.T97NOINT.ml 3 7 564.8574 589.9682 -275.4287 2 vs 3
70.09798 <.0001
lme.T972way.ml 4 12 507.0991 550.1461 -241.5496 3 vs 4
67.75832 <.0001
lme.T97fc.ml 5 9 502.5878 534.8730 -242.2939 4 vs 5
1.48869 0.6849
lme.T97ns.ml 6 8 504.8971 533.5951 -244.4485 5 vs 6
4.30928 0.0379
lme.T97min.ml 7 6 569.8188 591.3423 -278.9094 6 vs 7
68.92177 <.0001
Thanks
Ken
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
