24.02.2013 Views

Optimality

Optimality

Optimality

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

On the false discovery proportion 35<br />

Lehmann and Romano [10] are summarized to motivate our choice of critical values.<br />

Control of the FDP is then considered in Section 3. The main result is presented in<br />

Theorem 3.4 and generalized in Theorem 3.5. In Section 4, we prove that a certain<br />

stepdown procedure controls the FDR under a dependence assumption.<br />

2. A class of stepdown procedures<br />

A formal description of our setup is as follows. Suppose data X is available from<br />

some model P∈ Ω. A general hypothesis H can be viewed as a subset ω of Ω. For<br />

testing Hi : P ∈ ωi, i = 1, . . . , s, let I(P) denote the set of true null hypotheses<br />

when P is the true probability distribution; that is, i∈I(P) if and only if P∈ ωi.<br />

We assume that p-values ˆp1, . . . , ˆps are available for testing H1, . . . , Hs. Specifically,<br />

we mean that ˆpi must satisfy<br />

(2.1) P{ˆpi≤ u}≤u for any u∈(0,1) and any P∈ ωi,<br />

Note that we do not require ˆpi to be uniformly distributed on (0,1) if Hi is true,<br />

in order to accomodate discrete situations.<br />

In general, a p-value ˆpi will satisfy (2.1) if it is obtained from a nested set of<br />

rejection regions. In other words, suppose Si(α) is a rejection region for testing Hi;<br />

that is,<br />

(2.2) P{X∈ Si(α)}≤α for all 0 < α < 1, P∈ ωi<br />

and<br />

(2.3) Si(α)⊂Si(α ′ ) whenever α < α ′ .<br />

Then, the p-value ˆpi defined by<br />

(2.4) ˆpi = ˆpi(X) = inf{α : X∈ Si(α)}.<br />

satisfies (2.1).<br />

In this article, we will consider the following class of stepdown procedures. Let<br />

(2.5) α1≤ α2≤···≤αs<br />

be constants, and let ˆp (1)≤···≤ ˆp (s) denote the ordered p-values. If ˆp (1) > α1,<br />

reject no null hypotheses. Otherwise,<br />

(2.6) ˆp (1)≤ α1, . . . , ˆp (r)≤ αr,<br />

and hypotheses H (1), . . . , H (r) are rejected, where the largest r satisfying (2.6) is<br />

used. That is, a stepdown procedure starts with the most significant p-value and<br />

continues rejecting hypotheses as long as their corresponding p-values are small.<br />

The Holm [6] procedure uses αi = α/(s−i+1) and controls the FWER at level<br />

α under no assumptions on the joint distribution of the p-values. Lehmann and<br />

Romano [10] generalized the Holm procedure to control the k-FWER. Specifically,<br />

consider the stepdown procedure described in (2.6), where we now take<br />

(2.7) αi =<br />

�<br />

kα<br />

s<br />

kα<br />

s+k−i<br />

i≤k<br />

i > k<br />

Of course, the αi depend on s and k, but we suppress this dependence in the<br />

notation.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!