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
它是問題還是答案? – Alexan
抱歉不清楚。我建議如果將參數添加到bug命令中,則有效數字的數量將會增加。這絕對適用於參數估計,但我不能100%確定這是否也適用於DIC。 – amg