Publication de 6 numéros par an
ISSN Imprimer: 2152-5080
ISSN En ligne: 2152-5099
Indexed in
A MULTI-INDEX MARKOV CHAIN MONTE CARLO METHOD
RÉSUMÉ
In this paper, we consider computing expectations with respect to probability laws associated with a certain class of stochastic systems. In order to achieve such a task, one must not only resort to numerical approximation of the expectation but also to a biased discretization of the associated probability. We are concerned with the situation for which the discretization is required in multiple dimensions, for instance in space-time. In such contexts, it is known that the multi-index Monte Carlo (MIMC) method of Haji-Ali, Nobile, and Tempone, (Numer. Math., 132, pp. 767– 806, 2016) can improve on independent identically distributed (i.i.d.) sampling from the most accurate approximation of the probability law. Through a nontrivial modification of the multilevel Monte Carlo (MLMC) method, this method can reduce the work to obtain a given level of error, relative to i.i.d. sampling and even to MLMC. In this paper, we consider the case when such probability laws are too complex to be sampled independently, for example a Bayesian inverse problem where evaluation of the likelihood requires solution of a partial differential equation model, which needs to be approximated at finite resolution. We develop a modification of the MIMC method, which allows one to use standard Markov chain Monte Carlo (MCMC) algorithms to replace independent and coupled sampling, in certain contexts. We prove a variance theorem for a simplified estimator that shows that using our MIMCMC method is preferable, in the sense above, to i.i.d. sampling from the most accurate approximation, under appropriate assumptions. The method is numerically illustrated on a Bayesian inverse problem associated to a stochastic partial differential equation, where the path measure is conditioned on some observations.
-
Hoang Viet Ha, Quek Jia Hao, Schwab Christoph, Analysis of a multilevel Markov chain Monte Carlo finite element method for Bayesian inversion of log-normal diffusions, Inverse Problems, 36, 3, 2020. Crossref
-
Zhang Jiaxin, Modern Monte Carlo methods for efficient uncertainty quantification and propagation: A survey, WIREs Computational Statistics, 13, 5, 2021. Crossref
-
Jasra Ajay, Law Kody J. H., Xu Yaxian, Markov chain simulation for multilevel Monte Carlo, Foundations of Data Science, 3, 1, 2021. Crossref
-
Chada Neil K., Franks Jordan, Jasra Ajay, Law Kody J., Vihola Matti, Unbiased Inference for Discretely Observed Hidden Markov Model Diffusions, SIAM/ASA Journal on Uncertainty Quantification, 9, 2, 2021. Crossref
-
Jasra Ajay, Law Kody J. H., Tarakanov Alexander, Yu Fangyuan, Randomized Multilevel Monte Carlo for Embarrassingly Parallel Inference, in Driving Scientific and Engineering Discoveries Through the Integration of Experiment, Big Data, and Modeling and Simulation, 1512, 2022. Crossref
-
Yang Juntao, Hoang Viet Ha, Multilevel Markov Chain Monte Carlo for Bayesian inverse problem for Navier-Stokes equation, Inverse Problems and Imaging, 2022. Crossref
-
Davis Andrew D., Marzouk Youssef, Smith Aaron, Pillai Natesh, Rate-optimal refinement strategies for local approximation MCMC, Statistics and Computing, 32, 4, 2022. Crossref