Introducing Mexican needlets for CMB analysis: Issues for practical applications and comparison with standard needlets

Over the last few years, needlets have a emerged as a useful tool for the analysis of Cosmic Microwave Background (CMB) data. Our aim in this paper is first to introduce in the CMB literature a different form of needlets, known as Mexican needlets, first discussed in the mathematical literature by Geller and Mayeli (2009a,b). We then proceed with an extensive study of the properties of both standard and Mexican needlets; these properties depend on some parameters which can be tuned in order to optimize the performance for a given application. Our second aim in this paper is then to give practical advice on how to adjust these parameters in order to achieve the best properties for a given problem in CMB data analysis. In particular we investigate localization properties in real and harmonic spaces and propose a recipe on how to quantify the influence of galactic and point source masks on the needlet coefficients. We also show that for certain parameter values, the Mexican needlets provide a close approximation to the Spherical Mexican Hat Wavelets (whence their name), with some advantages concerning their numerical implementation and the derivation of their statistical properties.Comment: 40 pages, 11 figures, published version, main modification: added section on more realistic galactic and point source mask

Directional Variations of the Non-Gaussianity Parameter f_NL

We investigate local variations of the primordial non-Gaussianity parameter f_NL in the WMAP data, looking for possible influence of foreground contamination in the full-sky estimate of f_NL. We first improve the needlet bispectrum estimate in (Rudjord et al. 2009) on the full-sky to f_NL= 73 +/- 31 using the KQ75 mask on the co-added V+W channel. We find no particular values of f_NL estimates close to the galactic plane and conclude that foregrounds are unlikely to affect the estimate of f_NL in the V and W bands even for the smaller KQ85 mask. In the Q band however, we find unexpectedly high values of f_NL in local estimates close to the galactic mask, as well as significant discrepancies between Q band estimates and V/W band estimates. We therefore conclude that the Q band is too contaminated to be used for non-Gaussianity studies even with the larger KQ75 mask. We further noted that the local f_NL estimates on the V+W channel are positive on all equatorial bands from the north to the south pole. The probability for this to happen in a universe with f_NL = 0 is less than one percent.Comment: 6 pages, 2 figures, included section on non-Gaussian maps, references adde

Data Quality of Surgery for Carotid Artery Stenosis. Are the National Vascular Registries Reliable?

AbstractObjectivesTo study completeness of reporting carotid endarterectomies, including peri-operative stroke and mortality rate, in a national vascular registry, NorKar, and a national administrative registry, The Norwegian Patient Register (NPR).DesignComparative registry-based national study.MaterialsMember hospitals of NorKar, including 89% of carotid endarterectomies in Norway, were compared with relevant data in NPR for the years 2000–2002.MethodsWe compared procedure-codes, diagnosis-codes, in-hospital death and the occurrence of peri-operative stroke after treatment for carotid artery stenosis in the two registries to evaluate completeness.ResultsCompared with the NPR numbers, 16% of carotid endarterectomies were missing in the reports from member hospitals of NorKar. Further, during this three-year period, there was an under-reporting of seven strokes and two deaths. The discrepancy was most pronounced in 2001.ConclusionsThere is an under-reporting of patients operated on for carotid artery stenosis in NorKar according to NPR numbers as well as an under-reporting of early deaths and strokes. There is a need for better quality data in the NorKar Registry. Registry quality would be likely to improve if patient identifiable data were available in both registries

On the linear term correction for needlets/wavelets non-Gaussianity estimators

We derive the linear correction term for needlet and wavelet estimators of the bispectrum and the non-linearity parameter fNL on cosmic microwave background radiation data. We show that on masked WMAP-like data with anisotropic noise, the error bars improve by 10-20% and almost reach the optimal error bars obtained with the KSW estimator (Komatsu et al 2005). In the limit of full-sky and isotropic noise, this term vanishes. We apply needlet and wavelet estimators to the WMAP 7-year data and obtain our best estimate fNL=37.5 \pm 21.8.Comment: 10 pages, submitted to Ap

Non-Gaussianities in the local curvature of the 5-year WMAP data

