nsims <- 11000
burnin <- 1000
#nsims <- 60000
#burnin <- 10000
model.file <- "pa_linear.bug"

require(R2WinBUGS)
wbugs_dir <- "/path/to/WinBUGS14"
wine <- "/opt/local/bin/wine"
winepath <- "/opt/local/bin/winepath"

# the data
source("PA.bugs.data")
# initial values
inits <- function() { list(hca=0.1, off=rep(0.2,data$P), def=rep(0.2,data$P), prec=0.3) }
# parameters to keep track of
parameters <- c("hca")
for (i in 1:data$P) { parameters <- c(parameters,paste("off[",i,"]",sep="")) }
for (i in 1:data$P) { parameters <- c(parameters,paste("def[",i,"]",sep="")) }

# run the sims
model.sim <- bugs(data,inits,parameters,model.file,n.iter=nsims,n.burnin=burnin,bugs.directory=wbugs_dir,working.directory=".",useWINE=TRUE,WINE=wine,WINEPATH=winepath,program="winbugs")
print(model.sim)