ipumsr Example - NHGIS

Minnesota Population Center

2017-12-15

IPUMS - NHGIS Extraction and Analysis

Exercise 1

OBJECTIVE: Gain an understanding of how the NHGIS datasets are structured and how they can be leveraged to explore your research interests. This exercise will use an NHGIS dataset to explore slavery in the United States in 1830.

This vignette is adapted from the NHGIS Data Training Exercise available here: https://pop.umn.edu/sites/pop.umn.edu/files/nhgis_training_ex1_2017-01.pdf

Research Question

What was the state‐level distribution of slavery in 1830?

Objectives

  • Create and download an NHGIS data extract
  • Unzip data file and open in R
  • Analyze the data using R
  • Validate data analysis work using answer key

Download Extract from IPUMS Website

  1. Log in to NHGIS
  1. Find Tables

Q1) How many tables are available from the 1830 Census?

Q2) Other than slave status, what are some other topics could we learn about for 1830?

Q3) Click the table name to see additional information. How many variables does this table contain?

Q4) For which geographic levels is the table available?

Q5) Close the table pop‐up window and inspect the Select Data table… What is the universe for this table?

Q6) What differentiates this table from the other available slavery tables from 1830?

Q7) Name a percentage or ratio this table would allow us to calculate that the other tables would not, based on the counts available in each table?

  1. Create a Data Extract Creating a data extract requires the user to select the table(s), specify a geographic level, and select the data layout structure…
  1. Download the Data Extract From the Extracts History page, you will be able to download your data extract once it has finished processing, typically within a few minutes. You may leave this page and return once you have received the email alerting you to your finished extract.

If you refresh your browser window (click on the loop icon at top, or press F5), you will see the extract status change from ‘queued’ to ‘in progress’ to ‘complete’, at which time you will be able to click the ‘tables’ link to download the data.

Getting the data into R

You will need to change the filepaths noted below to the place where you have saved the extracts.

Getting the data into R

You will need to change the filepaths noted below to the place where you have saved the extracts.

#> Could not find NHGIS data and so could not run vignette.
#> 
#> If you tried to download the data following the instructions above, please makesure that the filenames are correct: 
#> csv - nhgis0001_csv.zip
#> shape - nhgis0001_shape.zip
#> And that you are in the correct directory if you are using a relative path:
#> Current directory - C:/Users/umn-burkx031/AppData/Local/Temp/Rtmpm6X7mY/Rbuilde8418a62dc1/ipumsr/vignettes
#> 
#> The data is also available on github. You can install it using the following commands: 
#>   if (!require(devtools)) install.packages('devtools')
#>   devtools::install_github('mnpopcenter/ipumsr/ipumsexamples')
#> After installation, the data should be available for this vignette.

Note that read_nhgis_sf relies on package sf. You can also read NHGIS data into the format used by package sp with function read_nhgis_sp.

Exercises

These exercises include example code written in the “tidyverse” style, meaning that they use the dplyr package. This package provides easy to use functions for data analysis, including mutate(), select(), arrange(), slice() and the pipe (%>%). There a numerous other ways you could solve these answers, including using the base R, the data.table package and others.

Analyze the Data

Q8) How many states/territories are included in this table?

Q9) Why do you think other states are missing?

Q10) Create a new variable called total_pop, with the total population for each state, by summing the counts in columns ABO001 to ABO006. Which state had the largest population?

Q11) Create a variable called slave_pop, with the total slave population by summing the variables ABO003 and ABO004. Which state had the largest slave population?

Q12) Create a variable called pct_slave with the Slave Population divided by the Total Population. Which states had the highest and lowest Percent Slave Population?

Q13) Are there any surprises, or is it as you expected?

Inspect the Codebook

Open the .txt codebook file that is in the same folder as the comma delimited file you have already analyzed. The codebook file is a valuable reference containing information about the table or tables you’ve downloaded.

Some of the information provided in the codebook can be read into R, using the function read_ipums_codebook().

Q14) What is the proper citation to provide when using NHGIS data in publications or researcher reports?

Q15) What is the email address for NHGIS to share any research you have published? (You can also send questions you may have about the site. We’re happy to help!)

Bonus - Make Maps using R

One of the reasons we are excited about bringing IPUMS data to R is the GIS capabilities available for free in R.

Q16) Make a map of the percent of the population that are slaves.