Re: [R] standard dev in glmmPQL

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Fri 04 Aug 2006 - 14:10:43 EST

          Let's look at the first example on the help page for "glmmQPL":

library(MASS)
library(nlme) # will be loaded automatically if omitted fitPQL <- glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,

                      family = binomial, data = bacteria)
fitPQL
VarCorr(fitPQL)
VarCorr(fitPQL)[2,2]

          You might also be interested in some other functions that might help you solve problems like this in the future:

	  'str(fitPQL)' provides a compact summary of 'fitPQL'.
	  'class(fitPQL)' tells us it is of class 'glmmPQL' and 'lme'
	  'methods(class="lme")' identifies all the functions with special 
methods written for 'lme' objects.
	  Hope this helps.
	  Spencer Graves

Maria Salomé Esteves Cabral wrote:
> Hi!
>
> Can anyone let me know how can I get the stdDev of the random intercept from the output of glmmPQL?
>
>
> Thanks
>
> Salomé
>
> ______________________________________________
> 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
> and provide commented, minimal, self-contained, reproducible code.



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 and provide commented, minimal, self-contained, reproducible code. Received on Fri Aug 04 15:21:35 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 04 Aug 2006 - 16:16:59 EST.

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