21.06.2014 Views

Subsampling estimates of the Lasso distribution.

Subsampling estimates of the Lasso distribution.

Subsampling estimates of the Lasso distribution.

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.

58 Numerical results<br />

Figure 6.1: <strong>Subsampling</strong> confidence intervals for single scenarios <strong>of</strong> <strong>the</strong> models A and A’<br />

(n= 250). Red triangles stand for <strong>the</strong> true parameters.<br />

Hence, for each j ∈ {1, . . . , p} we define <strong>the</strong> test φ (j)<br />

k<br />

, corresponding to <strong>the</strong> intervals I(j)<br />

k<br />

(k = 1, . . . , 4) defined in <strong>the</strong> previous section.<br />

Rates <strong>of</strong> coverage and <strong>of</strong> false rejection based on 500 replications <strong>of</strong> <strong>the</strong> scenario are given<br />

in tables 6.1, 6.2 and 6.3 for <strong>the</strong> two sided interval I 2 (with corresponding test φ 2 ) and a<br />

subsample size b = n 0.85 . Results for <strong>the</strong> o<strong>the</strong>r intervals are omitted but it was noted that<br />

<strong>the</strong> one sided confidence interval I 1 has rates <strong>of</strong> coverage/false rejection similar to I 2 and<br />

that I 3 and I 4 are extremely conservative.<br />

Fluctuations from <strong>the</strong> nominal level for <strong>the</strong> interval I 2 can be attributed to <strong>the</strong> stochastic<br />

approximation to <strong>the</strong> true subsample quantile and to <strong>the</strong> fact that <strong>the</strong> penalization parameter<br />

is chosen by cross-validation for each scenario. Thus <strong>the</strong> coverage/false positive<br />

rates can be considered correct.<br />

6.1.3 F.W.E.R<br />

Instead <strong>of</strong> testing for individual hypo<strong>the</strong>ses, one can also be interrested in testing for a<br />

whole family <strong>of</strong> null hypo<strong>the</strong>ses. The goal is <strong>the</strong>n to control <strong>the</strong> probability to make

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

Saved successfully!

Ooh no, something went wrong!