[R] p-values < 2.2e-16 not reported
gunter.berton at gene.com
Wed May 19 17:16:07 CEST 2010
If you believe P-values that small have any meaning at all, I have a bridge
to sell you...
Genentech Nonclinical Statistics
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of Will Eagle
Sent: Wednesday, May 19, 2010 1:53 AM
To: r-help at r-project.org
Subject: [R] p-values < 2.2e-16 not reported
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,
R-help at r-project.org mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
More information about the R-help