R/ccir_stan_run_binomial.r

In AMCOOK/bio.ccir:

#' @export
ccir_stan_run_binomial <- function(dat,inits=c(A0=0.5, B0 = 0.5),save =T,season=NULL, test=NA){

#inital value sampler for stan

   ins <- ccir_inits(dat,A0=inits[1],B0=inits[2])

dat$Ap = as.numeric(ins[1])
dat$Bp = as.numeric(ins[2])
dat$phip = 0
dat$App = 5
dat$Bpp = 5
dat$phipp = 100
#Model choice

    ccir_stan <- ccir_stan_models_binomial()
      pars.to.monitor = c('A','B','ERp')
      #test = stan(model_code = ccir_stan, model_name = 'binomial', pars = pars.to.monitor, fit=test, data = dat, iter = 2000,  warmup=200, thin=1, chains = 1, verbose = F)

      fit = stan(model_code = ccir_stan, model_name = 'bino', fit=test, data = dat, iter = 3500, warmup=200, thin=20, cores = 4, verbose = F)



  fdir = file.path(project.datadirectory('bio.lobster'),'outputs','ccir')
  if(!dir.exists(fdir)) dir.create(fdir,recursive=T,showWarnings=F)
  ti = 1
  if(!is.null(season)) ti = season
  fname = paste('LFA',dat$LFA,'Year',dat$Yr,ti,'Grid',paste(min(dat$Grid),max(dat$Grid),sep='-'),'Season',paste(min(dat$Seas),max(dat$Seas),sep='-'),'Sex',dat$Sex,dat$method,sep='.')
  out = list(fit= fit,dat = dat,file=fname, test=test)
  if(save) save(out,file=file.path(fdir,paste(fname,'rdata',sep='.')))
return(out)
    }