2.1_Run_Fit.R

Run_Fit <- function(model,variable.names,n.iter,...){
  jags.out <- coda.samples(model,
                      variable.names = variable.names,
                      n.iter = n.iter,
                      ...)
  
  ## split output
  out <- list(params = NULL, predict = NULL)
  mfit <- as.matrix(jags.out, chains = TRUE)
  pred.cols <- grep("x[", colnames(mfit), fixed = TRUE)
  chain.col <- which(colnames(mfit) == "CHAIN")
  out$predict <- ecoforecastR::mat2mcmc.list(mfit[, c(chain.col, pred.cols)])
  out$params <- ecoforecastR::mat2mcmc.list(mfit[, -pred.cols])
  
  return(out)
}