Find and download a shape file

Note: this step is described in detail here. Read it if you are not familiar with geospatial data management in R.

The region boundaries required to make maps are usually stored in geospatial objects. Those objects can come from shapefiles, geojson files or provided in a R package. See the map section for possibilities.

Let’s get a geospatial object from a shape file available here.

# Download the shapefile. (note that I store it in a folder called DATA. You have to change that if needed.)
download.file("http://thematicmapping.org/downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip", destfile = "DATA/world_shape_file.zip")
# You now have it in your current working directory, have a look!

# Unzip this file. You can do it with R (as below), or clicking on the object you downloaded.
system("unzip DATA/world_shape_file.zip")
#  -- > You now have 4 files. One of these files is a .shp file! (TM_WORLD_BORDERS_SIMPL-0.3.shp)

And let’s load it in R

# Read this shape file with the sf library.
library(sf)
my_sf <- read_sf(
  paste0(getwd(), "/DATA/world_shape_file/TM_WORLD_BORDERS_SIMPL-0.3.shp")
)

# Select Africa only
africa <- my_sf[my_sf$REGION == 2, ]

# -- > Now you have a sf object (simple feature data frame). You can start doing maps!

Plot the map

The hardest part is to build the color palette that will be used on the map. The RColorBrewer package makes it a breeze to get a nice set of colors thanks to its brewer.pal() function.

Once the palette is created, POP2005 is cut in several bins using the cut() function, and the appropriate color is attributed to each bin.

Finally, build the map with plot() and give the color vector to the col argument to get the choropleth.

# Plot
plot(st_geometry(africa), xlim = c(-20, 60), ylim = c(-40, 40))

Visualize the numeric variable

In a choropleth map, each region has a color that represents the value of a numeric variable.

In this example, this variable is already available in the africa object. You can have a look to it typing africa in the terminal.

It is thus possible to visualize the distribution of POP2005 with a classic histogram made in ggplot2

# library
library(dplyr)
library(ggplot2)

# Make sure the variable you are studying is numeric
africa$POP2005 <- as.numeric(africa$POP2005)

# Distribution of the population per country?
africa %>%
  ggplot(aes(x = POP2005)) +
  geom_histogram(bins = 20, fill = "#69b3a2", color = "white")

Build the choropleth