diff --git a/Analysis/R/postprocess_results.R b/Analysis/R/postprocess_results.R index e82d95192..6520a8fc3 100644 --- a/Analysis/R/postprocess_results.R +++ b/Analysis/R/postprocess_results.R @@ -486,10 +486,10 @@ mai_grid_cases_draws <- run_all( # G. Risk categories --------------------------------------------------------- # Get the risk category by location at all admin levels -risk_categories <- run_all( +risk_categories_95 <- run_all( config_dir = opt$config_dir, fun = postprocess_risk_category, - fun_name = "risk_categories", + fun_name = "risk_categories_95", fun_opts = list(cum_prob_thresh = 0.95), prefix = opt$prefix, suffix = opt$suffix, @@ -503,6 +503,25 @@ risk_categories <- run_all( output_file_type = "rds", verbose = opt$verbose) + +# Get the risk category by location at all admin levels +risk_categories_50 <- run_all( + config_dir = opt$config_dir, + fun = postprocess_risk_category, + fun_name = "risk_categories_50", + fun_opts = list(cum_prob_thresh = 0.50), + prefix = opt$prefix, + suffix = opt$suffix, + error_handling = opt$error_handling, + redo = opt$redo, + redo_interm = opt$redo_interm, + redo_aux = opt$redo_auxilliary, + output_dir = opt$output_dir, + interm_dir = opt$interm_dir, + data_dir = opt$data_dir, + output_file_type = "rds", + verbose = opt$verbose) + # Get the population at risk in each risk category by country pop_at_risk <- run_all( config_dir = opt$config_dir,