#----------------------------------------------
# linear fit including inference on sigma and extrapolation
# (same as done in OpenBUGS)
#
# GdA 27 marzo 2015
#----------------------------------------------
# function pausa
pausa <- function() {
  cat ("\n >> Guarda il plot e dai enter per continuare\n")
  scan()
}

library(rjags)

data <- source("fit_data_OB.txt")$value
inits <- source("fit_inits_OB.txt")$value  
model <- "fit_model.txt"    

print(data)
print(inits)
print(model)

jm <- jags.model(model, data, inits) 

update(jm, 100)                  # burn in 

to.sample <- c('m', 'c', 'sigma', 'muy1', 'y1')      # variables of interest
catena <- coda.samples(jm, to.sample, n.iter=20000)  # sampling 

print(summary(catena))

pausa()

# transform the chain in a dataframe
catena.df <- as.data.frame( as.mcmc(catena) )

hist(catena.df$m, nc=100, col='blue')
# etc. etc.

pausa()

plot(catena.df$m, ty='l', col='blue')

pausa()

plot(catena.df)