library(coda)

gauss.msy <-read.coda("bugs.out","bugs.ind",quiet=TRUE)
str(gauss.msy)

sample<-as.matrix(gauss.msy)
mu <- sample[,1]
sigma <- sample[,2]
y <- sample[,3]

layout(matrix(c(1:3),c(1,3)))  
hist(mu,nc=100,freq=F)
hist(sigma,nc=100,fr=F)
hist(y,nc=100,fr=F)

summary(mu)
mean(mu) ; sd(mu)
summary(sigma)
mean(sigma) ; sd(sigma)
summary(y)
mean(y) ; sd(y)

cor(mu,sigma)