analysis/1_run_prioritization_exp_1.R

In LorenzoCattarino/priorityeffort: Prioritization of conservation management effort.

# set up context --------------------------------------------------------------


my_pkgs <- c("priorityeffort")

root <- "context"

context::context_log_start()
ctx <- context::context_save(path = root,
                             packages = my_pkgs)


# define parameters -----------------------------------------------------------


parameters <- list(
  Exp = 1,
  Replicates = 10,
  no_ITER = 1000000,
  Temp_zero = 1,
  cooling_factor = 0.99999,
  fixed_targets = FALSE,
  occurrence_limits = c(500, 10000),
  target_limits = c(1, 0.1),
  TargetLevel = c(seq(50, 1000, 50), seq(1250, 10000, 250)),
  spf = 10,
  print_every_iter = 5,
  start_prop = 0.2)


# load context ----------------------------------------------------------------


context::context_load(ctx)


# create combinations of factor -----------------------------------------------


exp_des <- create_exp_des(parameters)


# run -------------------------------------------------------------------------


solution <- wrapper(parms = parameters,
                    exp_des = exp_des,
                    site_threat_array = site_threat_array,
                    planning_unit = planning_unit,
                    cons_feat_array = cons_feat_array,
                    site_species_array = site_species_array,
                    species_responses = species_responses,
                    parallel = TRUE)


# plot ------------------------------------------------------------------------


# temperature, cost and species penalty
lapply(solution, plot_trackers)

# map of priority effort
lapply(solution, plot_effort_map, daly_prj, rivers_prj, outline_prj)


# save solution ---------------------------------------------------------------


write_out_rds(solution,
              file.path("output", paste0("exp_", parameters$Exp)),
              "solution.rds")