library(coda); args <- commandArgs(TRUE); print(args); chains <- list() for (i in args) { x <- read.table(i, header=T); x <- mcmc(x); chains[[length(chains) + 1]] <- x; print(summary(x)); plot(x); autocorr.plot(x); } chains <- mcmc.list(chains) gelman.diag(chains);