[R] AIC and anova, lme

2008-02-26 Thread Patrick Giraudoux
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,

Re: [R] AIC and anova, lme

2008-02-26 Thread ian white
Patrick, The likelihoods of two models fitted using REML cannot be compared unless the fixed effects are the same in the two models. On Tue, 2008-02-26 at 14:38 +0100, Patrick Giraudoux wrote: Dear listers, Here we have a strange result we can hardly cope with. We want to compare a null

Re: [R] AIC and anova, lme

2008-02-26 Thread Patrick Giraudoux
ian white a écrit : Patrick, The likelihoods of two models fitted using REML cannot be compared unless the fixed effects are the same in the two models. Many thanks for this reminder. Shame on me: it recalls me that this subject may have been already largely discussed on this list. Now, I

Re: [R] AIC and anova, lme

2008-02-26 Thread Dieter Menne
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,