5 个拟合值的贝叶斯区间

Question

我通过从 R 运行 Winbugs 进行了贝叶斯分析，并导出了拟合值及其贝叶斯区间。这是相关的 Winbugs 输出，其中 mu[i] 是第 i 个拟合值。

node     mean   2.5%    97.5%   
mu[1]   0.7699  0.6661  0.94    
mu[2]   0.8293  0.4727  1.022   
mu[3]   0.7768  0.4252  0.9707  
mu[4]   0.6369  0.4199  0.8254  
mu[5]   0.7704  0.5054  1.023   

我想要做的是找到这 5 个拟合值的均值的贝叶斯区间。知道怎么做吗？

Answer 1

Chris Jackson 的回答是正确的，但是，如果您的模型已经运行了几个小时，您将不会高兴，因为这意味着修改模型并再次运行它。但是您可以在后期处理中在 R 中实现您的目标，而无需再次运行模型 - 通过取后验样本的平均值:

out <- bugs(...)
sapply(out$sims.list$mu, mean, ...) # I'm not sure exactly about the structure of
                                    # out$sims.list$mu, so it might be slightly 
                                    # different