CLB <- read_excel("ALLCOUNTY/Civil Labor Ballance.xlsx")
Ballance_CLB<-CLB%>%
group_by(FIPS)%>%
filter(!any(is.na(Civil.Labor)))
Ballance_CLB<-Ballance_CLB%>%
filter(DATE<"2010-09-01")
Ballance_CLB$Treat<-ifelse(Ballance_CLB$DATE>="2005-08-01"&
Ballance_CLB$STATE=="LA",1,0)
Ballance_CLB<-data.frame(Ballance_CLB)
Ballance_CLB$DATE<-as.Date(Ballance_CLB$DATE)
str(Ballance_CLB)
## 'data.frame': 351360 obs. of 10 variables:
## $ CODE : chr "AKALEU0LFN" "AKALEU0LFN" "AKALEU0LFN" "AKALEU0LFN" ...
## $ DATE : Date, format: "2000-01-01" "2000-02-01" ...
## $ STATE : chr "AK" "AK" "AK" "AK" ...
## $ FIPS : num 2016 2016 2016 2016 2016 ...
## $ COUNTY : chr "aleutians west census area" "aleutians west census area" "aleutians west census area" "aleutians west census area" ...
## $ Civil.Labor: num 3186 3713 3511 2845 2246 ...
## $ state_name : chr "Alaska" "Alaska" "Alaska" "Alaska" ...
## $ lat : num 52.8 52.8 52.8 52.8 52.8 ...
## $ lng : num -107 -107 -107 -107 -107 ...
## $ Treat : num 0 0 0 0 0 0 0 0 0 0 ...
Ballance_CLB$STATE<-factor(Ballance_CLB$STATE)
matx<-panel.matrices(Ballance_CLB,unit = 4,time = 2,outcome = 6,treatment = 10)
sydid<-synthdid_estimate(matx$Y, matx$N0, matx$T0)
plot(sydid,overlay=1)+theme(legend.position = "bottom",
legend.direction = "horizontal",
legend.box.background = element_blank())+
ylab("Civil Labor")+xlab("Date")+
scale_x_continuous(labels = c("2002-02-01","2002-11-01",
"2005-08-01","2008-05-01",
"2011-01-01"))
Unemployee Rate
# Import Dataset
UR<-read_xlsx("UnemRate/Merge UnemRate ballance V2 FULL-Location.xlsx")
# Remove NA Value of Unemployee Rate
Ballance_UR<-UR%>%
group_by(FIPS)%>%
filter(!any(is.na(Unemployee.Rate)))
Ballance_UR<-Ballance_UR%>%
filter(DATE<"2010-09-01")
Ballance_UR$Treat<-ifelse(Ballance_UR$DATE>="2005-08-01"&
Ballance_UR$STATE=="LA",1,0)
Ballance_UR<-data.frame(Ballance_UR)
Ballance_UR$DATE<-as.Date(Ballance_UR$DATE)
Ballance_UR$STATE<-factor(Ballance_UR$STATE)
matx<-panel.matrices(Ballance_UR,unit = 6,time = 2,outcome = 5,treatment = 9)
sydid<-synthdid_estimate(matx$Y, matx$N0, matx$T0)
plot(sydid,overlay=1)+theme(legend.position = "bottom",
legend.direction = "horizontal",
legend.box.background = element_blank())+
ylab("Unemployee Rate")+xlab("Date")+
scale_x_continuous(labels = c("2002-02-01","2002-11-01",
"2005-08-01","2008-05-01",
"2011-01-01"))
