The stepping-stone sampling algorithm for calculating the evidence of gravitational wave models

Abstract

Bayesian statistical inference has become increasingly important for the analysis of observations from the Advanced LIGO and Advanced Virgo gravitational-wave detectors. To this end, iterative simulation techniques, in particular nested sampling and parallel tempering, have been implemented in the software library LALInference to sample from the posterior distribution of waveform parameters of compact binary coalescence events. Nested sampling was mainly developed to calculate the marginal likelihood of a model but can produce posterior samples as a by-product. Thermodynamic integration is employed to calculate the evidence using samples generated by parallel tempering but has been found to be computationally demanding. Here we propose the stepping-stone sampling algorithm, originally proposed by Xie et al. (2011) in phylogenetics and a special case of path sampling, as an alternative to thermodynamic integration. The stepping-stone sampling algorithm is also based on samples from the power posteriors of parallel tempering but has superior performance as fewer temperature steps and thus computational resources are needed to achieve the same accuracy. We demonstrate its performance and computational costs in comparison to thermodynamic integration and nested sampling in a simulation study and a case study of computing the marginal likelihood of a binary black hole signal model applied to simulated data from the Advanced LIGO and Advanced Virgo gravitational wave detectors. To deal with the inadequate methods currently employed to estimate the standard errors of evidence estimates based on power posterior techniques, we propose a novel block bootstrap approach and show its potential in our simulation study and LIGO application.

10 pages, 5 figures, 2 tables

References in corpus (15)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Advanced Virgo: a 2nd generation interferometric gravitational wave detector
• Multi-messenger Observations of a Binary Neutron Star Merger
• GW170817: Measurements of Neutron Star Radii and Equation of State
• GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence
• An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi-GBM Detection of GRB 170817A
• Properties of the binary neutron star merger GW170817
• Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library
• BayesWave: Bayesian Inference for Gravitational Wave Bursts and Instrument Glitches
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• A Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background
• The optimal search for an astrophysical gravitational-wave background
• Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data
• Constraining the p-mode--g-mode tidal instability with GW170817
• Bayesian support for Evolution: detecting phylogenetic signal in a subset of the primate family

Cited by in corpus (11)

• Gravitational Wave Denoising of Binary Black Hole Mergers with Deep Learning
• Parameter estimation with gravitational waves
• Bilby-MCMC: An MCMC sampler for gravitational-wave inference
• GPU-accelerated massive black hole binary parameter estimation with LISA
• Eryn : A multi-purpose sampler for Bayesian inference
• Variational Inference as an alternative to MCMC for parameter estimation and model selection
• Computational Techniques for Parameter Estimation of Gravitational Wave Signals
• Exploring the Potential for Detecting Rotational Instabilities in Binary Neutron Star Merger Remnants with Gravitational Wave Detectors
• Closing the Evidence Gap: reddemcee, a Fast Adaptive Parallel Tempering Sampler
• The NANOGrav 15 yr Data Set: Customized Chromatic Noise Models
• A fully-automated end-to-end pipeline for massive black hole binary signal extraction from LISA data