RESEARCH
<< back
COSMOABC: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation

E. E. O. Ishida, S. D. P. Vitenti, M. Penna-Lima, J. Cisewski-kehe, R. S. de Souza, A. M. M. Trindade, E. Cameron, V. C. Busti, G. Östlin

Abstract
Approximate Bayesian Computation (ABC) enables parameter inference for complex physical systems in cases where the true likelihood function is unknown, unavailable, or computationally too expensive. It relies on the forward simulation of mock data and comparison between observed and synthetic catalogues. Here we present COSMOABC, a Python ABC sampler featuring a Population Monte Carlo variation of the original ABC algorithm, which uses an adaptive importance sampling scheme. The code is very flexible and can be easily coupled to an external simulator, while allowing to incorporate arbitrary distance and prior functions. As an example of practical application, we coupled COSMOABC with the NUMCOSMO library and demonstrate how it can be used to estimate posterior probability distributions over cosmological parameters based on measurements of galaxy clusters number counts without computing the likelihood function. COSMOABC is published under the GPLv3 license on PyPI and GitHub and documentation is available at http://goo.gl/SmB8EX.

Keywords
Galaxies: statistics - (cosmology:) large-scale structure of universe

Astronomy and Computing
Volume 13, Page 1
2015 November

>> ADS>> DOI

Faculdade de Ciências da Universidade de Lisboa Universidade do Porto Faculdade de Ciências e Tecnologia da Universidade de Coimbra
Fundação para a Ciência e a Tecnologia COMPETE 2020 PORTUGAL 2020 União Europeia