[R] Multilevel Modeling using glmmPQL

[Chuck Cleland]

On 3/11/2009 11:29 AM, Howard Alper wrote:
>   Hi,
>
>   I'm trying to perform a power simulation for a simple multilevel
> model, using the function glmmPQL in R version 2.8.1.  I want to extract
> the p-value for the fixed-effects portion of the regression, but I'm
> having trouble doing that.  I can extract the coefficients
> (summary(fit)$coeff), and the covariance matrix (summary(fit)$varFix),
> but I can't grab the p-value (or the t-statistic.)  Could someone
> explain how to do this?
> 
>   Please send responses to halper at health.nyc.gov.

library(MASS)

fm1 <- glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
               family = binomial, data = bacteria)

summary(fm1)$tTable[,5]

  Look at str(summary(fm1)) .

>   Thanks in advance,
> 
>   Howard
> 
> DISCLAIMER:\ Sample Disclaimer added in a VBScript.\ \ 	...{{dropped:6}}
> 
> ______________________________________________
> 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. 

-- 
Chuck Cleland, Ph.D.
NDRI, Inc. (www.ndri.org)
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212) 845-4495 (Tu, Th)
tel: (732) 512-0171 (M, W, F)
fax: (917) 438-0894