From 6187a3bffa91fdf6bee42f2a63e1e52c01a87e1e Mon Sep 17 00:00:00 2001
From: Alex Gebben Work <agebben@uwyo.edu>
Date: Wed, 4 Mar 2026 17:20:31 -0700
Subject: [PATCH] Added more data, made graph and reg

---
 Employment_Est.r |  4 +++
 Scale.r          | 66 ++++++++++++++++++++++++++++++++++++++----------
 2 files changed, 56 insertions(+), 14 deletions(-)
 create mode 100644 Employment_Est.r

diff --git a/Employment_Est.r b/Employment_Est.r
new file mode 100644
index 0000000..a97d423
--- /dev/null
+++ b/Employment_Est.r
@@ -0,0 +1,4 @@
+AVG_VALUE <- (0.15 +0.20)/2
+YEARLY_CAPACITY <- c(256,750,1000,1500)
+NUM_OPERATING_STAFF <- round(AVG_VALUE*YEARLY_CAPACITY)
+NUM_OPERATING_STAFF 
diff --git a/Scale.r b/Scale.r
index 1af8aa5..45615b9 100644
--- a/Scale.r
+++ b/Scale.r
@@ -1,5 +1,7 @@
 library(tidyverse)
 library(fixest)
+
+NAMES <- c("Trumbull Energy Center","Orange County Advanced Power Station","Shady Hills Combined Cycle Facility","Cumberland (TN)","Kingston","Lincoln Land Energy Center","Delta Blues Advanced Power Station","Homer","Vicksburg","Cheyenne Prairie","Legend and Lone Star","Wolf Summit Energy","Franklin Farms","Viola Generating Station","CPV Basin Ranch Energy Center","Smarr Combined Cycle Energy Facility","Waterford 5 & 6","Lake Charles","Jefferson Power Station")
 CAP <- c(308.7*2+360,
 	 453*2+400,
 	 612,
@@ -8,23 +10,59 @@ CAP <- c(308.7*2+360,
 	 638.4*2,
 	 477+341.7,
 	 4678,
-	 819.4)
+	 819.4,
+	 132,
+	754+453,
+	600,
+	1560,
+	705,
+	1490,
+	1425,
+	1640,
+	994,
+	754
+)
 
-NAMES <- c("Trumbull Energy Center","Orange County Advanced Power Station","Shady Hills Combined Cycle Facility","Cumberland (TN)","Kingston","Lincoln Land Energy Center","Delta Blues Advanced Power Station","Homer","Vicksburg")
-EMP <- c(25,27,12.5,30,NA,34,21,250,21)
-PEAK_CONST <- c(920,NA,600,NA,300,500,300,2500,560)
-DATE <- c(2017,2018,2018,2022,2015,2020,2024,2025,NA)
-DATA <- cbind(CAP,EMP,PEAK_CONST,DATE)  %>% as.matrix %>% as_tibble 
+DETAILED_REPORT <- c(1,1,1,1,0,1,0,0,0,1,0,0,0,0,0,0,0,0,0)
+EMP <- c(25,27,12.5,30,NA,34,21,250,21,12,73,30,44,30,45,30,30,30,22)
+PEAK_CONST <- c(920,NA,600,NA,300,500,300,2500,560,400,NA,400,(1500+1800)/2,500,1000,1200,700,1100,NA)
+DATE <- c(2017,2018,2018,2022,2015,2020,2024,2025,2025,2012,2024,2025,2024,2024,2026,2026,2025,2020,2025)
+DATA <- cbind(CAP,EMP,PEAK_CONST,DETAILED_REPORT,DATE)  %>% as.matrix %>% as_tibble 
 DATA$NAMES <- NAMES
+DATA
 
-colnames(DATA ) <- c("Capacity","Emp","Peak_Const","Report_Year","Plant")
-feols(log(Emp)~log(Capacity),DATA %>% filter(Plant!='Homer'))
-feols(log(Emp)~log(Capacity),DATA %>% filter(Plant!='Homer'))
-
-feols(log(Emp)~log(Capacity),DATA)
-
-feols(log(Peak_Const)~log(Capacity)+Report_Year,DATA)
-plot(CAP,PEAK_CONST)
+colnames(DATA ) <- c("Capacity","Emp","Peak_Const","Report_Quality","Report_Year","Plant")
+REG_DATA <-DATA %>% select(Emp,Capacity,Report_Quality,Report_Year)
+REG_DATA_TRIMMED <-DATA %>% filter(Plant!='Homer') %>% select(Emp,Capacity,Report_Quality,Report_Year)
+DATA
+MOD1 <- feols(log(Emp)~Capacity+Report_Quality+Report_Year,REG_DATA,vcov='hetero')
+feols(log(Emp)~Capacity+Report_Quality+Report_Year,REG_DATA,vcov='hetero')
+PRED_DATA <- DATA[1:4,]
+PRED_DATA$Capacity <- c(250,750,1000,1500)
+PRED_DATA$Report_Year <- 2026
+PRED_DATA$Plant <- 'Enyo'
+PRED_DATA$Emp<- NA 
+round(exp(predict(MOD1,newdata=PRED_DATA)))
+#PRED_DATA$Capacity<- c(640,1280,1280,3680) #Larger number includes the single cyle turbines that look to be turned off most of the time
+round(exp(predict(MOD1,newdata=PRED_DATA)))
 
 
 
+
+
+
+ggplot(REG_DATA,aes(x=Capacity,y=log(Emp)))+geom_point()+geom_smooth(method='lm')
+#feols(log(Emp)~Capacity+Report_Quality+Report_Year,REG_DATA_TRIMMED,vcov='hetero')
+
+feols(Peak_Const~Capacity+Report_Year ,DATA,vcov='hetero')
+DATA
+
+feols(Peak_Const~Capacity+Report_Year+Report_Quality,DATA,vcov='hetero')
+
+MOD1 <- feols(log(Emp)~Capacity,REG_DATA,vcov='hetero')
+MOD1
+PRED_DATA <- DATA[1:4,]
+
+
+
+b