Using the 5 year WMAP data, we re-investigate claims of non-Gaussianities and asymmetries detected in local curvature statistics of the 1 year WMAP data. In Hansen et al 2004, it was found that the northern ecliptic hemisphere was non-Gaussian at the ~1% level testing the densities of hill-, lake and saddle points based on the second derivatives of the CMB temperature map. The 5 year WMAP data has a much lower noise level and better control of systematics. Using these, we find that the anomalies are still present at a consistent level. Also the direction of maximum non-Gaussianity remains. Due to limited availability of computer resources, Hansen et al. 2004 were unable to calculate the full covariance matrix for the chi^2 test used. Here we apply the full covariance matrix instead of the diagonal approximation and find that the non-Gaussianities disappear and there is no preferred non-Gaussian direction. We compare with simulations of weak lensing to see if this may cause the observed non-Gaussianity when using diagonal covariance matrix. We conclude that weak lensing does not produce non-Gaussianity in the local curvature statistics at the scales investigated in this paper. The cause of the non-Gaussian detection in the case of a diagonal matrix remains unclear.Comment: 7 pages, 6 figures, included test on weak lensing simulation

Marginal distributions for cosmic variance limited CMB polarization data

We provide computationally convenient expressions for all marginal distributions of the polarization CMB power spectrum distribution P(C_l|sigma_l), where C_l = {C_l^TT, C_l^TE, C_l^EE, C_l^BB} denotes the set of ensemble averaged polarization CMB power spectra, and sigma_l = {sigma_l^TT, sigma_l^TE, sigma_l^EE, sigma_l^BB} the set of the realization specific polarization CMB power spectra. This distribution describes the CMB power spectrum posterior for cosmic variance limited data. The expressions derived here are general, and may be useful in a wide range of applications. Two specific applications are described in this paper. First, we employ the derived distributions within the CMB Gibbs sampling framework, and demonstrate a new conditional CMB power spectrum sampling algorithm that allows for different binning schemes for each power spectrum. This is useful because most CMB experiments have very different signal-to-noise ratios for temperature and polarization. Second, we provide new Blackwell-Rao estimators for each of the marginal polarization distributions, which are relevant to power spectrum and likelihood estimation. Because these estimators represent marginals, they are not affected by the exponential behaviour of the corresponding joint expression, but converge quickly.Comment: 8 pages, 3 figures; minor adjustment, accepted for publication in ApJ

CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator

We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{\ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher l's than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s.Comment: 10 pages, 7 figures, submitted to Ap

Non-Gaussian bubbles in the sky

We point out a possible generation mechanism of non-Gaussian bubbles in the sky due to bubble nucleation in the early universe. We consider a curvaton scenario for inflation and assume that the curvaton field phi, whose energy density is subdominant during inflation but which is responsible for the curvature perturbation of the universe, is coupled to another field sigma which undergoes false vacuum decay through quantum tunneling. For this model, we compute the skewness of the curvaton fluctuations due to its interaction with sigma during tunneling, that is, on the background of an instanton solution that describes false vacuum decay. We find that the resulting skewness of the curvaton can become large in the spacetime region inside the bubble. We then compute the corresponding skewness in the statistical distribution of the cosmic microwave background (CMB) temperature fluctuations. We find a non-vanishing skewness in a bubble-shaped region in the sky. It can be large enough to be detected in the near future, and if detected it will bring us invaluable information about the physics in the early universe.Comment: 6 pages, 6 figure

Inhomogeneous non-Gaussianity

We propose a method to probe higher-order correlators of the primordial density field through the inhomogeneity of local non-Gaussian parameters, such as f_NL, measured within smaller patches of the sky. Correlators between n-point functions measured in one patch of the sky and k-point functions measured in another patch depend upon the (n+k)-point functions over the entire sky. The inhomogeneity of non-Gaussian parameters may be a feasible way to detect or constrain higher-order correlators in local models of non-Gaussianity, as well as to distinguish between single and multiple-source scenarios for generating the primordial density perturbation, and more generally to probe the details of inflationary physics.Comment: 16 pages, 2 figures; v2: Minor changes and references added. Matches the published versio

TEASING: a fast and accurate approximation for the low multipole likelihood of the Cosmic Microwave Background temperature

We explore the low-l likelihood of the angular spectrum C(l) of masked CMB temperature maps using an adaptive importance sampler. We find that, in spite of a partial sky coverage, the likelihood distribution of each C(l) closely follows an inverse gamma distribution. Our exploration is accurate enough to measure the inverse gamma parameters along with the correlation between multipoles. Those quantities are used to build an approximation of the joint posterior distribution of the low-l likelihood. The accuracy of the proposed approximation is established using both statistical criteria and a mock cosmological parameter fit. When applied to the WMAP5 data set, this approximation yields cosmological parameter estimates at the same level of accuracy as the best current techniques but with very significant speed gains.Comment: 10 pages, 10 figures, submitted to MNRA