R2OpenBUGS将DIC四舍五入为4位有效数字

Question

我正在使用R2OpenBUGS来查看许多不同的模型，我注意到每个报告的DIC最多有4位有效数字。这看起来很可疑，所以我设置了DEBUG = TRUE，果然，DIC = pD + Dbar正在被舍入。为什么是这样？我可以告诉BUGS不绕DIC，甚至更好，告诉它如何绕圆？我希望DIC至少是那个数字。

这是错误调用，下面是偏差输出：

out＆lt; - 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

Answer 1

在bugs命令中，如果添加参数digits = 7会发生什么？