# analisi R del risultato della mcmc 

library(coda)
pfit.mcmc <- read.coda("bugs.out","bugs.ind",quiet=TRUE)
plot(pfit.mcmc)
summary(pfit.mcmc)

pfit.m <- as.matrix(pfit.mcmc)
for (i in 1:4) print(mean(pfit.m[,i]))
for (i in 1:4) print(sd(pfit.m[,i]))