r - 如何将 rjags 中的多个链组合成 R 中的一条链？

Question

我通常从带有多个链的 rjags 调用 JAGS 以用于诊断目的(例如 4 个链)。之后，我经常想对后验参数估计进行一些后处理(例如，使用预测值、计算附加统计数据)。然而，此时将链存储在列表中是很麻烦的。

将链组合成单个参数列表的好方法是什么？

Answer 1

runjags 包具有函数 combine.mcmc 。它的默认设置是组合一个或多个链并返回单个链。例如，

library(runjags)
fit <- combine.mcmc(multichainfit)

它还有其他组合链的选项。