Quantifying uncertainties on excursion sets under a Gaussian random field prior

01/15/2015
by   Dario Azzimonti, et al.
0

We focus on the problem of estimating and quantifying uncertainties on the excursion set of a function under a limited evaluation budget. We adopt a Bayesian approach where the objective function is assumed to be a realization of a Gaussian random field. In this setting, the posterior distribution on the objective function gives rise to a posterior distribution on excursion sets. Several approaches exist to summarize the distribution of such sets based on random closed set theory. While the recently proposed Vorob'ev approach exploits analytical formulae, further notions of variability require Monte Carlo estimators relying on Gaussian random field conditional simulations. In the present work we propose a method to choose Monte Carlo simulation points and obtain quasi-realizations of the conditional field at fine designs through affine predictors. The points are chosen optimally in the sense that they minimize the posterior expected distance in measure between the excursion set and its reconstruction. The proposed method reduces the computational costs due to Monte Carlo simulations and enables the computation of quasi-realizations on fine designs in large dimensions. We apply this reconstruction approach to obtain realizations of an excursion set on a fine grid which allow us to give a new measure of uncertainty based on the distance transform of the excursion set. Finally we present a safety engineering test case where the simulation method is employed to compute a Monte Carlo estimate of a contour line.

READ FULL TEXT
research
12/15/2020

On the Selection of Random Field Evaluation Points in the p-MLQMC Method

Engineering problems are often characterized by significant uncertainty ...
research
11/01/2018

Stochastic turbulence modeling in RANS simulations via Multilevel Monte Carlo

A multilevel Monte Carlo (MLMC) method for quantifying model-form uncert...
research
06/25/2019

h- and p-refined Multilevel Monte Carlo Methods for Uncertainty Quantification in Structural Engineering

Practical structural engineering problems are often characterized by sig...
research
07/15/2022

Quasi-Monte Carlo and discontinuous Galerkin

In this study, we design and develop Quasi-Monte Carlo (QMC) cubatures f...
research
01/19/2020

Distributionally Robust Bayesian Quadrature Optimization

Bayesian quadrature optimization (BQO) maximizes the expectation of an e...
research
04/01/2015

Bayesian model comparison with un-normalised likelihoods

Models for which the likelihood function can be evaluated only up to a p...
research
02/19/2023

Identifying Heterogeneity in Regression Compositional Data Integration with Many Categories

In compositional data, detecting which part of the whole delineates hete...

Please sign up or login with your details

Forgot password? Click here to reset