我通过从R运行Winbugs进行了贝叶斯分析,并得出了拟合值及其贝叶斯间隔。下面是相关的Winbugs输出,其中mui是第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个拟合值的平均值的贝叶斯区间。知道怎么回事吗?
发布于 2013-01-28 19:06:46
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发布于 2013-01-14 20:44:23
在WinBUGS模型代码中定义另一个节点
均值<- mu.mean (mu[])
并对其进行监控?
https://stackoverflow.com/questions/14257482
复制相似问题