Post‐selection point and interval estimation of signal sizes in Gaussian samples. Issue 2 (13th April 2017)

• Record Type: Journal Article
• Authors: Reid, Stephen; Taylor, Jonathan; Tibshirani, Robert
• Is Part Of: Canadian journal of statistics. Volume 45:Issue 2(2017)
• Page Start: 128
• Page End: 148
• Publication Date: 2017-04-13
• Subjects: Gaussian sequence model -- mean estimation -- minimax estimation -- selective inference -- shrinkage estimation -- MSC 2010: Primary 62Fxx -- secondary 62F07 -- 62F10 -- 62F12 -- 62F25; Mathematical statistics -- Periodicals; 519.5
• DOI: 10.1002/cjs.11320
• Languages: English
• ISSNs: 0319-5724

Abstract

Abstract: We tackle the problem of the estimation of a vector of underlying means (signal sizes) from a single vector‐valued observation y . Often one is interested in estimating only a subvector of signals corresponding to a set of selected, "interesting" sample elements. These "interesting" sample elements tend to have the largest absolute size, gleaned by applying some selection procedure like that of Benjamini & Hochberg (2015). Previous work on this estimation task proposes the reduction in size of the largest (absolute) sample elements either via shrinkage (like James–Stein) or by subtracting biases estimated using empirical Bayes methodology. We take a novel approach and adapt recent developments by Lee et al. (2016) in post‐selection inference. Adapting and applying their distributional results to our problem post‐selection point and interval estimators for underlying signal sizes are proposed. Simulations suggest that our estimator seems to perform quite well against competitors. Furthermore we prove an upper bound to the so‐called "worst case risk" of our estimator—when combined with the Benjamini–Hochberg selection procedure—and show that it is within a constant multiple of the minimax risk over a rich set of parameter spaces meant to evoke sparsity. The Canadian Journal of Statistics 45: 128–148; 2017 © 2017 Statistical Society of Canada