[edit]
GPS-ABC: Gaussian Process Surrogate Approximate Bayesian Computation
Proceedings of the 30th Conference on Uncertainty in Artificial Intelligence, PMLR R12:598-607, 2014.
Abstract
Scientists often express their understanding of the world through a computationally demand- ing simulation program. Analyzing the posterior distribution of the parameters given observations (the inverse problem) can be extremely chal- lenging. The Approximate Bayesian Computa- tion (ABC) framework is the standard statisti- cal tool to handle these likelihood free problems, but they require a very large number of simula- tions. In this work we develop two new ABC sampling algorithms that significantly reduce the number of simulations necessary for posterior in- ference. Both algorithms use confidence esti- mates for the accept probability in the Metropo- lis Hastings step to adaptively choose the number of necessary simulations. Our GPS-ABC algo- rithm 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 biologi- cal problem illustrate the potential of these algo- rithms.