Fast sampling with Gaussian scale-mixture priors in high-dimensional regression. Academic Article uri icon


  • We propose an efficient way to sample from a class of structured multivariate Gaussian distributions. The proposed algorithm only requires matrix multiplications and linear system solutions. Its computational complexity grows linearly with the dimension, unlike existing algorithms that rely on Cholesky factorizations with cubic complexity. The algorithm is broadly applicable in settings where Gaussian scale mixture priors are used on high-dimensional parameters. Its effectiveness is illustrated through a high-dimensional regression problem with a horseshoe prior on the regression coefficients. Other potential applications are outlined.

published proceedings

  • Biometrika

altmetric score

  • 0.5

author list (cited authors)

  • Bhattacharya, A., Chakraborty, A., & Mallick, B. K.

citation count

  • 59

complete list of authors

  • Bhattacharya, Anirban||Chakraborty, Antik||Mallick, Bani K

publication date

  • December 2016