inst/simulation/find_true_population_size.R

In statguy/STREM: Winter track counts

# library(devtools); install_github("statguy/Winter-Track-Counts")

library(parallel)
library(doMC)
registerDoMC(cores=detectCores())
library(STREM)
source("~/git/STREM/setup/WTC-Boot.R")
library(plyr)

scenarios <- c("A","B","C","D","E","F")
allScenarios <- c(scenarios, paste0(scenarios, "combined"), paste0(scenarios, "10days"))
maxIterations <- 50

l_ply(allScenarios, function(scenario) {

  mss <- getMSS(scenario=scenario)
  study <- mss$study
  
  tracks <- SimulatedTracks(study=study)
  iterations <- tracks$getTracksFileIterations()
  iterations <- iterations[iterations<=maxIterations]
  
  for (iteration in iterations) {
    tracks <- study$loadTracks(iteration=iteration, addColumns=FALSE)
    truePopulationSize <- tracks$truePopulationSize
    fileName <- study$context$getLongFileName(dir=tracks$getTracksDirectory(), name="TruePopulationSize", response=study$response, region=study$studyArea$region, tag=iteration)
    message("Saving to ", fileName)
    save(truePopulationSize, file=fileName)
  }

}, .parallel=T)