Import
RESULT<- read_excel("New RESULT PAER Permutattion.xlsx")
AER_Per<- read_excel("New Employee Gap PAER Permutation.xlsx")
LA_AER.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_AER.Result$MSE_pre[1])
AER_Per$Chan<-AER_Per$PAER-AER_Per$sdid_intecept
LA_AER.i<-filter(AER_Per,STATE=="LA")
hashV<-LA_AER.i$DATE[31]
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.AER<-left_join(STATE.LOW.MSE.15[,1:2],AER_Per,by=c("State"="STATE"))
#SMALL.STATE.AER<-filter(SMALL.STATE.AER,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.AER,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_AER.i,colour="black")+
ggtitle(label = "Art, Entertain and Recreation")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.AER.15
SMALL.LA.AER.15
#@####################################
RESULT<- read_excel("New RESULT PAF Permutattion.xlsx")
AF_Per<- read_excel("New Employee Gap PAF Permutation.xlsx")
LA_AF.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_AF.Result$MSE_pre[1])
AF_Per$Chan<-AF_Per$PAF-AF_Per$sdid_intecept
LA_AF.i<-filter(AF_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.AF<-left_join(STATE.LOW.MSE.15[,1:2],AF_Per,by=c("State"="STATE"))
#SMALL.STATE.AF<-filter(SMALL.STATE.AF,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.AF,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_AF.i,colour="black")+
ggtitle(label = "Accomodation and \n Food Services")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.AF.15
SMALL.LA.AF.15
RESULT<- read_excel("New RESULT PASW Permutattion.xlsx")
ASW_Per<- read_excel("New Employee Gap PASW Permutation.xlsx")
LA_ASW.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_ASW.Result$MSE_pre[1])
ASW_Per$Chan<-ASW_Per$PASW-ASW_Per$sdid_intecept
LA_ASW.i<-filter(ASW_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.ASW<-left_join(STATE.LOW.MSE.15[,1:2],ASW_Per,by=c("State"="STATE"))
#SMALL.STATE.ASW<-filter(SMALL.STATE.ASW,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.ASW,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_ASW.i,colour="black")+
ggtitle(label = "ASW")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.ASW.15
SMALL.LA.ASW.15
RESULT<- read_excel("New RESULT PCON Permutattion.xlsx")
CON_Per<- read_excel("New Employee Gap PCON Permutation.xlsx")
LA_CON.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_CON.Result$MSE_pre[1])
CON_Per$Chan<-CON_Per$PCON-CON_Per$sdid_intecept
LA_CON.i<-filter(CON_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.CON<-left_join(STATE.LOW.MSE.15[,1:2],CON_Per,by=c("State"="STATE"))
#SMALL.STATE.CON<-filter(SMALL.STATE.CON,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.CON,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_CON.i,colour="black")+
ggtitle(label = "Construction")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.CON.15
SMALL.LA.CON.15
RESULT<- read_excel("New RESULT PES Permutattion.xlsx")
ES_Per<- read_excel("New Employee Gap PES Permutation.xlsx")
LA_ES.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_ES.Result$MSE_pre[1])
ES_Per$Chan<-ES_Per$PES-ES_Per$sdid_intecept
LA_ES.i<-filter(ES_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.ES<-left_join(STATE.LOW.MSE.15[,1:2],ES_Per,by=c("State"="STATE"))
#SMALL.STATE.ES<-filter(SMALL.STATE.ES,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.ES,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_ES.i,colour="black")+
ggtitle(label = "Educational Services")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.ES.15
SMALL.LA.ES.15
RESULT<- read_excel("New RESULT PFI Permutattion.xlsx")
FI_Per<- read_excel("New Employee Gap PFI Permutation.xlsx")
LA_FI.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_FI.Result$MSE_pre[1])
FI_Per$Chan<-FI_Per$PFI-FI_Per$sdid_intecept
LA_FI.i<-filter(FI_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.FI<-left_join(STATE.LOW.MSE.15[,1:2],FI_Per,by=c("State"="STATE"))
#SMALL.STATE.FI<-filter(SMALL.STATE.FI,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.FI,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_FI.i,colour="black")+
ggtitle(label = "Financial Insurance")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.FI.15
SMALL.LA.FI.15
RESULT<- read_excel("New RESULT PGOV Permutattion.xlsx")
GOV_Per<- read_excel("New Employee Gap PGOV Permutation.xlsx")
LA_GOV.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_GOV.Result$MSE_pre[1])
GOV_Per$Chan<-GOV_Per$PGOV-GOV_Per$sdid_intecept
LA_GOV.i<-filter(GOV_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.GOV<-left_join(STATE.LOW.MSE.15[,1:2],GOV_Per,by=c("State"="STATE"))
#SMALL.STATE.GOV<-filter(SMALL.STATE.GOV,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.GOV,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_GOV.i,colour="black")+
ggtitle(label = "Government")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.GOV.15
SMALL.LA.GOV.15
RESULT<- read_excel("New RESULT PHCS Permutattion.xlsx")
HCS_Per<- read_excel("New Employee Gap PHCS Permutation.xlsx")
LA_HCS.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_HCS.Result$MSE_pre[1])
HCS_Per$Chan<-HCS_Per$PHCS-HCS_Per$sdid_intecept
LA_HCS.i<-filter(HCS_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.HCS<-left_join(STATE.LOW.MSE.15[,1:2],HCS_Per,by=c("State"="STATE"))
#SMALL.STATE.HCS<-filter(SMALL.STATE.HCS,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.HCS,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_HCS.i,colour="black")+
ggtitle(label = "Health Care and \n Social Assistance")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.HCS.15
SMALL.LA.HCS.15
#########
RESULT<- read_excel("New RESULT PINF Permutattion.xlsx")
INF_Per<- read_excel("New Employee Gap PINF Permutation.xlsx")
LA_INF.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_INF.Result$MSE_pre[1])
INF_Per$Chan<-INF_Per$PINF-INF_Per$sdid_intecept
LA_INF.i<-filter(INF_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.INF<-left_join(STATE.LOW.MSE.15[,1:2],INF_Per,by=c("State"="STATE"))
SMALL.STATE.INF<-filter(SMALL.STATE.INF,State!="TX")
ggplot()+geom_line(data=SMALL.STATE.INF,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_INF.i,colour="black")+
ggtitle(label = "Information")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.INF.15
SMALL.LA.INF.15
RESULT<- read_excel("New RESULT PMF Permutattion.xlsx")
MAN_Per<- read_excel("New Employee Gap PMF Permutation.xlsx")
LA_MAN.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_MAN.Result$MSE_pre[1])
MAN_Per$Chan<-MAN_Per$PMF-MAN_Per$sdid_intecept
LA_MAN.i<-filter(MAN_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.MAN<-left_join(STATE.LOW.MSE.15[,1:2],MAN_Per,by=c("State"="STATE"))
#SMALL.STATE.MAN<-filter(SMALL.STATE.MAN,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.MAN,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_MAN.i,colour="black")+
ggtitle(label = "Manufacturing")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.MAN.15
SMALL.LA.MAN.15
RESULT<- read_excel("New RESULT PMCE Permutattion.xlsx")
MCE_Per<- read_excel("New Employee Gap PMCE Permutation.xlsx")
LA_MCE.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_MCE.Result$MSE_pre[1])
MCE_Per$Chan<-MCE_Per$PMCE-MCE_Per$sdid_intecept
LA_MCE.i<-filter(MCE_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.MCE<-left_join(STATE.LOW.MSE.15[,1:2],MCE_Per,by=c("State"="STATE"))
#SMALL.STATE.MCE<-filter(SMALL.STATE.MCE,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.MCE,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_MCE.i,colour="black")+
ggtitle(label = "Management of Companies \n and Enterprises")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.MCE.15
SMALL.LA.MCE.15
RESULT<- read_excel("New RESULT PRE Permutattion.xlsx")
RE_Per<- read_excel("New Employee Gap PRE Permutation.xlsx")
LA_RE.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_RE.Result$MSE_pre[1])
RE_Per$Chan<-RE_Per$PRE-RE_Per$sdid_intecept
LA_RE.i<-filter(RE_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.RE<-left_join(STATE.LOW.MSE.15[,1:2],RE_Per,by=c("State"="STATE"))
#SMALL.STATE.RE<-filter(SMALL.STATE.RE,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.RE,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_RE.i,colour="black")+
ggtitle(label = "Real Estate and \n Rental and Leasing")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.RE.15
SMALL.LA.RE.15
RESULT<- read_excel("New RESULT PRT Permutattion.xlsx")
RT_Per<- read_excel("New Employee Gap PRT Permutation.xlsx")
LA_RT.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_RT.Result$MSE_pre[1])
RT_Per$Chan<-RT_Per$PRT-RT_Per$sdid_intecept
LA_RT.i<-filter(RT_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.RT<-left_join(STATE.LOW.MSE.15[,1:2],RT_Per,by=c("State"="STATE"))
#SMALL.STATE.RT<-filter(SMALL.STATE.RT,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.RT,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_RT.i,colour="black")+
ggtitle(label = "Retail Trade")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.RT.15
SMALL.LA.RT.15
RESULT<- read_excel("New RESULT POS Permutattion.xlsx")
OS_Per<- read_excel("New Employee Gap POS Permutation.xlsx")
LA_OS.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_OS.Result$MSE_pre[1])
OS_Per$Chan<-OS_Per$POS-OS_Per$sdid_intecept
LA_OS.i<-filter(OS_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.OS<-left_join(STATE.LOW.MSE.15[,1:2],OS_Per,by=c("State"="STATE"))
SMALL.STATE.OS<-SMALL.STATE.OS%>%filter(State!="GA")
#SMALL.STATE.OS<-filter(SMALL.STATE.OS,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.OS,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_OS.i,colour="black")+
ggtitle(label = "Other Services")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.OS.15
SMALL.LA.OS.15
RESULT<- read_excel("New RESULT PSPS Permutattion.xlsx")
SPS_Per<- read_excel("New Employee Gap PSPS Permutation.xlsx")
LA_SPS.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_SPS.Result$MSE_pre[1])
SPS_Per$Chan<-SPS_Per$PSPS-SPS_Per$sdid_intecept
LA_SPS.i<-filter(SPS_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.SPS<-left_join(STATE.LOW.MSE.15[,1:2],SPS_Per,by=c("State"="STATE"))
#SMALL.STATE.SPS<-filter(SMALL.STATE.SPS,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.SPS,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_SPS.i,colour="black")+
ggtitle(label = "Professional, Scientific, \n and Technical Services")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.SPS.15
SMALL.LA.SPS.15
RESULT<- read_excel("New RESULT PTWU Permutattion.xlsx")
TWU_Per<- read_excel("New Employee Gap PTWU Permutation.xlsx")
LA_TWU.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_TWU.Result$MSE_pre[1])
TWU_Per$Chan<-TWU_Per$PTWU-TWU_Per$sdid_intecept
LA_TWU.i<-filter(TWU_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.TWU<-left_join(STATE.LOW.MSE.15[,1:2],TWU_Per,by=c("State"="STATE"))
#SMALL.STATE.TWU<-filter(SMALL.STATE.TWU,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.TWU,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_TWU.i,colour="black")+
ggtitle(label = "Transportation, Warehousing, \n and Utilities")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.TWU.15
SMALL.LA.TWU.15
RESULT<- read_excel("New RESULT PWT Permutattion.xlsx")
WT_Per<- read_excel("New Employee Gap PWT Permutation.xlsx")
LA_WT.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_WT.Result$MSE_pre[1])
WT_Per$Chan<-WT_Per$PWT-WT_Per$sdid_intecept
LA_WT.i<-filter(WT_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.WT<-left_join(STATE.LOW.MSE.15[,1:2],WT_Per,by=c("State"="STATE"))
#SMALL.STATE.WT<-filter(SMALL.STATE.WT,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.WT,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_WT.i,colour="black")+
ggtitle(label = "Wholesale Trade")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.WT.15
SMALL.LA.WT.15
RESULT<- read_excel("New RESULT PTNF Permutattion.xlsx")
TNF_Per<- read_excel("New Employee Gap PTNF Permutation.xlsx")
LA_TNF.Result<-filter(RESULT,State=="LA")
STATE.LOW.MSE.15<-filter(RESULT,MSE_pre<2*LA_TNF.Result$MSE_pre[1])
TNF_Per$Chan<-TNF_Per$PTNF-TNF_Per$sdid_intecept
LA_TNF.i<-filter(TNF_Per,STATE=="LA")
#############LESS THAN 2*LA_MSE####################
SMALL.STATE.TNF<-left_join(STATE.LOW.MSE.15[,1:2],TNF_Per,by=c("State"="STATE"))
#SMALL.STATE.TNF<-filter(SMALL.STATE.TNF,STATE!="VA")
ggplot()+geom_line(data=SMALL.STATE.TNF,aes(x=DATE,y=Chan,group=State),colour="grey")+
geom_line(aes(x=DATE,y=Chan),data = LA_TNF.i,colour="black")+
ggtitle(label = "Total Non-farm")+
theme_bw()+geom_vline(xintercept = hashV,linetype=2)+xlab("")+
theme(plot.title = element_text(h=0.5,size = 9),
axis.text = element_text(size =7))+ylab("")->SMALL.LA.TNF.15
SMALL.LA.TNF.15
