[R] generalized linear mixed model by ML

From: Abderrahim Oulhaj <abderrahim.oulhaj_at_pharmacology.oxford.ac.uk>
Date: Thu 15 Dec 2005 - 22:48:47 EST

Dear All,

I wonder if there is a way to fit a generalized linear mixed models (for repeated binomial data) via a direct Maximum Likelihood Approach. The "glmm" in the "repeated" package (Lindsey), the "glmmPQL" in the "MASS" package (Ripley) and "glmmGIBBS" (Myle and Calyton) are not using the full maximum likelihood as I understand. The "glmmML" of Brostrom uses the "full maximum likelihood" by approximating the integral via Gauss- Hermite quadrature. However, glmmML is only valid for the random intercept model and the binomial family must be represented only as binary data. Does the lmer do the work?

The reason for asking that is:

McCulloch (1997): "Unfortunately, PQL methods can perform poorly relative to Maximum likelihood" Agresti (2002) p. 524 "........but where possible we recommend using ML rather than PQL".

Thanks a lot,


        [[alternative HTML version deleted]]

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Dec 15 23:05:27 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:41:39 EST