obtaining exact p-values in mixed effects model
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