skateZESSdecZ <- function(n.iter =15000, n.burnin = 5000, n.chains=2, n.thin = 30, debug = T){
require(R2WinBUGS)

WINE <-'/usr/bin/winepath'
WINEPATH <- '/usr/bin/winepath'
bugs.seed <- 14324

data1 <- read.csv('/svn/non-linear-modeling/projects/skate/DATA/EasternScotianShelf_WinterSkate_Zdata-1.csv', header=T)

chain1 <- list(q=rep(.9,3),sdobs=rep(.65,3),sdpro=rep(1,3),theta=c(0.2,0.25),recrate=2,N69=c(1,0.3),Nad=rep(4,5),N= (t(data1[,-1])))
chain2 <- list(q=c(0.8,1,1),sdobs=rep(1,3),sdpro=rep(0.4,3),theta=c(0.25,0.20),recrate=5,N69=c(1,0.5),Nad=rep(3,5),N= (t(data1[,-1])))
inits <- list(chain1,chain2)
           
bugs.file <- "/svn/non-linear-modeling/projects/skate/BUGS/winbugs/bugs/skateZESSdecZ.bug"
wk.dir <- '/svn/non-linear-modeling/projects/skate/BUGS/winbugs/junk/'

data <-as.matrix(data1[,-1])
data <- list(Iobs=t(data))

parameters <- c('q','sdobs','sdpro','theta','recrate','N69','Nad','N','z')

browser()
 tmp <- bugs(data =data, inits, parameters.to.save = parameters, model.file = bugs.file, n.chains = n.chains,n.iter = n.iter, n.burnin = n.burnin, n.thin = n.thin, bugs.directory =  "/home/tdavies/.wine/drive_c/Program Files/WinBUGS14/", working.directory = wk.dir, useWINE = T, newWINE = T, WINEPATH = WINEPATH, clearWD=T,debug = debug,bugs.seed=bugs.seed)

assign('skateZESSdecZ.out', tmp, pos = 1, immediate = T)
}