We study different variants of the Gibbs sampler algorithm from the perspective of their applicability to the estimation of power spectra of the cosmic microwave background (CMB) anisotropies. These include approaches studied earlier in the CMB literature as well as new ones which are proposed in this work. We demonstrate all these variants on full and cut sky simulations and compare their performance, assessing both their computational and statistical efficiency. For this we employ a consistent comparison metric, an effective sample size (ESS) per second, commonly used in this context in the statistical literature. We show that one of the proposed approaches, referred to as Centered overrelax, which capitalizes on additional, auxiliary variables to minimize computational time needed per sample, and uses overrelaxation to decorrelate subsequent samples, performs better than the standard Gibbs sampler by a factor between one and two orders of magnitude in the nearly full-sky, satellite-like cases. It therefore potentially provides an interesting alternative to the currently favored approaches.
翻译:我们从适用于估计宇宙微波背景的电源光谱的角度研究Gibbs采样器算法的不同变方,这些变方包括以前在CMB文献中研究过的方法以及这项工作中提议的新方法。我们在全空和切开的天空模拟中展示所有这些变方,并比较其性能,同时评估其计算和统计效率。为此,我们采用了一致的比较指标,即每秒有效抽样规模(ESS),在统计文献中通常在此背景下使用。我们表明,其中一种拟议办法,即 " 中枢过度松动 ",即利用额外辅助变量来尽量减少每个样本所需的计算时间,并使用过度松动来调整以后的样品,在几乎全天空、卫星类的案例中,以一至两个数量级的系数比标准Gibs取样器更好。因此,它有可能为目前偏好的办法提供一个有趣的替代办法。