Sparse Deterministic Approximation of Bayesian Inverse Problems

Abstract

We present a parametric deterministic formulation of Bayesian inverse problems with input parameter from infinite dimensional, separable Banach spaces. In this formulation, the forward problems are parametric, deterministic elliptic partial differential equations, and the inverse problem is to determine the unknown, parametric deterministic coefficients from noisy observations comprising linear functionals of the solution. We prove a generalized polynomial chaos representation of the posterior density with respect to the prior measure, given noisy observational data. We analyze the sparsity of the posterior density in terms of the summability of the input data's coefficient sequence. To this end, we estimate the fluctuations in the prior. We exhibit sufficient conditions on the prior model in order for approximations of the posterior density to converge at a given algebraic rate, in terms of the number $N$ of unknowns appearing in the parameteric representation of the prior measure. Similar sparsity and approximation results are also exhibited for the solution and covariance of the elliptic partial differential equation under the posterior. These results then form the basis for efficient uncertainty quantification, in the presence of data with noise.

Cited by in corpus (30)

• Spectral gaps for a Metropolis-Hastings algorithm in infinite dimensions
• Optimal scaling and diffusion limits for the Langevin algorithm in high dimensions
• Scalable posterior approximations for large-scale Bayesian inverse problems via likelihood-informed parameter and state reduction
• On the well-posedness of Bayesian inverse problems
• Posterior Consistency for Bayesian Inverse Problems through Stability and Regression Results
• Hessian-based adaptive sparse quadrature for infinite-dimensional Bayesian inverse problems
• Spectral likelihood expansions for Bayesian inference
• Random fields representations for stochastic elliptic boundary value problems and statistical inverse problems
• Bayesian parameter identification in Cahn-Hilliard models for biological growth
• Projected Stein Variational Gradient Descent
• Higher order Quasi-Monte Carlo integration for Bayesian Estimation
• Generalized modes in Bayesian inverse problems
• Inverse problem for the wave equation with a white noise source
• Bayesian Analysis of ODE's: solver optimal accuracy and Bayes factors
• Adaptive sampling-based quadrature rules for efficient Bayesian prediction
• Cell Detection by Functional Inverse Diffusion and Non-negative Group Sparsity$-$Part I: Modeling and Inverse Problems
• Multilevel adaptive sparse Leja approximations for Bayesian inverse problems
• Infinite-dimensional Bayesian approach for inverse scattering problems of a fractional Helmholtz equation
• A path-integral approach to Bayesian inference for inverse problems using the semiclassical approximation
• Variational Inference for Nonlinear Inverse Problems via Neural Net Kernels: Comparison to Bayesian Neural Networks, Application to Topology Optimization
• Efficient approximation of high-dimensional exponentials by tensornetworks
• Bregman Cost for Non-Gaussian Noise
• Posterior Integration on a Riemannian Manifold
• Stein variational reduced basis Bayesian inversion
• Sequential Monte Carlo Methods for Bayesian Elliptic Inverse Problems
• Higher order Quasi-Monte Carlo integration for holomorphic, parametric operator equations
• Quasi-Monte Carlo and Multilevel Monte Carlo Methods for Computing Posterior Expectations in Elliptic Inverse Problems
• Complexity Analysis of Accelerated MCMC Methods for Bayesian Inversion
• Stochastic Galerkin finite element method with local conductivity basis for electrical impedance tomography
• Polynomial collocation for handling an inaccurately known measurement configuration in electrical impedance tomography