Multilevel Markov Chain Monte Carlo for Bayesian Elliptic Inverse Problems with Besov Random Tree Priors

02/01/2023
by   Andreas Stein, et al.
0

We propose a multilevel Monte Carlo-FEM algorithm to solve elliptic Bayesian inverse problems with "Besov random tree prior". These priors are given by a wavelet series with stochastic coefficients, and certain terms in the expansion vanishing at random, according to the law of so-called Galton-Watson trees. This allows to incorporate random fractal structures and large deviations in the log-diffusion, which occur naturally in many applications from geophysics or medical imaging. This framework entails two main difficulties: First, the associated diffusion coefficient does not satisfy a uniform ellipticity condition, which leads to non-integrable terms and thus divergence of standard multilevel estimators. Secondly, the associated space of parameters is Polish, but not a normed linear space. We address the first point by introducing cut-off functions in the estimator to compensate for the non-integrable terms, while the second issue is resolved by employing an independence Metropolis-Hastings sampler. The resulting algorithm converges in the mean-square sense with essentially optimal asymptotic complexity, and dimension-independent acceptance probabilities.

READ FULL TEXT
02/01/2023

Multilevel Monte Carlo FEM for Elliptic PDEs with Besov Random Tree Priors

We develop a multilevel Monte Carlo (MLMC)-FEM algorithm for linear, ell...
11/05/2021

Sampling Methods for Bayesian Inference Involving Convergent Noisy Approximations of Forward Maps

We present Bayesian techniques for solving inverse problems which involv...
11/22/2021

Bayesian Inversion of Log-normal Eikonal Equations

We study the Bayesian inverse problem for inferring the log-normal slown...
12/07/2022

Monte Carlo convergence rates for kth moments in Banach spaces

We formulate standard and multilevel Monte Carlo methods for the kth mom...
10/02/2020

Multilevel quasi-Monte Carlo for random elliptic eigenvalue problems I: Regularity and error analysis

Random eigenvalue problems are useful models for quantifying the uncerta...
02/08/2022

Multilevel Delayed Acceptance MCMC

We develop a novel Markov chain Monte Carlo (MCMC) method that exploits ...
03/10/2022

Multi-index Sequential Monte Carlo ratio estimators for Bayesian Inverse problems

We consider the problem of estimating expectations with respect to a tar...

Please sign up or login with your details

Forgot password? Click here to reset