library(coda) gauss.1m <-read.coda("bugs.out","bugs.ind",quiet=TRUE) str(gauss.1m) densplot(gauss.1m) s<-as.vector(gauss.1m) plot(s[1:1000],ty="l") hist(s,nc=100) mean(s) sd(s) summary(s) quantile(s,c(0.01,0.025,0.05,0.1,0.16,0.84,0.9,0.95,0.99))