Robust sampling for weak lensing and clustering analyses with the Dark Energy Survey

Abstract

Recent cosmological analyses rely on the ability to accurately sample from high-dimensional posterior distributions. A variety of algorithms have been applied in the field, but justification of the particular sampler choice and settings is often lacking. Here we investigate three such samplers to motivate and validate the algorithm and settings used for the Dark Energy Survey (DES) analyses of the first 3 years (Y3) of data from combined measurements of weak lensing and galaxy clustering. We employ the full DES Year 1 likelihood alongside a much faster approximate likelihood, which enables us to assess the outcomes from each sampler choice and demonstrate the robustness of our full results. We find that the ellipsoidal nested sampling algorithm $\text{MultiNest}$ reports inconsistent estimates of the Bayesian evidence and somewhat narrower parameter credible intervals than the sliced nested sampling implemented in $\text{PolyChord}$. We compare the findings from $\text{MultiNest}$ and $\text{PolyChord}$ with parameter inference from the Metropolis-Hastings algorithm, finding good agreement. We determine that $\text{PolyChord}$ provides a good balance of speed and robustness, and recommend different settings for testing purposes and final chains for analyses with DES Y3 data. Our methodology can readily be reproduced to obtain suitable sampler settings for future surveys.

arXiv comment16 pages, 11 figures

References in corpus (51)

• emcee: The MCMC Hammer
• Efficient Computation of CMB anisotropies in closed FRW models
• Cosmological parameters from SDSS and WMAP
• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• A Comprehensive Measurement of the Local Value of the Hubble Constant with 1 km/s/Mpc Uncertainty from the Hubble Space Telescope and the SH0ES Team
• Baryonic Features in the Matter Transfer Function
• dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• Dark Energy Survey Year 1 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing
• Dark Energy Survey Year 3 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing
• Revising the Halofit Model for the Nonlinear Matter Power Spectrum
• Importance Nested Sampling and the MultiNest Algorithm
• Cobaya: Code for Bayesian Analysis of hierarchical physical models
• Efficient sampling of fast and slow cosmological parameters
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Parameter Estimation Methodology
• Cosmology Intertwined II: The Hubble Constant Tension
• CMB power spectrum parameter degeneracies in the era of precision cosmology
• PolyChord: nested sampling for cosmology
• PolyChord: next-generation nested sampling
• Cosmology Intertwined III: $f{\sigma }_8$ and $S_8$
• A Nested Sampling Algorithm for Cosmological Model Selection
• Massive Neutrinos and the Non-linear Matter Power Spectrum
• Data analysis recipes: Using Markov Chain Monte Carlo
• CosmoSIS: modular cosmological parameter estimation
• Fast and reliable MCMC for cosmological parameter estimation
• CosmoLike - Cosmological Likelihood Analyses for Photometric Galaxy Surveys
• A buyer's guide to the Hubble Constant
• II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements
• FAST-PT: a novel algorithm to calculate convolution integrals in cosmological perturbation theory
• Efficient Bayesian inference for multimodal problems in cosmology
• Combining cosmological datasets: hyperparameters and Bayesian evidence
• A Bayesian model selection analysis of WMAP3
• The Age of the Universe and the Cosmological Constant Determined from Cosmic Microwave Background Anisotropy Measurements
• Quantifying tensions in cosmological parameters: Interpreting the DES evidence ratio
• Bayesian evidence as a tool for comparing datasets
• anesthetic: nested sampling visualisation
• Quantifying dimensionality: Bayesian cosmological model complexities
• zeus: A Python implementation of Ensemble Slice Sampling for efficient Bayesian parameter inference
• Comparison of sampling techniques for Bayesian parameter estimation
• Model selection as a science driver for dark energy surveys
• Ensemble Slice Sampling: Parallel, black-box and gradient-free inference for correlated & multimodal distributions
• Bayesian `Hyper-Parameters' Approach to Joint Estimation: The Hubble Constant from CMB Measurements
• nestcheck: diagnostic tests for nested sampling calculations
• Recursive Pathways to Marginal Likelihood Estimation with Prior-Sensitivity Analysis
• Constraints on dark matter to dark radiation conversion in the late universe with DES-Y1 and external data
• On statistical uncertainty in nested sampling
• Conservative cosmology: combining data with allowance for unknown systematics
• Model selection forecasts for the spectral index from the Planck satellite
• DES Y1 results: Splitting growth and geometry to test $\Lambda$CDM
• Dark energy constraints and correlations with systematics from CFHTLS weak lensing, SNLS supernovae Ia and WMAP5

Cited by in corpus (18)

• Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Power Spectra
• DES Y3 + KiDS-1000: Consistent cosmology combining cosmic shear surveys
• DES Y3 cosmic shear down to small scales: constraints on cosmology and baryons
• KiDS-1000: Cosmology with improved cosmic shear measurements
• A Unified Catalog-level Reanalysis of Stage-III Cosmic Shear Surveys
• Bayesian model comparison for simulation-based inference
• Vacuum Energy Density Measured from Cosmological Data
• Cluster cosmology with anisotropic boosts: Validation of a novel forward modeling analysis and application on SDSS redMaPPer clusters
• An empirical approach to model selection: weak lensing and intrinsic alignments
• Euclid preparation. LIII. LensMC, weak lensing cosmic shear measurement with forward modelling and Markov Chain Monte Carlo sampling
• Dark Energy Survey Year 3: Blue Shear
• Euclid and KiDS-1000: Quantifying the impact of source-lens clustering on cosmic shear analyses
• The Radius of PSR J0437-4715 from NICER Data
• CombineHarvesterFlow: Joint Probe Analysis Made Easy with Normalizing Flows
• Small scale clustering of BOSS galaxies: dependence on luminosity, color, age, stellar mass, specific star formation rate and other properties
• What We Have Learned from NICER About Neutron Star Radii
• Bayesian evidence and model selection approach for time-dependent dark energy
• Cosmological constraints on viable $f(R)$ models using weak lensing