我使用了OpenBUGS,它生成了MCMC输出的coda文件。要计算和绘制Gelman Rubin和Geweke诊断,我需要将此coda.odc文件转换为R中的mcmc对象?有没有办法做到这一点?或者,您是否建议我采用其他方式进行此分析?
由于
答案 0 :(得分:1)
coda包中的read.coda
,...类似于:
library(coda)
my.coda <- read.coda("Chain1.txt", "CODAIndex.txt")
可以使用.odc文件。如果您将CODA文件保存为OpenBUGS中的.txt,则可以正常工作
my.coda
将是mcmc
类型对象。