Recognition: unknown
Efficient sampling of fast and slow cosmological parameters
read the original abstract
Physical parameters are often constrained from the data likelihoods using sampling methods. Changing some parameters can be much more computationally expensive (`slow') than changing other parameters (`fast parameters'). I describe a method for decorrelating fast and slow parameters so that parameter sampling in the full space becomes almost as efficient as sampling in the slow subspace when the covariance is well known and the distributions are simple. This gives a large reduction in computational cost when there are many fast parameters. The method can also be combined with a fast 'dragging' method proposed by Neal (2005) that can be more robust and efficient when parameters cannot be fully decorrelated a priori or have more complicated dependencies. I illustrate these methods for the case of cosmological parameter estimation using data likelihoods from the Planck satellite observations with dozens of fast nuisance parameters, and demonstrate a speed up by a factor of five or more. In more complicated cases, especially where the fast subspace is very fast but complex or highly correlated, the fast-slow sampling methods can in principle give arbitrarily large performance gains. The new samplers are implemented in the latest version of the publicly available CosmoMC code.
This paper has not been read by Pith yet.
Forward citations
Cited by 7 Pith papers
-
DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations
First-year DESI BAO data are consistent with flat LambdaCDM and, when combined with CMB, show a 2.5-3.9 sigma preference for evolving dark energy (w0 > -1, wa < 0) that strengthens with certain supernova datasets.
-
Cobaya: Code for Bayesian Analysis of hierarchical physical models
Cobaya is a modular Bayesian analysis code that exploits model interdependencies via automatic caching and a novel parameter-blocking algorithm to minimize sampling cost.
-
Into the Gompverse: A robust Gompertzian reionization model for CMB analyses
A Gompertzian reionization model with three nuisance parameters demotes optical depth to a derived quantity, reducing its uncertainty by a factor of three and revealing potential neutrino mass tension in CMB analyses.
-
SWIM: Stochastic Warm Inflation Module to generate and analyse Warm Inflationary power spectrum
SWIM is a publicly available numerical platform that generates full Warm Inflation scalar power spectra and enables MCMC parameter estimation with CMB data for arbitrary potentials and dissipation forms.
-
GetDist: a Python package for analysing Monte Carlo samples
GetDist implements boundary-corrected KDE with automatic smoothing for analyzing weighted and correlated Monte Carlo samples, plus plotting and diagnostic tools.
-
Planck 2018 results. VI. Cosmological parameters
Final Planck CMB data confirms the flat 6-parameter ΛCDM model with Ω_c h² = 0.120 ± 0.001, Ω_b h² = 0.0224 ± 0.0001, n_s = 0.965 ± 0.004, τ = 0.054 ± 0.007, H_0 = 67.4 ± 0.5 km/s/Mpc, and no strong evidence for extensions.
-
Extended Dark Energy analysis using DESI DR2 BAO measurements
Extended analysis of DESI DR2 data confirms robust evidence for dynamical dark energy with phantom crossing preference, stable under parametric and non-parametric modeling.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.