{"title":"Variance-based sensitivity of Bayesian inverse problems to the prior distribution","authors":"John Darges, Alen Alexanderian, Pierre Gremaud","doi":"10.1615/int.j.uncertaintyquantification.2024051475","DOIUrl":null,"url":null,"abstract":"The formulation of Bayesian inverse problems involves choosing prior\ndistributions; choices that seem equally reasonable may lead to significantly\ndifferent conclusions. We develop a computational approach to better\nunderstand the impact of the hyperparameters defining the prior on the\nposterior statistics of the quantities of interest. Our approach relies on\nglobal sensitivity analysis (GSA) of Bayesian inverse problems with respect to\nthe hyperparameters defining the prior. This, however, is a challenging\nproblem---a naive double loop sampling approach would require running a prohibitive\nnumber of Markov chain Monte Carlo (MCMC) sampling procedures. The present\nwork takes a foundational step in making such a sensitivity analysis practical\nthrough (i) a judicious combination of efficient surrogate models and (ii) a\ntailored importance sampling method. In particular, we can perform accurate\nGSA of posterior prediction statistics with respect to prior hyperparameters\nwithout having to repeat MCMC runs. We demonstrate the effectiveness of the\napproach on a simple Bayesian linear inverse problem and a nonlinear inverse\nproblem governed by an epidemiological model.","PeriodicalId":48814,"journal":{"name":"International Journal for Uncertainty Quantification","volume":"9 1","pages":""},"PeriodicalIF":1.5000,"publicationDate":"2024-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"International Journal for Uncertainty Quantification","FirstCategoryId":"5","ListUrlMain":"https://doi.org/10.1615/int.j.uncertaintyquantification.2024051475","RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"ENGINEERING, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
The formulation of Bayesian inverse problems involves choosing prior
distributions; choices that seem equally reasonable may lead to significantly
different conclusions. We develop a computational approach to better
understand the impact of the hyperparameters defining the prior on the
posterior statistics of the quantities of interest. Our approach relies on
global sensitivity analysis (GSA) of Bayesian inverse problems with respect to
the hyperparameters defining the prior. This, however, is a challenging
problem---a naive double loop sampling approach would require running a prohibitive
number of Markov chain Monte Carlo (MCMC) sampling procedures. The present
work takes a foundational step in making such a sensitivity analysis practical
through (i) a judicious combination of efficient surrogate models and (ii) a
tailored importance sampling method. In particular, we can perform accurate
GSA of posterior prediction statistics with respect to prior hyperparameters
without having to repeat MCMC runs. We demonstrate the effectiveness of the
approach on a simple Bayesian linear inverse problem and a nonlinear inverse
problem governed by an epidemiological model.
期刊介绍:
The International Journal for Uncertainty Quantification disseminates information of permanent interest in the areas of analysis, modeling, design and control of complex systems in the presence of uncertainty. The journal seeks to emphasize methods that cross stochastic analysis, statistical modeling and scientific computing. Systems of interest are governed by differential equations possibly with multiscale features. Topics of particular interest include representation of uncertainty, propagation of uncertainty across scales, resolving the curse of dimensionality, long-time integration for stochastic PDEs, data-driven approaches for constructing stochastic models, validation, verification and uncertainty quantification for predictive computational science, and visualization of uncertainty in high-dimensional spaces. Bayesian computation and machine learning techniques are also of interest for example in the context of stochastic multiscale systems, for model selection/classification, and decision making. Reports addressing the dynamic coupling of modern experiments and modeling approaches towards predictive science are particularly encouraged. Applications of uncertainty quantification in all areas of physical and biological sciences are appropriate.