library(tidyverse)
library(fixest)
CAP <- c(308.7*2+360,
	 453*2+400,
	 612,
	 454.8+323,
	323,
	 638.4*2,
	 477+341.7,
	 4678,
	 819.4)

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 
DATA$NAMES <- NAMES

colnames(DATA ) <- c("Capacity","Emp","Peak_Const","Report_Year","Plant")
DATA
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)