setwd("~/mailman/qmssviz/hw2/_posts/")

library(plyr)
library(dplyr)

lgas <- read.csv("lgas.csv")
niger <- read.csv("sample_health_facilities.csv")

names(lgas)
names(niger)

## select facilities located in Southern
## zones of Nigeria
unique(niger$zone) # look at different zones
nigerSouth <- niger[niger$zone %in% 
                      c("South-South", "Southeast", "Southwest"), ]
## get pop_2006 column
pop2006 <- lgas[, c("lga_id","pop_2006")]
## merge pop_2006 with southern facilities 
lgasSouth <- join(nigerSouth, pop2006, by = "lga_id")

## caculate total number of full time nurses and doctors
## for all health facilities in each state, 
## sorted by state population
lgasSouth %>%
    group_by(state) %>%
    summarise(num_nurses = sum(num_nurses_fulltime),
              num_doctors = sum(num_doctors_fulltime),
              population = sum(pop_2006)) %>%
    arrange(desc(population))