Fwd: Re: [Rd] p.adjust(<NA>s), was "Re: [BioC] limma and p-values"

From: Gordon Smyth <smyth_at_wehi.edu.au>
Date: Tue 18 Jan 2005 - 11:32:10 EST


On Sun, 16 Jan 2005 19:44:26 +1100, I wrote
>I can only think of one situation in which the NAs might represent unknown
>but existing p-values. This would be when a large experiment has been
>conducted leading to many p-values. Instead of inputing all the p-values
>to the p.adjust() function, you decide to enter only the smallest p-values
>and represent the others using NAs. The trouble with this approach is that
>it can only be used with method="bonferroni".

Actually this could be done with "holm" as well as "bonferroni". Holm's method enforces monotonicity of the adjusted p-values by working up from the low end rather than working down from the high end, so just knowing the small ones is still ok.

Gordon

> All the other adjustment methods are step-up or step-down methods or
> involve closure like Hommel's method. For these methods, you simply have
> to know all the p-values before you can adjust any of them.



R-devel@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Tue Jan 18 10:48:07 2005

This archive was generated by hypermail 2.1.8 : Tue 18 Jan 2005 - 11:22:12 EST