[R] AIC and anova, lme
Dieter Menne
dieter.menne at menne-biomed.de
Tue Feb 26 15:05:28 CET 2008
Patrick Giraudoux <patrick.giraudoux <at> univ-fcomte.fr> writes:
>
> Dear listers,
>
> Here we have a strange result we can hardly cope with. We want to
> compare a null mixed model with a mixed model with one independent
> variable.
>
> > lmmedt1<-lme(mediane~1, random=~1|site, na.action=na.omit, data=bdd2)
> > lmmedt9<-lme(mediane~log(0.0001+transat), random=~1|site,
> na.action=na.omit, data=bdd2)
...
> The usual conclusion would be that the two models are equivalent and to
> keep the null model for parsimony (!).
>
> However, an anova shows that the variable 'log(1e-04 + transat)' is
> significantly different from 0 in model 2 (lmmedt9)
>
> > anova(lmmedt9)
> numDF denDF F-value p-value
> (Intercept) 1 20 289.43109 <.0001
> log(1e-04 + transat) 1 20 31.18446 <.0001
>
Ask the author of pgirmess to add some checks for the model as anova and
stepAIC do:
Dieter
-----
library(MASS)
library(nlme)
fm1 <- lme(distance ~ age, data = Orthodont)
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)
>>In anova.lme(fm1, fm2) :
<< Fitted objects with different fixed effects. REML comparisons are not<<
meaningful.
stepAIC(fm2)
>>Error in extractAIC.lme(fit, scale, k = k, ...) :
>> AIC undefined for REML fit
More information about the R-help
mailing list