Maximum a posteriori probability estimates in infinite-dimensional Bayesian inverse problems

被引:50
作者
Helin, T. [1 ]
Burger, M. [2 ,3 ]
机构
[1] Univ Helsinki, FI-00014 Helsinki, Finland
[2] Univ Munster, Inst Computat & Appl Math, D-48149 Munster, Germany
[3] Cells Mot Cluster Excellence, D-48149 Munster, Germany
基金
芬兰科学院; 欧洲研究理事会;
关键词
inverse problems; Bayesian; maximum a posteriori; non-Gaussian; Besov; FUNCTIONAL BREGMAN DIVERGENCE; PRIORS; DISTRIBUTIONS;
D O I
10.1088/0266-5611/31/8/085009
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
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.
引用
收藏
页数:22
相关论文
empty
未找到相关数据