[R] obtaining exact p-values in mixed effects model
Yves Magliulo
ym at climpact.com
Wed Sep 1 12:38:41 CEST 2004
hi,
> Is there a way we can obtain
> the exact p-values from lme without rounding?
use summary instead.
> used commands:
>
> library(nlme)
> g<-lme(value~factor(fac1)+factor(fac2)+factor(fac1):factor(fac2),data=mydataframe,random=~1|factor(fac3))
> ag<-anova(g)
summary(g)$tTable[,5] will provide the exact p-value
you can always get the exact value of summary or predict etc... remember
it's not only printing, it's also a list
see ?summary.lme for more info
hope it's help
> kind regards, R. Alberts
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
------
Yves Magliulo <ym at climpact.com>
R&D Engineer, CLIMPACT
Tel. : +33 (0) 1 44 27 34 31
Fax. : +33 (0) 1 44 27 49 96
Universite Pierre et Marie Curie
Boite 101 - Tour 45 - 5eme etage - Couloir 45/46
4 place Jussieu, 75252 Paris CEDEX 05, France
More information about the R-help
mailing list