2

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つの近似値の平均のベイズ区間を見つけることです。どのようにアイデアはありますか？

4

2 に答える 2

2

クリスジャクソンの答えは正しいですが、モデルが何時間も実行されている場合、それはモデルを変更して再度実行することを意味するため、満足できません。ただし、後のサンプルを平均することで、モデルを再度実行しなくても、後処理で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-28T11:06:46.413 に答える
1

WinBUGSモデルコードで別のノードを定義します

mu.mean <-mean（mu []）

そしてそれを監視しますか？

于 2013-01-14T12:44:22.890 に答える