Multilevel sequential Monte Carlo samplers
Abstract
Here, we study the approximation of expectations w.r.t. probability distributions associated to the solution of partial differential equations (PDEs); this scenario appears routinely in Bayesian inverse problems. In practice, one often has to solve the associated PDE numerically, using, for instance finite element methods and leading to a discretisation bias, with the step-size level hL. In addition, the expectation cannot be computed analytically and one often resorts to Monte Carlo methods. In the context of this problem, it is known that the introduction of the multilevel Monte Carlo (MLMC) method can reduce the amount of computational effort to estimate expectations, for a given level of error. This is achieved via a telescoping identity associated to a Monte Carlo approximation of a sequence of probability distributions with discretisation levels $${\infty}$$ >h0>h1 ...>hL. In many practical problems of interest, one cannot achieve an i.i.d. sampling of the associated sequence of probability distributions. A sequential Monte Carlo (SMC) version of the MLMC method is introduced to deal with this problem. In conclusion, it is shown that under appropriate assumptions, the attractive property of a reduction of the amount of computational effort to estimate expectations, for a given level of error, can be maintainedmore »
- Authors:
-
- Univ. College London, London (United Kingdom). Dept. of Statistical Science
- National Univ. of Singapore (Singapore). Dept. of Statistics & Applied Probability
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Computer Science and Mathematics Division
- King Abdullah Univ. of Science and Technology, Thuwal (Saudi Arabia)
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Org.:
- USDOE Laboratory Directed Research and Development (LDRD) Program
- OSTI Identifier:
- 1302922
- Grant/Contract Number:
- AC05-00OR22725; R-155-000-143-112
- Resource Type:
- Journal Article: Accepted Manuscript
- Journal Name:
- Stochastic Processes and Their Applications
- Additional Journal Information:
- Journal Name: Stochastic Processes and Their Applications; Journal ID: ISSN 0304-4149
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; 77 NANOSCIENCE AND NANOTECHNOLOGY; multilevel Monte Carlo; sequential Monte Carlo; Bayesian inverse problems
Citation Formats
Beskos, Alexandros, Jasra, Ajay, Law, Kody, Tempone, Raul, and Zhou, Yan. Multilevel sequential Monte Carlo samplers. United States: N. p., 2016.
Web. doi:10.1016/j.spa.2016.08.004.
Beskos, Alexandros, Jasra, Ajay, Law, Kody, Tempone, Raul, & Zhou, Yan. Multilevel sequential Monte Carlo samplers. United States. https://doi.org/10.1016/j.spa.2016.08.004
Beskos, Alexandros, Jasra, Ajay, Law, Kody, Tempone, Raul, and Zhou, Yan. 2016.
"Multilevel sequential Monte Carlo samplers". United States. https://doi.org/10.1016/j.spa.2016.08.004. https://www.osti.gov/servlets/purl/1302922.
@article{osti_1302922,
title = {Multilevel sequential Monte Carlo samplers},
author = {Beskos, Alexandros and Jasra, Ajay and Law, Kody and Tempone, Raul and Zhou, Yan},
abstractNote = {Here, we study the approximation of expectations w.r.t. probability distributions associated to the solution of partial differential equations (PDEs); this scenario appears routinely in Bayesian inverse problems. In practice, one often has to solve the associated PDE numerically, using, for instance finite element methods and leading to a discretisation bias, with the step-size level hL. In addition, the expectation cannot be computed analytically and one often resorts to Monte Carlo methods. In the context of this problem, it is known that the introduction of the multilevel Monte Carlo (MLMC) method can reduce the amount of computational effort to estimate expectations, for a given level of error. This is achieved via a telescoping identity associated to a Monte Carlo approximation of a sequence of probability distributions with discretisation levels ${\infty}$ >h0>h1 ...>hL. In many practical problems of interest, one cannot achieve an i.i.d. sampling of the associated sequence of probability distributions. A sequential Monte Carlo (SMC) version of the MLMC method is introduced to deal with this problem. In conclusion, it is shown that under appropriate assumptions, the attractive property of a reduction of the amount of computational effort to estimate expectations, for a given level of error, can be maintained within the SMC context.},
doi = {10.1016/j.spa.2016.08.004},
url = {https://www.osti.gov/biblio/1302922},
journal = {Stochastic Processes and Their Applications},
issn = {0304-4149},
number = ,
volume = ,
place = {United States},
year = {Wed Aug 24 00:00:00 EDT 2016},
month = {Wed Aug 24 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
A general theory of particle filters in hidden Markov models and some applications
journal, December 2013
- Chan, Hock Peng; Lai, Tze Leung
- The Annals of Statistics, Vol. 41, Issue 6
Central limit theorem for sequential Monte Carlo methods and its application to Bayesian inference
journal, December 2004
- Chopin, Nicolas
- The Annals of Statistics, Vol. 32, Issue 6
Multilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients
journal, January 2011
- Cliffe, K. A.; Giles, M. B.; Scheichl, R.
- Computing and Visualization in Science, Vol. 14, Issue 1
Sequential Monte Carlo samplers
journal, June 2006
- Del Moral, Pierre; Doucet, Arnaud; Jasra, Ajay
- Journal of the Royal Statistical Society: Series B (Statistical Methodology), Vol. 68, Issue 3
On adaptive resampling strategies for sequential Monte Carlo methods
journal, February 2012
- Del Moral, Pierre; Doucet, Arnaud; Jasra, Ajay
- Bernoulli, Vol. 18, Issue 1
Limit theorems for weighted samples with applications to sequential Monte Carlo methods
journal, October 2008
- Douc, Randal; Moulines, Eric
- The Annals of Statistics, Vol. 36, Issue 5
Multilevel Monte Carlo Path Simulation
journal, June 2008
- Giles, Michael B.
- Operations Research, Vol. 56, Issue 3
Monte Carlo Complexity of Global Solution of Integral Equations
journal, June 1998
- Heinrich, S.
- Journal of Complexity, Vol. 14, Issue 2
Complexity analysis of accelerated MCMC methods for Bayesian inversion
journal, July 2013
- Hoang, Viet Ha; Schwab, Christoph; Stuart, Andrew M.
- Inverse Problems, Vol. 29, Issue 8
Forward and Inverse Uncertainty Quantification Using Multilevel Monte Carlo Algorithms for an Elliptic Nonlocal Equation
journal, January 2016
- Jasra, Ajay; Law, Kody J. H.; Zhou, Yan
- International Journal for Uncertainty Quantification, Vol. 6, Issue 6
A Hierarchical Multilevel Markov Chain Monte Carlo Algorithm with Applications to Uncertainty Quantification in Subsurface Flow
journal, January 2015
- Dodwell, T. J.; Ketelsen, C.; Scheichl, R.
- SIAM/ASA Journal on Uncertainty Quantification, Vol. 3, Issue 1
Inference for a class of partially observed point process models
journal, August 2012
- Martin, James S.; Jasra, Ajay; McCoy, Emma
- Annals of the Institute of Statistical Mathematics, Vol. 65, Issue 3
Unbiased Estimation with Square Root Convergence for SDE Models
journal, October 2015
- Rhee, Chang-Han; Glynn, Peter W.
- Operations Research, Vol. 63, Issue 5
The Finite Element Method for Elliptic Problems
journal, December 1978
- Ciarlet, Philippe G.; Oden, J. T.
- Journal of Applied Mechanics, Vol. 45, Issue 4
On the stability of sequential Monte Carlo methods in high dimensions
journal, August 2014
- Beskos, Alexandros; Crisan, Dan; Jasra, Ajay
- The Annals of Applied Probability, Vol. 24, Issue 4
On adaptive resampling strategies for sequential Monte Carlo methods
text, January 2012
- Del Moral, Pierre; Doucet, Arnaud; Jasra, Ajay
- arXiv
A Hierarchical Multilevel Markov Chain Monte Carlo Algorithm with Applications to Uncertainty Quantification in Subsurface Flow
preprint, January 2013
- Dodwell, T. J.; Ketelsen, C.; Scheichl, R.
- arXiv
Works referencing / citing this record:
Unbiased multi-index Monte Carlo
journal, December 2017
- Crisan, Dan; Moral, Pierre Del; Houssineau, Jeremie
- Stochastic Analysis and Applications, Vol. 36, Issue 2
Multilevel sequential Monte Carlo: Mean square error bounds under verifiable conditions
journal, December 2016
- Del Moral, Pierre; Jasra, Ajay; Law, Kody J. H.
- Stochastic Analysis and Applications, Vol. 35, Issue 3
Multilevel Monte Carlo in approximate Bayesian computation
journal, January 2019
- Jasra, Ajay; Jo, Seongil; Nott, David
- Stochastic Analysis and Applications, Vol. 37, Issue 3
A transport-based multifidelity preconditioner for Markov chain Monte Carlo
journal, November 2019
- Peherstorfer, Benjamin; Marzouk, Youssef
- Advances in Computational Mathematics
Correction of coarse-graining errors by a two-level method: Application to the Asakura-Oosawa model
journal, October 2019
- Kobayashi, Hideki; Rohrbach, Paul B.; Scheichl, Robert
- The Journal of Chemical Physics, Vol. 151, Issue 14
Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
journal, February 2019
- Warne, David J.; Baker, Ruth E.; Simpson, Matthew J.
- Journal of The Royal Society Interface, Vol. 16, Issue 151
Correction of coarse-graining errors by a two-level method: Application to the Asakura-Oosawa model.
text, January 2019
- Kobayashi, Hideki; Rohrbach, Paul; Scheichl, Robert
- Apollo - University of Cambridge Repository
Unbiased Multi-index Monte Carlo
preprint, January 2017
- Crisan, Dan; Del Moral, Pierre; Houssineau, Jeremie
- arXiv
Multilevel Monte Carlo in Approximate Bayesian Computation
preprint, January 2017
- Jasra, Ajay; Jo, Seongil; Nott, David
- arXiv
A transport-based multifidelity preconditioner for Markov chain Monte Carlo
preprint, January 2018
- Peherstorfer, Benjamin; Marzouk, Youssef
- arXiv
Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
text, January 2018
- Warne, David J.; Baker, Ruth E.; Simpson, Matthew J.
- arXiv
Correction of coarse-graining errors by a two-level method: application to the Asakura-Oosawa model
text, January 2019
- Kobayashi, Hideki; Rohrbach, Paul B.; Scheichl, Robert
- arXiv
On coupling particle filter trajectories
journal, March 2017
- Sen, Deborshee; Thiery, Alexandre H.; Jasra, Ajay
- Statistics and Computing, Vol. 28, Issue 2
Unbiased estimation of the gradient of the log-likelihood in inverse problems
journal, March 2021
- Jasra, Ajay; Law, Kody J. H.; Lu, Deng
- Statistics and Computing, Vol. 31, Issue 3
Correction of coarse-graining errors by a two-level method: Application to the Asakura-Oosawa model
journal, October 2019
- Kobayashi, Hideki; Rohrbach, Paul B.; Scheichl, Robert
- The Journal of Chemical Physics, Vol. 151, Issue 14
A Wasserstein coupled particle filter for multilevel estimation
journal, June 2022
- Ballesio, Marco; Jasra, Ajay; von Schwerin, Erik
- Stochastic Analysis and Applications
Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
journal, February 2019
- Warne, David J.; Baker, Ruth E.; Simpson, Matthew J.
- Journal of The Royal Society Interface, Vol. 16, Issue 151
Multilevel Sequential Monte Carlo Samplers for Normalizing Constants
journal, July 2017
- Moral, Pierre Del; Jasra, Ajay; Law, Kody J. H.
- ACM Transactions on Modeling and Computer Simulation, Vol. 27, Issue 3
Multi-Index Sequential Monte Carlo Methods for Partially Observed Stochastic Partial Differential Equations
journal, January 2021
- Jasra, Ajay; Law, Kody J. H.; Xu, Yaxian
- International Journal for Uncertainty Quantification, Vol. 11, Issue 3
Multilevel ensemble Kalman filtering for spatio-temporal processes
text, January 2021
- Chernov, Alexey; Hoel, Hakon Andreas; Law, Kody J. H.
- RWTH Aachen University
Error bounds for sequential Monte Carlo samplers for multimodal distributions
journal, February 2019
- Paulin, Daniel; Jasra, Ajay; Thiery, Alexandre
- Bernoulli, Vol. 25, Issue 1
Unbiased Multi-index Monte Carlo
preprint, January 2017
- Crisan, Dan; Del Moral, Pierre; Houssineau, Jeremie
- arXiv
Multilevel Sequential Monte Carlo with Dimension-Independent Likelihood-Informed Proposals
preprint, January 2017
- Beskos, Alexandros; Jasra, Ajay; Law, Kody
- arXiv
Advanced Multilevel Monte Carlo Methods
preprint, January 2017
- Jasra, Ajay; Law, Kody; Suciu, Carina
- arXiv
Operator-Based Uncertainty Quantification of Stochastic Fractional PDEs
preprint, January 2018
- Kharazmi, Ehsan; Zayernouri, Mohsen
- arXiv
A practical example for the non-linear Bayesian filtering of model parameters
preprint, January 2018
- Bulté, Matthieu; Latz, Jonas; Ullmann, Elisabeth
- arXiv
Tree-based Particle Smoothing Algorithms in a Hidden Markov Model
preprint, January 2018
- Ding, Dong; Gandy, Axel
- arXiv
Vector operations for accelerating expensive Bayesian computations -- a tutorial guide
text, January 2019
- Warne, David J.; Sisson, Scott A.; Drovandi, Christopher
- arXiv
Multilevel adaptive sparse Leja approximations for Bayesian inverse problems
text, January 2019
- Farcas, Ionut-Gabriel; Latz, Jonas; Ullmann, Elisabeth
- arXiv
Multilevel Sequential Importance Sampling for Rare Event Estimation
text, January 2019
- Wagner, Fabian; Latz, Jonas; Papaioannou, Iason
- arXiv
On Unbiased Estimation for Discretized Models
preprint, January 2021
- Heng, Jeremy; Jasra, Ajay; Law, Kody J. H.
- arXiv