abc: an R package for approximate Bayesian computation (ABC)

Katalin Csilléry , Olivier François (1) , Michael G B Blum (1)

1 TIMC-IMAG-BCM - Biologie Computationnelle et Mathématique

1. Many recent statistical applications involve inference under complex models, where it is computationally prohibitive to calculate likelihoods but possible to simulate data. Approximate Bayesian computation (ABC) is devoted to these complex models because it bypasses the evaluation of the likelihood function by comparing observed and simulated data. 2. We introduce the R package ‘abc’ that implements several ABC algorithms for performing parameter estimation and model selection. In particular, the recently developed nonlinear heteroscedastic regression methods for ABC are implemented. The ‘abc’ package also includes a cross-validation tool for measuring the accuracy of ABC estimates and to calculate the misclassification probabilities when performing model selection. The main functions are accompanied by appropriate summary and plotting tools. 3. R is already widely used in bioinformatics and several fields of biology. The R package ‘abc’ will make the ABC algorithms available to a large number of R users. ‘abc’ is a freely available R package under the GPL license, and it can be downloaded at http://cran.r-project.org/web/packages/ abc/index.html.

Domaines: Génétique des populations [q-bio.PE]

Katalin Csilléry, Olivier François, Michael G B Blum. abc: an R package for approximate Bayesian computation (ABC). Methods in Ecology and Evolution, 2012, 3 (3), pp.475-479. ⟨10.1111/j.2041-210X.2011.00179.x⟩. ⟨hal-00861843⟩