[R] p-values < 2.2e-16 not reported

Will Eagle will.eagle at gmx.net
Wed May 19 10:53:20 CEST 2010


Dear all,

how can I get the exact p-value of a statistical test like cor.test() if 
the p-value is below the default machine epsilon value of 
.Machine$double.eps =  2.220446e-16?

At the moment smaller p-values are reported as "p-value < 2.2e-16".
.Machine$double.eps <- 1E-100 does not solve this issue, although this 
value should be used by the format.pval() function.

To know the exact p-values down to 1E-200 is very important since I have 
multiple tests which require a alpha error-threshold below 2.2E-16.

Thanks in advance,

Will



More information about the R-help mailing list