2014-10-15 6 views
1

Я использовал OpenBUGS, и он выпустил файлы кодов выхода MCMC. Чтобы рассчитать и построить диагностику Gelman Rubin и Geweke, мне нужно преобразовать этот файл coda.odc в объект mcmc в R? Есть какой-либо способ сделать это? Или вы порекомендовали мне другие способы сделать этот анализ?Как преобразовать файл OpenBUGS Coda в объект mcmc в R?

Благодаря

ответ

1

Ues read.coda в coda упаковке, ... что-то вроде:

library(coda) 
my.coda <- read.coda("Chain1.txt", "CODAIndex.txt") 

Может работать с ODC-файлов. Непосредственно работает, если вы сохраняете файлы CODA как .txt из OpenBUGS . my.coda будет объектом типа mcmc.