Bayesian computation and model selection without likelihoods.
Until recently, the use of Bayesian inference was limited to a few cases because for many realistic probability models the likelihood function cannot be calculated analytically. The situation changed with the advent of likelihood-free inference algorithms, often subsumed under the term approximate Bayesian computation (ABC). A key innovation was the use of a postsampling regression adjustment, allowing larger tolerance values and as such shifting computation time to realistic orders of magnitude. Here we propose a reformulation of the regression adjustment in terms of a general linear model (GLM). This allows the integration into the sound theoretical framework of Bayesian statistics and the use of its methods, including model selection via Bayes factors. We then apply the proposed methodology to the question of population subdivision among western chimpanzees, Pan troglodytes verus.
Computational and Molecular Population Genetics Laboratory, Institute of Ecology and Evolution, University of Bern, 3012 Bern, Switzerland. christoph.leuenberger@unifr.ch
Christoph
C
Leuenberger
Daniel
D
Wegmann
D000465:Algorithms; D000818:Animals; D001499:Bayes Theorem; D005828:Genetics, Population; D006801:Humans; D016013:Likelihood Functions; D016014:Linear Models; D008957:Models, Genetic; D002679:Pan troglodytes / Q000235:genetics