#-------------------------------------------------------------------
# inferenza con vincolo: angoli di triangolo
#   
#  G. D'Agostini, marzo 2011
#-------------------------------------------------------------------
library(rjags)

dati <- NULL
dati$a <- 58  # gradi
dati$b <- 73  # gradi
dati$g <- 54  # gradi
dati$tau.alpha <- 1/(2^2)  # sigma = 2 gradi  -> tau=1/sigma^2 
dati$tau.beta  <- 1/(2^2)  #   idem
dati$tau.gamma <- 1/(2^2)  #   idem

# usa rjags -----------------------------------------------------
jm <- jags.model("triangolo.bug", dati)

update(jm, 1000)

samples <- coda.samples(jm, c("alpha","beta","gamma"), n.iter=10000)

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

# cambiamo formato
samples.df <- as.data.frame(as.mcmc(samples))
# correlazioni
print(cor(samples.df))

# etc. etc. a piacere