R help - Sep 2004 - 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

Hello,

Using a fixed effects linear model (with lm), I can get exact p-values
out of the AVOVA table, even if they are very small, eg. 1.0e-200.

Using lme (linear mixed effects) from the nlme library,
it appears that there is rounding of the p-values to zero, if
the p-value is less than about 1.0e-16. Is there a way we can obtain 
the exact p-values from lme without rounding?


used commands:

library(nlme)
g<-lme(value~factor(fac1)+factor(fac2)+factor(fac1):factor(fac2),data=mydataframe,random=~1|factor(fac3))
ag<-anova(g)


kind regards, R. Alberts