[R] Coda, HPDinterval and multiple chains

Dear all, I'm not 100% sure if this question is best directed at the r-list, or a mailing list concerned with Bayesian analysis, so please accept my apologies if another audience may be more appropriate. I have been using the rjags package to run Jags models with multiple chains and store the

Re: [R] Coda, HPDinterval and multiple chains

How about require(coda) data(line) str(line) cline - as.mcmc(do.call(rbind,line)) str(cline) # Thus HPDinterval(cline) # (or any FUN.mcmc) sum(cline[,'alpha'] = median(cline[,'alpha']))/nrow(cline) Cheers On Thu, Oct 4, 2012 at 12:22 PM, NORRIS Paul p.nor...@ed.ac.uk wrote: Dear all, I'm