From 2a666eefda8dd51973d7cbda944b4b2c793b0369 Mon Sep 17 00:00:00 2001 From: alex Date: Thu, 27 Aug 2026 12:40:03 -0600 Subject: [PATCH] Updating data --- IMPLAN_Visuals.r | 26 ++++++++++++++++++++------ 1 file changed, 20 insertions(+), 6 deletions(-) diff --git a/IMPLAN_Visuals.r b/IMPLAN_Visuals.r index 3c535bd..12e71a7 100644 --- a/IMPLAN_Visuals.r +++ b/IMPLAN_Visuals.r @@ -1,8 +1,22 @@ library(tidyverse) -GET_FILE <- function(YEAR){ -read_csv(paste0("Model_Outputs/IMPLAN/US/US_",YEAR,"/economic_indicators_by_impact.csv")) %>% mutate(year=YEAR) %>% filter(!is.na(Impact)) -} -DATA <- lapply(2028:2035,GET_FILE) %>% bind_rows()%>% rename(Income=`Labor Income`) %>% mutate(Income=parse_number(Income),`Value Added`=parse_number(`Value Added`),Output=parse_number(Output)) -DATA -DATA$Impact <- gsub("3 - ","",gsub("2 - ","",gsub("1 - ","",DATA$Impact))) +GET_DATA <- function(){ + GET_FILE_US <- function(YEAR){read_csv(paste0("Model_Outputs/IMPLAN/US/US_",YEAR,"/economic_indicators_by_impact.csv")) %>% mutate(year=YEAR) %>% filter(!is.na(Impact))%>% mutate(Region="US") } + GET_FILE_WY <- function(YEAR){read_csv(paste0("Model_Outputs/IMPLAN/Wyoming/WY_",YEAR,"/economic_indicators_by_impact.csv")) %>% mutate(year=YEAR) %>% filter(!is.na(Impact))%>% mutate(Region="WY") } +DATA <-rbind(lapply(2028:2035,GET_FILE_US) %>% bind_rows()%>% rename(Income=`Labor Income`) %>% mutate(Income=parse_number(Income),`Value Added`=parse_number(`Value Added`),Output=parse_number(Output)) , +lapply(2028:2035,GET_FILE_WY) %>% bind_rows()%>% rename(Income=`Labor Income`) %>% mutate(Income=parse_number(Income),`Value Added`=parse_number(`Value Added`),Output=parse_number(Output))) + + DATA$Impact <- gsub("3 - ","",gsub("2 - ","",gsub("1 - ","",DATA$Impact))) + DATA <- DATA %>% pivot_longer(-c(Impact,year,Region)) %>% rename(Type=Impact,Impact=name) + DATA$Impact <- gsub("3 - ","",gsub("2 - ","",gsub("1 - ","",DATA$Impact))) +return(DATA) +} + +DATA <- GET_DATA() +EMP_WY <- DATA %>% filter(Impact=="Employment",Region=="WY") +DOL_WY <- DATA %>% filter(Impact!="Employment",Region=="WY") + +ggplot(EMP_WY, aes(x = year, y = value, fill = Type)) + geom_bar(stat = "identity", position = "stack") + theme_minimal() +ggplot(DOL_WY, aes(x = year, y = value/10^9, fill = Type)) + geom_bar(stat = "identity", position = "stack") + theme_minimal()+facet_wrap(~Impact) + +