Data Science Project Revision

newcdata <- co2_monthly %>%
  group_by(decade) %>%
  mutate(meanc = mean(int_mean_co2))


ggplot(co2_monthly, aes(x = year, y = int_mean_co2))+
  #geom_line(aes(x = date, y = int_mean_co2), col = "blue", alpha =0.75) + 
  geom_point(aes(x = date, y = int_mean_co2), col = "black", shape = 20, size =0.5) + 
  geom_smooth(method = "loess", aes(col=decade))+
  geom_vline(xintercept = 1963, color = "black")+
  annotate("label", x =1963, y = 360, label = "JFK Assassinated")+
  geom_vline(xintercept = 1980, color = "black")+
  annotate("label", x =1980, y = 370, label = "John Lennon Killed")+
  geom_vline(xintercept = 1994, color = "black")+
  annotate("label", x =1994, y = 380, label = "Kurt Cobain Suicide")+
  geom_vline(xintercept = 2018, color = "black")+
  annotate("label", x =2010, y = 340, label = "XXXtentacion Murdered")+
  geom_hline(yintercept=400, color="red")+
    annotate("label", x =1990, y = 400, label = "Crisis Threshold")+
  ggtitle(expression("Yearly CO"[2]*" AT MAUNA LOA and Important Historical Deaths")) +
  ylab(expression("CO"[2]*"(ppm)")) +
  xlab("Year") +
  theme_linedraw()+
  theme(axis.line.y.left = element_line(),
      axis.line.x.bottom = element_line(),
        panel.grid.major = element_line(colour = "gray", size = .5))

  #      panel.grid.minor = element_line(colour = "light gray", size = .25, linetype = "solid"),
  #     panel.background = element_rect(NULL))