A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data

Abstract

We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes.

Submitted to ApJ

References in corpus (14)

• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Data Processing, Sky Maps, and Basic Results
• Three Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Polarization Analysis
• Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Temperature Analysis
• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Galactic Foreground Emission
• Joint Bayesian component separation and CMB power spectrum estimation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Galactic Foreground Emission
• Bayesian analysis of sparse anisotropic universe models and application to the 5-yr WMAP data
• A re-analysis of the three-year WMAP temperature power spectrum and likelihood
• Estimation of Polarized Power Spectra by Gibbs sampling
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• The joint large-scale foreground-CMB posteriors of the 3-year WMAP data
• Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps
• Marginal distributions for cosmic variance limited CMB polarization data

Cited by in corpus (28)

• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Power Spectra and WMAP-Derived Parameters
• Bayesian power-spectrum inference for Large Scale Structure data
• Polarized CMB power spectrum estimation using the pure pseudo-cross-spectrum approach
• Planck Intermediate Results. IX. Detection of the Galactic haze with Planck
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• Bayesian delensing delight: sampling-based inference of the primordial CMB and gravitational lensing
• Methods for Bayesian power spectrum inference with galaxy surveys
• Bayesian component separation and CMB estimation for the 5-year WMAP temperature data
• Matrix-free Large Scale Bayesian inference in cosmology
• The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field
• Methods for pixel domain correction of EB leakage
• BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
• Bayesian power-spectrum inference with foreground and target contamination treatment
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• Bayesian angular power spectrum analysis of interferometric data
• Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze
• Marginal distributions for cosmic variance limited CMB polarization data
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• Almanac: MCMC-based signal extraction of power spectra and maps on the sphere
• A new method for obtaining the star formation law in galaxies
• General solutions of the leakage in integral transforms and applications to the EB-leakage and detection of the cosmological gravitational wave background
• A CMB Gibbs sampler for localized secondary anisotropies
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations
• Pixel domain multi-resolution minimum variance painting of CMB maps
• Accurate and efficient likelihood modeling for large-scale CMB data