abctoolbox a versatile toolkit for approximate bayesian computationsabctoolbox近似贝叶斯计算的通用的工具包.pdfVIP

Wegmann et al. BMC Bioinformatics 2010, 11:116 /1471-2105/11/116 SOFTWARE Open Access ABCtoolbox: a versatile toolkit for approximate Bayesian computations 1,2* 3 4 1,2 Daniel Wegmann , Christoph Leuenberger , Samuel Neuenschwander , Laurent Excoffier Abstract Background: The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results: Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher le