Can anyone explain to me why the REML option is not available when using lme4::glmer?  I believe the recommendation for fitting linear mixed effects models using lmer is to first determine the random effects structure using REML (with the full fixed effect structure), then to conduct LR tests (or model selection) for fixed effects structure using ML, and finally to make inferences from the final model selected fit with REML.  However, when fitting a GLMM in lem4 there is no longer any REML option.  I am confused about why this is and would really appreciate knowing some more of the theory behind it.

Eric Stolen




	[[alternative HTML version deleted]]

