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

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 afternoo

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

Hi Jorge, and thanks for your answer, it looks promising. However, I have a question. First of all, I am a lowly biologist so please excuse any horrible errors of understanding I may make. So, the BH correction depends on, among other things, sorting the vector of P-values from the smallest to th

[R] Multiple testing corrections on very large vector

Hello all, I am calculating probabilities of association for terms in the GeneOntology database. I have ~4e7 probabilities and I would like to correct for multiple testing. At the moment I have something like the following: pv is a vector containing my 4e7 probabilities. To run the multipl