plotting_code.R

library(data.table)
library(ggplot2)
library(ggsci)

dat <- fread("~/Downloads/neoplasm_vs_circulatory_mortality_parsed.txt")
colnames(dat) <- c("V0", "Gender", "GenderCode", "CauseOfDeath", "ICD", "Year", "Year1", "Deaths", "Population", "CrudeRate","AgeAdjustedRate")
dat$Cause <- plyr::mapvalues(dat$CauseOfDeath, from=c("Diseases of the circulatory system","Neoplasms"), to=c("Circulatory\nSystem","Neoplasm"))

ggplot(dat, aes(x=Year, y=CrudeRate, color=Gender, linetype=Cause)) + 
  geom_line() + 
  theme_bw() + 
  scale_color_aaas() +
  labs(x="Year", y="Death Rate", title="Crude Death Rate") 

ggplot(dat, aes(x=Year, y=AgeAdjustedRate, color=Gender, linetype=Cause)) + 
  geom_line() + 
  theme_bw() + 
  scale_color_aaas() +
  labs(x="Year", y="Death Rate", title="Age-adjusted Death Rate")