Consistent High-Dimensional Bayesian Variable Selection via Penalized Credible Regions

• Published in: Journal of the American Statistical Association 2012-12, Vol.107 (500), p.1610-1624
• Main Authors: Bondell, Howard D, Reich, Brian J
• Format: Article
• Language: English
• Subjects: Algorithms; Bayesian analysis; Bayesian method; Consistency; Covariance; Covariance matrices; Credible region; data analysis; Datasets; Eigenvalues; Feature selection; LASSO; Markov analysis; Markov chain; Markovian processes; Mathematical vectors; Modeling; Monte Carlo simulation; Multidimensional analysis; Probabilities; Regression analysis; Sample size; screening; Standard error; Statistical analysis; Statistical variance; Statistics; Stochastic search; Theory and Methods
• ISSN: 1537-274X; 0162-1459; 1537-274X
• DOI: 10.1080/01621459.2012.716344

For high-dimensional data, particularly when the number of predictors greatly exceeds the sample size, selection of relevant predictors for regression is a challenging problem. Methods such as sure screening, forward selection, or penalized regressions are commonly used. Bayesian variable selection methods place prior distributions on the parameters along with a prior over model space, or equivalently, a mixture prior on the parameters having mass at zero. Since exhaustive enumeration is not feasible, posterior model probabilities are often obtained via long Markov chain Monte Carlo (MCMC) runs. The chosen model can depend heavily on various choices for priors and also posterior thresholds. Alternatively, we propose a conjugate prior only on the full model parameters and use sparse solutions within posterior credible regions to perform selection. These posterior credible regions often have closed-form representations, and it is shown that these sparse solutions can be computed via existing algorithms. The approach is shown to outperform common methods in the high-dimensional setting, particularly under correlation. By searching for a sparse solution within a joint credible region, consistent model selection is established. Furthermore, it is shown that, under certain conditions, the use of marginal credible intervals can give consistent selection up to the case where the dimension grows exponentially in the sample size. The proposed approach successfully accomplishes variable selection in the high-dimensional setting, while avoiding pitfalls that plague typical Bayesian variable selection methods.