Page 1 of 1

[astro-ph/0606088] A re-analysis of the three-year WMAP temp

Posted: October 04 2006
by Hans Kristian Eriksen
I just wanted to note that we have now published the Gibbs sampling data from the WMAP3 re-analysis at

http://www.astro.uio.no/~hke

under the "Research" tab. (Couldn't figure out how to link to a particular frame, so sorry for not having a direct link.. However, this will eventually all appear on Lambda, I think.)

These data form the basis of the low-l likelihood we use, and we also provide an F90 module to compute that likelihood easily. If you want to use any of this, please feel free. However, I would really appreciate it if you would send me an email if you can't make it work – I'm not sure if the README I wrote is sufficiently complete.. In fact, it would also be great if you would send me an email even if you can make it work, so that I know somebody is actually using this – if not, I'm not sure how much time I'll spend on making updates available.. :-)

Anyway, as described in the revised paper, the low-l WMAP likelihood part should be replaced by either this code or by an exact likelihood evaluation at Nside=16 and lmax=30, as the MASTER-based likelihood is not quite good enough in the l=13−30 range. Personally, I prefer the Gibbs+BR approach, since it's much faster, and significantly reduces the time spent on likelihood evaluations in the MCMC sampling process. But both are of course equivalent in terms of numerical results.

It is also worth noting that we make Gibbs sampled sky maps available at the above link. These may be useful for analyses that requires phase information. However, if you do use these for publications, I strongly recommend that 1) you make sure you really understand the properties of the Gibbs sampled sky maps properly, and 2) *really* check the convergence properties of your particular statistic! While the samples we provide give very good convergence characteristics for parameter estimation purposes up to lmax=30 or 50 or so, the story may be very different for some complicated non-Gaussianity statistic. So *please* check this carefully before going public with some funny effect :-)

Good luck!