29 lines
1.0 KiB
R
29 lines
1.0 KiB
R
library(tidyverse)
|
|
getwd()
|
|
#Extracted data from EIA see https://www.eia.gov/dnav/ng/ng_cons_sum_a_EPG0_veu_mmcf_m.htm
|
|
Gas_Cons <- read_csv("./Natural_Gas_Consumption_by_End_Use.csv",skip=6)
|
|
colnames(Gas_Cons) <- gsub(' Natural Gas Deliveries to Electric Power Consumers MMcf',"",colnames(Gas_Cons ) )
|
|
colnames(Gas_Cons) <- gsub(' ',"_",colnames(Gas_Cons ) )
|
|
|
|
colnames(Gas_Cons)[1] <- "Obs_Date"
|
|
State_Order <- sort(colnames(Gas_Cons[-1]))
|
|
Gas_Cons <- Gas_Cons[,c("Obs_Date",State_Order)]
|
|
|
|
Months <- substr(Gas_Cons$Obs_Date,1,3)
|
|
Years <- substr(Gas_Cons$Obs_Date,5,8)
|
|
Full_Dates <- paste("01", Gas_Cons$Obs_Date)
|
|
|
|
Gas_Cons$Month<- factor(Months,ordered = FALSE)
|
|
Gas_Cons$Year <- factor(Years,ordered = TRUE)
|
|
Gas_Cons$Obs_Date <- as.Date(Full_Dates, format = "%d %b %Y")
|
|
|
|
|
|
Col_Names <- colnames(Gas_Cons)
|
|
Gas_Cons <- Gas_Cons[,c("Year","Month",Col_Names)]
|
|
Gas_Cons[is.na(Gas_Cons )] <- 0
|
|
|
|
Gas_Cons <- Gas_Cons[Gas_Cons$Yeara>2014,]
|
|
Wyoming_Gas_Cons <- Gas_Cons[,c("Year","Month","Obs_Date","Wyoming")]
|
|
write_csv(Gas_Cons,
|
|
|