[R-sig-ME] ML vs. REML to find a parsimonious mixed model
Maarten.Jung at mailbox.tu-dresden.de
Sun Apr 15 13:00:08 CEST 2018
I want to use LRTs via anova() on fitted linear mixed models (merMod
objects) to find a parsimonious mixed model containing only variance
components supported by the data (e.g. Matuschek et al. 2017 , Bates et
al. 2015 ).
In this situation my focus is *only on the reduction of the random effects
part* of the models.
The aforementioned papers use ML instead of REML estimation within this
process. Douglas Bates seems to prefer ML model comparison due to the
skewed nature of the distribution of variance estimators  and the user
Wolfgang states that "the ML estimator usually has lower mean-squared error
(MSE) than the REML estimator" . However, literally every textbook I
know suggests using REML estimation when comparing mixed models that differ
only in their random effect parts.
What would you suggest in this particular situation? ML or REML?
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models