MAP estimators and their consistency in Bayesian nonparametric inverse problems
Abstract
We consider the inverse problem of estimating an unknown function u from noisy measurements y of a known, possibly nonlinear, map {G} applied to u. We adopt a Bayesian approach to the problem and work in a setting where the prior measure is specified as a Gaussian random field μ0. We work under a natural set of conditions on the likelihood which implies the existence of a well-posed posterior measure, μy. Under these conditions, we show that the maximum a posteriori (MAP) estimator is well defined as the minimizer of an Onsager-Machlup functional defined on the Cameron-Martin space of the prior; thus, we link a problem in probability with a problem in the calculus of variations. We then consider the case where the observational noise vanishes and establish a form of Bayesian posterior consistency for the MAP estimator. We also prove a similar result for the case where the observation of {G}(u) can be repeated as many times as desired with independent identically distributed noise. The theory is illustrated with examples from an inverse problem for the Navier-Stokes equation, motivated by problems arising in weather forecasting, and from the theory of conditioned diffusions, motivated by problems arising in molecular dynamics.
- Publication:
-
Inverse Problems
- Pub Date:
- September 2013
- DOI:
- 10.1088/0266-5611/29/9/095017
- arXiv:
- arXiv:1303.4795
- Bibcode:
- 2013InvPr..29i5017D
- Keywords:
-
- Mathematics - Probability;
- 49N45 60G35 62F15
- E-Print:
- changed title, minor fixes