Feature request: adding a log.p= option to stats::p.adjust?
Is there any interest in adding a log.p= option to p.adjust() so that it
can accept log-transformed p-values and return log-transformed adjusted
I have some functions that, on occasion, return very low p-values. To
avoid underflow in such cases, I allow my users to set log.p=TRUE, which
is passed onto pt(), pchisq(), etc. to return log-transformed p-values.
However, currently, these log-p-values need to be run through exp()
before being used in p.adjust(). This re-introduces the possibility of
underflow errors, sometimes leading to adjusted p-values of zero.