2014/03/31 by Ismaël Castillo, Johannes Schmidt-Hieber, Aad van der Vaart
Mathematics · #math.ST #stat.ME #stat.TH
paper · pdf · doi:10.1214/15-aos1334
published as Annals of Statistics 2015, Vol. 43, No. 5, 1986-2018 · Published at http://dx.doi.org/10.1214/15-AOS1334 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org)
arxiv created 2015/10/14 · arxiv updated 2015/10/15
We study full Bayesian procedures for high-dimensional linear regression under sparsity constraints. The prior is a mixture of point masses at zero and continuous distributions. Under compatibility conditions on the design matrix, the posterior distribution is shown to contract at the optimal rate for recovery of the unknown sparse vector, and to give optimal prediction of the response vector. It is also shown to select the correct sparse model, or at least the coefficients that are significantly different from zero. The asymptotic shape of the posterior distribution is characterized and employed to the construction and study of credible sets for uncertainty quantification.