We report on the results from two independent but complementary statistical analyses of the WMAP first-year data, based on the power spectrum and N-point correlation functions. We focus on large and intermediate scales (larger than about 3 degrees) and compare the observed data against Monte Carlo ensembles with WMAP-like properties. In both analyses, we measure the amplitudes of the large-scale fluctuations on opposing hemispheres and study the ratio of the two amplitudes. The power-spectrum analysis shows that this ratio for WMAP, as measured along the axis of maximum asymmetry, is high at the 95%-99% level (depending on the particular multipole range included). The axis of maximum asymmetry of the WMAP data is weakly dependent on the multipole range under consideration but tends to lie close to the ecliptic axis. In the N-point correlation function analysis we focus on the northern and southern hemispheres defined in ecliptic coordinates, and we find that the ratio of the large-scale fluctuation amplitudes is high at the 98%-99% level. Furthermore, the results are stable with respect to choice of Galactic cut and also with respect to frequency band. A similar asymmetry is found in the COBE-DMR map, and the axis of maximum asymmetry is close to the one found in the WMAP data.

We perform a metastudy of recently published redshift space distortion (RSD) measurements of the cosmological growth rate, $f(z){\sigma }_8(z)$. We analyze the latest results from the 6dFGS, BOSS, LRG, WiggleZ, and VIPERS galaxy redshift surveys, and compare the measurements to expectations from Planck. In this Letter we point out that the RSD measurements are consistently lower than the values expected from Planck, and the relative scatter between the RSD measurements is lower than expected. A full resolution of this issue may require a more robust treatment of nonlinear effects in RSD models, although the trend for a low ${\sigma }_8$ agrees with recent constraints on ${\sigma }_8$ and ${\Omega }_m$ from Sunyaev-Zeldovich cluster counts identified in Planck.

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional sampling of foreground spectral parameters, and 2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground-CMB posterior distribution, and therefore all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multi-resolution observations. To verify the method, we analyse simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3-yr WMAP data, downgraded to a common resolution of 3 degree FWHM. The results from the actual 3-yr WMAP temperature analysis are presented in a companion Letter.Comment: 23 pages, 16 figures; version accepted for publication in ApJ -- only minor changes, all clarifications. More information about the WMAP3 analysis available at http://www.astro.uio.no/~hke under the Research ta

We study the Internal Linear Combination (ILC) method presented by the Wilkinson Microwave Anisotropy Probe (WMAP) science team, with the goal of determining whether it may be used for cosmological purposes, as a template-free alternative to existing foreground correction methods. We conclude that the method does have the potential to do just that, but great care must be taken both in implementation, and in a detailed understanding of limitations caused by residual foregrounds which can still affect cosmological results. As a first step we demonstrate how to compute the ILC weights both accurately and efficiently by means of Lagrange multipliers, and apply this method to the observed data to produce a new version of the ILC map. Next we describe how to generate Monte Carlo simulations of the ILC map, and find that these agree well with the observed map on angular scales up to l~200. Finally we make two comments to the on-going debates concerning the large-scale properties of the WMAP data. First, we note that the Galactic south-eastern quadrant is associated with notably different ILC weights than the other three quadrants, possibly indicating a foreground related anisotropy. Second, we study the properties of the quadrupole and octopole (amplitude, alignment and planarity), and confirm the previously reported anomalies. Even more interestingly, we find that the l=5 mode is spherically symmetric at about 3 sigma, and that the l=6 mode is planar at the 2 sigma level. However, we also assess the impact of residual foregrounds on these statistics, and find that the ILC map is not clean enough to allow for cosmological conclusions. Alternative methods must be developed to study these issues further.Comment: 15 pages, 9 figures, accepted for publication in Ap

Motivated by the large-scale asymmetry observed in the cosmic microwave background sky, we consider a specific class of anisotropic cosmological models -- Bianchi type VII_h -- and compare them to the WMAP first-year data on large angular scales. Remarkably, we find evidence of a correlation which is ruled out as a chance alignment at the 3sigma level. The best fit Bianchi model corresponds to x=0.55, Omega_0=0.5, a rotation axis in the direction (l,b)=(222degr,-62degr), shear (sigma/H)_0=2.4e-10 and a right--handed vorticity (omega/H)_0=6.1e-10. Correcting for this component greatly reduces the significance of the large-scale power asymmetry, resolves several anomalies detected on large angular scales (ie. the low quadrupole amplitude and quadrupole/octopole planarity and alignment), and can account for a non--Gaussian "cold spot" on the sky. Despite the apparent inconsistency with the best-fit parameters required in inflationary models to account for the acoustic peaks, we consider the results sufficiently provocative to merit further consideration.

(Abridged)Motivated by the recent results of Hansen et al. (2008) concerning a noticeable hemispherical power asymmetry in the WMAP data on small angular scales, we revisit the dipole modulated signal model introduced by Gordon et al. (2005). This model assumes that the true CMB signal consists of a Gaussian isotropic random field modulated by a dipole, and is characterized by an overall modulation amplitude, A, and a preferred direction, p. Previous analyses of this model has been restricted to very low resolution due to computational cost. In this paper, we double the angular resolution, and compute the full corresponding posterior distribution for the 5-year WMAP data. The results from our analysis are the following: The best-fit modulation amplitude for l <= 64 and the ILC data with the WMAP KQ85 sky cut is A=0.072 +/- 0.022, non-zero at 3.3sigma, and the preferred direction points toward Galactic coordinates (l,b) = (224 degree, -22 degree) +/- 24 degree. The corresponding results for l <~ 40 from earlier analyses was A = 0.11 +/- 0.04 and (l,b) = (225 degree,-27 degree). The statistical significance of a non-zero amplitude thus increases from 2.8sigma to 3.3sigma when increasing l_max from 40 to 64, and all results are consistent to within 1sigma. Similarly, the Bayesian log-evidence difference with respect to the isotropic model increases from Delta ln E = 1.8 to Delta ln E = 2.6, ranking as "strong evidence" on the Jeffreys' scale. The raw best-fit log-likelihood difference increases from Delta ln L = 6.1 to Delta ln L = 7.3. Similar, and often slightly stronger, results are found for other data combinations. Thus, we find that the evidence for a dipole power distribution in the WMAP data increases with l in the 5-year WMAP data set, in agreement with the reports of Hansen et al. (2008).