Cargando…

Generating MCMC proposals by randomly rotating the regular simplex

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....

Descripción completa

Detalles Bibliográficos
Autor principal: Holbrook, Andrew J.
Formato: Online Artículo Texto
Lenguaje:English
Publicado: 2023
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10553364/
https://www.ncbi.nlm.nih.gov/pubmed/37799825
http://dx.doi.org/10.1016/j.jmva.2022.105106
_version_ 1785116152083513344
author Holbrook, Andrew J.
author_facet Holbrook, Andrew J.
author_sort Holbrook, Andrew J.
collection PubMed
description 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.
format Online
Article
Text
id pubmed-10553364
institution National Center for Biotechnology Information
language English
publishDate 2023
record_format MEDLINE/PubMed
spelling pubmed-105533642023-10-05 Generating MCMC proposals by randomly rotating the regular simplex Holbrook, Andrew J. J Multivar Anal Article 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. 2023-03 2022-09-23 /pmc/articles/PMC10553364/ /pubmed/37799825 http://dx.doi.org/10.1016/j.jmva.2022.105106 Text en https://creativecommons.org/licenses/by-nc-nd/4.0/This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/ (https://creativecommons.org/licenses/by-nc-nd/4.0/) ).
spellingShingle Article
Holbrook, Andrew J.
Generating MCMC proposals by randomly rotating the regular simplex
title Generating MCMC proposals by randomly rotating the regular simplex
title_full Generating MCMC proposals by randomly rotating the regular simplex
title_fullStr Generating MCMC proposals by randomly rotating the regular simplex
title_full_unstemmed Generating MCMC proposals by randomly rotating the regular simplex
title_short Generating MCMC proposals by randomly rotating the regular simplex
title_sort generating mcmc proposals by randomly rotating the regular simplex
topic Article
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10553364/
https://www.ncbi.nlm.nih.gov/pubmed/37799825
http://dx.doi.org/10.1016/j.jmva.2022.105106
work_keys_str_mv AT holbrookandrewj generatingmcmcproposalsbyrandomlyrotatingtheregularsimplex