# Program to analyse coda files generated by OpenBugs
#
# GdA 23/10/2010


library(coda)

#- i file bugs_out.txt e bugs_ind.txt
#  sono stati salvati dall'interno di OpenBUGS
angoli <-read.coda("bugs_out.txt","bugs_ind.txt",quiet=TRUE)
str(angoli)
plot(angoli)
summary(angoli)

storia <- as.matrix(angoli)

somma <- storia[,1] + storia[,2] +  storia[,3]

summary(somma)


cor(storia[,1],storia[,2])
cor(storia[,1],storia[,3])
cor(storia[,2],storia[,3])