emcee: The MCMC Hammer
Pith reviewed 2026-05-13 21:55 UTC · model grok-4.3
The pith
The emcee package implements an affine-invariant ensemble sampler for MCMC that requires tuning only one or two parameters.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
emcee is a Python implementation of the Goodman and Weare affine-invariant ensemble sampler for MCMC. It achieves excellent performance as measured by autocorrelation time while requiring hand-tuning of only one or two parameters instead of approximately N squared for an N-dimensional space. The code exploits the parallelism of the ensemble method to allow easy use of multiple processors.
What carries the argument
The affine-invariant ensemble sampler, which proposes moves based on the current positions of other walkers in the ensemble to maintain invariance under affine transformations.
If this is right
- Users can apply MCMC methods to complex models in astrophysics with minimal parameter tuning.
- The sampler can be run in parallel without additional coding effort.
- Published projects have already used it, demonstrating its stability.
Where Pith is reading between the lines
- Similar ensemble methods could be adapted for other sampling algorithms beyond MCMC.
- Integration with modern hardware like GPUs might further improve performance.
Load-bearing premise
That the affine-invariant ensemble method converges reliably to the true posterior distribution for arbitrary likelihood functions.
What would settle it
Demonstrating a likelihood function where the emcee sampler fails to converge or has very long autocorrelation times despite using the recommended settings.
read the original abstract
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 $\sim 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 and API. 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 MIT License.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents emcee, a stable and well-tested open-source Python implementation of the affine-invariant ensemble MCMC sampler introduced by Goodman & Weare (2010). It describes the algorithm, implementation details including parallelism via multiprocessing, the user API, and performance measured by autocorrelation time. The central claims are that the method requires hand-tuning of only 1-2 parameters (versus ~N² for traditional Metropolis-Hastings in N dimensions) and delivers efficient sampling, with the code already used in multiple published astrophysics projects.
Significance. If the reported stability, performance, and ease of use hold, this work supplies a practical, accessible tool that lowers the barrier for high-dimensional Bayesian inference in astrophysics. The open MIT-licensed code, documented usage in the literature, and exploitation of ensemble parallelism for multi-core performance constitute concrete strengths that directly benefit the community.
minor comments (1)
- Abstract: the URL for code availability should be supplemented with a permanent identifier (e.g., Zenodo DOI) to ensure long-term accessibility.
Simulated Author's Rebuttal
We thank the referee for their positive review of the manuscript and their recommendation to accept. No major comments were raised in the report.
Circularity Check
No significant circularity; software implementation of prior algorithm
full rationale
This is a software description paper that implements the affine-invariant ensemble sampler proposed in the independent prior reference Goodman & Weare (2010). The central claims about reduced hand-tuning (1-2 parameters vs ~N²) and performance are directly inherited from that cited work without new derivation, fitting, or self-referential steps in the present manuscript. Implementation details, API, parallelism, and autocorrelation reporting are practical contributions backed by the open MIT-licensed code. No equations or results reduce by construction to the paper's own inputs; the derivation chain is external and the paper is self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The affine-invariant ensemble sampler proposed by Goodman & Weare (2010) produces correct samples from the target distribution when properly implemented.
Lean theorems connected to this paper
-
IndisputableMonolith.Foundation.DAlembert.Inevitabilitybilinear_family_forced unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
One major advantage of the algorithm is that it requires hand-tuning of only 1 or 2 parameters compared to ∼N² for a traditional algorithm in an N-dimensional parameter space.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 60 Pith papers
-
Dispersion Measure Distribution of Unlocalized Fast Radio Bursts as a Probe of the Hubble Constant
The DM distribution of unlocalized FRBs yields H0 = 73.8 +14.0/-12.3 km/s/Mpc with 18% uncertainty.
-
New Way to Date Globular Clusters: Brown Dwarf Cooling Sequences
A new histogram-free likelihood method applied to simulated JWST observations of brown dwarfs shows that globular cluster ages can be determined with formal errors under 0.2 Gyr.
-
Prospect of Measuring the Cosmic Dipole by Strongly Lensed Gravitational Waves Associated with Galaxy Surveys
Simulations forecast that 10 years of Einstein Telescope and Cosmic Explorer data could detect the cosmic dipole magnitude using strongly lensed GW events, with tighter bounds from combining double, triple, and quadru...
-
VLTI/PIONIER imaging of post-AGB binaries. An INSPIRING hunt for inner rim substructures in circumbinary discs
High-resolution interferometric imaging of eight post-AGB circumbinary discs reveals diverse inner-rim substructures including azimuthal brightness enhancements and arc-like features not explained by inclination alone.
-
Hidden Monsters with SPHEREx I: A goldmine for heavily reddened quasars at cosmic noon
SPHEREx data confirm 77 new luminous heavily reddened quasars at 1.5<z<3.9 that are hot-dust poor relative to unobscured quasars, supporting a blow-out feedback phase.
-
Precision hyperfine spectroscopy of an individual nuclear-spin-9/2
An Er3+ spin sensor in CaWO4 measures the hyperfine NMR spectrum of a single 93Nb nuclear spin-9/2 with Hertz resolution, determining its site, position, quadrupolar tensor, and two new terms in the spin Hamiltonian.
-
Measuring neutrino mass and asymmetry through galaxy pairwise peculiar velocity
Galaxy pairwise peculiar velocities from Cosmicflows-4 yield M_ν = 0.24^{+0.34}_{-0.18} eV and η² = 2.14^{+0.30}_{-0.32} (7σ non-zero asymmetry) in the CMB framework, consistent with prior Planck results.
-
Cosmologically viable non-polynomial quasi-topological gravity: explicit models, $\Lambda$CDM limit and observational constraints
Non-polynomial quasi-topological gravity models reproduce the standard thermal history, generate dynamical dark energy of geometric origin, and fit supernova, cosmic chronometer, and BAO data competitively with ΛCDM.
-
Measurement of the galaxy-velocity power spectrum of DESI tracers with the kinematic Sunyaev-Zeldovich effect using DESI DR2 and ACT DR6
DESI DR2 and ACT DR6 data yield 17σ LRG-velocity, 8.3σ ELG-velocity, and 6.8σ QSO-velocity detections plus a 3.1σ velocity-velocity signal, producing f_NL^loc = 15.9_{-34.4}^{+34.6} from the velocity field.
-
Correlation between nuclear isospin asymmetry and $\alpha$-particle preformation probability for superheavy nuclei from a Bayesian inference
Bayesian MCMC fitting of a phenomenological model reveals that isospin asymmetry strongly suppresses alpha preformation probability in superheavy nuclei while reproducing shell effects and experimental half-lives.
-
Accelerated Time-domain Analysis for Gravitational Wave Astronomy
Presents a practical fully time-domain end-to-end likelihood for gravitational-wave inference with structured linear algebra and GPU acceleration.
-
Cosmological analysis of the DESI DR1 Lyman alpha 1D power spectrum
DESI DR1 Lyman-alpha data yields Δ²★=0.379±0.032 and n★=-2.309±0.019 at k★=0.009 km⁻¹s and z=3, sharpening N_eff, α_s, and β_s constraints by factors of 1.18-1.90 when combined with other probes.
-
The velocity coherence scale: a novel probe of cosmic homogeneity and a potential standard ruler
The velocity coherence scale R_v marks the onset of statistical homogeneity, is redshift-independent in comoving coordinates, and connects directly to the matter-radiation equality scale k_eq in standard cosmology.
-
High CO/H2 ratios supports an exocometary origin for a CO-rich debris disk
First direct H2 measurements in two CO-rich exocometary belts yield CO/H2 lower limits of >1.35e-3 and >3.09e-5, showing the gas is H2-poor and supporting secondary exocometary origin.
-
First astrometric constraints on parity-violation in the gravitational wave background
First astrometric constraints on parity-violating SGWB amplitude are reported as h70²ΩV = -0.020 ± 0.025 (Gaia) and -0.004 ± 0.010 (VLBA) at 2σ, consistent with zero, over 4.2e-18 Hz to 1.1e-8 Hz.
-
Tomographic Alcock-Paczynski Test with Marked Correlation Functions
First integration of tomographic AP tests with MCFs and PCA compression yields 48% and 45% tighter errors on Ω_m and w versus standard two-point functions.
-
DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars
DESI measures BAO scales in six redshift bins with 0.52% combined precision using 5.7 million objects, detecting the signal at up to 9.1 sigma and finding larger scales than Planck LCDM at z<0.8.
-
JWST Observations of Asteroid 2024 YR4 Rule Out a 2032 Lunar Impact and Demonstrate a New Regime for Planetary Defense Follow-up
JWST/NIRCam observations of 2024 YR4 extend the orbital arc by eight months, reduce 2032 lunar encounter uncertainty by >30x, and rule out impact with a 22,900 ± 800 km miss distance.
-
Latitudinal variations in Neptune's temperature profile observed with ALMA
ALMA CO(3-2) observations with MCMC retrievals show latitudinal temperature trends on Neptune, including a cold southern polar layer at 300-600 mbar, consistent with Voyager at 10-100 mbar but new at deeper levels.
-
Milky Way Dynamics Favor Dark Matter over Modified Gravity Models
Milky Way radial and vertical dynamics are inconsistent with MOND and STVG but consistent with dark matter halos, disfavoring the former at high significance.
-
A single power law for the TRAPPIST-1 flare distribution across four orders of magnitude in energy
TRAPPIST-1 flares follow a single power law N(≥E_TESS) ∝ E_TESS^{-0.753} from 10^{29} to 10^{33} erg after sensitivity corrections and bandpass conversion.
-
Confronting Color Glass Condensate at next-to-leading order with HERA data
A Bayesian global fit at full NLO+NLL accuracy extracts the posterior distribution for the non-perturbative initial condition of the NLO Balitsky-Kovchegov equation from HERA inclusive and charm data.
-
If at First You Don't Succeed, Trispectrum: I. Estimating the Matter Power Spectrum Covariance with Higher-Order Statistics
Estimators from squeezed bispectrum and collapsed trispectrum recover unbiased small-scale matter power spectrum covariance at the percent level using 25 Quijote simulations.
-
AI-assisted modeling and Bayesian inference of unpolarized quark transverse momentum distributions from Drell-Yan data
An AI-assisted Bayesian framework extracts TMD PDFs from global Drell-Yan data using surrogate models for scalable MCMC sampling.
-
Probing Collapsed Dark Matter Halos with Fast Radio Bursts
Core-collapsed SIDM halos produce longer FRB image time delays than CDM halos, enabling future surveys to constrain self-interaction cross sections above roughly 18-40 cm²/g depending on collapse timing.
-
Posterior Predictive Checks for Gravitational-wave Populations: Limitations and Improvements
Maximum-likelihood-based posterior predictive checks detect model misspecification better than event-level versions for uncertain spin tilts, but current detector sensitivity limits their power; the Gaussian Component...
-
$\boldsymbol{B_c}$ Meson Spectroscopy from Bayesian MCMC: Probing Confinement and State Mixing
Bayesian MCMC sampling of Cornell and log-modified Cornell potentials reproduces known B_c states and supplies mass predictions for higher excitations with propagated uncertainties.
-
Model-independent test of the cosmic distance duality relation with recent observational data
Two model-independent methods applied to latest SN and BAO data find the cosmic distance duality relation consistent with observations within 1 sigma and no evidence of violation.
-
Simultaneous Estimation of Seabed and Its Roughness With Longitudinal Waves
An infinite-dimensional Bayesian framework estimates seabed topography and roughness simultaneously from acoustic data by assuming statistical isotropy and using fractional differentiability.
-
Impact of projection-induced optical selection bias on the weak lensing mass calibration of galaxy clusters
Projection-induced selection bias causes 20-50% overestimation of weak lensing masses for optically selected galaxy clusters, larger on scales >3 Mpc.
-
The universal growth of magnetic energy during the nonlinear phase of subsonic and supersonic small-scale dynamos
Ensemble simulations find linear magnetic energy growth (p_nl=1) in subsonic small-scale dynamos and quadratic growth (p_nl=2) in supersonic ones, with universal conversion efficiency ~1/100 of turbulent kinetic energ...
-
Bring the noise: exact inference from noisy simulations in collider physics
Introduces pseudo-marginal MCMC with unbiased Poisson likelihood estimator for exact inference despite noisy collider Monte Carlo simulations.
-
Reconstructing the Thermal Sunyaev Zeldovich Power Spectrum from Planck using the ABS Method
The ABS method applied to Planck PR3 data yields a tSZ power spectrum amplitude 34% lower than the Planck 2015 best-fit when trispectrum is included.
-
Modeling and measuring the anisotropic halo 3-point correlation function: a coordinated study
Novel implementation of anisotropic 3PCF model and estimator tested on 298 halo catalogs shows degeneracy breaking between f and b1 in 3PCF-only analysis but limited added value in joint 2PCF+3PCF due to tree-level mo...
-
Gapless neutron superfluidity can explain the late time cooling of transiently accreting neutron stars
Gapless neutron superfluidity from vortex pinning explains late-time cooling of KS 1731-260 and MXB 1659-29 without requiring suppression of superfluidity.
-
A closer look at the electromagnetic signatures of Bethe-Heitler pair production process in blazars
An empirical kernel for Bethe-Heitler pair spectra is introduced and applied to argue that synchrotron emission from these pairs can explain gamma-ray emission in low- and intermediate-peaked blazars for jet regions o...
-
Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing
DES Y3 3x2pt analysis constrains S8=0.776±0.017 and Ωm=0.339±0.032 in flat ΛCDM, consistent with Planck CMB results at p=0.13-0.48.
-
Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method
Bayesian data augmentation reintroduces missing LISA data segments as auxiliary variables during posterior sampling to enable consistent parameter estimation for galactic binaries despite gaps.
-
Cosmological constraints from neighbor-density-weighted marked correlation functions
Neighbor-density-weighted marked correlation functions improve FoM for Ωm–σ8 by 1.7–2.5× over standard 2PCF using Gaussian-process emulators on 129 w0waCDM+∑mν simulations.
-
Constraints on Schwarzschild Black Hole in a Generalized Dehnen-Type $(1,4,\gamma)$ Dark Matter Halo via the S2 Star Orbit around Sgr A$^\star$
Derives perihelion shift equations for S2 star in generalized Schwarzschild-Dehnen BH-DM spacetime and constrains gamma, rho_s, rs via MCMC on Do et al. and Gillessen et al. datasets.
-
JWST Advanced Deep Extragalactic Survey (JADES) Data Release 5: stellar population catalogue for galaxies in GOODS-N and GOODS-S
JADES DR5 delivers a public catalog of Bayesian-inferred stellar masses, SFRs, SFHs, dust, metallicities, and AGN contributions for ~500k galaxies via Prospector with an evolving SFMS prior.
-
Euclid preparation: Testing multi-field inflation with galaxy power spectrum and bispectrum
Validates redshift-space power spectrum and bispectrum analysis on Abacus-PNG mocks to recover unbiased f_NL constraints for Euclid spectroscopic sample.
-
Detection of a dark matter subhalo in the strongly lensed system PJ011646
A subhalo of M200 = 2.78e10 solar masses and concentration 30 is detected at 5.8 sigma in PJ011646 via ALMA imaging and grid-based NFW search after fitting an elliptical power-law plus multipole macromodel.
-
Connecting CGM enrichment with Lyman alpha emitters at 2.9 < z < 6.7
Blind MUSE search detects 156 LAEs at 2.9<z<6.7 with 34 C IV and 14 Mg II absorber associations, indicating low-mass galaxies enrich the CGM.
-
Mare versus highland lunar impact flash light curve dichotomy
Lunar impact flashes on highlands exhibit shallower, longer-lasting decay than on mare terrain, driven by fine ejecta in dual-size cooling models and indicating lithology-dependent early cratering.
-
Mare versus highland lunar impact flash light curve dichotomy
Lunar impact flashes on highlands exhibit longer and shallower decay than on maria, with the dual-size ejecta model attributing the difference mainly to fine droplets and indicating lithology dependence in early cratering.
-
Bulgeless Evolution And the Rise of Discs (BEARD) III. A numerical simulation view of satellites around Milky-Way analogues
Simulation comparison finds bulgeless galaxies host more centrally concentrated, disc-aligned satellites with steeper faint-end luminosity functions than bulge-dominated controls, reflecting co-evolution and quieter m...
-
Modeling the probability distribution for cosmological analysis with photometrically classified samples
A simplified likelihood treating photometric contamination as a redshift-dependent Gaussian mean shift is strongly favored by Bayes factor over the BEAMS two-component model and improves cosmological constraints on th...
-
On the quenching of LRD X-ray emission by both Compton-thick gas and high accretion rates
LRDs require Compton-thick gas at moderate metallicity plus high accretion rates producing weak X-rays to explain their non-detection, implying they are not chemically pristine.
-
A single power law for the TRAPPIST-1 flare distribution across four orders of magnitude in energy
TRAPPIST-1 flares obey a single power law N(≥E) ∝ E^{-0.753} from 10^29 to 10^33 erg in TESS energy after sensitivity corrections.
-
Perturbation Dynamics and Structure Formation in Extended Proca-Nuevo Gravity
Extended Proca-Nuevo gravity modifies the background expansion via a vector field algebraic constraint but leaves the matter growth equation identical to general relativity.
-
The January 2010 flare of Mrk421: Insights from a stochastic acceleration model
The January 2010 flare of Mrk421 shows spectral variability consistent with stochastic acceleration, including development of a Maxwellian pile-up component in the particle distribution on three nights.
-
Dynamical dark energy from Kretschmann scalar at low redshifts
Replacing the cosmological constant with the Kretschmann scalar yields a dynamical dark energy model that fits supernova and cosmic chronometer data and produces a phantom-crossing equation-of-state parameter w(z) sim...
-
Directional Tests of the Cosmic Distance Duality Relation using Pantheon+ and BAO
No evidence for directional anisotropy in the cosmic distance duality relation is found, yielding a robust 95% upper bound of 0.025 on the dipole amplitude after accounting for survey selection effects.
-
Probing soft signals of gravitational-wave memory with space-based interferometers
Space-based detectors can measure soft displacement-memory signals from gravitational waves at SNR greater than or equal to 10.
-
Testing the Spacetime Geometry of Sgr A* with the Relativistic Orbit of S2 star
Multiple black-hole and naked-singularity spacetimes are tested against S2 star astrometric and spectroscopic data, showing statistical degeneracy among several models under current observations.
-
Revisiting the Reported Period of FRB 20201124A Using MCMC Methods
An MCMC-based period search method recovers previously reported candidate periods in FRB 20201124A observations.
-
Probing the sensitivity of dark energy dynamics to equation of state parametrization flexibility
Flexible dark energy equation of state parametrizations mildly favor dynamical phantom-like behavior at z~1-2, but the preference is only ~2 sigma and sensitive to the model form and extrapolation.
-
GW250114: testing Hawking's area law and the Kerr nature of black holes
GW250114 data confirm the remnant black hole ringdown frequencies lie within 30% of Kerr predictions and that the final horizon area is larger than the sum of the progenitors' areas to high credibility.
-
Combined tracer analysis for DESI 2024 BAO
Combining LRG and ELG tracers with bias weighting improves BAO constraints by 11% on alpha_iso and 7% on alpha_AP in DESI DR1 data for the 0.8<z<1.1 bin.
Reference graph
Works this paper leans on
-
[1]
Akeret, J., Seehars, S., Amara, A., Refregier, A., & Csillaghy, A. 201 2, arXiv:1212.1721
- [2]
- [3]
- [4]
-
[5]
Brammer, G. B., S´ anchez-Janssen, R., Labb´ e, I., et al. 2012, ApJ, 758, L17 Brewer B. J., P´ artay L. B., Cs´ anyi G., 2011, Statistics and Computing, 21, 4, 649-656, arXiv:0912.2380
work page Pith review arXiv 2012
-
[6]
J., Foreman-Mackey, D., & Hogg, D
Brewer, B. J., Foreman-Mackey, D., & Hogg, D. W. 2012, arXiv:1211.5805
-
[7]
Brown, J. M., Rosenfeld, K. A., Andrews, S. M., Wilner, D. J., & van Dish oeck, E. F. 2012, ApJ, 758, L30
work page 2012
- [8]
-
[9]
Christen, J., A general purpose scale-independent MCMC algorithm , technical report I-07-16, CIMAT, Guanajuato, 2007
work page 2007
-
[10]
Cieza, L. A., Olofsson, J., Harvey, P. M., et al. 2013, ApJ, 762, 100
work page 2013
-
[11]
Cowles, M. K. & Carlin, B. P., 1996, Journal of the American Statistic al Association, 91, 883 Crossfield, I. J. M., Barman, T., Hansen, B. M. S., Tanaka, I., & Koda ma, T. 2012, ApJ, 760, 140
work page 1996
-
[12]
E., Guhathakurta, P., Fardal, M
Dorman, C. E., Guhathakurta, P., Fardal, M. A., et al. 2012, ApJ, 75 2, 147
work page 2012
-
[13]
Goodman, J. & Weare, J., 2010, Comm. App. Math. Comp. Sci., 5, 65
work page 2010
-
[14]
C., Bayesian Logical Data Analysis for the Physical Sciences , Cambridge University Press, 2005
Gregory, P. C., Bayesian Logical Data Analysis for the Physical Sciences , Cambridge University Press, 2005
work page 2005
-
[15]
Data analysis recipes: Fitting a model to data
Hogg, D. W., Bovy, J., & Lang, D., 2010, arXiv:1008.4686 [astro-ph.IM]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[16]
Hou, F., Goodman, J., Hogg, D. W., Weare, J., & Schwab, C. 2012, ApJ , 745, 198
work page 2012
-
[17]
Quasi-Periodic Oscillations and broadband variability in short magnetar bursts
Huppenkothen, D., Watts, A. L., Uttley, P., et al. 2012, arXiv:1212.1011
work page Pith review arXiv 2012
-
[18]
Lang, D., & Hogg, D. W. 2012, AJ, 144, 46 – 14 –
work page 2012
-
[19]
MacKay, D., Information Theory, Inference, and Learning Algorithms , Cambridge University Press, 2003
work page 2003
-
[20]
Monnier, J. D., Che, X., Zhao, M., et al. 2012, ApJ, 761, L3
work page 2012
-
[21]
Morton, T. D. 2012, ApJ, 761, 6
work page 2012
-
[22]
Olofsson, J., Juh´ asz, A., Henning, T., et al. 2012, A&A, 542, A90
work page 2012
-
[23]
& Gelman, A., 2010, Statistica Sinica, 20, 343–364
Pasarica, C. & Gelman, A., 2010, Statistica Sinica, 20, 343–364
work page 2010
-
[24]
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P., Numerical Recipes: The Art of Scientific Computing , Cambridge University Press, 2007
work page 2007
-
[25]
Reis, R. C., Miller, J. M., Reynolds, M. T., et al. 2013, ApJ, 763, 48 Roˇ skar, R., Debattista, V. P., & Loebman, S. R. 2012, arXiv:1211.1982
-
[26]
Sanders, J. S., & Fabian, A. C. 2013, MNRAS, 453 van Haasteren, R., Mingarelli, C. M. F., Vecchio, A., & Lassus, A. 2013 , arXiv:1301.6673
-
[27]
Weisz, D. R., Fouesneau, M., Hogg, D. W., et al. 2013, ApJ, 762, 123
work page 2013
-
[28]
Widrow, L. M. and Pym, B. and Dubinski, J., 2008, ApJ, 679, 1239This preprint was prepared with the AAS L ATEX macros v5.2. – 15 – A. Installation The easiest way to install emcee is using pip4. Running the command % pip i n s t a l l emcee at the command line of a UNIX-based system will install the package in your Python path. If you would like to instal...
work page 2008
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.