[R] p.adjust( , fdr)

2010-08-08 Thread qing chen
Hello, I am not sure about the p.adjust( , fdr). How do these adjusted p-values get? I have read papers of BH method. For independent case, we compare the ordered p-values with the alfa*i/m, where m is the number of tests. But I have checked that result based on the adjusted p-values is different

Re: [R] p.adjust( , fdr)

2010-08-08 Thread Tal Galili
Hi Qing. I believe what you are referring to is: p.adjust(..., BH) The fdr option uses the method of Benjamini, Hochberg, and Yekutieli. (Not the original Benjamini, Hochberg article) It might be the method described here: