.png){kind=icon}
Return
Efficient hyperparameter estimation in Bayesian inverse problems using sample average approximation
.png){kind=icon}
PRE
AI
DOI：10.1098/rsta.2024.0046.png)
.png){kind=icon}
Original request for help.png){kind=icon}
Share.png){kind=icon}
Save
Abstract
En 中文
In Bayesian inverse problems, it is common to consider several hyperparameters that define the prior and the noise model that must be estimated from the data. In particular, we are interested in linear inverse problems with additive Gaussian noise and Gaussian priors defined using Mat & eacute;rn covariance models. In this case, we estimate the hyperparameters using the maximum a posteriori (MAP) estimate of the marginalized posterior distribution. However, this is a computationally intensive task since it involves computing log determinants. To address this challenge, we consider a stochastic average approximation (SAA) of the objective function and use the preconditioned Lanczos method to compute efficient approximations of the function and gradient evaluations. We propose a new preconditioner that can be updated cheaply for new values of the hyperparameters and an approach to compute approximations of the gradient evaluations, by reutilizing information from the function evaluations. We demonstrate the performance of our approach on static and dynamic seismic tomography problems.This article is part of the theme issue 'Frontiers of applied inverse problems in science and engineering'.
Keywords:
Monte Carlo
preconditioning
inverse problems
Journal
P
IF:
3.7
Papers:
7.8K
Citations:
2.8W
Organization
View more organizations
Cited Papers
.png){kind=icon}
ShareSave
ShareSaveShareSaveShareSaveEfficient iterative methods for hyperparameter estimation in large-scale linear inverse problems
PREAI
ShareSaveShareSave View more