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: http://www.math.tau.ac.il/~ybenja/MyPapers/benjamini_yekutieli_ANNSTAT2001.pdf ----------------Contact Details:------------------------------------------------------- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English) ---------------------------------------------------------------------------------------------- On Sun, Aug 8, 2010 at 11:06 PM, qing chen <qingche...@gmail.com> wrote: > 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 with > that by using the independent case method. > Then how do the result of p.adjust( , fdr) come? > And how can we interpret them? > Thank you! > > > > Your truly, Qing Chen > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.