Abstract:In various industrial contexts, estimating the distribution of unobserved random vectors Xi from some noisy indirect observations H(Xi) + Ui is required. If the relation between Xi and the quantity H(Xi), measured with the error Ui, is implemented by a CPU-consuming computer model H, a major practical difficulty is to perform the statistical inference with a relatively small number of runs of H. Following Fu et al. (2014), a Bayesian statistical framework is considered to make use of possible prior knowledge on the parameters of the distribution of the Xi, which is assumed Gaussian. Moreover, a Markov Chain Monte Carlo (MCMC) algorithm is carried out to estimate their posterior distribution by replacing H by a kriging metamodel build from a limited number of simulated experiments. Two heuristics, involving two different criteria to be optimized, are proposed to sequentially design these computer experiments in the limits of a given computational budget. The first criterion is a Weighted Integrated Mean Square Error (WIMSE). The second one, called Expected Conditional Divergence (ECD), developed in the spirit of the Stepwise Uncertainty Reduction (SUR) criterion, is based on the discrepancy between two consecutive approximations of the target posterior distribution. Several numerical comparisons conducted over a toy example then a motivating real case-study show that such adaptive designs can significantly outperform the classical choice of a maximin Latin Hypercube Design (LHD) of experiments. Dealing with a major concern in hydraulic engineering, a particular emphasis is placed upon the prior elicitation of the case-study, highlighting the overall feasibility of the methodology. Faster convergences and manageability considerations lead to recommend the use of the ECD criterion in practical applications.

Rational Kriging

Rational Kriging

Rational Quadratic Approximation to Real Algebraic Curves

Least Squares Rational Approximation

Rational function approximation with normalized positive denominators

A geostatistical model based on random walks to krige regions with irregular boundaries and holes

Rational Krylov methods for functions of matrices with applications to fractional partial differential equations

Non-Homogeneity Estimation and Universal Kriging on the Sphere

Error bounds for the approximation of matrix functions with rational Krylov methods

Rational inner functions on the polydisk -- a survey

Rational Maximum Likelihood Estimators of Kronecker Covariance Matrices

Scalable Stochastic Kriging with Markovian Covariances

Robust Parameter Estimation for Rational Ordinary Differential Equations

An algorithm for two‐variable rational interpolation suitable for matrix manipulations with the evaluation–interpolation method

Kriging: a method of interpolation for geographical information systems

Polynomial-Chaos-based Kriging

Resolution Adaptive Fixed Rank Kriging

An adaptive kriging method for solving nonlinear inverse statistical problems

Enhanced Universal Kriging for Transformed Input Parameter Spaces

Kryging: Geostatistical analysis of large-scale datasets using Krylov subspace methods

Rational minimax approximation via adaptive barycentric representations