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

Fixed bug

Browse Source
This commit is contained in:
Alex Gebben Work 2026-10-05 16:53:33 -06:00
parent 6edbfa0df7
commit ec38b9b1bf
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
Data_Center_Visuals.r
View File

@ -120,24 +120,26 @@ TAX_PLOT <- ggplot(IMPLAN_TAX_FIG_DATA, aes(x=Year,y=value/10^6,fill=Impact,grou
#############Dollar values table
NPV <- IMPLAN_DOLLAR %>% mutate(value=value/10^9,Discount_0=value,Discount_2=value/(1+0.02)^(Year-2026),Discount_5=value/(1+0.05)^(Year-2026),Discount_10=value/(1+0.1)^(Year-2026)) %>% group_by(Impact,Region) %>% summarize(Discount_0=sum(Discount_0),Discount_2=sum(Discount_2),Discount_5=sum(Discount_5),Discount_10=sum(Discount_10)) %>% ungroup
NPV[,3:6] <-round(NPV[,3:6],2)
NPV
WY_NPV <- cbind(c("Discount","0%","2%","5%","10%"),NPV %>% filter(Region=='Wyoming') %>% select(-Region) %>% t) %>% as_tibble
colnames(WY_NPV ) <- WY_NPV[1,]
WY_NPV <- WY_NPV[-1,]
WY_NPV <- WY_NPV %>% mutate_if(is.character, str_trim)
WY_NPV[,-1] <- lapply( lapply(WY_NPV[,-1] ,as.numeric) %>% as_tibble,dollar) %>% as_tibble
WY_TAX <- WY_NPV[,c(1,5:7)]
WY_NPV <- WY_NPV[,c(1:4)]
WY_TAX <- WY_NPV %>% select(Discount,Federal,`Wyoming State Taxes`,`Wyoming County Taxes`,`Sub County General`)
WY_NPV <- WY_NPV %>% select(Discount,,Output,`Value Added`,Income)
#US
US_NPV <- cbind(c("Discount","0%","2%","5%","10%"),NPV %>% filter(Region!='Wyoming') %>% select(-Region) %>% t) %>% as_tibble
colnames(US_NPV ) <- US_NPV[1,]
US_NPV <- US_NPV[-1,]
US_NPV <- US_NPV %>% mutate_if(is.character, str_trim)
US_NPV[,-1] <- lapply( lapply(US_NPV[,-1] ,as.numeric) %>% as_tibble,dollar) %>% as_tibble
US_TAX <- US_NPV[,c(1,5:6)]
US_NPV <- US_NPV[,c(1:4)]
US_TAX <- US_NPV %>% select(Discount,Federal,`Other States`,`Sub County General`)
US_NPV <- US_NPV %>% select(Discount,Output,`Value Added`,Income)
ALL_TAX <- US_TAX
ALL_TAX[,2] <- dollar(parse_number(as.character(US_TAX[,2] %>% t))+ parse_number(as.character(WY_TAX[,2] %>% t)))
ALL_TAX <- cbind(ALL_TAX,WY_TAX[,-1:-2])
@ -227,7 +229,6 @@ RANK <- EMP_SUMMARY %>% select(Occupation,Period,Rank) %>% ungroup
PEAK_WY <- RES %>% mutate(Period=ifelse(Year>=2030,"Operations","Construction")) %>% left_join(RANK) %>% mutate(Occupation=ifelse(Rank==16,"Other",Occupation)) %>% group_by(Occupation,Year,Period,Rank) %>% summarize(Peak_Employment=max(Employment,na.rm=TRUE)) %>% group_by(Period,Occupation,Rank) %>% summarize(Wy_Peak_Employment=max(Peak_Employment,na.rm=TRUE)) %>% ungroup %>% select(Occupation,Period,Wy_Peak_Employment)
PEAK_US <- RES_US %>% select(Occupation,Year,US_Employment=Employment) %>% left_join(RES %>% mutate(Wy_Employment=Employment)) %>% mutate(Employment=US_Employment-Wy_Employment)%>% select(-US_Employment,-Wages) %>% mutate(Period=ifelse(Year>=2030,"Operations","Construction")) %>% left_join(RANK) %>% mutate(Occupation=ifelse(Rank==16,"Other",Occupation)) %>% group_by(Occupation,Year,Period,Rank) %>% summarize(Peak_Employment=max(Employment,na.rm=TRUE)) %>% group_by(Period,Occupation,Rank) %>% summarize(US_Peak_Employment=max(Peak_Employment,na.rm=TRUE)) %>% ungroup %>% select(Occupation,Period,US_Peak_Employment)
rbind(RES_US,RES) %>% group_by(Year) %>% summarize(Employment=max(Employment)) %>% plot
PEAK_ALL <- rbind(RES_US,RES) %>% group_by(Occupation,Year) %>% summarize(Employment=max(Employment)) %>% mutate(Period=ifelse(Year>=2030,"Operations","Construction")) %>% left_join(RANK) %>% mutate(Occupation=ifelse(Rank==16,"Other",Occupation)) %>% group_by(Occupation,Year,Period,Rank) %>% summarize(Peak_Employment=max(Employment,na.rm=TRUE)) %>% group_by(Occupation,Period) %>% summarize(Total_Peak_Employment=max(Peak_Employment)) %>% ungroup
#PEAK_ALL <- rbind(US ,WY) %>% group_by(Occupation,Year,Period) %>% summarize(Employment=sum(Employment,na.rm=TRUE)) %>% left_join(RANK) %>% mutate(Occupation=ifelse(Rank==16,"Other",Occupation)) %>% group_by(Occupation,Period) %>% summarize(Total_Peak_Employment=max(Employment,na.rm=TRUE)) %>% ungroup