alaska.R

# devtools::install_github("walkerke/tidycensus")
library(tidycensus)
library(tidyverse)
library(geofacet)
library(stringr)
library(extrafont)

# 2010 Census

age <- get_decennial(geography = "county", state = "AK", table = "P012", 
                     summary_var = "P0010001") %>%
  mutate(variable = str_replace(variable, "P01200", "")) %>%
  filter(!variable %in% c("01", "02", "26")) %>%
  arrange(NAME, variable)


# Wrangle it
agegroups <- c("0-4", "5-9", "10-14", "15-19", "15-19", "20-24", "20-24", 
               "20-24", "25-29", "30-34", "35-39", "40-44", "45-49", "50-54", 
               "55-59", "60-64", "60-64", "65-69", "65-69", "70-74", "75-79", 
               "80-84", "85+")

agesex <- c(paste("Male", agegroups), 
            paste("Female", agegroups))

agefull <- rep(agesex, length(unique(age$NAME)))

age$group <- agefull

age2 <- age %>%
  group_by(NAME, group) %>%
  mutate(group_est = sum(value)) %>%
  distinct(NAME, group, .keep_all = TRUE) %>%
  ungroup() %>%
  mutate(percent = 100 * (group_est / summary_value), 
         NAME = str_replace(NAME, ", Alaska", "")) %>%
  mutate(NAME = str_replace(NAME, " Borough| Census Area| Municipality| City and Borough", "")) %>%
  mutate(NAME = str_replace(NAME, "Wade Hampton", "Kusilvak")) %>%
  select(name = NAME, group, percent) %>%
  separate(group, into = c("sex", "age"), sep = " ") %>%
  mutate(age = factor(age, levels = unique(age)), 
         percent = ifelse(sex == "Female", percent, -percent)) 

xlabs = c("0-4" = "0-4", "5-9" = "", "10-14" = "", "15-19" = "", "20-24" = "", 
          "25-29" = "", "30-34" = "", "35-39" = "", "40-44" = "", "45-49" = "", 
          "50-54" = "", "55-59" = "", "60-64" = "", "65-69" = "", "70-74" = "", 
          "75-79" = "", "80-84" = "", "85+" = "85+")

mygrid <- data.frame(
  row = c(1, 1, 2, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 8, 8, 9, 9, 10, 10),
  col = c(4, 5, 3, 3, 4, 3, 5, 6, 3, 5, 6, 4, 5, 6, 7, 8, 3, 1, 2, 8, 9, 10, 4, 8, 9, 8, 9, 9, 10),
  code = c("Northwest Arctic", "North Slope", "Nome", "Kusilvak", "Yukon-Koyukuk", "Bethel", "Denali", "Fairbanks North Star", "Dillingham", "Matanuska-Susitna", "Southeast Fairbanks", "Lake and Peninsula", "Kenai Peninsula", "Anchorage", "Valdez-Cordova", "Yakutat", "Bristol Bay", "Aleutians West", "Aleutians East", "Hoonah-Angoon", "Haines", "Skagway", "Kodiak Island", "Sitka", "Juneau", "Prince of Wales-Hyder", "Petersburg", "Wrangell", "Ketchikan Gateway"),
  name = c("Northwest Arctic", "North Slope", "Nome", "Kusilvak", "Yukon-Koyukuk", "Bethel", "Denali", "Fairbanks North Star", "Dillingham", "Matanuska-Susitna", "Southeast Fairbanks", "Lake and Peninsula", "Kenai Peninsula", "Anchorage", "Valdez-Cordova", "Yakutat", "Bristol Bay", "Aleutians West", "Aleutians East", "Hoonah-Angoon", "Haines", "Skagway", "Kodiak Island", "Sitka", "Juneau", "Prince of Wales-Hyder", "Petersburg", "Wrangell", "Ketchikan Gateway"),
  stringsAsFactors = FALSE
)

ggplot(data = age2, aes(x = age, y = percent, fill = sex)) +
  geom_bar(stat = "identity", width = 1) + 
  scale_y_continuous(breaks=c(-5, 0, 5),labels=c("5%", "0%", "5%")) + 
  coord_flip() + 
  theme_minimal(base_family = "Tahoma") + 
  scale_x_discrete(labels = xlabs) + 
  scale_fill_manual(values = c("red", "navy")) + 
  theme(panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(), 
        strip.text.x = element_text(size = 7.75, hjust = 0.62)) + 
  labs(x = "", y = "", fill = "", 
       title = "Demographic structure of Alaska counties", 
       caption = "Data source: 2010 US Census, tidycensus R package.  Chart by @kyle_e_walker.") + 
  facet_geo(~ name, grid = mygrid) 

ggsave("plots/alaska.png", dpi = 300, width = 12, height = 10)

I'd like to suggest this arrangement which to me looks a little closer in my mind and has I think has more relative directions.

mygrid <- data.frame(
row = c(1, 1, 2, 3, 2, 4, 3, 2, 5, 4, 2, 6, 6, 5, 5, 5, 6, 8, 7, 7, 6, 6, 7, 8, 7, 9, 8, 8, 9),
col = c(4, 5, 3, 3, 4, 3, 4, 5, 3, 5, 6, 4, 6, 5, 7, 8, 3, 1, 2, 9, 9, 10, 4, 8, 10, 9, 9, 10, 10),
code = c("Northwest Arctic", "North Slope", "Nome", "Kusilvak", "Yukon-Koyukuk", "Bethel", "Denali", "Fairbanks North Star", "Dillingham", "Matanuska-Susitna", "Southeast Fairbanks", "Lake and Peninsula", "Kenai Peninsula", "Anchorage", "Valdez-Cordova", "Yakutat", "Bristol Bay", "Aleutians West", "Aleutians East", "Hoonah-Angoon", "Haines", "Skagway", "Kodiak Island", "Sitka", "Juneau", "Prince of Wales-Hyder", "Petersburg", "Wrangell", "Ketchikan Gateway"),
name = c("Northwest Arctic", "North Slope", "Nome", "Kusilvak", "Yukon-Koyukuk", "Bethel", "Denali", "Fairbanks North Star", "Dillingham", "Matanuska-Susitna", "Southeast Fairbanks", "Lake and Peninsula", "Kenai Peninsula", "Anchorage", "Valdez-Cordova", "Yakutat", "Bristol Bay", "Aleutians West", "Aleutians East", "Hoonah-Angoon", "Haines", "Skagway", "Kodiak Island", "Sitka", "Juneau", "Prince of Wales-Hyder", "Petersburg", "Wrangell", "Ketchikan Gateway"),
stringsAsFactors = FALSE
)

ggplot(mygrid) + geom_text(aes(x = col, y = -row, label = name))

Agreed, that looks much better! I'm going to make a website with all of the state plots and I'll use your grid for it.