[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