REMI_Talk/R_Work/make_figures.r
2026-05-11 17:03:33 -06:00

29 lines
1.5 KiB
R

library(tidyverse)
library(janitor)
library(scales)
library(paletteer)
CONSUMPTION <- read_csv('EIA_Commerical_Sector_Data.csv',skip=4) %>% clean_names()
CONSUMPTION[,c(2:14)] <- 33.43* CONSUMPTION[,c(2:14)]
colnames(CONSUMPTION) <- gsub("_quads","",colnames(CONSUMPTION))
CONSUMPTION <- CONSUMPTION %>% pivot_longer(-year,values_to='GW',names_to="Use")
TOTAL <- CONSUMPTION %>% filter(Use=='total_gross_end_use_consumption')
CONSUMPTION <-CONSUMPTION %>% filter(Use!='total_gross_end_use_consumption')
CONSUMPTION$Use <- ifelse(CONSUMPTION$Use!='data_center_servers',"Other","Data Centers")
CONSUMPTION <- CONSUMPTION %>% group_by(year,Use) %>% summarize(GW=sum(GW)) %>% ungroup
CONSUMPTION$Use <- factor(CONSUMPTION$Use,level=c("Other","Data Centers"))
CONSUMPTION <- CONSUMPTION %>% group_by(Use) %>% mutate(Change=GW-sum(ifelse(year==2025,GW,0))) %>% ungroup
STACK_PLOT <- ggplot(CONSUMPTION,aes(x=year,y=Change,fill=Use)) +geom_area()+theme_bw()+facet_wrap(~Use)+ theme(legend.position = "top",text = element_text(size = 20))+scale_fill_manual(values=paletteer_d("PNWColors::Shuksan2"))+ylab("Change GW-year")
STACK_TOTAL<- ggplot(CONSUMPTION,aes(x=year,y=GW,fill=Use)) +geom_area()+theme_bw()+ theme(legend.position = "top",text = element_text(size = 20))+scale_fill_manual(values=paletteer_d("PNWColors::Shuksan2"))+ylab("Total Consumpiotn (GW-year)")
ggsave("../Plot.png",STACK_PLOT,width=12,height=8,units="in",dpi=600)
ggsave("../Total.png",STACK_TOTAL,width=12,height=8,units="in",dpi=600)