[论文解读] Comment on "Bayesian evidence: can we beat MultiNest using traditional MCMC methods", by Rutger van Haasteren (arXiv:0911.2150)

F. Feroz, M. P. Hobson+1 more|arXiv (Cornell University)|Jan 5, 2010

Gaussian Processes and Bayesian Inference参考文献 30被引用 3

一句话总结

本文批判了 Rutger van Haasteren 声称基于 Voronoi 图划分的方法在贝叶斯证据评估中优于 MultiNest 的观点。作者认为，这种比较本质上是不公平的，因为 van Haasteren 的方法依赖于预先存在的后验样本，而 MultiNest 则执行端到端的贝叶斯推断——包括后验探索和证据计算——而无需此类假设，因此在复杂、多模态问题上更具鲁棒性和可推广性。

ABSTRACT

In arXiv:0911.2150, Rutger van Haasteren seeks to criticize the nested sampling algorithm for Bayesian data analysis in general and its MultiNest implementation in particular. He introduces a new method for evidence evaluation based on the idea of Voronoi tessellation and requiring samples from the posterior distribution obtained through MCMC based methods. He compares its accuracy and efficiency with MultiNest, concluding that it outperforms MultiNest in several cases. This comparison is completely unfair since the proposed method can not perform the complete Bayesian data analysis including posterior exploration and evidence evaluation on its own while MultiNest allows one to perform Bayesian data analysis end to end. Furthermore, their criticism of nested sampling (and in turn MultiNest) is based on a few conceptual misunderstandings of the algorithm. Here we seek to set the record straight.

研究动机与目标

• 挑战 van Haasteren 声称其基于 Voronoi 图划分的方法在贝叶斯证据评估中优于 MultiNest 的有效性。
• 澄清 van Haasteren 对嵌套采样和 MultiNest 批判中的概念性误解。
• 强调 MultiNest 提供了完整且自包含的贝叶斯推断解决方案，而其他方法则依赖于外部后验样本。
• 证明即使在高维、多模态问题中，MultiNest 也能准确计算证据，而其他方法会失效。
• 重申嵌套采样和 MultiNest 在宇宙学与天体物理学复杂推断任务中的优越性。

提出的方法

• 作者分析了 van Haasteren 的方法，该方法利用后验样本上的 Voronoi 图划分来估计贝叶斯证据。
• 他们指出，该方法并非自包含，且依赖于通过 MCMC 获得的先验后验样本，而 MultiNest 内部即可完成此过程。
• 批判的核心在于概念性错误：将一个后处理工具（van Haasteren 的方法）与一个完整的推断流程（MultiNest）进行比较。
• 作者使用具有均匀先验的高维椭球形高斯似然函数，测试 MultiNest 在维度 2 到 32 之间的证据计算准确性。
• 他们使用 MultiNest 的活跃点数和自适应效率设置，计算对数证据和信息量（H）。
• 解析的对数证据在所有维度下均为 0.0，作为准确性的基准。

实验结果

研究问题

• RQ1 基于 MCMC 获得的后验样本的方法能否在贝叶斯证据评估中优于 MultiNest？
• RQ2 将 van Haasteren 的基于 Voronoi 图的方法与 MultiNest 进行比较是否公平且有效？
• RQ3 MultiNest 是否能在高维、多模态及退化参数空间中准确计算贝叶斯证据？
• RQ4 与端到端推断工具相比，依赖外部后验样本的后处理方法存在哪些局限性？
• RQ5 在证据评估中，随着维度增加，MultiNest 的性能如何退化？

主要发现

• 在维度 2 到 32 的范围内，MultiNest 准确计算了贝叶斯证据，log(Z) ≈ 0.0，其结果在误差范围内与解析值一致。
• 对于 32 维问题，MultiNest 共执行了 1410 万次似然计算，得到 log(Z) = 0.40 ± 0.31，表明即使在高维情况下仍具有可接受的准确性。
• 在 32 维时，信息量 H 增加至 96.10，确认了复杂度呈指数级增长，但 MultiNest 依然有效。
• 当后验仅占据先验体积的 e^−96.10 时，MultiNest 仍能保持准确性，证明其在高维、低似然区域的鲁棒性。
• 作者确认，通过哈密顿采样，MultiNest 可在高达约 1000 维的条件下准确计算证据，表明其具有强大的可扩展性。
• van Haasteren (2009) 中的比较无效，因为其将一个依赖后验样本的工具与一个能内建生成后验样本的方法进行比较，导致评估从本质上就不公平。