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