output/script/example001.R

In Chang-Yu-Chang/MigrationCommunity: Microbial Consumer-Resource Model

library(tidyverse)
library(deSolve)
library(data.table)
library(micro.crm)

# Make regional species pool
parameter_load(parameter_text = "example001.txt")
pool_build(save_data = TRUE, data_name = "species_pool.Rdata")

# Make resource environment
resource <- setNames(c(1, rep(0, P-1)), paste0("R", sprintf("%03d", 1:P)))
assign("Rin", c(Rin1, rep(0, times = P-1)), envir = .GlobalEnv) # Supply for each resource. Set to 0

# Run
for (r in 1:20) {
  # Simulation
  result <- community_generate(pool = regional_pool, I = I, threshold = threshold, seed = r) %>%
    c(resource, .) %>%
    CR_model(time_limit = time_limit, time_step = 1)
  
  # Write result file to txt
  parameter_write(result = result, 
                  parameter_list = parameter_list, 
                  file_name = paste0("simulation_result/example001_", sprintf("%03d", r),".txt"))
  
  # Report progress
  print(r)
}