Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling

Abstract

We describe an efficient and exact method that enables global Bayesian analysis of cosmic microwave background (CMB) data. The method reveals the joint posterior density (or likelihood for flat priors) of the power spectrum $C_{\ell }$ and the CMB signal. Foregrounds and instrumental parameters can be simultaneously inferred from the data. The method allows the specification of a wide range of foreground priors. We explicitly show how to propagate the non-Gaussian dependency structure of the $C_{\ell }$ posterior through to the posterior density of the parameters. If desired, the analysis can be coupled to theoretical (cosmological) priors and can yield the posterior density of cosmological parameter estimates directly from the time-ordered data. The method does not hinge on special assumptions about the survey geometry or noise properties, etc. It is based on a Monte Carlo approach and hence parallelizes trivially. No trace or determinant evaluations are necessary. The feasibility of this approach rests on the ability to solve the systems of linear equations which arise. These are of the same size and computational complexity as the map-making equations. We describe a pre-conditioned conjugate gradient technique that solves this problem and demonstrate in a numerical example that the computational time required for each Monte Carlo sample scales as $n_p^{3/2}$ with the number of pixels $n_p$. We test our method using the COBE-DMR data and explore the non-Gaussian joint posterior density of the COBE-DMR $C_{\ell }$ in several projections.

arXiv comment12 pages, 6 EPS figures, LaTeX. Supporting material and high-res figures will be made available at http://www.astro.uiuc.edu/projects/wandelt_research/GEM/

References in corpus (30)

• MAXIMA-1: A Measurement of the Cosmic Microwave Background Anisotropy on angular scales of 10 arcminutes to 5 degrees
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Foreground Emission
• DASI First Results: A Measurement of the Cosmic Microwave Background Angular Power Spectrum
• Cosmic Microwave Background Anisotropies
• MASTER of the CMB Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex CMB Data Sets
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Parameter Estimation Methodology
• Determination of Inflationary Observables by Cosmic Microwave Background Anisotropy Experiments
• Estimating the Power Spectrum of the Cosmic Microwave Background
• A Measurement of the Angular Power Spectrum of the CMB from l = 100 to 400
• How to measure CMB power spectra without losing information
• High Resolution Observations of the CMB Power Spectrum with ACBAR
• Radical Compression of Cosmic Microwave Background Data
• Wiener Reconstruction of The Large Scale Structure
• Power Spectrum of Primordial Inhomogeneity Determined from the 4-Year COBE DMR Sky Maps
• Myths and Truths Concerning Estimation of Power Spectra
• The Pseudo-$C_l$ method: Cosmic microwave background anisotropy power spectrum statistics for high precision cosmology
• Cosmic microwave background anisotropy power spectrum statistics for high precision cosmology
• First Intrinsic Anisotropy Observations with the Cosmic Background Imager
• An Efficient Technique to Determine the Power Spectrum from Cosmic Microwave Background Sky Maps
• Application of Monte Carlo Algorithms to the Bayesian Analysis of the Cosmic Microwave Background
• Improved Measurement of the Angular Power Spectrum of Temperature Anisotropy in the CMB from Two New Analyses of BOOMERANG Observations
• The CMB power spectrum out to l=1400 measured by the VSA
• The Statistical Significance of the Low CMB Mulitipoles
• Wiener Filtering of the COBE DMR Data
• Fast, exact CMB power spectrum estimation for a certain class of observational strategies
• The Wiener-Filtered COBE DMR Data and Predictions for the Tenerife Experiment
• Power Spectrum Estimators For Large CMB Datasets
• Fast Power Spectrum Estimation
• Cosmic Microwave Background Anisotropy in the COBE DMR 4-yr Sky Maps
• Approximating the Likelihood Function of CMB Experiments

Cited by in corpus (167)

• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP data
• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Power Spectra and WMAP-Derived Parameters
• Bayes in the sky: Bayesian inference and model selection in cosmology
• Weak Gravitational Lensing of the CMB
• Planck 2013 results. XV. CMB power spectra and likelihood
• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Are There Cosmic Microwave Background Anomalies?
• Is the low-l microwave background cosmic?
• Planck 2018 results. IV. Diffuse component separation
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• Planck intermediate results. LVII. Joint Planck LFI and HFI data processing
• Joint Bayesian component separation and CMB power spectrum estimation
• Planck 2013 results. XII. Component separation
• On the large-angle anomalies of the microwave sky
• Planck 2015 results. X. Diffuse component separation: Foreground maps
• Likelihood Analysis of CMB Temperature and Polarization Power Spectra
• The Atacama Cosmology Telescope: likelihood for small-scale CMB data
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Bayesian analysis of sparse anisotropic universe models and application to the 5-yr WMAP data
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• CMB component separation by parameter estimation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Bayesian Estimation of CMB Polarization Maps
• Constraints on primordial non-Gaussianity from 800,000 photometric quasars
• Algorithms for bispectra: forecasting, optimal analysis, and simulation
• Planck pre-launch status: the Planck-LFI programme
• Bayesian power-spectrum inference for Large Scale Structure data
• CMB temperature and polarisation pseudo-Cl estimators and covariances
• Bayesian analysis of an anisotropic universe model: systematics and polarization
• Reconstruction of Gaussian and log-normal fields with spectral smoothness
• Dark Energy Survey Year 3 results: curved-sky weak lensing mass map reconstruction
• Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment
• Unmasking the Masked Universe: the 2M++ catalogue through Bayesian eyes
• Fast and Efficient Template Fitting of Deterministic Anisotropic Cosmological Models Applied to WMAP Data
• A re-analysis of the three-year WMAP temperature power spectrum and likelihood
• An Affine-Invariant Sampler for Exoplanet Fitting and Discovery in Radial Velocity Data
• Bayesian optimisation for likelihood-free cosmological inference
• MADmap: A Massively Parallel Maximum-Likelihood Cosmic Microwave Background Map-Maker
• Hierarchical Cosmic Shear Power Spectrum Inference
• Cosmic Cartography of the Large-Scale Structure with Sloan Digital Sky Survey Data Release 6
• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• Bayesian Power Spectrum Analysis of the First-Year WMAP data
• Exploring Cosmic Origins with CORE: B-mode Component Separation
• Estimation of Polarized Power Spectra by Gibbs sampling
• Recovering the nonlinear density field from the galaxy distribution with a Poisson-Lognormal filter
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• 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
• Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator
• Efficient Wiener filtering without preconditioning
• CMB power spectrum estimation using wavelets
• Large-Angle Correlations in the Cosmic Microwave Background
• Methods for Bayesian power spectrum inference with galaxy surveys
• Bayesian inference from photometric redshift surveys
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• First measurement of the bulk flow of nearby galaxies using the cosmic microwave background
• Inference with minimal Gibbs free energy in information field theory
• Bayesian 3d velocity field reconstruction with VIRBIuS
• New estimates of the CMB angular power spectra from the WMAP 5 yrs low resolution data
• Searching for a dipole modulation in the large-scale structure of the Universe
• Filling in CMB map missing data using constrained Gaussian realizations
• Bayesian component separation and CMB estimation for the 5-year WMAP temperature data
• The joint large-scale foreground-CMB posteriors of the 3-year WMAP data
• TEASING: a fast and accurate approximation for the low multipole likelihood of the Cosmic Microwave Background temperature
• Point-Source Power in 3 Year Wilkinson Microwave Anisotropy Probe Data
• A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data
• D${}^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations
• Properties and use of CMB power spectrum likelihoods
• The integrated bispectrum as a test of CMB non-Gaussianity: detection power and limits on f_NL with WMAP data
• Using rotation measure grids to detect cosmological magnetic fields -- a Bayesian approach
• Bayesian semi-blind component separation for foreground removal in interferometric 21-cm observations
• Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps
• Matrix-free Large Scale Bayesian inference in cosmology
• The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field
• Likelihood methods for CMB experiments
• Simulations and cosmological inference: A statistical model for power spectra means and covariances
• Optimal Image Reconstruction in Radio Interferometry
• $B$-mode forecast of CMB-Bh$\overline{a}$rat
• Bayesian estimation of our local motion from the Planck-2018 CMB temperature map
• Methods for pixel domain correction of EB leakage
• Foreground Separation and Constraints on Primordial Gravitational Waves with the PICO Space Mission
• A new approach for the statistical denoising of Planck interstellar dust polarization data
• BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data
• CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator
• Weak lensing mass reconstruction using sparsity and a Gaussian random field
• CMB Tomography: Reconstruction of Adiabatic Primordial Scalar Potential Using Temperature and Polarization Maps
• BeyondPlanck XIII. Intensity foreground sampling, degeneracies, and priors
• Cosmology with the largest galaxy cluster surveys: Going beyond Fisher matrix forecasts
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling
• Wiener filter reloaded: fast signal reconstruction without preconditioning
• Harmonic in-painting of CMB sky by constrained Gaussian realization
• A comparison of CMB Angular Power Spectrum Estimators at Large Scales: the TT case
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• Primordial power spectrum and cosmology from black-box galaxy surveys
• Pure E and B polarization maps via Wiener filtering
• Bayesian power-spectrum inference with foreground and target contamination treatment
• BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
• BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
• Harmonic Inpainting of the Cosmic Microwave Background Sky: Formulation and Error Estimate
• Hybrid-basis inference for large-scale galaxy clustering: combining spherical and Cartesian Fourier analyses
• BeyondPlanck VI. Noise characterization and modelling
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• Bayesian angular power spectrum analysis of interferometric data
• Bayesian astrostatistics: a backward look to the future
• Full-sky beam convolution for cosmic microwave background applications
• Semi-blind Bayesian inference of CMB map and power spectrum
• Likelihood, Fisher information, and systematics of cosmic microwave background experiments
• Data analysis methods for the cosmic microwave background
• Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze
• Luminous giants populate the dense Cosmic Web: The radio luminosity-environmental density relation for radio galaxies in action
• Joint Bayesian analysis of large angular scale CMB temperature anomalies
• Hierarchical Bayesian CMB Component Separation with the No-U-Turn Sampler
• Marginal distributions for cosmic variance limited CMB polarization data
• Hierarchical probabilistic inference of the color-magnitude diagram and shrinkage of stellar distance uncertainties
• Probing local non-Gaussianities within a Bayesian framework
• Improved large scales interstellar dust foreground model and CMB solar dipole measurement
• Large-Scale Polarized Foreground Component Separation for Planck
• Optimized Large-Scale CMB Likelihood And Quadratic Maximum Likelihood Power Spectrum Estimation
• Testing the ABS method with the simulated Planck temperature maps
• BeyondPlanck X. Planck LFI frequency maps with sample-based error propagation
• Optimal and fast E/B separation with a dual messenger field
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• Non-Circular beam correction to the CMB power spectrum
• Sparse Inpainting and Isotropy
• Filamentary Baryons and Where to Find Them: A forecast of synchrotron radiation from merger and accretion shocks in the local Cosmic Web
• BeyondPlanck IV. On end-to-end simulations in CMB analysis -- Bayesian versus frequentist statistics
• Joint Bayesian estimation of tensor and lensing B-modes in the power spectrum of CMB polarization data
• Bayesian inference on the sphere beyond statistical isotropy
• Almanac: MCMC-based signal extraction of power spectra and maps on the sphere
• Fast scalar quadratic maximum likelihood estimators for the CMB B-mode power spectrum
• Forecast of foreground cleaning strategies for AliCPT-1
• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• Preconditioner-free Wiener filtering with a dense noise matrix
• Cosmological Parameters from a re-analysis of the WMAP-7 low resolution maps
• Testing synchrotron models and frequency resolution in BINGO 21 cm simulated maps using GNILC
• Bayesian noise estimation for non-ideal CMB experiments
• Probabilistic Cosmological Mass Mapping from Weak Lensing Shear
• Field-level multiprobe analysis of the CMB, integrated Sachs-Wolfe effect, and the galaxy density maps
• Dynamic system classifier
• Improving self-calibration
• Bayesian large-scale structure inference and cosmic web analysis
• Enhancing CMB map reconstruction and power spectrum estimation with convolutional neural networks
• A Bayesian estimate of the CMB-large-scale structure cross-correlation
• Bayesian analysis of cosmic structures
• A Monte Carlo comparison between template-based and Wiener-filter CMB dipole estimators
• A CMB Gibbs sampler for localized secondary anisotropies
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Renormalization group computation of likelihood functions for cosmological data sets
• Approximate methods for the generation of dark matter halo catalogs in the age of precision cosmology
• Foreground removal and angular power spectrum estimation of 21 cm signal using harmonic space ILC method
• Improved renormalization group computation of likelihood functions for cosmological data sets
• An Optimal and Model-Independent Measurement of the Intracluster Pressure Profile I: Methodology and First Applications
• Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations
• Shrinking the Quadratic Estimator
• Pixel domain multi-resolution minimum variance painting of CMB maps
• A Foreground Model Independent Bayesian CMB Temperature and Polarization Signal Reconstruction and Cosmological Parameter Estimation over Large Angular Scales
• Pico: Parameters for the Impatient Cosmologist
• Diagnosing Systematic Effects Using the Inferred Initial Power Spectrum
• Using hybrid GPU/CPU kernel splitting to accelerate spherical convolutions
• Generalized Redundant Calibration of Radio Interferometers
• Non-Gaussian likelihood of weak lensing power spectra