期刊
SIAM-ASA JOURNAL ON UNCERTAINTY QUANTIFICATION
卷 9, 期 1, 页码 29-64出版社
SIAM PUBLICATIONS
DOI: 10.1137/20M1318365
关键词
inverse problems; hierarchical Bayes; Markov chain Monte Carlo; pseudomarginalization; Poisson likelihood; positron emission tomography
资金
- Gordon Preston Fellowship by Monash University
- Australian Research Council [CE140100049]
- Australian Research Council [CE140100049] Funding Source: Australian Research Council
This work introduces scalable optimization-based Markov chain Monte Carlo (MCMC) methods for solving hierarchical Bayesian inverse problems with nonlinear parameter-to-observable maps and a broader class of hyperparameters. The algorithmic development is based on the recently developed scalable randomize-then-optimize (RTO) method and integrates RTO into Metropolis-within-Gibbs updates or pseudomarginal MCMC for efficient sampling in hierarchical Bayesian inversion. The integration of RTO and pseudomarginal MCMC provides sampling performance robust to model parameter dimensions, as demonstrated in numerical examples of PDE-constrained inverse problems and positron emission tomography.
In many hierarchical inverse problems, not only do we want to estimate high- or infinite-dimensional model parameters in the parameter-to-observable maps, but we also have to estimate hyperparameters that represent critical assumptions in the statistical and mathematical modeling processes. As a joint effect of high-dimensionality, nonlinear dependence, and nonconcave structures in the joint posterior distribution over model parameters and hyperparameters, solving inverse problems in the hierarchical Bayesian setting poses a significant computational challenge. In this work, we develop scalable optimization-based Markov chain Monte Carlo (MCMC) methods for solving hierarchical Bayesian inverse problems with nonlinear parameter-to-observable maps and a broader class of hyperparameters. Our algorithmic development is based on the recently developed scalable randomize-then-optimize (RTO) method [J. M. Bardsley et al., SIAM J. Sci. Comput., 42 (2016), pp. A1317-A1347] for exploring the high- or infinite-dimensional parameter space. We first extend the RTO machinery to the Poisson likelihood and discuss the implementation of RTO in the hierarchical setting. Then, by using RTO either as a proposal distribution in a Metropolis-withinGibbs update or as a biasing distribution in the pseudomarginal MCMC [C. Andrieu and G. O. Roberts, Ann. Statist., 37 (2009), pp. 697-725], we present efficient sampling tools for hierarchical Bayesian inversion. In particular, the integration of RTO and the pseudomarginal MCMC has sampling performance robust to model parameter dimensions. Numerical examples in PDE-constrained inverse problems and positron emission tomography are used to demonstrate the performance of our methods.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据