vignettes/simulations/runSimsOsensei.R

In evanchildress/perform:

# library(dplyr)
# library(jagsUI)
# library(data.table)
# library(reshape2)
#
# load("~/perform/wbTemps.RData")
# load("~/perform/wbLengths.RData")
fullResults<-NULL
library(perform)

#grab the simulation number from the environment passed from the slurm call
# simNum<- Sys.getenv('SLURM_ARRAY_TASK_ID') %>% as.numeric()
for(s in 1:10){
simNum<-s
#set random seed based on simNum because all iterations were returning identical results
set.seed(simNum)
#
# source("perform/pSim.R")
# source("perform/fitModel.R")
# source("perform/predictPerformance.R")

results<-NULL
iter<-1
for(opt in c(10,15,20)){
 for(tMax in c(16,20,24,27)){
   for(r in c("wb jimmy")){
     for(e in c(0.00075,0.0015,0.003)){
        for(s in c(T)){
          if(tMax<=opt){next}
          print(paste0("starting simulation ",simNum,
                       " iteration ",iter," of 30 at ",
                       as.POSIXct(Sys.time())))
          re<-pSimGrowth(tOpt=opt,ctMax=tMax,sigma=4,eps=e,seasonal=s,
               nYoy=50,river=r)
          results<-rbind(results,re)
	  iter<-iter+1
        }
      }
    }
  }
}
fullResults<-rbind(results,fullResults)
}
# saveRDS(results,paste0("~/perform/output/pSim",simNum,".rds"))
# cat("sim ",simNum," complete")
saveRDS(fullResults,"simResult.rds")