Re: [R] Multiple testing corrections on very large vector

From: terdon <>
Date: Tue, 08 Mar 2011 09:33:30 -0800 (PST)

Hi Jorge,

    first of all THANKS! I just ran you suggestion and got blown away:

 system.time(res <- p.adjust(pv, method = 'fdr'))    user system elapsed
 55.052 3.100 62.560

I had tried the same thing using mt.rawp2adjp as per my original post, sent it to a cluster here on friday afternoon and it was still not finished (AND using ~16GB of memory) three days later. wow.

OK, so, great but what is the difference between mt.rawp2adjp and p.adjust? Do you know? Is it a simple difference in how the same algorithm is implemented? Is there any reason one would be more trustworthy than the other?

And again, thanks, wow...

View this message in context:
Sent from the R help mailing list archive at

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Tue 08 Mar 2011 - 18:34:51 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 09 Mar 2011 - 16:50:20 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive