Home > Research > Publications & Outputs > The Apogee to Apogee Path Sampler

Electronic data

  • The Apogee to Apogee Path Sampler

    Rights statement: 12m

    Accepted author manuscript, 1.07 MB, PDF document

    Available under license: CC BY-NC: Creative Commons Attribution-NonCommercial 4.0 International License

Links

Text available via DOI:

View graph of relations

The Apogee to Apogee Path Sampler

Research output: Contribution to Journal/MagazineJournal articlepeer-review

Published
Close
<mark>Journal publication date</mark>2/10/2023
<mark>Journal</mark>Journal of Computational and Graphical Statistics
Issue number4
Volume32
Number of pages11
Pages (from-to)1436-1446
Publication StatusPublished
Early online date24/04/23
<mark>Original language</mark>English

Abstract

Among Markov chain Monte Carlo algorithms, Hamiltonian Monte Carlo (HMC) is often the algorithm of choice for complex, high-dimensional target distributions; however, its efficiency is notoriously sensitive to the choice of the integration-time tuning parameter. When integrating both forward and backward in time using the same leapfrog integration step as HMC, the set of apogees, local maxima in the potential along a path, is the same whatever point (position and momentum) along the path is chosen to initialize the integration. We present the Apogee to Apogee Path Sampler (AAPS), which uses this invariance to create a simple yet generic methodology for constructing a path, proposing a point from it and accepting or rejecting that proposal so as to target the intended distribution. We demonstrate empirically that AAPS has a similar efficiency to HMC but is much more robust to the setting of its equivalent tuning parameter, the number of apogees that the path crosses. Supplementary materials for this article are available online.