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

by   Andreas Stein, et al.

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.


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

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

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

We present Bayesian techniques for solving inverse problems which involv...

Bayesian Inversion of Log-normal Eikonal Equations

We study the Bayesian inverse problem for inferring the log-normal slown...

Monte Carlo convergence rates for kth moments in Banach spaces

We formulate standard and multilevel Monte Carlo methods for the kth mom...

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

Random eigenvalue problems are useful models for quantifying the uncerta...

Multilevel Delayed Acceptance MCMC

We develop a novel Markov chain Monte Carlo (MCMC) method that exploits ...

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