scripts/rename_comsie_daisie_input.R

In TheoPannetier/hamr: Utility package for comrad

# DAISIE
args <- tidyr::expand_grid(
  "siga" = c(0.091, 0.369),
  "gamma" = c(1e-04),
  "i" = 1:100
) %>% pwalk(function(siga, gamma, i) {
  cat("siga =", siga, "gamma =", gamma, "i =", i, "\t")
  # Input data
  input_data <- readRDS(
    glue::glue("comsie_data/daisie/input/daisie_input_siga_{siga}_gamma_{gamma}_{i}.rds")
  )
  saveRDS(input_data, file = glue::glue(
    "comsie_data/daisie/input/daisie_input_siga_{siga}_gamma_{gamma}_{i}_f_1.rds"
  ))
  # initpars
  initpars_data <- readRDS(
    glue::glue("comsie_data/daisie/input/daisie_initpars_siga_{siga}_gamma_{gamma}_{i}.rds")
  )
  saveRDS(initpars_data, file = glue::glue(
    "comsie_data/daisie/input/daisie_initpars_siga_{siga}_gamma_{gamma}_{i}_f_1.rds"
  ))
})

# DDD
args <- tidyr::expand_grid(
  "siga" = c(0.091),
  "gamma" = c(0, 1e-04, 5e-04, 0.001, 0.01),
  "i" = 1:100
)
args[263:500,] %>% pwalk(function(siga, gamma, i) {
  cat("siga =", siga, "gamma =", gamma, "i =", i, "\t")
  # Input data
  input_data <- readRDS(
    glue::glue("comsie_data/ddd/input/ddd_input_siga_{siga}_gamma_{gamma}_{i}.rds")
  )
  saveRDS(input_data, file = glue::glue(
    "comsie_data/ddd/input/ddd_input_siga_{siga}_gamma_{gamma}_{i}_f_1.rds"
  ))
  # initpars
  initpars_data <- readRDS(
    glue::glue("comsie_data/ddd/input/ddd_initpars_siga_{siga}_gamma_{gamma}_{i}.rds")
  )
  saveRDS(initpars_data, file = glue::glue(
    "comsie_data/ddd/input/ddd_initpars_siga_{siga}_gamma_{gamma}_{i}_f_1.rds"
  ))
})

args <- tidyr::expand_grid(
  "siga" = c(0.091),
  "gamma" = c(0, 1e-04, 5e-04, 0.001, 0.01),
  "i" = 1:100,
  "params_i" = 1:4
)

args[1049:2000,] %>% pwalk(function(siga, gamma, i, params_i) {
  cat("siga =", siga, "gamma =", gamma, "i =", i, "\t")
  # Input data
  dd_ml <- readRDS(
    glue::glue("comsie_data/ddd/output/ddd_ml_siga_{siga}_gamma_{gamma}_rep_{i}_params_{params_i}.rds")
  )
  saveRDS(dd_ml, file = glue::glue(
    "comsie_data/ddd/output/ddd_ml_siga_{siga}_gamma_{gamma}_rep_{i}_f_1_params_{params_i}.rds"
  ))

  cr_ml <- readRDS(
    glue::glue("comsie_data/ddd/output/cr_ml_siga_{siga}_gamma_{gamma}_rep_{i}_params_{params_i}.rds")
  )
  saveRDS(cr_ml, file = glue::glue(
    "comsie_data/ddd/output/cr_ml_siga_{siga}_gamma_{gamma}_rep_{i}_f_1_params_{params_i}.rds"
  ))

})

#########

args <- tidyr::expand_grid(
  "siga_ibm" = c(0.091, 0.369),
  "gamma_ibm" = c(0.001),
  "replicate" = 1:100,
  "f" = 1
) %>% pwalk(create_comsie_initparsopt_daisie)