library(coda)
mcmc.list <- read.openbugs('coda_output')
plot(mcmc.list)
geweke.diag(mcmc.list)
Geweke.plot(mcmc.list)
heidel.diag(mcmc.list)
gelman.diag(mcmc.list)
gelman.plot(mcmc.list)
autocorr.diag(mcmc.list)
autorcorr.plot(mcmc.list)
raftery.diag(mcmc.list, q=0.5, r=0.05)
thinned <- window(mcmc.list, 2000, 25000, 50)
