Re: [R] {Spam?} RE: {Spam?} RE: {Spam?} Re: mgcv, testing gamm vs lme, which degrees of freedom?

2010-06-22 Thread Joris Meys
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

Re: [R] {Spam?} RE: {Spam?} RE: {Spam?} Re: mgcv, testing gamm vs lme, which degrees of freedom?

2010-06-21 Thread Christos Argyropoulos
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

Re: [R] {Spam?} RE: {Spam?} RE: {Spam?} Re: mgcv, testing gamm vs lme, which degrees of freedom?

2010-06-21 Thread Carlo Fezzi
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