The Rate of Convergence for Approximate Bayesian Computation

Approximate Bayesian Computation (ABC) is a popular computational method for likelihood-free Bayesian inference. The term "likelihood-free" refers to problems where the likelihood is intractable to compute or estimate directly, but where it is possible to generate simulated data $X$ relatively easily given a candidate set of parameters $\theta$ simulated from a prior distribution. Parameters which generate simulated data within some tolerance $\delta$ of the observed data $x^*$ are regarded as plausible, and a collection of such $\theta$ is used to estimate the posterior distribution $\theta\,|\,X\!=\!x^*$. Suitable choice of $\delta$ is vital for ABC methods to return good approximations to $\theta$ in reasonable computational time. While ABC methods are widely used in practice, particularly in population genetics, study of the mathematical properties of ABC estimators is still in its infancy. We prove that ABC estimates converge to the exact solution under very weak assumptions and, under slightly stronger assumptions, quantify the rate of this convergence. Our results can be used to guide the choice of the tolerance parameter $\delta$.

[1]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[2]  Jochen Voss,et al.  An Introduction to Statistical Computing: A Simulation-based Approach , 2013 .

[3]  S. Filippi,et al.  Optimizing threshold-schedules for sequential approximate Bayesian computation: applications to molecular systems , 2013, Statistical applications in genetics and molecular biology.

[4]  R. Wilkinson Approximate Bayesian computation (ABC) gives exact results under the assumption of model error , 2008, Statistical applications in genetics and molecular biology.

[5]  Arnaud Guyader,et al.  New insights into Approximate Bayesian Computation , 2012, 1207.6461.

[6]  Paul Fearnhead,et al.  Constructing summary statistics for approximate Bayesian computation: semi‐automatic approximate Bayesian computation , 2012 .

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

[8]  C C Drovandi,et al.  Estimation of Parameters for Macroparasite Population Evolution Using Approximate Bayesian Computation , 2011, Biometrics.

[9]  R. Cooke Real and Complex Analysis , 2011 .

[10]  S. Tavaré,et al.  Dating primate divergences through an integrated analysis of palaeontological and molecular data. , 2011, Systematic biology.

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

[12]  Laurent Excoffier,et al.  Bayesian inference of the demographic history of chimpanzees. , 2010, Molecular biology and evolution.

[13]  David M. Walker,et al.  Parameter inference in small world network disease models with approximate Bayesian Computational methods , 2009, Physica A: Statistical Mechanics and its Applications.

[14]  Viet Chi Tran,et al.  HIV with contact tracing: a case study in approximate Bayesian computation. , 2008, Biostatistics.

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

[16]  Simon Tavaré,et al.  Integrating Approximate Bayesian Computation with Complex Agent-Based Models for Cancer Research , 2010, COMPSTAT.

[17]  M. Beaumont,et al.  Inferring introduction routes of invasive species using approximate Bayesian computation on microsatellite data , 2010, Heredity.

[18]  J. Lopes,et al.  The use of approximate Bayesian computation in conservation genetics and its application in a case study on yellow-eyed penguins , 2010, Conservation Genetics.

[19]  Christophe Andrieu,et al.  Model criticism based on likelihood-free inference, with an application to protein network evolution , 2009, Proceedings of the National Academy of Sciences.

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

[21]  L. Excoffier,et al.  Statistical evaluation of alternative models of human evolution , 2007, Proceedings of the National Academy of Sciences.

[22]  S. Coles,et al.  Inference for Stereological Extremes , 2007 .

[23]  Mark M. Tanaka,et al.  Sequential Monte Carlo without likelihoods , 2007, Proceedings of the National Academy of Sciences.

[24]  Andrew R. Francis,et al.  Using Approximate Bayesian Computation to Estimate Tuberculosis Transmission Parameters From Genotype Data , 2006, Genetics.

[25]  Kevin R. Thornton,et al.  Approximate Bayesian Inference Reveals Evidence for a Recent, Severe Bottleneck in a Netherlands Population of Drosophila melanogaster , 2006, Genetics.

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

[27]  Amir Dembo,et al.  Large Deviations Techniques and Applications , 1998 .

[28]  P. Donnelly,et al.  Inferring coalescence times from DNA sequence data. , 1997, Genetics.

[29]  H. P. Annales de l'Institut Henri Poincaré , 1931, Nature.