Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy

Abstract

Markov Chain Monte Carlo based Bayesian data analysis has now become the method of choice for analyzing and interpreting data in almost all disciplines of science. In astronomy, over the last decade, we have also seen a steady increase in the number of papers that employ Monte Carlo based Bayesian analysis. New, efficient Monte Carlo based methods are continuously being developed and explored. In this review, we first explain the basics of Bayesian theory and discuss how to set up data analysis problems within this framework. Next, we provide an overview of various Monte Carlo based methods for performing Bayesian data analysis. Finally, we discuss advanced ideas that enable us to tackle complex problems and thus hold great promise for the future. We also distribute downloadable computer software (available at https://github.com/sanjibs/bmcmc/ ) that implements some of the algorithms and examples discussed here.

49 pages, draft version, to appear in Annual Review of Astronomy and Astrophysics

References in corpus (20)

• Bayes in the sky: Bayesian inference and model selection in cosmology
• The mass distribution and gravitational potential of the Milky Way
• Trigonometric Parallaxes of Massive Star Forming Regions: VI. Galactic Structure, Fundamental Parameters and Non-Circular Motions
• The pseudo-marginal approach for efficient Monte Carlo computations
• The Milky Way's circular velocity curve between 4 and 14 kpc from APOGEE data
• On the Shoulders of Giants: Properties of the Stellar Halo and the Milky Way Mass Distribution
• Inferring the eccentricity distribution
• Adaptive Importance Sampling in General Mixture Classes
• Constraining the Galaxy's dark halo with RAVE stars
• Kinematics of the stellar halo and the mass distribution of the Milky Way using BHB stars
• On Correlated-noise Analyses Applied To Exoplanet Light Curves
• Kinematic modelling of the Milky Way using the RAVE and GCS stellar surveys
• A Bayesian model selection analysis of WMAP3
• Extended distribution functions for our Galaxy
• Estimation of cosmological parameters using adaptive importance sampling
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• From EM to Data Augmentation: The Emergence of MCMC Bayesian Computation in the 1980s
• Bayesian Model Averaging in Astrophysics: A Review
• Marginal likelihoods of distances and extinctions to stars: computation and compact representation

Cited by in corpus (6)

• The TESS-HERMES survey Data Release 1: high-resolution spectroscopy of the TESS southern continuous viewing zone
• The Need for Speed: Escape velocity and dynamical mass measurements of the Andromeda galaxy
• Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method
• Bayesian Coronal Seismology
• A Spatially Resolved Radio Spectral Index Study of the Dwarf Irregular Galaxy NGC\,1569
• Double-peaked OI profile: a likely signature of the gaseous ring around KH 15D