library(tidyverse)
library(scales)
library(janitor)
library(paletteer)
source("Scripts/Return_Data_Center_Model_Inputs.r")
DATA_CENTER <- DATA_CENTER_INPUTS()
COMP_PURCHASE <- read_csv("Model_Outputs/IMPLAN/Data_Center/Wyoming/Computer_Purchase_Mill_Dollars/economic_indicators_by_impact.csv") %>% rename(Income=`Labor Income`)  %>% mutate(Income=parse_number(Income),`Value Added`=parse_number(`Value Added`),Output=parse_number(Output))  %>% mutate(Region='Wyoming',Event='COMPUTERS')
	COMP_PURCHASE$Impact <- gsub("3 - ","",gsub("2 - ","",gsub("1 - ","",COMP_PURCHASE$Impact)))
	COMP_PURCHASE <-  COMP_PURCHASE %>% pivot_longer(-c(Impact,Region,Event)) %>% rename(Type=Impact,Impact=name)
	COMP_PURCHASE$Impact <- gsub("3 - ","",gsub("2 - ","",gsub("1 - ","",COMP_PURCHASE$Impact)))
	COMP_PURCHASE <- COMP_PURCHASE %>% mutate(Source='IMPLAN') %>% select(Type,Region,Impact,Event,Source,value)
	COMP_PURCHASE <- COMP_PURCHASE %>% filter(!is.na(Type))
	COMP_PURCHASE$value <- COMP_PURCHASE$value/10^6 #shift to impact per dollar from impact per million dollar
COMP_PURCHASE <-	COMP_PURCHASE %>% rename('marginal'=value)
COMP_PURCHASE <- COMP_PURCHASE   %>% left_join(DATA_CENTER) %>% mutate(value=marginal*total)

COMP_PURCHASE <- COMP_PURCHASE %>% group_by(Type,Year,Region,Impact,Source) %>% summarize(value=sum(value)) %>% ungroup
COMP_PURCHASE 
EMP <-COMP_PURCHASE %>% filter(Impact=='Employment')
EMP 
ggplot(EMP ,aes(x=Year,y=value,fill=Type)) +geom_area()