我通常从带有几个链条的rjags中调用JAGS，以进行诊断(例如4条链条)。
之后，我经常想对后验参数估计值进行一些后期处理(例如，使用预测值，计算其他统计信息)。但是，此时将链存储在列表中是很麻烦的。
将链条组合成单个参数列表的好方法是什么？
最佳答案
runjags软件包具有功能 combine.mcmc 。它的默认设置是合并一个或多个链并返回一个链。例如。，
library(runjags)
fit <- combine.mcmc(multichainfit)
它还具有用于组合链的其他选项。