Concave-Convex PDMP-based sampling

Sutton, Matthew and Fearnhead, Paul (2023) Concave-Convex PDMP-based sampling. Journal of Computational and Graphical Statistics, 32 (4). ISSN 1061-8600

[thumbnail of Concave_Convex-2]
Text (Concave_Convex-2)
Concave_Convex_2.pdf - Accepted Version
Available under License Creative Commons Attribution.

Download (535kB)


Recently nonreversible samplers based on simulating piecewise deterministic Markov processes (PDMPs) have shown potential for efficient sampling in Bayesian inference problems. However, there remains a lack of guidance on how to best implement these algorithms. If implemented poorly, the computational costs of simulating event times can outweigh the statistical efficiency of the nonreversible dynamics. Drawing on the adaptive rejection literature, we propose the concave-convex adaptive thinning approach for simulating a piecewise deterministic Markov process, which we call CC-PDMP. This approach provides a general guide for constructing bounds that may be used to facilitate PDMP-based sampling. A key advantage of this method is its additive structure—adding concave-convex decompositions yields a concave-convex decomposition. This makes the construction of bounds modular, as given a concave-convex decomposition for a class of likelihoods and a family of priors, they can be combined to construct bounds for the posterior. We show that constructing our bounds is simple and leads to computationally efficient thinning. Our approach is well suited to local PDMP simulation where conditional independence of the target can be exploited for potentially huge computational gains. We provide an R package and compare with existing approaches to simulating events in the PDMP literature. Supplementary materials for this article are available online.

Item Type:
Journal Article
Journal or Publication Title:
Journal of Computational and Graphical Statistics
Uncontrolled Keywords:
Research Output Funding/yes_externally_funded
?? stat.meyes - externally fundeddiscrete mathematics and combinatoricsstatistics and probabilitystatistics, probability and uncertainty ??
ID Code:
Deposited By:
Deposited On:
15 May 2023 14:15
Last Modified:
27 Jun 2024 00:50