[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