GPS-ABC: Gaussian Process Surrogate Approximate Bayesian Computation

01/13/2014

∙

Scientists often express their understanding of the world through a computationally demanding simulation program. Analyzing the posterior distribution of the parameters given observations (the inverse problem) can be extremely challenging. The Approximate Bayesian Computation (ABC) framework is the standard statistical tool to handle these likelihood free problems, but they require a very large number of simulations. In this work we develop two new ABC sampling algorithms that significantly reduce the number of simulations necessary for posterior inference. Both algorithms use confidence estimates for the accept probability in the Metropolis Hastings step to adaptively choose the number of necessary simulations. Our GPS-ABC algorithm stores the information obtained from every simulation in a Gaussian process which acts as a surrogate function for the simulated statistics. Experiments on a challenging realistic biological problem illustrate the potential of these algorithms.

READ FULL TEXT

GPS-ABC: Gaussian Process Surrogate Approximate Bayesian Computation

Introduction To Gaussian Process Regression In Bayesian Inverse Problems, With New ResultsOn Experimental Design For Weighted Error Measures

POPE: Post Optimization Posterior Evaluation of Likelihood Free Models

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

Bayesian optimisation for likelihood-free cosmological inference

Accelerating Approximate Bayesian Computation with Quantile Regression: Application to Cosmological Redshift Distributions

Sampling of Bayesian posteriors with a non-Gaussian probabilistic learning on manifolds from a small dataset

Emulating satellite drag from large simulation experiments

GPS-ABC: Gaussian Process Surrogate Approximate Bayesian Computation

Related Research

Introduction To Gaussian Process Regression In Bayesian Inverse Problems, With New ResultsOn Experimental Design For Weighted Error Measures

POPE: Post Optimization Posterior Evaluation of Likelihood Free Models

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

Bayesian optimisation for likelihood-free cosmological inference

Accelerating Approximate Bayesian Computation with Quantile Regression: Application to Cosmological Redshift Distributions

Sampling of Bayesian posteriors with a non-Gaussian probabilistic learning on manifolds from a small dataset

Emulating satellite drag from large simulation experiments