
   Hello, list members:

   I am fitting a mixed effects model with a split-plot analysis. I have read that in order to be able to use anova to compare models with a different fixed-effects structure, it's not correct to use the default REML method, so it must be changed to ML. 

   - When using "lme", that can be done just writing: method="ML", whereas when using "lmer" you achieve that by writing "REML=false".
   - My questions are: 
      1. Is it necessary to specify the ML method when fitting the maximum model, or is it better when simplifying the first term (with the "update" function), so that the maximum model is analized by REML and the subsequent simplified models by ML?
      2. Once you change the method to "ML" either in the maximum model or in the first simplified model, is it necessary to specify "method="ML"" or or "REML=false" (depending on the function) on each of the subsequent model simplifications or is it enough with writing it once?

   Thank you very much!

   Iker   



      
	[[alternative HTML version deleted]]


