rm(list = ls())
# enter data size
n=300
iterates=80000
burn_in=5000
library(R2OpenBUGS)
BUGS_path <-
"c:\\Program Files (x86)\\OpenBUGS\\OpenBUGS323\\OpenBUGS.exe"
setwd("c:\\E-drive\\Books\\LOGLIN3\\BAYES\\")
working_dir <- "c:\\E-drive\\Books\\LOGLIN3\\BAYES\\"

# Enter data
Trauma <- read.table(
url("http://stat.unm.edu/~fletcher/LLM/DATA/TRAUMAa.DAT"),
#"C:\\E-drive\\Books\\LOGLIN3\\DATA\\TRAUMAa.DAT",
		  sep="",col.names=c("ID","Death","ISS","TI","RTS","AGE"))
attach(Trauma)
Death=1-Death

data <- list( "n", "Death","ISS","RTS","AGE","TI")

inits <- function() {
  list(beta=c(0.5,0.5,0.5,0.5,0.5,0.5))
}

parameters <- list( "beta" )


Trauma.sim <- bugs( data, inits, parameters,
                  model.file="trauma_flat.txt",
                   n.chains=1, n.iter=iterates+burn_in,
                   n.thin=1, n.burnin=burn_in,
                   OpenBUGS.pgm=BUGS_path,
                   working.directory=working_dir,
                   debug=F )
Trauma.sim$summary
# save p x t matrix of gamma posterior iterates
beta <- t(Trauma.sim$sims.list$beta)
save(beta,file="trauma_flat_samp.Rda")