Laplace priors and spatial inhomogeneity in Bayesian inverse problems

12/10/2021
by   Sergios Agapiou, et al.
0

Spatially inhomogeneous functions, which may be smooth in some regions and rough in other regions, are modelled naturally in a Bayesian manner using so-called Besov priors which are given by random wavelet expansions with Laplace-distributed coefficients. This paper studies theoretical guarantees for such prior measures - specifically, we examine their frequentist posterior contraction rates in the setting of non-linear inverse problems with Gaussian white noise. Our results are first derived under a general local Lipschitz assumption on the forward map. We then verify the assumption for two non-linear inverse problems arising from elliptic partial differential equations, the Darcy flow model from geophysics as well as a model for the Schrödinger equation appearing in tomography. In the course of the proofs, we also obtain novel concentration inequalities for penalized least squares estimators with ℓ^1 wavelet penalty, which have a natural interpretation as maximum a posteriori (MAP) estimators. The true parameter is assumed to belong to some spatially inhomogeneous Besov class B^α_11, α>0. In a setting with direct observations, we complement these upper bounds with a lower bound on the rate of contraction for arbitrary Gaussian priors. An immediate consequence of our results is that while Laplace priors can achieve minimax-optimal rates over B^α_11-classes, Gaussian priors are limited to a (by a polynomial factor) slower contraction rate. This gives information-theoretical justification for the intuition that Laplace priors are more compatible with ℓ^1 regularity structure in the underlying parameter.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/25/2018

Bayesian linear inverse problems in regularity scales

We obtain rates of contraction of posterior distributions in inverse pro...
research
09/13/2022

Adaptive inference over Besov spaces in the white noise model using p-exponential priors

In many scientific applications the aim is to infer a function which is ...
research
08/30/2022

Besov priors in density estimation: optimal posterior contraction rates and adaptation

Besov priors are nonparametric priors that model spatially inhomogeneous...
research
09/24/2018

Convergence rates for Penalised Least Squares Estimators in PDE-constrained regression problems

We consider PDE constrained nonparametric regression problems in which t...
research
10/04/2018

Posterior contraction for empirical Bayesian approach to inverse problems under non-diagonal assumption

We investigate an empirical Bayesian nonparametric approach to a family ...
research
06/28/2023

Spatiotemporal Besov Priors for Bayesian Inverse Problems

Fast development in science and technology has driven the need for prope...
research
07/20/2021

On some information-theoretic aspects of non-linear statistical inverse problems

Results by van der Vaart (1991) from semi-parametric statistics about th...

Please sign up or login with your details

Forgot password? Click here to reset