Script/run_edd_rtree.R

In EvoLandEco/eve: Evolution Emulator, phylogenetically dependent

args <- commandArgs(TRUE)

name <- args[1]
set <- as.numeric(args[2])
size <- as.numeric(args[3])
n <- as.numeric(args[4])

combo1 <- eve::edd_combo_maker(
  la = c(0.4, 0.5, 0.6),
  mu = c(0, 0.1, 0.2),
  beta_n = c(-0.04, -0.02, 0),
  beta_phi = c(-0.04, -0.02, 0, 0.001, 0.002),
  age = c(6),
  model = "dsce2",
  metric = c("pd"),
  offset = c("simtime")
)

combo2 <- eve::edd_combo_maker(
  la = c(0.4, 0.5, 0.6),
  mu = c(0, 0.1, 0.2),
  beta_n = c(-0.04, -0.02, 0),
  beta_phi = c(-0.04, -0.02, 0, 0.001, 0.002),
  age = c(6),
  model = "dsce2",
  metric = c("ed"),
  offset = c("none")
)

combo3 <- eve::edd_combo_maker(
  la = c(0.4, 0.5, 0.6),
  mu = c(0, 0.1, 0.2),
  beta_n = c(-0.04, -0.02, 0),
  beta_phi = c(-0.04, -0.02, 0, 0.001, 0.002),
  age = c(6),
  model = "dsce2",
  metric = c("nnd"),
  offset = c("none")
)

combo <- c(combo1, combo2, combo3)

names(combo) <- as.character(seq(1, length(combo)))

if (!file.exists("params")) {
  params <- eve:::combo_to_tibble(combo)
  write.table(params, "params")
}

out <-
  eve::search_tree(
    pars = combo[[set]]$pars,
    age = combo[[set]]$age,
    model = combo[[set]]$model,
    metric = combo[[set]]$metric,
    offset = combo[[set]]$offset,
    size = size,
    n = n,
    history = FALSE,
    verbose = FALSE
  )

print(paste0("Starting parameter set ", set))

eve:::check_folder(name)
eve:::save_result(out, name, set)