我有一个runjags对象,它有两个混合非常好的链(链1和链3),而另一个链没有混合(链2)。如何将runjags对象修剪为只包含链1和链3?这是一个使用runjags生成JAGS模型的可重复示例(尽管这里的链混合得很好)。library(runjags)
x <- seq(1,10, by = 0.1)trim.jags <- mcm
我把R中runjags模型的输出作为mcmc.list。下面是生成3条1,000个样本链的代码。我想把所有的12条链都剪到最后的400个样品上。下面是运行runjags模型并将输出转换为mcmc.list的一些数据model {for (i in 1:N){tau <- pow(sigma, -2)}jags.