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.

On gamma estimation via matrix kriging

Quasi-Monte Carlo methods for calculating derivatives sensitivities on the GPU

Chebyshev Greeks: Smoothing Gamma without Bias

Monte Carlo Greeks for financial products via approximative transition densities

Kernel Estimation of the Greeks for Options with Discontinuous Payoffs

ACCURATE AND EFFICIENT COMPUTATIONS FOR THE GREEKS OF EUROPEAN MULTI-ASSET OPTIONS

Real-Time Derivative Pricing and Hedging with Consistent Metamodels

Efficient global estimation of conditional-value-at-risk through stochastic kriging and extreme value theory

Derivatives Sensitivities Computation under Heston Model on GPU

Sliced gradient-enhanced Kriging for high-dimensional function approximation

Pricing and Risk Management with High-Dimensional Quasi Monte Carlo and Global Sensitivity Analysis

Efficient Computation of Option Prices and Greeks by Quasi-Monte Carlo Method with Smoothing and Dimension Reduction.

Scalable Stochastic Kriging with Markovian Covariances

Gamma Hedging and Rough Paths

An adaptive kriging method for solving nonlinear inverse statistical problems

Quasi-Monte Carlo-Based Conditional Malliavin Method for Continuous-Time Asian Option Greeks

Failure of Fourier pricing techniques to approximate the Greeks

Generalized Exponentiated Gradient Algorithms and Their Application to On-Line Portfolio Selection

Inf-convolution of g_Γ-solution and its applications

Quasi-Monte Carlo Simulation for American Option Sensitivities

Convergence of Milstein Brownian bridge Monte Carlo methods and stable Greeks calculation