Needles and straw in haystacks: Empirical Bayes estimates of possibly sparse sequences

Top Cited Papers
Open Access
Abstract
An empirical Bayes approach to the estimation of possibly sparse sequences observed in Gaussian white noise is set out and investigated. The prior considered is a mixture of an atom of probability at zero and a heavy-tailed density \gamma, with the mixing weight chosen by marginal maximum likelihood, in the hope of adapting between sparse and dense sequences. If estimation is then carried out using the posterior median, this is a random thresholding procedure. Other thresholding rules employing the same threshold can also be used. Probability bounds on the threshold chosen by the marginal maximum likelihood approach lead to overall risk bounds over classes of signal sequences of length n, allowing for sparsity of various kinds and degrees. The signal classes considered are ``nearly black'' sequences where only a proportion \eta is allowed to be nonzero, and sequences with normalized \ell_p norm bounded by \eta, for \eta >0 and 01. Simulations show excellent performance.Comment: Published by the Institute of Mathematical Statistics (http://www.imstat.org) in the Annals of Statistics (http://www.imstat.org/aos/) at http://dx.doi.org/10.1214/00905360400000003
All Related Versions

This publication has 12 references indexed in Scilit: