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