[R] nlme for exact binomial model

From: Nonyane, Bareng A.S. <bnonyane_at_jhsph.edu>
Date: Thu, 31 Mar 2011 16:13:19 -0400


Hi
I am conducting a meta-analysis on proportions and I'm interested in using the exact binomial likeliood for the within-study effects rather than the approximate approaches and I don't want to use the quasi-likelihood approach in R's glmmPQL. Does anyone know how I can do this using nlme? I can do it in SAS nlmixed but I would rather use R -this is my SAS code:

proc nlmixed data=work.temp ;
parms mup=2 vp= 08; /*initial values*/ rawprop=1/(1+exp(-truep)); /*unkown true p logit proportion*/ model numberofisolates~binomial(sumisolatestudy, rawprop) ; random truep~normal(mup,vp) subject=studyid; run;

Thanks

B. Aletta Nonyane, PhD
Assistant Scientist
Department of International Health
Johns Hopkins Bloomberg School of Public Health 615 N. Wolfe Street
Baltimore, MD, 21205

        [[alternative HTML version deleted]]



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 31 Mar 2011 - 20:33:10 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 31 Mar 2011 - 20:50:25 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive