[R] Multiple testing corrections on very large vector

terdon selenop at gmail.com
Tue Mar 8 18:33:30 CET 2011


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: http://r.789695.n4.nabble.com/Multiple-testing-corrections-on-very-large-vector-tp3341398p3341856.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list