Hi Carlos,
there is no possible way you can compare both models using a classical
statistical framework, be it ML, REML or otherwise. The assumptions
are violated. Regarding the df, see my previous mail.
In your case, I'd resort to the AIC/BIC criteria, and if prediction is
the main focus, compar
Hi Carlo,
One "quick" way to verify that the fixed effects parts are the same, is to
examine the design matrices that "gamm" passes on to lme. They should be the
same if the interpretation is correct. If you want to follow the mixed effect
relation to its end, why not specify "REML" as the fi
Hi Christos,
thanks for your kind reply, I agree entirely with your interpreation.
In the first model comparison, however, "anova" does seem to work
according to our interpretation, i.e. the "df" are equal in the two model.
My intuition is that the "anova" command does a fixed-effect test rather
3 matches
Mail list logo