[R] AIC and anova, lme

Patrick Giraudoux patrick.giraudoux at univ-fcomte.fr
Tue Feb 26 14:38:53 CET 2008


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)

Using the Akaike Criterion and selMod of the package pgirmess gives the 
following output:

 > selMod(list(lmmedt1,lmmedt9))
                 model       LL K  N2K       AIC  deltAIC  w_i      AICc 
deltAICc w_ic
2 log(1e-04 + transat) 44.63758 4  7.5 -81.27516 0.000000 0.65 -79.67516 
0.000000 0.57
1                    1 43.02205 3 10.0 -80.04410 1.231069 0.35 -79.12102 
0.554146 0.43

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

Has anyone an opinion about what looks like a paradox here ?

Patrick





More information about the R-help mailing list