#-------------------------------------------------------------------
# script per simulare un campione gaussiano e analizzarlo con rbugs
#   
#  G. D'Agostini, gennaio 2011
#-------------------------------------------------------------------

library(rjags)

# simula i dati --------------------------------------------------
mu    = ....  # sostituire i puntini con valori numerici
sigma = ....  # sostituire i puntini con valori numerici
N     = ....  # sostituire i puntini con valori numerici
data <- NULL
data$X <- rnorm(N, mu, sigma)

# usa rjags -----------------------------------------------------
jm <- jags.model("gauss_Nm_sigma_pred.bug", data)

update(jm, 3000)

samples <- coda.samples(jm, c("mu","sigma","y"), n.iter=30000)

# analizza la catena --------------------------------------------
summary(samples)
plot(samples)

samples.mcmc <- as.mcmc(samples)
samples.df <- as.data.frame(samples.mcmc)
cor(samples.df)
var(samples.df)
# etc. etc.