Sdid For Hourly Wage per Quartal

SDID

CountyWage <- read_excel("Merge County Hourly Earning.xlsx")
CountyWage$St<-substr(x = CountyWage$AreaCode,start = 1,stop = 2)
CountyWage$Treat<-ifelse(CountyWage$St==22& 
                         CountyWage$Date>="2005-07-01",1,0)

str(CountyWage)

## tibble [1,861,306 × 11] (S3: tbl_df/tbl/data.frame)
##  $ AreaCode           : chr [1:1861306] "01001" "01001" "01001" "01001" ...
##  $ Year               : chr [1:1861306] "2000" "2000" "2000" "2000" ...
##  $ Qtr                : chr [1:1861306] "1" "1" "1" "1" ...
##  $ Area               : chr [1:1861306] "Autauga County, Alabama" "Autauga County, Alabama" "Autauga County, Alabama" "Autauga County, Alabama" ...
##  $ Industry           : chr [1:1861306] "Goods-producing" "Natural resources and mining" "Construction" "Manufacturing" ...
##  $ Average Weekly Wage: num [1:1861306] 709 394 408 835 343 390 533 528 472 393 ...
##  $ AverageHourWage    : num [1:1861306] 17.73 9.85 10.2 20.88 8.57 ...
##  $ Qm                 : num [1:1861306] 1 1 1 1 1 1 1 1 1 1 ...
##  $ Date               : POSIXct[1:1861306], format: "2000-01-01" "2000-01-01" ...
##  $ St                 : chr [1:1861306] "01" "01" "01" "01" ...
##  $ Treat              : num [1:1861306] 0 0 0 0 0 0 0 0 0 0 ...

#Manufcaturing
sectors<-unique(CountyWage$Industry)

dumy_date=seq.Date(from = as.Date("2000-01-01"),to = as.Date("2008-01-01"),by = "quarter")
Fix_df<-data.frame(AreaCode=rep(unique(CountyWage$AreaCode),length(dumy_date)))
Fix_df<-arrange(Fix_df,AreaCode)
Fix_df$Date<-rep(dumy_date,length(unique(CountyWage$AreaCode)))

for (i in 1:(length(sectors)-1)){
cf<-CountyWage%>%
  filter(Industry==sectors[i])%>%
  select(AreaCode,Date,AverageHourWage,Treat)

cf<-left_join(Fix_df,cf,by=c("AreaCode","Date"))


cf<-cf%>%
  group_by(AreaCode)%>%
  filter(!any(is.na(AverageHourWage)))

cf<-data.frame(cf)
cf$Date<-as.Date(cf$Date)

pn<-panel.matrices(cf,unit = "AreaCode",time = "Date",
                   outcome = "AverageHourWage",treatment = "Treat")

sdid<-synthdid_estimate(pn$Y, pn$N0, pn$T0)

grph<-plot(sdid)+
  ggtitle(label = paste0(sectors[i]))+
  theme(legend.position="bottom", legend.direction='horizontal',
        legend.key=element_blank(), legend.background=element_blank(), 
        plot.title = element_text(hjust = 0.5,size = 9),
        axis.text = element_text(size=7),axis.title = element_text(size=7))

print(grph)
}