skate_ADMB_mcmc <- function() {  ## This is not quite working yet (at least the read_admn part and I don't know where the MCMC samples are going -- model does run)
require(R2admb)
  file.copy('../DATA/skate.dat', './skate.dat')

#params1 <- c("recrate","theta.1","theta.2","z.1.1","z.2.1","z.3.1","z.1.2","z.2.2","z.3.2","z.1.3","z.2.3","z.3.3")

  run_admb(fn="skate", verbose = TRUE,mcmc=TRUE,mcmc.opts = mcmc.control(mcmc=1000,mcmc2=0,mcsave=10,mcmcpars=('recrate')))

 #run_admb(fn="skate", verbose = TRUE,mcmc=TRUE,mcmc.opts = mcmc.control(mcmc=100,mcmc2=0,mcsave=10,mcmcpars=params1))


#  L.mcmc <- read_admb("skate",mcmc=TRUE,mcmc.opts = mcmc.control(mcmcpars=params1)) 
#  L.mcmc <- read_admb("skate",mcmc=TRUE,mcmc.opts = mcmc.control(mcmcpars='recrate')) 
 L.mcmc <- read_admb("skate",mcmc=TRUE)
#,mcmc.opts = mcmc.control(mcmcpars='recrate')) 

#clean_admb("skate")
  return(L.mcmc)
}