English

Generating MCMC proposals by randomly rotating the regular simplex

Computation 2022-09-15 v3 Other Statistics

Abstract

We present the simplicial sampler, a class of parallel MCMC methods that generate and choose from multiple proposals at each iteration. The algorithm's multiproposal randomly rotates a simplex connected to the current Markov chain state in a way that inherently preserves symmetry between proposals. As a result, the simplicial sampler leads to a simplified acceptance step: it simply chooses from among the simplex nodes with probability proportional to their target density values. We also investigate a multivariate Gaussian-based symmetric multiproposal mechanism and prove that it also enjoys the same simplified acceptance step. This insight leads to significant theoretical and practical speedups. While both algorithms enjoy natural parallelizability, we show that conventional implementations are sufficient to confer efficiency gains across an array of dimensions and a number of target distributions.

Keywords

Cite

@article{arxiv.2110.06445,
  title  = {Generating MCMC proposals by randomly rotating the regular simplex},
  author = {Andrew J. Holbrook},
  journal= {arXiv preprint arXiv:2110.06445},
  year   = {2022}
}

Comments

To appear in Journal of Multivariate Analysis. Code here: https://github.com/andrewjholbrook/simplicialSampler

R2 v1 2026-06-24T06:50:50.886Z