Show simple item record

dc.contributor.authorMadrigal-Cianci, Juan P.
dc.contributor.authorNobile, Fabio
dc.contributor.authorTempone, Raul
dc.identifier.citationMadrigal-Cianci, J. P., Nobile, F., & Tempone, R. (2023). Analysis of a Class of Multilevel Markov Chain Monte Carlo Algorithms Based on Independent Metropolis–Hastings. SIAM/ASA Journal on Uncertainty Quantification, 11(1), 91–138.
dc.description.abstractIn this work, we present, analyze, and implement a class of multilevel Markov chain Monte Carlo (ML-MCMC) algorithms based on independent Metropolis–Hastings proposals for Bayesian inverse problems. In this context, the likelihood function involves solving a complex differential model, which is then approximated on a sequence of increasingly accurate discretizations. The key point of this algorithm is to construct highly coupled Markov chains together with the standard multilevel Monte Carlo argument to obtain a better cost-tolerance complexity than a single-level MCMC algorithm. Our method extends the ideas of Dodwell et al., [SIAM/ASA J. Uncertain. Quantif., 3 (2015), pp. 1075–1108] to a wider range of proposal distributions. We present a thorough convergence analysis of the ML-MCMC method proposed, and show, in particular, that (i) under some mild conditions on the (independent) proposals and the family of posteriors, there exists a unique invariant probability measure for the coupled chains generated by our method, and (ii) that such coupled chains are uniformly ergodic. We also generalize the cost-tolerance theorem of Dodwell et al. to our wider class of ML-MCMC algorithms. Finally, we propose a self-tuning continuation-type ML-MCMC algorithm. The presented method is tested on an array of academic examples, where some of our theoretical results are numerically verified. These numerical experiments evidence how our extended ML-MCMC method is robust when targeting some pathological posteriors, for which some of the previously proposed ML-MCMC algorithms fail.
dc.description.sponsorshipThis project was partially funded by the King Abdullah University of Science and Technology through grant OSR-2015-CRG4-2584-01:   ``Advanced Multi-Level sampling techniques for Bayesian Inverse Problems with applications to subsurface." The first and second authors  also acknowledge support from the Swiss Data Science Center(SDSC) through grant P18-09. The third author acknowledges support from  the  Alexander von Humboldt foundation. The authors would like to extend their gratitude to the anonymous reviewers whose comments significantly improved the quality of this manuscript. The authors would also like to thank Dr. Sebastian Krumscheid, Dr. Panagiotis Tsilifis, and Dr. Anamika Pandey for fruitful discussions during the very early stages of this work.
dc.publisherSociety for Industrial & Applied Mathematics (SIAM)
dc.rightsPublished version archived with thanks to SIAM/ASA Journal on Uncertainty Quantification.
dc.titleAnalysis of a Class of Multilevel Markov Chain Monte Carlo Algorithms Based on Independent Metropolis–Hastings
dc.contributor.departmentComputer, Electrical and Mathematical Sciences and Engineering, KAUST, Thuwal, Saudi Arabia, and Alexander von Humboldt professor in Mathematics of Uncertainty Quantification, RWTH Aachen University, Aachen 52062, Germany.
dc.contributor.departmentApplied Mathematics and Computational Science Program
dc.contributor.departmentComputer, Electrical and Mathematical Science and Engineering (CEMSE) Division
dc.identifier.journalSIAM/ASA Journal on Uncertainty Quantification
dc.eprint.versionPublisher's Version/PDF
dc.contributor.institutionProtocol Labs.
dc.contributor.institutionSB-MATH-CSQI, École Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland.
kaust.personTempone, Raul

Files in this item

Published Version

This item appears in the following Collection(s)

Show simple item record


*Selected version