library(tidyverse)
Nobel Laureates
In January 2017, Buzzfeed published an article on why Nobel laureates show immigration is so important for American science. You can read the article here. In the article they show that while most living Nobel laureates in the sciences are based in the US, many of them were born in other countries. This is one reason why scientific leaders say that immigration is vital for progress. In this lab we will work with the data from this article to recreate some of their visualizations as well as explore new questions.
Learning goals
- Replicating published results
- Data wrangling and visualization
Read the Buzzfeed article titled These Nobel Prize Winners Show Why Immigration Is So Important For American Science. We will replicate this analysis in the workshop so it’s crucial that you’re familiar with it ahead of time.
Packages
We’ll use the tidyverse package for much of the data wrangling. This package is already installed for you. You can load them by running the following in your Console:
Data
The dataset for this assignment can be found as a CSV (comma separated values) file. You can read it in using the following.
<- read_csv("data04/nobel.csv") nobel
The variable descriptions are as follows:
id
: ID numberfirstname
: First name of laureatesurname
: Surnameyear
: Year prize woncategory
: Category of prizeaffiliation
: Affiliation of laureatecity
: City of laureate in prize yearcountry
: Country of laureate in prize yearborn_date
: Birth date of laureatedied_date
: Death date of laureategender
: Gender of laureateborn_city
: City where laureate was bornborn_country
: Country where laureate was bornborn_country_code
: Code of country where laureate was borndied_city
: City where laureate dieddied_country
: Country where laureate dieddied_country_code
: Code of country where laureate diedoverall_motivation
: Overall motivation for recognitionshare
: Number of other winners award is shared withmotivation
: Motivation for recognition
In a few cases the name of the city/country changed after laureate was given (e.g. in 1975 Bosnia and Herzegovina was called the Socialist Federative Republic of Yugoslavia). In these cases the variables below reflect a different name than their counterparts without the suffix `_original`.
born_country_original
: Original country where laureate was bornborn_city_original
: Original city where laureate was borndied_country_original
: Original country where laureate dieddied_city_original
: Original city where laureate diedcity_original
: Original city where laureate lived at the time of winning the awardcountry_original
: Original country where laureate lived at the time of winning the award
Exercises
Get to know your data
- How many observations and how many variables are in the dataset? Use inline code to answer this question. What does each row represent?
There are some observations in this dataset that we will exclude from our analysis to match the Buzzfeed results.
- Create a new data frame called
nobel_living
that filters for
- laureates for whom
country
is available - laureates who are people as opposed to organizations (organizations are denoted with
"org"
as theirgender
) - laureates who are still alive (their
died_date
isNA
)
Confirm that once you have filtered for these characteristics you are left with a data frame with 228 observations, once again using inline code.
Most living Nobel laureates were based in the US when they won their prizes
… says the Buzzfeed article. Let’s see if that’s true.
First, we’ll create a new variable to identify whether the laureate was in the US when they won their prize. We’ll use the mutate()
function for this. The following pipeline mutates the nobel_living
data frame by adding a new variable called country_us
. We use an if statement to create this variable. The first argument in the if_else()
function we’re using to write this if statement is the condition we’re testing for. If country
is equal to "USA"
, we set country_us
to "USA"
. If not, we set the country_us
to "Other"
.
Note: we can achieve the same result using the fct_other()
function we’ve seen before (i.e. with country_us = fct_other(country, "USA")
). We decided to use the if_else()
here to show you one example of an if statement in R.
<- nobel_living %>%
nobel_living mutate(
country_us = if_else(country == "USA", "USA", "Other")
)
Next, we will limit our analysis to only the following categories: Physics, Medicine, Chemistry, and Economics.
<- nobel_living %>%
nobel_living_science filter(category %in% c("Physics", "Medicine", "Chemistry", "Economics"))
For the next exercise work with the nobel_living_science
data frame you created above.
Create a faceted bar plot visualizing the relationship between the category of prize and whether the laureate was in the US when they won the nobel prize. Interpret your visualization, and say a few words about whether the Buzzfeed headline is supported by the data.
- Your visualization should be faceted by category.
- For each facet you should have two bars, one for winners in the US and one for Other.
- Flip the coordinates so the bars are horizontal, not vertical.
But of those US-based Nobel laureates, many were born in other countries
Hint: You should be able to cheat borrow from code you used earlier to create the country_us
variable.
Create a new variable called
born_country_us
that has the value"USA"
if the laureate is born in the US, and"Other"
otherwise. How many of the winners are born in the US?Add a second variable to your visualization from Exercise 3 based on whether the laureate was born in the US or not. Based on your visualization, do the data appear to support Buzzfeed’s claim? Explain your reasoning in 1-2 sentences.
- Your final visualization should contain a facet for each category.
- Within each facet, there should be a bar for whether the laureate won the award in the US or not.
- Each bar should have segments for whether the laureate was born in the US or not.
Here’s where those immigrant Nobelists were born
Note: your bar plot won’t exactly match the one from the Buzzfeed article. This is likely because the data has been updated since the article was published.
- In a single pipeline, filter for laureates who won their prize in the US, but were born outside of the US, and then create a frequency table (with the
count()
function) for their birth country (born_country
) and arrange the resulting data frame in descending order of number of observations for each country. Which country is the most common?
Interested in how Buzzfeed made their visualizations?
The plots in the Buzzfeed article are called waffle plots. You can find the code used for making these plots in Buzzfeed’s GitHub repo (yes, they have one!) here. You’re not expected to recreate them as part of your assignment, but you’re welcomed to do so for fun!