[R] obtaining exact p-values in mixed effects model

Rudi Alberts r.alberts at cs.rug.nl
Wed Sep 1 20:15:26 CEST 2004


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:


kind regards, R. Alberts

More information about the R-help mailing list