Maximum a posteriori probability estimates in infinite-dimensional Bayesian inverse problems. (17th July 2015)

• Record Type: Journal Article
• Authors: Helin, T; Burger, M
• Is Part Of: Inverse problems. Volume 31:Number 8(2015:Aug.)
• Publication Date: 2015-07-17
• Subjects: inverse problems -- Bayesian -- maximum a posteriori -- non-Gaussian -- Besov; 65J22 -- 62F15 -- 62C10 -- 65C60 -- 62F10; Inverse problems (Differential equations) -- Periodicals; 515.357
• DOI: 10.1088/0266-5611/31/8/085009
• Languages: English
• ISSNs: 0266-5611

Abstract

Abstract: A demanding challenge in Bayesian inversion is to efficiently characterize the posterior distribution. This task is problematic especially in high-dimensional non-Gaussian problems, where the structure of the posterior can be very chaotic and difficult to analyse. Current inverse problem literature often approaches the problem by considering suitable point estimators for the task. Typically the choice is made between the maximum a posteriori (MAP) or the conditional mean (CM) estimate. The benefits of either choice are not well-understood from the perspective of infinite-dimensional theory. Most importantly, there exists no general scheme regarding how to connect the topological description of a MAP estimate to a variational problem. The recent results by Dashti and others (Dashti et al 2013 Inverse Problems 29 095017 ) resolve this issue for nonlinear inverse problems in Gaussian framework. In this work we improve the current understanding by introducing a novel concept called the weak MAP (wMAP) estimate. We show that any MAP estimate in the sense of Dashti et al (2013 Inverse Problems 29 095017 ) is a wMAP estimate and, moreover, how the wMAP estimate connects to a variational formulation in general infinite-dimensional non-Gaussian problems. The variational formulation enables to study many properties of the infinite-dimensional MAP estimate that were earlier impossible to study. In a recent work by the authors (Burger and Lucka 2014 Maximum a posteriori estimates in linear inverse problems with logconcave priors are proper bayes estimators preprint) the MAP estimator was studied in the context of the Bayes cost method. Using Bregman distances, proper convex Bayes cost functions were introduced for which the MAP estimator is the Bayes estimator. Here, we generalize these results to the infinite-dimensional setting. Moreover, we discuss the implications of our results for some examples of prior models such as the Besov prior and hierarchical prior.