Zigzag path connects two Monte Carlo samplers: Hamiltonian counterpart to a piecewise deterministic Markov process
Zigzag and other piecewise deterministic Markov process samplers have attracted significant interest for their non-reversibility and other appealing properties for Bayesian posterior computation. Hamiltonian Monte Carlo is another state-of-the-art sampler, exploiting fictitious momentum to guide Mar...
Saved in:
Published in | arXiv.org |
---|---|
Main Authors | , , |
Format | Paper |
Language | English |
Published |
Ithaca
Cornell University Library, arXiv.org
29.07.2024
|
Subjects | |
Online Access | Get full text |
Cover
Loading…
Abstract | Zigzag and other piecewise deterministic Markov process samplers have attracted significant interest for their non-reversibility and other appealing properties for Bayesian posterior computation. Hamiltonian Monte Carlo is another state-of-the-art sampler, exploiting fictitious momentum to guide Markov chains through complex target distributions. We establish an important connection between the zigzag sampler and a variant of Hamiltonian Monte Carlo based on Laplace-distributed momentum. The position and velocity component of the corresponding Hamiltonian dynamics travels along a zigzag path paralleling the Markovian zigzag process; however, the dynamics is non-Markovian in this position-velocity space as the momentum component encodes non-immediate pasts. This information is partially lost during a momentum refreshment step, in which we preserve its direction but re-sample magnitude. In the limit of increasingly frequent momentum refreshments, we prove that Hamiltonian zigzag converges strongly to its Markovian counterpart. This theoretical insight suggests that, when retaining full momentum information, Hamiltonian zigzag can better explore target distributions with highly correlated parameters by suppressing the diffusive behavior of Markovian zigzag. We corroborate this intuition by comparing performance of the two zigzag cousins on high-dimensional truncated multivariate Gaussians, including a 11,235-dimensional target arising from a Bayesian phylogenetic multivariate probit modeling of HIV virus data. |
---|---|
AbstractList | Zigzag and other piecewise deterministic Markov process samplers have attracted significant interest for their non-reversibility and other appealing properties for Bayesian posterior computation. Hamiltonian Monte Carlo is another state-of-the-art sampler, exploiting fictitious momentum to guide Markov chains through complex target distributions. We establish an important connection between the zigzag sampler and a variant of Hamiltonian Monte Carlo based on Laplace-distributed momentum. The position and velocity component of the corresponding Hamiltonian dynamics travels along a zigzag path paralleling the Markovian zigzag process; however, the dynamics is non-Markovian in this position-velocity space as the momentum component encodes non-immediate pasts. This information is partially lost during a momentum refreshment step, in which we preserve its direction but re-sample magnitude. In the limit of increasingly frequent momentum refreshments, we prove that Hamiltonian zigzag converges strongly to its Markovian counterpart. This theoretical insight suggests that, when retaining full momentum information, Hamiltonian zigzag can better explore target distributions with highly correlated parameters by suppressing the diffusive behavior of Markovian zigzag. We corroborate this intuition by comparing performance of the two zigzag cousins on high-dimensional truncated multivariate Gaussians, including a 11,235-dimensional target arising from a Bayesian phylogenetic multivariate probit modeling of HIV virus data. |
Author | Zhang, Zhenyu Nishimura, Akihiko Suchard, Marc A |
Author_xml | – sequence: 1 givenname: Akihiko surname: Nishimura fullname: Nishimura, Akihiko – sequence: 2 givenname: Zhenyu surname: Zhang fullname: Zhang, Zhenyu – sequence: 3 givenname: Marc surname: Suchard middlename: A fullname: Suchard, Marc A |
BookMark | eNqNyrFOAkEQgOEN0QRE3mESa5K7XU7Qlmho6KxsyGQZcPFuZt2Zk8SndwsfwOovvv_O3bAwTdzMh9AuNyvvp26hemmaxj-ufdeFmSvv6fyDZ8hoHxCFmaIp2FVgL2wEWyy9gOKQeyr6DDscUm_CCbnuY11KxmJgAgg5UaRrUoIjVRgSJ7UUYY_lU74hF4mkeu9uT9grLf46dw-vL2_b3bL610hqh4uMhSsdfNeuNk9-3Ybwv-sXkYlOPg |
ContentType | Paper |
Copyright | 2024. This work is published under http://arxiv.org/licenses/nonexclusive-distrib/1.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License. |
Copyright_xml | – notice: 2024. This work is published under http://arxiv.org/licenses/nonexclusive-distrib/1.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License. |
DBID | 8FE 8FG ABJCF ABUWG AFKRA AZQEC BENPR BGLVJ CCPQU DWQXO HCIFZ L6V M7S PIMPY PQEST PQQKQ PQUKI PRINS PTHSS |
DatabaseName | ProQuest SciTech Collection ProQuest Technology Collection Materials Science & Engineering Collection ProQuest Central (Alumni) ProQuest Central ProQuest Central Essentials ProQuest Central Technology Collection ProQuest One Community College ProQuest Central Korea SciTech Premium Collection ProQuest Engineering Collection Engineering Database Publicly Available Content Database ProQuest One Academic Eastern Edition (DO NOT USE) ProQuest One Academic ProQuest One Academic UKI Edition ProQuest Central China Engineering Collection |
DatabaseTitle | Publicly Available Content Database Engineering Database Technology Collection ProQuest Central Essentials ProQuest One Academic Eastern Edition ProQuest Central (Alumni Edition) SciTech Premium Collection ProQuest One Community College ProQuest Technology Collection ProQuest SciTech Collection ProQuest Central China ProQuest Central ProQuest Engineering Collection ProQuest One Academic UKI Edition ProQuest Central Korea Materials Science & Engineering Collection ProQuest One Academic Engineering Collection |
DatabaseTitleList | Publicly Available Content Database |
Database_xml | – sequence: 1 dbid: 8FG name: ProQuest Technology Collection url: https://search.proquest.com/technologycollection1 sourceTypes: Aggregation Database |
DeliveryMethod | fulltext_linktorsrc |
Discipline | Physics |
EISSN | 2331-8422 |
Genre | Working Paper/Pre-Print |
GroupedDBID | 8FE 8FG ABJCF ABUWG AFKRA ALMA_UNASSIGNED_HOLDINGS AZQEC BENPR BGLVJ CCPQU DWQXO FRJ HCIFZ L6V M7S M~E PIMPY PQEST PQQKQ PQUKI PRINS PTHSS |
ID | FETCH-proquest_journals_25148927133 |
IEDL.DBID | 8FG |
IngestDate | Thu Oct 10 22:32:06 EDT 2024 |
IsOpenAccess | true |
IsPeerReviewed | false |
IsScholarly | false |
Language | English |
LinkModel | DirectLink |
MergedId | FETCHMERGED-proquest_journals_25148927133 |
OpenAccessLink | https://www.proquest.com/docview/2514892713?pq-origsite=%requestingapplication% |
PQID | 2514892713 |
PQPubID | 2050157 |
ParticipantIDs | proquest_journals_2514892713 |
PublicationCentury | 2000 |
PublicationDate | 20240729 |
PublicationDateYYYYMMDD | 2024-07-29 |
PublicationDate_xml | – month: 07 year: 2024 text: 20240729 day: 29 |
PublicationDecade | 2020 |
PublicationPlace | Ithaca |
PublicationPlace_xml | – name: Ithaca |
PublicationTitle | arXiv.org |
PublicationYear | 2024 |
Publisher | Cornell University Library, arXiv.org |
Publisher_xml | – name: Cornell University Library, arXiv.org |
SSID | ssj0002672553 |
Score | 3.5617158 |
SecondaryResourceType | preprint |
Snippet | Zigzag and other piecewise deterministic Markov process samplers have attracted significant interest for their non-reversibility and other appealing properties... |
SourceID | proquest |
SourceType | Aggregation Database |
SubjectTerms | Bayesian analysis Markov analysis Markov chains Momentum Multivariate analysis Samplers |
Title | Zigzag path connects two Monte Carlo samplers: Hamiltonian counterpart to a piecewise deterministic Markov process |
URI | https://www.proquest.com/docview/2514892713 |
hasFullText | 1 |
inHoldings | 1 |
isFullTextHit | |
isPrint | |
link | http://utb.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwfV3dS8MwED90Q_DNT_yY40Bfi21a280XwdE6hI0hCsOXkaa3MZC1a6oDH_zbzcVOH4Q9hkBIcuE-frn7HcBVJ1IqVZKZj93ICaYknG4gpk7mGRlLV5DKLNvnMOy_BI_jm3ENuOk6rXKtE62iznLFGPm1scNBpytMTHVXLB3uGsW_q3ULjW1oesyEx5XiycMvxiLCyHjM_j81a21HsgfNkSyo3IctWhzAjk25VPoQytf57FPOkHsCo-J8E1VprFY5DpgxCnuyfMtRS6bvLfUt9hmKMJ6akSfaDg9UFkbuWOUosZiTotVcE2Z1foslYEYuxck_sPgpBziCyyR-7vWd9UYn9VPSk7-D-8fQWOQLOgEMMo98n2M_475NQzJeh09pRK6XSRUFnVNobVrpbPP0OewKY7sZwhTdFjSq8p0ujO2t0ra94DY07-Ph6MmMBl_xN0cnkaA |
link.rule.ids | 786,790,12792,21416,33406,33777,43633,43838 |
linkProvider | ProQuest |
linkToHtml | http://utb.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwfV3fS8MwED50Q_TNn_hj6oG-Fre0tp0vgmOj6jaGTBi-lDS9jcJYa1Md-Nebq50-CHsOhCQX7rv7cvkO4Nr3lIqUZOXjpmc5UxJW2xFTK24ZG8umIBWXap9DN3h1nia3k4pw01VZ5conlo46ThVz5DcGhx2_LUxOdZ-9W9w1il9XqxYam1BnyU2_BvWH7nD08suyCNczMbP9z9GW6NHbhfpIZpTvwQYt9mGrLLpU-gDyt2T2JWfIXYFRccWJKjQWyxQHrBmFHZnPU9SSBXxzfYcBkxEmVjMWxbLHA-WZsTwWKUrMElK0TDRhXFW4lBLMyJ9x0k_Mfj4EHMJVrzvuBNZqoWF1mXT4t3X7CGqLdEHHgE7cItvm7M8EcFOXTNxhU-RRsxVL5Tn-CTTWzXS6fvgStoPxoB_2H4fPZ7AjDJIzoSnaDagV-QedGyQuoovquL8Bss2TLA |
openUrl | ctx_ver=Z39.88-2004&ctx_enc=info%3Aofi%2Fenc%3AUTF-8&rfr_id=info%3Asid%2Fsummon.serialssolutions.com&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=Zigzag+path+connects+two+Monte+Carlo+samplers%3A+Hamiltonian+counterpart+to+a+piecewise+deterministic+Markov+process&rft.jtitle=arXiv.org&rft.au=Nishimura%2C+Akihiko&rft.au=Zhang%2C+Zhenyu&rft.au=Suchard%2C+Marc+A&rft.date=2024-07-29&rft.pub=Cornell+University+Library%2C+arXiv.org&rft.eissn=2331-8422 |