在R中,有用于 MCMC 的包和函数,例如jags.model允许我们控制要运行的链数,即在 内jags.model,我们可以设置n.chains = 44链。生成的跟踪图是多色的,因为它跟踪了四个独立链的后验估计。我的问题是,当我们从对象中检索后绘制或估计时,使用了链中的哪一个?是平均值吗?谢谢。
当我们在 MCMC 模型中一次运行多个链时,它们如何组合在一起以进行后向抽签?
机器算法验证
r
马尔可夫链蒙特卡罗
2022-03-23 10:34:24
1个回答
一般来说,每条链的老化迭代都被丢弃,其余的迭代被简单地组合起来,就像来自一个长链一样。拥有独立链的主要原因是为了更好地诊断收敛问题。
您可能还会发现这个问题很有帮助。