Page 1 of 1

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 08 2006
by Antony Lewis
This paper gives a detailed analysis of 2002 ACBAR data, including a nice treatment of the SZ effect, which they claim is the most likely explanation of the excess power seen by CBI on small scales.

1. What's the status of the rest of the data? (out soon, or another three years?)

2. Parameters constraints seem to have \tau = 0.097\pm 0.014, an error less than half the size of that usually obtained from WMAP3. This doesn't look consistent with Fig 5.

3. Are CosmoMC data files available with recommended settings for a non-SZ analysis?

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 08 2006
by Hans Kristian Eriksen
Antony Lewis wrote: 2. Parameters constraints seem to have \tau = 0.097\pm 0.014, an error less than half the size of that usually obtained from WMAP3. This doesn't look consistent with Fig 5.
I may be wrong here, but several parameters look a little funny, I think. First, n_s for the updated WMAP likelihood is listed to be 0.966, while Jostein Kristiansen found 0.960 in his test run, I believe. At least a part of this could be due to inclusion of lensing, though.. Second, the marginalized likelihoods in Figure 5 have too many bumps and wiggles on them for my taste.

All in all, I have a slightly uneasy feeling about convergence in these results. The authors say in the text that they run CosmoMC until "the largest eigenvalue of the GR test is less than 0.1". Not sure what eigenvalue means in this respect, but if it's simply the GR statistic as outputted from CosmoMC, then I think 1.1 is a little too liberal. For publications, I generally prefer R < 1.01. But of course, I may be entirely wrong here, and perhaps eigenvalue means something different in this respect.

I also noted that there seems to be a trend for high values in Figure 3 -- except for the fifth Acbar data point, which is spot on the theoretical curve, the first nine points are high. Any idea what this could be due to? Statistical fluctuation seems unlikely..?

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 08 2006
by Carlo Contaldi
I may be wrong here, but several parameters look a little funny, I think. First, ns for the updated WMAP likelihood is listed to be 0.966, while Jostein Kristiansen found 0.960 in his test run, I believe. At least a part of this could be due to inclusion of lensing, though.. Second, the marginalized likelihoods in Figure 5 have too many bumps and wiggles on them for my taste.


The value of [tex]n_s[/tex] is due to lensing as you mentioned. We get 0.960 without lensing for the WMAP3 (updated likelihood) only case.
Second, the marginalized likelihoods in Figure 5 have too many bumps and wiggles on them for my taste...All in all, I have a slightly uneasy feeling about convergence in these results. The authors say in the text that they run CosmoMC until "the largest eigenvalue of the GR test is less than 0.1".
The largest GR eigenvalue should read 0.01 in the text, this was our criterion for convergence. The current version's 0.1 is a typo. e.g. WMAP3 only has 0.008.
I also noted that there seems to be a trend for high values in Figure 3 – except for the fifth Acbar data point, which is spot on the theoretical curve, the first nine points are high. Any idea what this could be due to? Statistical fluctuation seems unlikely..?
In power the calibration error is 12% so you should bear this in mind when doing a [tex]\chi^2[/tex]-by-eye since we haven't rescaled the power spectrum in any of the plots.

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 08 2006
by Carlo Contaldi
Antony Lewis wrote: 2. Parameters constraints seem to have \tau = 0.097\pm 0.014, an error less than half the size of that usually obtained from WMAP3. This doesn't look consistent with Fig 5.
Thanks for spotting that Antony! The errors in the table for [tex]\tau[/tex] are a typo, they should be around 0.03 as you said. Figure 5 (and median values of [tex]\tau[/tex]) are correct. [tex]\tau[/tex] line of table 4 should read

[tex]\tau=0.097^{+0.033}_{-0.030} \ \ 0.092^{+0.029}_{-0.030}\ \ 0.092^{+0.028}_{-0.029} \ \ 0.090^{+0.027}_{-0.027}[/tex].

All other numbers are unaffected.

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 09 2006
by William Holzapfel
Antony Lewis wrote:This paper gives a detailed analysis of 2002 ACBAR data, including a nice treatment of the SZ effect, which they claim is the most likely explanation of the excess power seen by CBI on small scales.

1. What's the status of the rest of the data? (out soon, or another three years?)
I suppose there is no good excuse for taking three years to publish this. As it turned out, we spent a long time sorting out the details of the calibration. In the complete data set, there is sufficient sky coverage to calibrate directly from WMAP->ACBAR and we anticipate that this will go more smoothly. It seems reasonable that we will have the final power spectrum out in ~ six months.

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 09 2006
by Hans Kristian Eriksen
Carlo Contaldi wrote:
I may be wrong here, but several parameters look a little funny, I think. First, ns for the updated WMAP likelihood is listed to be 0.966, while Jostein Kristiansen found 0.960 in his test run, I believe. At least a part of this could be due to inclusion of lensing, though.. Second, the marginalized likelihoods in Figure 5 have too many bumps and wiggles on them for my taste.


