Hi,
We are trying to calculate ML-parameter estimates of a mixed effects models
where the observations are weighted. The "weights" option in lmer() with
RELM=FALSE seems attractive. Does anyone know the mechanism it uses to
calculate weighted ML estimates? (Is there a paper?). Also, our model
includes a serial correlation structure(e.g. AR(1)) among the residuals. As
far as I know, no such "cor" option is available in lmer() unlike in nlme.
Is this correct and if so, any suggestions on how to proceed?
thanks,
Harlan Campbell
McGill University
[[alternative HTML version deleted]]