在 R2OpenBUGS 中打印 DIC

Question

我在 R 的“R2OpenBUGS”中工作。我有一些 mcmc 链：

mcmc <- bugs(data = ... DIC=TRUE ... codaPkg=TRUE)

我的问题是当 codaPkg = TRUE 时如何打印 DIC。如果 codaPkg = FaLSE，那么如果我只执行 print(mcmc)，它会在最后打印 DIC 值。

但是当 "codaPkf=TRUE" 时，它只会打印如下内容：

[1] "/tmp/Rtmp3OlPh5/CODAchain1.txt" "/tmp/Rtmp3OlPh5/CODAchain2.txt"

那么，我现在如何打印 DIC 值？

Answer 1

这应该可行：

bugs.log(paste(substr(mcmc[1],1,nchar(mcmc[1])-14),"log.txt", sep = ""))$DIC