1

我使用 Runjags 运行贝叶斯模型,然后使用 coda 包转换 MCMC.list 中的输出。我使用 Gelman-Rubin 诊断(单变量)检查收敛性。

有时,PSRF 很大只是因为一个链在某个点采样了一个很大的值(PsRF 从 ~1 变为 1,2)。另一方面,PSRF 与 Runjags 接近 1。有时 Runjags 的 PSRF 比 Coda 大。

我没有发现计算上有什么区别。你知道吗?即使 PSRF 为 1,2,是否可以认为参数从图中收敛?

4

0 回答 0