Import Dataset
sec<-c("AER","AF","ASW","CON","ES","FI",
"HCS","INF","MAN","MCE","OS","PSP",
"RE","RT","TW","U","WT")
sect<-data.frame("sec"=sec,
"names"=c("Arts, Entertainment, and Recreation",
"Accommodation and Food Services",
"Administrative and Support and Waste Management and Remediation Services",
"Construction","Educational Services",
"Finance Insurance","Health Care and Social Assistance",
"Information","Manufacturing","Management of Companies and Enterprises",
"Other Services","Professional, Scientific, and Technical Services",
"Real Estate and Rental and Leasing","Retail Trade","Transportation and Warehousing",
"Utilities","Wholesale Trade"))
##Imporing Data
for( i in sec){
assign(i,value = read_xlsx(paste0("MSA ",i," PLOT.xlsx")))
}
plot.theme <- theme(legend.position="bottom",
legend.direction = 'horizontal',
legend.key = element_blank(),axis.text.x = element_text(angle = 45, hjust = 1),
legend.background = element_blank(),
plot.title = element_text(hjust = 0.5, size = 12,face = "bold"),
axis.text = element_text(size = 8),axis.title = element_text(size = 9),
)
AER
AER$Date<-as.Date(AER$Date)
pl.AER<- ggplot(AER, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="AER",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.AER

AF
AF$Date<-as.Date(AF$Date)
pl.AF<- ggplot(AF, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="AF",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.AF

ASW
ASW$Date<-as.Date(ASW$Date)
pl.ASW<- ggplot(ASW, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="ASW",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.ASW

CON
CON$Date<-as.Date(CON$Date)
pl.CON<- ggplot(CON, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="CON",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.CON

ES
ES$Date<-as.Date(ES$Date)
pl.ES<- ggplot(ES, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="ES",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.ES

FI
FI$Date<-as.Date(FI$Date)
pl.FI<- ggplot(FI, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="FI",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.FI

HCS
HCS$Date<-as.Date(HCS$Date)
pl.HCS<- ggplot(HCS, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="HCS",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.HCS

INF
INF$Date<-as.Date(INF$Date)
pl.INF<- ggplot(INF, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="INF",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.INF

MAN
MAN$Date<-as.Date(MAN$Date)
pl.MAN<- ggplot(MAN, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="MAN",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.MAN

MCE
MCE$Date<-as.Date(MCE$Date)
pl.MCE<- ggplot(MCE, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="MCE",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.MCE

OS
OS$Date<-as.Date(OS$Date)
pl.OS<- ggplot(OS, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="OS",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.OS

PSP
PSP$Date<-as.Date(PSP$Date)
pl.PSP<- ggplot(PSP, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="PSP",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.PSP

RE
RE$Date<-as.Date(RE$Date)
pl.RE<- ggplot(RE, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="RE",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.RE

RT
RT$Date<-as.Date(RT$Date)
pl.RT<- ggplot(RT, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="RT",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.RT

TW
TW$Date<-as.Date(TW$Date)
pl.TW<- ggplot(TW, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="TW",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.TW

U
U$Date<-as.Date(U$Date)
pl.U<- ggplot(U, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="U",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.U

WT
WT$Date<-as.Date(WT$Date)
pl.WT<- ggplot(WT, aes(x = Date, y = average)) +
geom_point()+
geom_line() + theme_bw()+ylab(sect[sec=="WT",]$names)+xlab("Year")+
geom_vline(xintercept = as.Date("2005-07-01"),linetype=2)+
geom_ribbon(aes(ymin = lo, ymax = up), alpha = 0.2)+plot.theme
pl.WT
