#-----------------------------------------------------------------
#  other plots of the MC results
#
#   GdA 26/1/2016
#----------------------------------------------------------------

# to be 'sourced' after simpleMC_unif_pois_binom.R

mc <- as.matrix(chain)

head(mc)

par(mfrow=c(2,2))

hist(mc[,1], nc=50, col='green', xlab='n')
hist(mc[,2], nc=50, col='blue', xlab='p')
hist(mc[,3], nc=50, col='magenta', xlab='x', main='')
hist(mc[,4], nc=50, col='red', xlab='y', main="Distribuzione di x^p")

par(mfrow=c(1,1))