mort <- read.csv('~/Desktop/Data Incubator/Proposed Project/mortality/2005_data.csv')
#summary(mort)
#names(mort)

library(ggplot2)

mort_by_brain <- subset(mort, (mort$X358_cause_recode >= 174) & (mort$X358_cause_recode <= 194) & 
                          (mort$detail_age >= 20 & mort$detail_age <= 80))

p1 <- ggplot(aes(x = mort_by_brain$X358_cause_recode), data = mort_by_brain) +
  geom_histogram() +
  facet_wrap(~sex) +
  xlab("Deaths by cause recode")
p1

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.