# Program to analyse coda files generated by OpenBugs
#
# GdA 8/04/09
#-------------

library(coda)

#- i file bugs_out.txt e bugs_ind.txt
#  sono stati salvati dall'interno di OpenBUGS
gauss.msy <-read.coda("bugs_out.txt","bugs_ind.txt",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)
layout(1)

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

cor(mu,sigma)
cor(mu,y)