Bayesian inference for inverse problems occurring in uncertainty analysis

Shuai Fu 1 Gilles Celeux 1 Nicolas Bousquet 2 Mathieu Couplet 2
1 SELECT - Model selection in statistical learning
LMO - Laboratoire de Mathématiques d'Orsay, INRIA Saclay - Ile de France, CNRS - Centre National de la Recherche Scientifique : UMR
Abstract : The inverse problem considered here is to estimate the distribution of a non-observed random variable $X$ from some noisy observed data $Y$ linked to $X$ through a time-consuming physical model $H$. Bayesian inference is considered to take into account prior expert knowledge on $X$ in a small sample size setting. A Metropolis-Hastings within Gibbs algorithm is proposed to compute the posterior distribution of the parameters of $X$ through a data augmentation process. Since calls to $H$ are quite expensive, this inference is achieved by replacing $H$ with a kriging emulator interpolating $H$ from a numerical design of experiments. This approach involves several errors of different nature and, in this paper, we pay effort to measure and reduce the possible impact of those errors. In particular, we propose to use the so-called DAC criterion to assess in the same exercise the relevance of the numerical design and the prior distributions. After describing how computing this criterion for the emulator at hand, its behavior is illustrated on numerical experiments.
Document type :
[Research Report] RR-7995, INRIA. 2012
Contributor : Shuai Fu <>
Submitted on : Friday, June 15, 2012 - 6:16:51 PM
Last modification on : Saturday, June 16, 2012 - 2:53:46 PM
Document(s) archivé(s) le : Sunday, September 16, 2012 - 3:05:23 AM


Files produced by the author(s)


  • HAL Id : hal-00708814, version 1



Shuai Fu, Gilles Celeux, Nicolas Bousquet, Mathieu Couplet. Bayesian inference for inverse problems occurring in uncertainty analysis. [Research Report] RR-7995, INRIA. 2012. <hal-00708814>




Record views


Document downloads