-
Estimation of Polarized Power Spectra by Gibbs sampling
Authors:
D. L. Larson,
H. K. Eriksen,
B. D. Wandelt,
K. M. Gorski,
Greg Huey,
J. B. Jewell,
I. J. O'Dwyer
Abstract:
Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measuremen…
▽ More
Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measurements than for temperature measurements. While approximate methods can alias power from the larger E-mode spectrum into the weaker B-mode spectrum, the Gibbs sampler (or equivalently, exact likelihood evaluations) allows for a statistically optimal separation of these modes in terms of power spectra. To demonstrate the method, we analyze two simulated data sets: 1) a hypothetical future CMBPol mission, with the focus on B-mode estimation; and 2) a Planck-like mission, to highlight the computational feasibility of the method.
△ Less
Submitted 2 October, 2006; v1 submitted 31 July, 2006;
originally announced August 2006.
-
A re-analysis of the three-year WMAP temperature power spectrum and likelihood
Authors:
H. K. Eriksen,
Greg Huey,
R. Saha,
F. K. Hansen,
J. Dick,
A. J. Banday,
K. M. Gorski,
P. Jain,
J. B. Jewell,
L. Knox,
D. L. Larson,
I. J. O'Dwyer,
T. Souradeep,
B. D. Wandelt
Abstract:
We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do…
▽ More
We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions.
△ Less
Submitted 2 October, 2006; v1 submitted 5 June, 2006;
originally announced June 2006.
-
A Statistically Robust 3-Sigma Detection of Non-Gaussianity in the WMAP Data Using Hot and Cold Spots
Authors:
David L. Larson,
Benjamin D. Wandelt
Abstract:
We present a careful frequentist analysis of one- and two-point statistics of the hot and cold spots in the cosmic microwave background (CMB) data obtained by the Wilkinson Microwave Anisotropy Probe (WMAP). Our main result is the detection of a new anomaly at the 3-sigma level using temperature-weighted extrema correlation functions. We obtain this result using a simple hypothesis test which re…
▽ More
We present a careful frequentist analysis of one- and two-point statistics of the hot and cold spots in the cosmic microwave background (CMB) data obtained by the Wilkinson Microwave Anisotropy Probe (WMAP). Our main result is the detection of a new anomaly at the 3-sigma level using temperature-weighted extrema correlation functions. We obtain this result using a simple hypothesis test which reduces the maximum risk of a false detection to the same level as the claimed significance of the test. We further present a detailed study of the robustness of our earlier claim (Larson and Wandelt 2004) under variations in the noise model and in the resolution of the map. Free software which implements our test is available online.
△ Less
Submitted 3 May, 2005;
originally announced May 2005.
-
Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator
Authors:
M. Chu,
H. K. Eriksen,
L. Knox,
K. M. Gorski,
J. B. Jewell,
D. L. Larson,
I. J. O'Dwyer,
B. D. Wandelt
Abstract:
We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence…
▽ More
We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence of the BR estimate rises rapidly with increasing l, at least at low l. Our existing sample chains are only long enough to achieve convergence at l less than about 40. In comparison with P(C_l|d) as reported by the WMAP team we find significant differences at these low l values. These differences lead to up to approximately 0.5 sigma shifts in the estimates of parameters in a 7-parameter Lambda CDM model with non-zero dn_s/dln(k). Fixing dn_s/dln(k)= 0 makes these shifts much less significant.
Unlike existing analytic approximations, the BR estimator can be straightforwardly extended for the case of power spectra from correlated fields, such as temperature and polarization. We discuss challenges to extending the procedure to higher l and provide some solutions.
△ Less
Submitted 26 November, 2004;
originally announced November 2004.
-
Power spectrum estimation from high-resolution maps by Gibbs sampling
Authors:
H. K. Eriksen,
I. J. O'Dwyer,
J. B. Jewell,
B. D. Wandelt,
D. L. Larson,
K. M. Gorski,
S. Levin,
A. J. Banday,
P. B. Lilje
Abstract:
We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify tha…
▽ More
We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify that a particular solution to any given problem does not affect the scientific results. We then apply these programs to simulated data with known properties at intermediate (N_side = 128) and high (N_side = 512) resolutions, to study effects such as incomplete sky coverage and white vs. correlated noise. From these simulations we also establish the Markov chain correlation length as a function of signal-to-noise ratio, and give a few comments on the properties of the correlation matrices involved. Parallelization issues are also discussed, with emphasis on real-world limitations imposed by current super-computer facilities. The scientific results from the analysis of the first-year WMAP data are presented in a companion letter.
△ Less
Submitted 23 August, 2004; v1 submitted 1 July, 2004;
originally announced July 2004.
-
Bayesian Power Spectrum Analysis of the First-Year WMAP data
Authors:
I. J. O'Dwyer,
H. K. Eriksen,
B. D. Wandelt,
J. B. Jewell,
D. L. Larson,
K. M. Gorski,
A. J. Banday,
S. Levin,
P. B. Lilje
Abstract:
We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotrop…
▽ More
We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation.
△ Less
Submitted 7 December, 2004; v1 submitted 1 July, 2004;
originally announced July 2004.
-
The Hot and Cold Spots in the WMAP Data are Not Hot and Cold Enough
Authors:
David L. Larson,
Benjamin D. Wandelt
Abstract:
This paper presents a frequentist analysis of the hot and cold spots of the cosmic microwave background data collected by the Wilkinson Microwave Anisotropy Probe (WMAP). We compare the WMAP temperature statistics of extrema (number of extrema, mean excursion, variance, skewness and kurtosis of the excursion) to Monte-Carlo simulations. We find that, on average, the local maxima (high temperatur…
▽ More
This paper presents a frequentist analysis of the hot and cold spots of the cosmic microwave background data collected by the Wilkinson Microwave Anisotropy Probe (WMAP). We compare the WMAP temperature statistics of extrema (number of extrema, mean excursion, variance, skewness and kurtosis of the excursion) to Monte-Carlo simulations. We find that, on average, the local maxima (high temperatures in the anisotropy) are too cold and the local minima are too warm. In order to quantify this claim we describe a two-sided statistical hypothesis test which we advocate for other investigations of the Gaussianity hypothesis. Using this test we reject the isotropic Gaussian hypothesis at more than 99% confidence in a well-defined way. Our claims are based only on regions that are outside the most conservative WMAP foreground mask. We perform our test separately on maxima and minima, and on the north and south ecliptic and Galactic hemispheres and reject Gaussianity at above 95% confidence for almost all tests of the mean excursions. The same test also shows the variance of the maxima and minima to be low in the ecliptic north (99% confidence), but consistent in the south; this effect is not as pronounced in the Galactic north and south hemispheres.
△ Less
Submitted 25 October, 2004; v1 submitted 2 April, 2004;
originally announced April 2004.
-
Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
Authors:
Benjamin D. Wandelt,
David L. Larson,
Arun Lakshminarayanan
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 fore…
▽ More
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.
△ Less
Submitted 5 October, 2003; v1 submitted 2 October, 2003;
originally announced October 2003.
