Abstract:We consider a novel Bayesian approach to estimation, uncertainty quantification, and variable selection for a high-dimensional linear regression model under sparsity. The number of predictors can be nearly exponentially large relative to the sample size. We put a conjugate normal prior initially disregarding sparsity, but for making an inference, instead of the original multivariate normal posterior, we use the posterior distribution induced by a map transforming the vector of regression coefficients to a sparse vector obtained by minimizing the sum of squares of deviations plus a suitably scaled $\ell_1$-penalty on the vector. We show that the resulting sparse projection-posterior distribution contracts around the true value of the parameter at the optimal rate adapted to the sparsity of the vector. We show that the true sparsity structure gets a large sparse projection-posterior probability. We further show that an appropriately recentred credible ball has the correct asymptotic frequentist coverage. Finally, we describe how the computational burden can be distributed to many machines, each dealing with only a small fraction of the whole dataset. We conduct a comprehensive simulation study under a variety of settings and found that the proposed method performs well for finite sample sizes. We also apply the method to several real datasets, including the ADNI data, and compare its performance with the state-of-the-art methods. We implemented the method in the \texttt{R} package called \texttt{sparseProj}, and all computations have been carried out using this package.

Bayesian nonparametric regression with varying residual density

Bayesian nonparametric mean residual life regression

Bayesian Nonparametric Covariance Regression

A generalized Bayesian approach for high-dimensional robust regression with serially correlated errors and predictors

Adaptive Bayesian Regression on Data with Low Intrinsic Dimensionality

Robust Bayesian Inference for Berkson and Classical Measurement Error Models

Bayesian inference for high-dimensional nonstationary Gaussian processes

A Fully Nonparametric Modelling Approach to Binary Regression

Functional regression approximate Bayesian computation for Gaussian process density estimation

Bayesian Manifold Regression

A Latent-Variable Bayesian Nonparametric Regression Model

Adaptive Bayesian density estimation using Pitman-Yor or normalized inverse-Gaussian process kernel mixtures

Bayesian Nonparametric Regression and Density Estimation Using Integrated Nested Laplace Approximations

Bayesian Nonparametric Regression Analysis of Data with Random Effects Covariates from Longitudinal Measurements

Latent Factor Models for Density Estimation

Nonparametric Bayesian inference for mean residual life functions in survival analysis

Bayesian High-dimensional Linear Regression with Sparse Projection-posterior

Robust Bayesian Model Averaging for Linear Regression Models With Heavy-Tailed Errors

Quasi-Bayesian Nonparametric Density Estimation via Autoregressive Predictive Updates

Marginally Specified Priors for Nonparametric Bayesian Estimation

Tail-adaptive Bayesian shrinkage