For each species, the following order to run the code should be followed:
i. The first code to run is the pbrasiliensis_1_dataGeneration.qmd html. Four new files will be created:
PA_species_time1_blobs and PA_species_time2_blobs).PO_species_time1_raster and PO_species_time2_raster).ii. After, run pbrasiliensis_2_variableSelection.qmd html to get the best predictors for the modeled species. The predictors will be: cov.1, cov.2, cov.3, and cov.4.
iii. Then run the pbrasiliensis_3_dataPreparation.qmd html code with the selected predictors to prepare the data for the model. Four new files will be created:
PO_pbrasiliensis_time1.model.data and PO_pbrasiliensis_time2.model.data).PA_pbrasiliensis_time1.model.data and PA_pbrasiliensis_time2.model.data).iv. Run the model with pbrasiliensis_4_modelRun.qmd html. Two new files will be created:
pbrasiliensis_model.txt)pbrasiliensis_model.rds) - caution this is a >5GB filev. And, finally, create the output figures with pbrasiliensis_5_modelOutputs.qmd html. The html output will show all the results, plus a new file will be created:
docs/pbrasiliensis_model_diagnostics.pdf).