The value of [tex]n_s[/tex] is due to lensing as you mentioned. We get 0.960 without lensing for the WMAP3 (updated likelihood) only case.
Second, the marginalized likelihoods in Figure 5 have too many bumps and wiggles on them for my taste...All in all, I have a slightly uneasy feeling about convergence in these results. The authors say in the text that they run CosmoMC until "the largest eigenvalue of the GR test is less than 0.1".
The largest GR eigenvalue should read 0.01 in the text, this was our criterion for convergence. The current version's 0.1 is a typo. e.g. WMAP3 only has 0.008.
The marginalized distributions for WMAP-only still looked funny to me (take a look at the bumps at high values for tau and A, and the funny shape of n_s), so I started a run including lensing last night. The parameters I got were these:

Omega_b h^2 0.0224 +/- 0.0007
Omega_c h^2 0.107 +/- 0.008
theta 1.04 +/- 0.004
tau 0.0923 +/- 0.03
n_s 0.963 +/- 0.016
log(10^10 A_s) 3.04 +/- 0.07

That is, the shift in n_s is smaller than what's listed in Kuo et al., where the value is 0.966, and tau is also smaller. Also, the distributions I get all look smooth and nice, without any bumps and wiggles. For reference, R-1 is 0.011 for this run.

Then again, I don't have a lot of experience with lensing, so I may have forgotten some switch -- all I did, was to set CMB_lensing=T in params.ini.

Did you use the WMAP likelihood out of the box, so to speak, or did you make some other changes as well?

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 09 2006
by Carlo Contaldi
Hans Kristian Eriksen wrote: Did you use the WMAP likelihood out of the box, so to speak, or did you make some other changes as well?
The version we used may have some small differences from the lambda version as it comes directly from the development of the updated likelihood code that was done a couple of weeks ago at CITA by Mike Nolta and Jon Sievers. I will check against the official version to see if there are any differences.

There are some lensing accuracy tweaks that you can set in the CAMB params but I think that should not affect results for WMAP's relatively low [tex] l[/tex] limit.

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 10 2006
by Håvard Alnes
A stuipd question (probably): Why is the errorbar for the WMAP quadrupole extending below zero in figure 3?

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 12 2006
by Kazuhide Ichikawa
I read the paper with interest but I could not follow some of the details.

1) When they say "WMAP3+ACBAR" in sec 6 and 7, do they use all the bands in table 3 ? I thought highest (4?) [tex]\ell[/tex] bands should not be used when SZ is not included. Or does it make no difference?

2) In sec 7.3 (or table 5), I infer "CMBall" does not include the highest [tex]\ell[/tex] band of CBI but "CMBall+BIMA" does. Am I correct?
In connection with this, I would like to know whether the constraint on [tex]\alpha^{SZ}[/tex] for "CMBall+BIMA" changes when BIMA is not used.

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 14 2006
by Massimiliano Lattanzi
Hi,

I was willing to run some cosmoMC chains using the new ACBAR data, and so I was trying to organize the data in a "cosmoMC-friendly" form, i.e. a form similar to the one of the datafiles included in the cosmoMc distribution like ACBAR_lge800_.dataset. However I have some questions on this:

1. According to the paper the calibration uncertainity is 6%, and it is said that this value is nearly unchanged from the 2002 data release. However the 2002 ACBAR data file in cosmoMC reads

calib_uncertainity = 0.2

so am I missing something or this is just a more conservative estimate?

2. I cannot find, neither in the website nor in the paper, the values of the beam uncertainities. Is this ok (i.e., we have to wait the collaboration to make them available) or again am I missing something? If the values are for the moment unavailable, can I simply put

beam_uncertainity = F

in the .dataset file?

Thanks

Massimiliano

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 16 2006
by Chao-Lin Kuo
Hi Massimiliano,

The central value didn't change much, but we improve its uncertainty.
The calibration uncertainty was 10% in the previous release,
in temperature unit. The new calibration uncertainty (6%) would
make calib_uncertainity = 0.12.

The beam uncertainty is 3%, as stated in the first paper.

Thank you.

Chao-Lin

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 16 2006
by Massimiliano Lattanzi
Hi Chao-Lin

Thank you very much for your answer. Things are clearer now.

However I still cannot exactly figure out where the numbers in the fourth column (beam uncertainity) of the ACBAR_lge800_.dataset file provided with cosmomc, come from. They all amount to a few percent of the measured value, but some are large as 7-8%.

Probably this will not make so large a difference in the final MCMC results, but I was wondering anyway. Can you or someone else clarify on this?

Thanks

Massimiliano

Re: [astro-ph/0611198] Improved Measurements of the CMB Powe

Posted: November 21 2006
by Jason Dick
Håvard Alnes wrote:A stuipd question (probably): Why is the errorbar for the WMAP quadrupole extending below zero in figure 3?
It looks like those errors are symmetric about the data point, i.e. they are Gaussian errors. A more careful analysis would lead to a lower bound on the quadrupole that is closer to the data point than the upper bound, and above zero.

[astro-ph/0611198] Improved Measurements of the CMB Power Sp

Posted: November 21 2006
by Joe Zuntz
A more careful analysis would lead to a lower bound on the quadrupole that is closer to the data point than the upper bound,
Yes.

and above zero.
Ideally, but not necessarily. Some methods can have error bars that legitimately extend below zero in high noise situations (essentially where you're finding the difference between two large numbers, data and noise, where you haven't perfectly estimated the noise). Bayesian approaches can deal with these situations with a sensible prior.