2016-02-29 60 views
1

我使用R2OpenBUGS查看了許多不同的模型,並且我注意到每個報告的DIC至多有4位有效數字。這看起來很可疑,所以我設置了DEBUG = TRUE,果然,DIC = pD + Dbar正在四捨五入。爲什麼是這樣?我可以告訴BUGS不圍繞DIC,或者甚至更好地告訴它如何輪迴?我想DIC至少是個位數字。R2OpenBUGS將DIC舍入爲4位有效數字

這裏是錯誤調用和下方是偏差輸出:

出< - 錯誤(DATA1,inits = initial.values,parameters.to.save = C( 「ERR」, 「B1」,」 B0" , 「pred365」),model.file =爲MyModel,n.chains = 1,n.iter = 25000,DIC = TRUE,調試= TRUE)

Dbar Dhat DIC  pD 

ALS 31430.0 29700.0 33150.0 1723.0

合計31430.0 29700.0 33150.0 1723.0

回答

0

在bugs命令中,如果添加參數digits = 7,會發生什麼情況?

+1

它是問題還是答案? – Alexan

+0

抱歉不清楚。我建議如果將參數添加到bug命令中,則有效數字的數量將會增加。這絕對適用於參數估計,但我不能100%確定這是否也適用於DIC。 – amg