Multilevel Modeling using glmmPQL

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090311/b543a014/attachment-0002.pl>
  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