一种改进的马尔可夫链蒙特卡抽样方法,用于CMB数据中的宇宙学参数估计
天体物理学
2007-05-23 v1
摘要
马尔可夫链蒙特卡抽样(MCMC)已成为探索由CMB各向异性观测约束的宇宙学参数后验分布的标准技术。给定足够的时间,任何MCMC抽样器最终都会收敛,其平稳分布即为感兴趣的后验分布。然而,实际中,朴素的抽样器需要相当长的时间才能充分探索后验分布。在最佳情况下,这仅导致浪费的CPU时间,但在更糟的情况下,可能导致对宇宙学参数值的置信区间被低估。即使对于当前的CMB数据集,广泛使用的cosmomc软件包中所使用的抽样器也未能高效采样。对于预期将用于Planck任务的数据集,这种困难更为明显。因此,我们提出了一种用于分析总能量CMB观测的新型MCMC抽样器,可轻松集成到cosmomc软件中,但具有快速收敛特性并能产生可靠的置信区间。通过使用用于提议分布的动态宽度、动态协方差矩阵抽样以及针对沿已知退化方向移动的专用提议分布,实现了上述目标。
关键词
引用
@article{arxiv.astro-ph/0307219,
title = {An improved Markov-chain Monte Carlo sampler for the estimation of cosmological parameters from CMB data},
author = {Anze Slosar and Michael Hobson},
journal= {arXiv preprint arXiv:astro-ph/0307219},
year = {2007}
}
备注
8 pages, 6 figures, submitted to MNRAS