rm(list = ls())
library(R2jags)

setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
working_dir <- "c:\\E-drive\\Books\\LOGLIN3\\BAYES\\"
ORing_data <- read.csv("O-ring_data.csv", header=T)

iterates <- 10000
burn_in <- 1000

n <- dim(ORing_data)[1]
y <- ORing_data$Failure
temp <- ORing_data$Temperature
# The next two commands are only needed
# for running the centered model
centered_temp <- scale(temp, scale=F)
c.temp <- centered_temp[,1]

a <- c( 1, .577 )
b <- c( .577, 1 )
inits <- function() {
  list( ptilde = c( 0.5, 0.5 ) )
}

sim.dat.jags<-list("y","n","a","b","temp")
ORing.params<-c("beta")
ORing.fit <- jags(data = sim.dat.jags, inits = inits,
  parameters.to.save = ORing.params, n.chains = 1, n.thin=1,
n.iter = iterates+burn_in, n.burnin = burn_in,
model.file = "O-ring_model_a.txt")
beta <- t(ORing.fit$BUGSoutput$sims.list$beta)
#save(beta,file="post-samp.Rda")
summary(ORing.fit$BUGSoutput$sims.list$beta)