{"id":11445,"date":"2018-08-27T13:45:04","date_gmt":"2018-08-27T18:45:04","guid":{"rendered":"https:\/\/rud.is\/b\/?p=11445"},"modified":"2018-08-27T13:45:04","modified_gmt":"2018-08-27T18:45:04","slug":"simplifying-world-tile-grid-creation-with-geom_wtg","status":"publish","type":"post","link":"https:\/\/rud.is\/b\/2018\/08\/27\/simplifying-world-tile-grid-creation-with-geom_wtg\/","title":{"rendered":"Simplifying World Tile Grid Creation with geom_wtg()"},"content":{"rendered":"

Nowadays (I’ve seen that word used so<\/em> much in journal articles lately that I could not resist using it) I’m using world tile grids more frequently as the need arises to convey the state of exposure of various services at a global (country) scale. Given that necessity fosters invention it seemed that having a ggplot2 geom for world tile grids would make my life easier and also make work more efficient.<\/p>\n

To that end, there’s a nascent ggplot2 extension package for making world tile grids called (uncreatively) worldtilegrid<\/code>?<\/a>. It’s also at GitHub if you’re more comfortable working from code from there.<\/p>\n

The README<\/a> has examples but we’ll walk through another one here and work with life expectancy data curated by Our World in Data<\/a>. Rather than draw out this post to a less tenable length with an explanation of how to pull the XHR JSON data into R, just grab the CSV linked with the visualization on that page and substitute the path in the code for where you stored it.<\/p>\n

We do need to clean up this data a bit since it has some issues. Let’s do that and carve out some slices into two new data frames so we can work with the most recent curated year and some historical data:<\/p>\n

library(hrbrthemes) # gitlab\/github :: hrbrmstr\/hrbrthemes (use devtools installers)\nlibrary(worldtilegrid) # gitlab\/github :: hrbrmstr\/worldtilegrid (use devtools installers)\nlibrary(tidyverse)\n\ncols(\n  Entity = col_character(),\n  Code = col_character(),\n  Year = col_integer(),\n  `Life expectancy (Clio-Infra up to 1949; UN Population Division for 1950 to 2015)` = col_double()\n) -> le_cols\n\nread_csv(\"~\/Downloads\/life-expectancy.csv\", col_types = le_cols) %>%\n  set_names(c(\"country\", \"iso3c\", \"year\", \"life_expectancy\")) -> lifexp\n\n# clean it up a bit since it's not great data and bust it into groups\n# that match the vis on the site vs use continuous raw data, esp since\n# the data is really just estimates\nfilter(lifexp, !is.na(iso3c)) %>%\n  filter(nchar(iso3c) == 3) %>%\n  mutate(\n    grp = cut(\n      x = life_expectancy,\n      breaks = c(10, 20, 30, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85),\n      include.lowest = TRUE\n    )\n  ) -> lifexp\n\nlast_year <-  filter(lifexp, year == last(year))\na_few_years <- filter(lifexp, year %in% c(seq(1900, 2015, 20), 2015))<\/code><\/pre>\n
Now we can use it to make some cartograms.<\/p>\n
Using the World Tile Grid geom<\/h3>\n
You can reference a previous post<\/a> for what it would normally take to create these tile grids. The idea was to simplify the API down to the caller just needing to specify some country names\/ISO[23]Cs (to a country<\/code> aesthetic) and a value (to a fill<\/code> aesthetic) and let the package do the rest. To that end, here's how to display the life expectancy data for 2015:<\/p>\n
ggplot(last_year) +\n  geom_wtg(aes(country = iso3c, fill = grp), border_col = \"#2b2b2b\") + # yep, you can sum up the entire blog post to this one line\n  viridis::scale_fill_viridis(\n    name = \"Life Expectancy\", discrete = TRUE, na.value=alpha(ft_cols$gray, 1\/10),\n    drop = FALSE\n  ) +\n  coord_equal() +\n  labs(\n    x=NULL, y=NULL,\n    title = \"Life expectancy, 2015\",\n    subtitle = \"Shown is period life expectancy at birth. This corresponds to an estimate of the average number\\nof years a newborn infant would live if prevailing patterns of mortality at the time of its birth were to\\nstay the same throughout its life\",\n    caption = \"Data source: \"\n\n  ) +\n  theme_ft_rc(grid=\"\") +\n  theme(axis.text=element_blank()) +\n  theme(legend.position = \"bottom\")<\/code><\/pre>\n
\"\"<\/a><\/p>\n
You can add labels via geom_text()<\/code> and use the wtg<\/code> stat for it as it provides a number of computed variables you can work with:<\/p>\n