[R] singular convergence in glmmPQL

[dave fournier]

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