распечатать DIC в R2OpenBUGS

Я работаю в "R2OpenBUGS" в R. У меня есть цепочка mcmc:

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

Мой вопрос заключается в том, как распечатать DIC, когда codaPkg = TRUE. Если codaPkg = FaLSE, то если я просто выполню print(mcmc), в конце будет напечатано значение DIC.

Но когда «codaPkf=TRUE», он просто печатает что-то вроде:

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

Итак, как я могу распечатать значение DIC сейчас?


person Babak    schedule 02.04.2017    source источник


Ответы (1)


Это должно работать:

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

person Ashkan Mirzaee    schedule 16.03.2018