5 个拟合值的贝叶斯区间答案

【问题标题】：Bayesian interval of 5 fitted values5 个拟合值的贝叶斯区间
【发布时间】：2012-12-24 19:32:46
【问题描述】：

我通过从 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 个拟合值的平均值的贝叶斯区间。知道怎么做吗？

【问题讨论】：

标签： bayesian confidence-interval winbugs credible-interval

【解决方案1】：

在 WinBUGS 模型代码中定义另一个节点

mu.mean

并监控它？

【讨论】：

【解决方案2】：

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

【讨论】：