singular convergence in glmmPQL

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

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

______________________________________________
[hidden email] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html