Updates

expl.R

 str(dat)
 
 yeartab <- table(dat$CRASH_YEAR)
-mpy <- mean(yeartab)
+mpy <- sum(yeartab)/(18.75)
 
 # Obviously, data missing from 2018
 ggplot(dat) + geom_bar(aes(x=CRASH_YEAR), fill=NA, col="black") +
   theme_classic() + labs(x = "Crash severity", y = "Number of crashes") +
   scale_y_continuous(labels = scales::comma)
 
+knitr::kable(table(dat$MULTI_VEH, dat$CRASH_SEV))
+
 # Not all openroad crashes are on state highways
 table(dat$URBAN, dat$CRASH_SH_DESC)