alex/Population_Study
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Updated key year from 2029 to 2027

Browse Source
This commit is contained in:
Alex Gebben Work 2025-12-15 12:27:39 -07:00
parent 80b23d6b3a
commit a5b30bfde4
3 changed files with 8 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
2C_Upper_Bound_Result_Analysis.r
View File

@ -20,7 +20,7 @@ png(paste0(SAVE_RES_LOC,"Population_Fan_Chart_Upper_Bound.png"), width = 12, hei
POP_PLOT POP_PLOT
dev.off() dev.off()
#Population Summary table for the key years #Population Summary table for the key years
KEY_YEARS <- c(2029,2030,2035,2045,2055,2065) KEY_YEARS <- c(2027,2030,2035,2045,2055,2065)
KEY_YEAR_SUMMARY <- POP_DATA %>% filter(Year %in% KEY_YEARS,Group==5) %>% select(-Group) %>% rename(CI_95_Lower=MIN,CI_95_Upper=MAX) %>% left_join(POP_DATA %>% filter(Year %in% KEY_YEARS,Group==0) %>% select(-Group) %>% select(Year,Median=MIN)) %>% select(Year,CI_95_Lower,Median,CI_95_Upper) KEY_YEAR_SUMMARY <- POP_DATA %>% filter(Year %in% KEY_YEARS,Group==5) %>% select(-Group) %>% rename(CI_95_Lower=MIN,CI_95_Upper=MAX) %>% left_join(POP_DATA %>% filter(Year %in% KEY_YEARS,Group==0) %>% select(-Group) %>% select(Year,Median=MIN)) %>% select(Year,CI_95_Lower,Median,CI_95_Upper)
KEY_YEAR_SUMMARY <- round(KEY_YEAR_SUMMARY ) KEY_YEAR_SUMMARY <- round(KEY_YEAR_SUMMARY )
KEY_YEAR_SUMMARY <- t(KEY_YEAR_SUMMARY ) KEY_YEAR_SUMMARY <- t(KEY_YEAR_SUMMARY )

2
2D_Lower_Bound_Result_Analysis.r
View File

@ -21,7 +21,7 @@ POP_PLOT
dev.off() dev.off()
#Population Summary table for the key years #Population Summary table for the key years
KEY_YEARS <- c(2029,2030,2035,2045,2055,2065) KEY_YEARS <- c(2027,2030,2035,2045,2055,2065)
KEY_YEAR_SUMMARY <- POP_DATA %>% filter(Year %in% KEY_YEARS,Group==5) %>% select(-Group) %>% rename(CI_95_Lower=MIN,CI_95_Upper=MAX) %>% left_join(POP_DATA %>% filter(Year %in% KEY_YEARS,Group==0) %>% select(-Group) %>% select(Year,Median=MIN)) %>% select(Year,CI_95_Lower,Median,CI_95_Upper) KEY_YEAR_SUMMARY <- POP_DATA %>% filter(Year %in% KEY_YEARS,Group==5) %>% select(-Group) %>% rename(CI_95_Lower=MIN,CI_95_Upper=MAX) %>% left_join(POP_DATA %>% filter(Year %in% KEY_YEARS,Group==0) %>% select(-Group) %>% select(Year,Median=MIN)) %>% select(Year,CI_95_Lower,Median,CI_95_Upper)
KEY_YEAR_SUMMARY <- round(KEY_YEAR_SUMMARY ) KEY_YEAR_SUMMARY <- round(KEY_YEAR_SUMMARY )
KEY_YEAR_SUMMARY <- t(KEY_YEAR_SUMMARY ) KEY_YEAR_SUMMARY <- t(KEY_YEAR_SUMMARY )

10
Check_Simulation_Counts.r
View File

@ -1,5 +1,7 @@
library(tidyverse) library(tidyverse)
read_csv("Kemmerer_2016_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year)) %>% pull(Sim_UUID) %>% unique %>% length()/10^5 PRELIM <- read_csv("./Results/Simulations/Kemmerer_2016_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year),is.numeric(Population)) %>% pull(Sim_UUID) %>% unique %>% length()
read_csv("Kemmerer_2024_Upper_Bound_With_Data_Center_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year)) %>% pull(Sim_UUID) %>% unique %>% length()/10^5 UPPER <- read_csv("./Results/Simulations/Kemmerer_2024_Upper_Bound_With_Data_Center_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year),!is.na(Year),is.numeric(Population)) %>% pull(Sim_UUID) %>% unique %>% length()
read_csv("Kemmerer_2024_With_Lower_Bound_With_Data_Center_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year)) %>% pull(Sim_UUID) %>% unique %>% length()/10^5 LOWER <- read_csv("./Results/Simulations/Kemmerer_2024_Lower_Bound_With_Data_Center_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year),is.numeric(Population)) %>% pull(Sim_UUID) %>% unique %>% length()
read_csv("Kemmerer_2024_With_Data_Center_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year)) %>% pull(Sim_UUID) %>% unique %>% length()/10^5 MAIN <- read_csv("./Results/Simulations/Kemmerer_2024_With_Data_Center_Simulation.csv") %>% filter(!is.na(Population),!is.na(Sim_UUID),!is.na(Year),!is.na(Year),is.numeric(Population)) %>% pull(Sim_UUID) %>% unique %>% length()
rbind(PRELIM,UPPER,LOWER,MAIN)/10^5