多参数似然函数评估通过二维切片的迭代方法
Evaluation of multi-parameter likelihoods through iteration of two-dimensional slices
摘要 Abstract
本文介绍了一种通过固定除两个参数外的所有参数值来解决多参数似然函数的方法。对这两个变量的评估随后通过对每个参数依次循环迭代,直至收敛。我们通过透镜化宇宙微波背景(CMB)的温度功率谱测试了该技术。由于定义功率谱的六个参数之一——功率谱振幅\(A_s\)在小偏差下接近线性,这使得计算从二维网格上的增量计算简化为一维增量计算,因此这一演示特别有效。在每次迭代步骤中,\(A_s\)与其他参数配对。迭代过程得出的参数值与由“普朗克”得出的结果一致,并且对于谱的计算在几百次调用内完成。我们进一步计算了最大多极矩\(\ell_{\text{max}}\)作为参数的函数,范围从\(\ell_{\text{max}}=959\)到2500,并发现在该范围的低端存在双峰行为。在一般情况下,即当没有变量是线性的情况下,我们识别出一些调节因素,例如在每次迭代步骤中同时改变两个参数,从而减少每一步迭代的次数。马尔可夫链蒙特卡洛(MCMC)计算一直是评估多参数函数的主要工具。对于具有准线性变量(如\(A_s\))的应用,二维迭代方法比MCMC效率高出几个数量级。
In this paper we introduce a method for resolving multi-parameter likelihoods by fixing all parameter values, but two. Evaluation of those two variables is followed by iteratively cycling through each of the parameters in turn until convergence. We test the technique on the temperature power spectrum of the lensed cosmic microwave background (CMB). That demonstration is particularly effective since one of the six parameters that define the power spectra, the power spectrum amplitude, $A_{s}$, nears linearity at small deviations, reducing computation to incrementation in one-dimension, rather than over a 2D grid. At each iterative step $A_{s}$ is paired with a different parameter. The iterative process yields parameter values in agreement with those derived by \textit{Planck}, and results are obtained within a few hundred calls for spectra. We further compute parameter values as a function of maximum multipole, $\ell_{\text{max}}$, spanning a range from $\ell_{\text{max}}$=959 to 2500, and uncover bi-modal behavior at the lower end of that range. In the general case, in which neither variable is linear, we identify moderating factors, such as changing both parameters each iterative step, reducing the number of steps per iteration. Markov chain Monte Carlo (MCMC) computation has been the dominant instrument for evaluating multi-parameter functions. For applications with a quasi-linear variable such as, $A_s$, the 2D iterative method is orders of magnitude more efficient than MCMC.