library(ggplot2)
library(hrbrthemes)

p_cases_vs_pop <- ggplot(cases_pop, aes(x = population, y = confirmed_cases)) + 
                  geom_point() +
                  labs(x="Population", y="Confirmed cases",
                  title="Population versus COVID-19 Cases, by ZIP Code",
                  subtitle="Restricted to San Diego County",
                  caption="Updated April 4, 2020") + 
                  theme_ipsum_rc()