UCL Discovery Stage
UCL home » Library Services » Electronic resources » UCL Discovery Stage

Manifold lifting: scaling Markov chain Monte Carlo to the vanishing noise regime

Au, Khai Xiang; Graham, Matthew M; Thiery, Alexandre H; (2023) Manifold lifting: scaling Markov chain Monte Carlo to the vanishing noise regime. Journal of the Royal Statistical Society Series B: Statistical Methodology , 85 (3) pp. 757-782. 10.1093/jrsssb/qkad023. Green open access

[thumbnail of paper_and_supplementary_material.pdf]
Preview
Text
paper_and_supplementary_material.pdf - Accepted Version

Download (12MB) | Preview

Abstract

Standard Markov chain Monte Carlo methods struggle to explore distributions that concentrate in the neighbourhood of low-dimensional submanifolds. This pathology naturally occurs in Bayesian inference settings when there is a high signal-to-noise ratio in the observational data but the model is inherently over-parametrised or nonidentifiable. In this paper, we propose a strategy that transforms the original sampling problem into the task of exploring a distribution supported on a manifold embedded in a higher-dimensional space; in contrast to the original posterior this lifted distribution remains diffuse in the limit of vanishing observation noise. We employ a constrained Hamiltonian Monte Carlo method, which exploits the geometry of this lifted distribution, to perform efficient approximate inference. We demonstrate in numerical experiments that, contrarily to competing approaches, the sampling efficiency of our proposed methodology does not degenerate as the target distribution to be explored concentrates near low-dimensional submanifolds. Python code reproducing the results is available at https://doi.org/10.5281/zenodo.6551654.

Type: Article
Title: Manifold lifting: scaling Markov chain Monte Carlo to the vanishing noise regime
Open access status: An open access version is available from UCL Discovery
DOI: 10.1093/jrsssb/qkad023
Publisher version: https://doi.org/10.1093/jrsssb/qkad023
Language: English
Additional information: This version is the author accepted manuscript. For information on re-use, please refer to the publisher’s terms and conditions.
Keywords: Hamiltonian Monte Carlo, Bayesian inverse problems, nonidentifiability
UCL classification: UCL
URI: https://discovery-pp.ucl.ac.uk/id/eprint/10169057
Downloads since deposit
250Downloads
Download activity - last month
Download activity - last 12 months
Downloads by country - last 12 months

Archive Staff Only

View Item View Item