24.02.2013 Views

Optimality

Optimality

Optimality

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

46 J. P. Romano and A. M. Shaikh<br />

Using this notation, observe that<br />

t<br />

E(FDP|ˆr1, . . . , ˆr s−|I|) = E(<br />

t + f {t + f > 0}|ˆr1, . . . , ˆr s−|I|)<br />

t<br />

≤ E(<br />

t + j {t > 0}|ˆr1, . . . , ˆr s−|I|)<br />

≤ |I|<br />

|I| + j E({t > 0}|ˆr1, . . . , ˆr s−|I|)<br />

≤ |I|<br />

|I| + j P{ˆq (1)≤ α ∗ j+1|ˆr1, . . . , ˆr s−|I|)<br />

≤ |I|<br />

|I| + j<br />

|I|<br />

�<br />

P{ˆqi≤ α ∗ j+1|ˆr1, . . . , ˆr s−|I|}<br />

i=1<br />

(4.2) ≤ |I|<br />

|I| + j |I|α∗ j+1<br />

≤ |I|2<br />

|I| + j min{<br />

sα<br />

,1}<br />

(s−j) 2<br />

(4.3) ≤ |I|α |I|s<br />

(s−j) (|I| + j)(s−j) .<br />

The inequality (4.2) follows from the assumption (3.5) on the joint distribution<br />

of p-values. To complete the proof, note that|I| + j≤ s. It follows that |I|α<br />

(s−j) ≤ α<br />

and (|I| + j)(s−j)−|I|s = j(s−|I|)−j 2 = j(s−|I|−j)≥0. Combining these<br />

two inequalities, we have that the expression in (4.3) is bounded above by α. The<br />

desired bound for the FDR follows immediately.<br />

The following simple example illustrates the fact that the FDR is not controlled<br />

by the stepdown procedure with constants α∗ i absent the restriction (3.5) on the<br />

dependence structure of the p-values.<br />

Example 4.1. Suppose there are s = 3 hypotheses, two of which are true. In this<br />

case, α∗ 1 = α<br />

3 , α∗ 2 = 3α<br />

4 , and α∗ 3 = min{3α,1}. Define the joint distribution of the<br />

two true p-values q1 and q2 as follows: Denote by Ii the half open interval [ i−1 i<br />

3 , 3 )<br />

and let (q1, q2)∼U(Ii×Ij) with probability 1<br />

6 for all (i, j) such that i�= j, 1≤i≤3<br />

and 1≤j≤ 3. It is easy to see that (q (1), q (2))∼U(Ii× Ij) with probability 1<br />

3 for<br />

all (i, j) such that i < j, 1≤i≤3and 1≤j≤ 3. Now define the distribution<br />

of the false p-value r1 conditional on (q1, q2) by the following rule: If q (1)≤ α/3,<br />

then let r1 = 1; otherwise, let r1 = 0. For such a joint distribution of (q1, q2, r1), we<br />

1<br />

and is at least 2 whenever<br />

have that the FDP is identically one whenever q (1)≤ α<br />

3<br />

α<br />

3 < q (1)≤ 3α<br />

4 . Hence,<br />

For α < 4<br />

9 , we therefore have that<br />

FDR≥P{q (1)≤ α 1<br />

} +<br />

3 2 P{α<br />

3 < q (1)≤ 3α<br />

4 }.<br />

FDR≥ 2α<br />

3<br />

+ (3α<br />

4<br />

α 13α<br />

− ) =<br />

3 12<br />

> α.

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

Saved successfully!

Ooh no, something went wrong!