Batch simulations and uncertainty quantification in Gaussian process surrogate approximate Bayesian computation

by   Marko Järvenpää, et al.

Surrogate models such as Gaussian processes (GP) have been proposed to accelerate approximate Bayesian computation (ABC) when the statistical model of interest is expensive-to-simulate. In one such promising framework the discrepancy between simulated and observed data is modelled with a GP which is further used to form a model-based estimator for the intractable posterior. In this article we improve this approach in several ways. We develop batch-sequential Bayesian experimental design strategies to parallellise the expensive simulations. In earlier work only sequential strategies have been used. Current surrogate-based ABC methods also do not fully account the uncertainty due to the limited budget of simulations as they output only a point estimate of the ABC posterior. We propose a numerical method to fully quantify the uncertainty in, for example, ABC posterior moments. We also provide some new analysis on the GP modelling assumptions in the resulting improved framework called Bayesian ABC and on its connection to Bayesian quadrature (BQ) and Bayesian optimisation (BO). Experiments with several toy and real-world simulation models demonstrate advantages of the proposed techniques.


Batch simulations and uncertainty quantification in Gaussian process surrogate-based approximate Bayesian computation

Surrogate models such as Gaussian processes (GP) have been proposed to a...

Leveraging Bayesian Analysis To Improve Accuracy of Approximate Models

We focus on improving the accuracy of an approximate model of a multisca...

Parallel Gaussian process surrogate method to accelerate likelihood-free inference

We consider Bayesian inference when only a limited number of noisy log-l...

Gaussian process modeling in approximate Bayesian computation to estimate horizontal gene transfer in bacteria

Approximate Bayesian computation (ABC) can be used for model fitting whe...

Hierarchical shrinkage Gaussian processes: applications to computer code emulation and dynamical system recovery

In many areas of science and engineering, computer simulations are widel...

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

Several applications such as nuclear forensics, nuclear fuel cycle simul...

Please sign up or login with your details

Forgot password? Click here to reset