[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 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.


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:
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.