贝叶斯稀疏线性回归的神经元化先验

Neuronized Priors for Bayesian Sparse Linear Regression

Journal of the American Statistical Association · 2021
被引 16
ABS 4

中文导读

提出神经元化先验统一并扩展了拉普拉斯、柯西、马蹄铁和尖峰-板等收缩先验,通过高斯权重变量和激活函数变换实现高效贝叶斯变量选择,无需潜在指示变量,并提供了R包NPrior。

Abstract

Although Bayesian variable selection methods have been intensively studied, their routine use in practice has not caught up with their non-Bayesian counterparts such as Lasso, likely due to difficulties in both computations and flexibilities of prior choices. To ease these challenges, we propose the neuronized priors to unify and extend some popular shrinkage priors, such as Laplace, Cauchy, horseshoe, and spike-and-slab priors. A neuronized prior can be written as the product of a Gaussian weight variable and a scale variable transformed from Gaussian via an activation function. Compared with classic spike-and-slab priors, the neuronized priors achieve the same explicit variable selection without employing any latent indicator variables, which results in both more efficient and flexible posterior sampling and more effective posterior modal estimation. Theoretically, we provide specific conditions on the neuronized formulation to achieve the optimal posterior contraction rate, and show that a broadly applicable MCMC algorithm achieves an exponentially fast convergence rate under the neuronized formulation. We also examine various simulated and real data examples and demonstrate that using the neuronization representation is computationally more or comparably efficient than its standard counterpart in all well-known cases. An R package NPrior is provided for using neuronized priors in Bayesian linear regression.

贝叶斯统计变量选择机器学习计算统计