emcee : The MCMC Hammer

作者：Daniel Foreman-Mackey, David W. Hogg, Dustin Lang, Jonathan Goodman · 发表于：Publications of the Astronomical Society of the Pacific · 年份：2013 · DOI：10.1086/670067 · 被引用次数：12279 · 研究领域：Markov Chains and Monte Carlo Methods、Statistical Mechanics and Entropy、Galaxies: Formation, Evolution, Phenomena

We introduce a stable, well tested Python implementation of the affine-invariant ensemble sampler for Markov chain Monte Carlo (MCMC) proposed by Goodman & Weare ( 2010 ). The code is open source and has already been used in several published projects in the astrophysics literature. The algorithm behind emcee has several advantages over traditional MCMC sampling methods and it has excellent performance as measured by the autocorrelation time (or function calls per independent sample). One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ∼ N 2 for a traditional algorithm in an N -dimensional parameter space. In this document, we describe the algorithm and the details of our implementation. Exploiting the parallelism of the ensemble method, emcee permits any user to take advantage of multiple CPU cores without extra effort. The code is available online at http://dan.iel.fm/emcee under the GNU General Public License v2.