by_common_jobs <- family %>%
  group_by(Household.Head.Occupation) %>% filter(!Household.Head.Occupation == "NA") %>%
  dplyr::summarise_at(vars(Total.Household.Income), funs(mean(., na.rm = TRUE)))  %>%
  head(20) %>% ungroup()
by_common_jobs
ggplot(by_common_jobs, aes(x=reorder(Household.Head.Occupation, Total.Household.Income), Total.Household.Income)) +
  geom_bar(fill = "#3288BD", stat = "identity") + geom_text(aes(x = Household.Head.Occupation,
  y = 1, label = paste0("(",Total.Household.Income,")",sep="")), hjust=0, vjust=.5, size = 4, colour = 'black', fontface = 'italic') + bluetheme + coord_flip() + xlab("")