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

From: Martin Maechler <maechler_at_stat.math.ethz.ch>
Date: Tue 18 Jan 2005 - 07:45:40 EST

>>>>> "GS" == Gordon Smyth <smyth@wehi.edu.au>
>>>>> on Sun, 16 Jan 2005 19:44:26 +1100 writes:

    GS> The new committed version of p.adjust() contains some     GS> problems:
    >> p.adjust(c(0.05,0.5),method="hommel")     GS> [1] 0.05 0.50

    GS> No adjustment!

yes, but that's still better than what the current version of R 2.0.1 does, namely to give NA NA + two warnings ..

    GS> I can't see how the new treatment of NAs can be
    GS> justified. One needs to distinguish between NAs which
    GS> represent missing p-values and NAs which represent
    GS> unknown p-values. In virtually all applications giving
    GS> rise to NAs, the NAs represent missing p-values which
    GS> could not be computed because of missing data. In such
    GS> cases, the observed p-values should definitely be
    GS> adjusted as if the NAs weren't there, because NAs
    GS> represent p-values which genuinely don't exist.

hmm, "definitely" being a bit strong. One could argue that ooonoe should use multiple imputation of the underlying missing data, or .. other scenarios.

I'll reply to your other, later, more detailed message separately and take the liberty to drop the other points here...

Martin



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

This archive was generated by hypermail 2.1.8 : Fri 18 Mar 2005 - 09:02:36 EST