[R-sig-ME] Extracting p-values from an lme object (any help would be appreciated)

Ben Bolker bolker at ufl.edu
Wed Nov 25 16:40:00 CET 2009


Ben Zuckerberg wrote:
> Greetings,
> 
> I know that p-value estimation in mixed models are somewhat 
> controversial at the moment, but I am simply curious how one extracts 
> p-values for the fixed factors from an lme object (using nlme)? The 
> summary() gives the desired p-values; all I want to do is access them.  
> The names() command shows the returned values (e.g., fitted, residuals) 
> for many components of the model, but I can't seem to access the 
> p-values.  Thank you in advance!
> 

library(nlme)
example(lme)
summary(fm1)$tTable[,"p-value"]



-- 
Ben Bolker
Associate professor, Biology Dep't, Univ. of Florida
bolker at ufl.edu / www.zoology.ufl.edu/bolker
GPG key: www.zoology.ufl.edu/bolker/benbolker-publickey.asc




More information about the R-sig-mixed-models mailing list