我正在使用 R2OpenBUGS 查看许多不同的模型，我注意到每个报告的 DIC 最多有 4 个有效数字。这看起来很可疑，所以我设置了 DEBUG=TRUE，果然，DIC = pD + Dbar 正在四舍五入。为什么是这样？我可以告诉 BUGS 不要舍入 DIC，或者更好的是，告诉它如何舍入？我希望 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=TRUE,debug=TRUE)
Dbar Dhat DIC pD
ALS 31430.0 29700.0 33150.0 1723.0
总计 31430.0 29700.0 33150.0 1723.0