ABC-EP: Expectation Propagation for Likelihoodfree Bayesian Computation

Many statistical models of interest to the natural and social sciences have no tractable likelihood function. Until recently, Bayesian inference for such models was thought in-feasible. Pritchard et al. (1999) introduced an algorithm known as ABC, for Approximate Bayesian Computation, that enables Bayesian computation in such models. Despite steady progress since this first breakthrough, such as the adaptation of MCMC and Sequential Monte Carlo techniques to likelihood-free inference, state-of-the art methods remain hard to use and require enormous computation times. Among other issues, one faces the difficult task of finding appropriate summary statistics for the model, and tuning the algorithm can be time-consuming when little prior information is available. We show that Expectation Propagation, a widely successful approximate inference technique, can be adapted to the likelihood-free context. The resulting algorithm does not require summary statistics, is an order of magnitude faster than existing techniques, and remains usable when prior information is vague.

[1]  Ole Winther,et al.  Expectation Consistent Approximate Inference , 2005, J. Mach. Learn. Res..

[2]  Olivier François,et al.  Non-linear regression models for Approximate Bayesian Computation , 2008, Stat. Comput..

[3]  R. Plevin,et al.  Approximate Bayesian Computation in Evolution and Ecology , 2011 .

[4]  M. Feldman,et al.  Population growth of human Y chromosomes: a study of Y chromosome microsatellites. , 1999, Molecular biology and evolution.

[5]  Matthias Bethge,et al.  Bayesian Inference for Sparse Generalized Linear Models , 2007, ECML.

[6]  J. Rosenthal,et al.  Optimal scaling for various Metropolis-Hastings algorithms , 2001 .

[7]  Paul Marjoram,et al.  Markov chain Monte Carlo without likelihoods , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Nicolas Chopin,et al.  Expectation Propagation for Likelihood-Free Inference , 2011, 1107.5959.

[9]  D. Balding,et al.  Approximate Bayesian computation in population genetics. , 2002, Genetics.

[10]  Yuan Qi,et al.  Sparse-posterior Gaussian Processes for general likelihoods , 2010, UAI.

[11]  M. Beaumont Approximate Bayesian Computation in Evolution and Ecology , 2010 .

[12]  David Welch,et al.  Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems , 2009, Journal of The Royal Society Interface.

[13]  Jean-Michel Marin,et al.  Approximate Bayesian computational methods , 2011, Statistics and Computing.

[14]  Roman Holenstein,et al.  Particle Markov chain Monte Carlo , 2009 .

[15]  Ole Winther,et al.  Perturbation Corrections in Approximate Inference: Mixture Modelling Applications , 2009, J. Mach. Learn. Res..

[16]  M. Seeger Expectation Propagation for Exponential Families , 2005 .

[17]  D. Balding,et al.  Statistical Applications in Genetics and Molecular Biology On Optimal Selection of Summary Statistics for Approximate Bayesian Computation , 2011 .

[18]  M. Blum Approximate Bayesian Computation: A Nonparametric Perspective , 2009, 0904.0635.

[19]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[20]  Tom Minka,et al.  Expectation Propagation for approximate Bayesian inference , 2001, UAI.

[21]  Arnaud Doucet,et al.  An adaptive sequential Monte Carlo method for approximate Bayesian computation , 2011, Statistics and Computing.

[22]  C. Rasmussen,et al.  Approximations for Binary Gaussian Process Classification , 2008 .

[23]  A. Doucet,et al.  Particle Markov chain Monte Carlo methods , 2010 .

[24]  P. Moral,et al.  Sequential Monte Carlo samplers , 2002, cond-mat/0212648.

[25]  Darren J. Wilkinson,et al.  Bayesian inference for a discretely observed stochastic kinetic model , 2008, Stat. Comput..