[R] singular convergence in glmmPQL

From: dave fournier <otter_at_otter-rsch.com>
Date: Tue 28 Feb 2006 - 08:38:40 EST

Hi,

If you are having trouble with convergence in glmmPQL you might want to try our glmmADMB package for R at

http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html

It is more stable according to some users.

I believe the corresponding commands for your model are

# paste this after your other commands

  df<-data.frame(y.long,arm.long,id.long)
# load glmmADMB package

 library("glmmADMB")
# invoke glmmADMB routine

 re.1 <- glmm.admb(y.long ~ arm.long,
random=~1,group='id.long',family='binomial',link='logit',data=df)

and that should do it.

   Cheers,

   Dave

-- 
David A. Fournier
P.O. Box 2040,
Sidney, B.C. V8l 3s3
Canada
http://otter-rsch.com

______________________________________________
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 Tue Feb 28 08:35:18 2006

This archive was generated by hypermail 2.1.8 : Tue 28 Feb 2006 - 09:08:49 EST