私は R2OpenBUGS を使用してさまざまなモデルを調べていますが、報告されたすべての DIC の有効数字が最大で 4 桁であることに気付きました。これは疑わしいと思われたので、DEBUG=TRUE に設定すると、DIC = pD + Dbar が丸められます。どうしてこれなの？DIC を丸めないように BUGS に指示することはできますか? DICを少なくとも一桁にしたい。
これがバグコールであり、以下は逸脱出力です。
out <- bugs(data1, inits=initial.values,parameters.to.save = c("err","b1","b0","pred365"), model.file = mymodel, n.chains = 1, n.iter = 25000、DIC=真、デバッグ=真)
Dbar Dhat DIC pD
ALS 31430.0 29700.0 33150.0 1723.0
合計 31430.0 29700.0 33150.0 1723.0