From fce74d225f96fcf4e3f89491c89a673b1d87dd3e Mon Sep 17 00:00:00 2001 From: Alex Gebben Work Date: Tue, 2 Dec 2025 17:23:47 -0700 Subject: [PATCH] Fixed small bug --- 2_Result_Analysis.r | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/2_Result_Analysis.r b/2_Result_Analysis.r index 2b0804f..2a2623c 100644 --- a/2_Result_Analysis.r +++ b/2_Result_Analysis.r @@ -1,7 +1,7 @@ library(tidyverse) - ###Process the simulations and save the main percentile results by year RES <- read_csv("Results/Simulations/Kemmerer_2016_Simulation.csv") + HIST <- readRDS("Data/Cleaned_Data/Population_Data/RDS/Kemmerer_Diamondville_Population_Data.Rds") %>% filter(County=='Lincoln') %>% mutate(Percentile="Actual Population") %>% filter(Year>=1940) YEARS <- min(RES$Year):max(RES$Year) ######Population GRAPH_DATA <- do.call(rbind,lapply(YEARS,function(x){quantile(RES %>% filter(Year==x) %>% pull(Population),c(0.025,0.05,0.1,0.25,0.4,0.5,0.6,0.75,0.9,0.95,0.975))})) %>% as_tibble %>% mutate(Year=YEARS) @@ -25,7 +25,6 @@ library(tidyverse) #Add historic MEDIAN_PRED <- GRAPH_DATA %>% filter(Percentile=='50%') GRAPH_DATA <- GRAPH_DATA %>% filter(Percentile!='50%') - HIST <- readRDS("Data/Cleaned_Data/Population_Data/RDS/Kemmerer_Diamondville_Population_Data.Rds") %>% filter(County=='Lincoln') %>% mutate(Percentile="Actual Population") %>% filter(Year>=1940) ALPHA=0.2 COLOR <- 'black'