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))