1

我在 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 值呢？

4

1 回答 1

0

这应该有效：

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

于 2018-03-16T20:26:23.673 回答