analysis/01_model_runs.R

In mrc-ide/rtss_impact_cea_2021: Update of analysis from Penny et al 2015

library(dplyr)
library(ggplot2)
library(tidyr)
library(patchwork)

source("R/set_up.R")
source("R/run.R")
source("R/post_processing.R")

# Specify suite of runs:
model_runs <- expand_grid(pfpr = pfpr_levels(), season = "perennial", rtss_coverage = c(0, 0.5, 0.8, 0.9), draw = 1:50)
dim(model_runs)

# Run in parallel (locally)
library(furrr)
future::plan(multisession)
model_output_raw <- future_pmap(model_runs, model_run, .options = furrr_options(seed = TRUE))
saveRDS(model_output_raw, "analysis/data/derived_data/model_output_raw.RDS")