r4ds/joins.qmd

957 lines
38 KiB
Plaintext
Raw Normal View History

# Joins {#sec-joins}
```{r}
#| results: "asis"
#| echo: false
source("_common.R")
2022-09-10 21:38:12 +08:00
status("polishing")
```
2016-07-19 21:01:50 +08:00
## Introduction
It's rare that a data analysis involves only a single data frame.
2022-06-20 23:08:33 +08:00
Typically you have many data frames, and you must **join** them together to answer the questions that you're interested in.
This chapter will introduce you to two important types of joins:
2016-01-15 21:05:45 +08:00
- Mutating joins, add new variables to one data frame from matching observations in another.
- Filtering joins, filter observations from one data frame based on whether or not they match an observation in another.
We'll begin by discussing keys, the variables used to connect a pair of data frames in a join.
You'll then see how to use joins to tackle a variety of challenges from the nycflights13 dataset.
Next we'll discuss how joins work, focusing on their action on the rows.
We'll finish up with a discussion of non-equi-joins, a family of joins that provide a more flexible way of matching keys than the default equality relationship.
If you're familiar with SQL, you should find the ideas in this chapter familiar, as their realization in dplyr is very similar.
2016-07-19 21:01:50 +08:00
### Prerequisites
2022-09-19 05:23:42 +08:00
::: callout-important
2022-09-26 21:38:18 +08:00
This chapter relies on features only found in dplyr 1.1.0, which is still in development.
2022-09-19 05:23:42 +08:00
If you want to live life on the edge you can get the dev version with `devtools::install_github("tidyverse/dplyr")`.
:::
2022-09-06 23:25:59 +08:00
We'll explore the five related datasets from nycflights13 using the join functions from dplyr.
2016-07-19 21:01:50 +08:00
```{r}
#| label: setup
#| message: false
2016-10-04 01:30:24 +08:00
library(tidyverse)
2016-07-19 21:01:50 +08:00
library(nycflights13)
```
2022-09-06 23:25:59 +08:00
## Keys
To understand joins, you need to first understand how two tables might be connected.
The connection between a pair of tables is defined by a pair of keys, which each consist of one or more variables.
In this section, you'll learn about the two types of key and their realization in the datasets of the nycflights13 package.
2022-09-06 23:25:59 +08:00
You'll also learn how to check that your keys are valid, and what to do if your table lacks a key.
2022-09-06 23:25:59 +08:00
### Primary and foreign keys
Every join involves a pair of keys: a primary key and a foreign key.
2022-09-16 21:00:35 +08:00
A **primary key** is a variable that uniquely identifies an observation.
A **foreign key** is the corresponding variable in another table.
Both primary and foreign keys can consist of more than one variable, which we'll call a **compound key**.
Let's make those terms concrete by looking more of the data in nycfights13:
2022-09-06 23:25:59 +08:00
- `airlines` lets you look up the full carrier name from its abbreviated code.
Its primary key is the two letter `carrier` code.
2016-01-30 21:13:46 +08:00
```{r}
airlines
```
2022-09-06 23:25:59 +08:00
- `airports` gives information about each airport.
Its primary key is the three letter `faa` airport code.
2016-01-30 21:13:46 +08:00
```{r}
airports
```
2016-01-13 23:06:43 +08:00
2022-09-06 23:25:59 +08:00
- `planes` gives information about each plane.
It's primary key is the `tailnum`.
2016-01-13 23:06:43 +08:00
```{r}
planes
```
2016-01-30 21:13:46 +08:00
2022-09-06 23:25:59 +08:00
- `weather` gives the weather at each NYC airport for each hour.
It has a compound primary key; to uniquely identify each observation you need to know both `origin` (the location) and `time_hour` (the time).
```{r}
weather
```
These datasets are all connected via the `flights` data frame because the `tailnum`, `carrier`, `origin`, `dest`, and `time_hour` variables are all foreign keys:
2022-09-06 23:25:59 +08:00
- `flights$tailnum` connects to primary key `planes$tailnum`.
- `flights$carrier` connects to primary key `airlines$carrier`.
2022-09-06 23:25:59 +08:00
- `flights$origin` connects to primary key `airports$faa`.
- `flights$dest` connects to primary key `airports$faa` .
- `flights$origin`-`flights$time_hour` connects to primary key `weather$origin`-`weather$time_hour`.
2022-09-16 21:00:35 +08:00
You'll notice a nice feature in the design of these keys: they almost all have the same name in both tables, which, as you'll see shortly, will make your joining life much easier.
It's also worth noting the opposite relationship: almost every variable name used in multiple tables has the same meaning in each place.
There's only one exception: `year` means year of departure in `flights` and year of manufacturer in `planes`.
This will become important when we start actually joining tables together.
2022-09-06 23:25:59 +08:00
We can also draw these relationships, as in @fig-flights-relationships.
This diagram is a little overwhelming, but it's simple compared to some you'll see in the wild!
The key to understanding diagrams like this is that you'll solve real problems by working with pairs of data frames.
You don't need to understand the whole thing; you just need to understand the chain of connections between the two data frames that you're interested in.
```{r}
#| label: fig-flights-relationships
#| echo: false
#| out-width: ~
#| fig-cap: >
#| Connections between all five data frames in the nycflights13 package.
2022-09-06 23:25:59 +08:00
#| Variables making up a primary key are coloured grey, and are connected
2022-09-20 12:09:54 +08:00
#| to their corresponding foreign keys with arrows.
#| fig-alt: >
#| Diagram showing the relationships between airports, planes, flights,
#| weather, and airlines datasets from the nycflights13 package. The faa
#| variable in the airports data frame is connected to the origin and dest
#| variables in the flights data frame. The tailnum variable in the planes
#| data frame is connected to the tailnum variable in flights. The
#| time_hour and origin variables in the weather data frame are connected
#| to the variables with the same name in the flights data frame. And
#| finally the carrier variables in the airlines data frame is connected
#| to the carrier variable in the flights data frame. There are no direct
#| connections between airports, planes, airlines, and weather data
#| frames.
knitr::include_graphics("diagrams/relational.png", dpi = 270)
```
2022-09-06 23:25:59 +08:00
### Checking primary keys
2016-01-15 21:05:45 +08:00
Now that that we've identified the primary keys in each table, it's good practice to verify that they do indeed uniquely identify each observation.
2022-09-06 23:25:59 +08:00
One way to do that is to `count()` the primary keys and look for entries where `n` is greater than one.
This reveals that `planes` and `weather` both look good:
2016-01-15 21:05:45 +08:00
```{r}
2022-02-24 03:15:52 +08:00
planes |>
count(tailnum) |>
2016-08-01 00:32:16 +08:00
filter(n > 1)
2022-02-24 03:15:52 +08:00
weather |>
count(time_hour, origin) |>
2016-08-01 00:32:16 +08:00
filter(n > 1)
2016-01-15 21:05:45 +08:00
```
2022-09-06 23:25:59 +08:00
You should also check for missing values in your primary keys --- if a value is missing then it can't identify an observation!
```{r}
planes |>
filter(is.na(tailnum))
weather |>
filter(is.na(time_hour) | is.na(origin))
```
### Surrogate keys
So far we haven't talked about the primary key for `flights`.
It's not super important here, because there are no data frames that use it as a foreign key, but it's still useful to consider because it's easier to work with observations if have some way to describe them to others.
2022-09-06 23:25:59 +08:00
After a little thinking and experimentation, we determined that there are three variables that together uniquely identify each flight:
2016-01-15 21:05:45 +08:00
```{r}
2022-02-24 03:15:52 +08:00
flights |>
count(time_hour, carrier, flight) |>
2016-08-01 00:32:16 +08:00
filter(n > 1)
2016-01-15 21:05:45 +08:00
```
Does the absence of duplicates automatically make `time_hour`-`carrier`-`flight` a primary key?
2022-09-06 23:25:59 +08:00
It's certainly a good start, but it doesn't guarantee it.
For example, are altitude and longitude a good primary key for `airports`?
2022-09-06 23:25:59 +08:00
```{r}
airports |>
count(alt, lat) |>
filter(n > 1)
```
Identifying an airport by it's altitude and latitude is clearly a bad idea, and in general it's not possible to know from the data alone whether or not a combination of variables makes a good a primary key.
But for flights, the combination of `time_hour`, `carrier`, and `flight` seems reasonable because it would be really confusing for an airline and its customers if there were multiple flights with the same number in the air at the same time.
2022-09-06 23:25:59 +08:00
That said, we might be better off introducing a simple numeric surrogate key using the row number:
2022-09-06 23:25:59 +08:00
```{r}
flights2 <- flights |>
mutate(id = row_number(), .before = 1)
flights2
```
2022-09-16 21:00:35 +08:00
Surrogate keys can be particular useful when communicating to other humans: it's much easier to tell someone to take a look at flight 2001 than to say look at UA430 which departed 9am 2013-01-03.
2016-01-15 21:05:45 +08:00
### Exercises
2022-09-06 23:25:59 +08:00
1. We forgot to draw the relationship between `weather` and `airports` in @fig-flights-relationships.
What is the relationship and how should it appear in the diagram?
2016-08-12 21:58:32 +08:00
2. `weather` only contains information for the three origin airports in NYC.
If it contained weather records for all airports in the USA, what additional connection would it make to `flights`?
2022-09-06 23:25:59 +08:00
3. The `year`, `month`, `day`, `hour`, and `origin` variables almost form a compound key for `weather`, but there's one hour that has duplicate observations.
Can you figure out what's special about that hour?
4. We know that some days of the year are special and fewer people than usual fly on them.
How might you represent that data as a data frame?
What would be the primary key?
How would it connect to the existing data frames?
5. Draw a diagram illustrating the connections between the `Batting`, `People`, and `Salaries` data frames in the Lahman package.
Draw another diagram that shows the relationship between `People`, `Managers`, `AwardsManagers`.
How would you characterise the relationship between the `Batting`, `Pitching`, and `Fielding` data frames?
2016-01-15 21:05:45 +08:00
2022-09-06 23:25:59 +08:00
## Basic joins {#sec-mutating-joins}
Now that you understand how data frames are connected via keys, we can start using joins to better understand the `flights` dataset.
dplyr provides six join functions: `left_join()`, `inner_join()`, `right_join()`, `semi_join()`, and `anti_join()`.
They all have the same interface: they take a pair of data frames `x` and `y` and return a data frame.
The order of the rows and columns in the output is primarily determined by `x`.
2022-09-06 23:25:59 +08:00
In this section, you'll learn how to use one mutating join, `left_join()`, and two filtering joins, `semi_join()` and `anti_join()`.
In the next section, you'll learn exactly how these functions work, and about the remaining `inner_join()`, `right_join()` and `full_join()`.
2022-09-06 23:25:59 +08:00
### Mutating joins
A **mutating join** allows you to combine variables from two data frames: it first matches observations by their keys, then copies across variables from one data frame to the other.
Like `mutate()`, the join functions add variables to the right, so if your dataset has many variables, you won't see the new ones.
For these examples, we'll make it easier to see what's going on by creating a narrower dataset:
2022-09-06 23:25:59 +08:00
```{r}
flights2 <- flights |>
select(year, time_hour, origin, dest, tailnum, carrier)
flights2
```
(Remember that in RStudio you can also use `View()` to avoid this problem.)
There are four types of mutating join, but there's one that you'll use almost all of the time: `left_join()`.
It's special because the output will always have the same rows as `x`[^joins-1].
2022-09-06 23:25:59 +08:00
The primary use of `left_join()` is to add in additional metadata.
For example, we can use `left_join()` to add the full airline name to the `flights2` data:
[^joins-1]: That's not 100% true, but you'll get a warning whenever it isn't.
2022-09-06 23:25:59 +08:00
```{r}
flights2 |>
left_join(airlines)
```
Or we could find out the temperature and wind speed when each plane departed:
```{r}
flights2 |>
left_join(weather |> select(origin, time_hour, temp, wind_speed))
```
Or what size of plane was flying:
2022-09-06 23:25:59 +08:00
```{r}
flights2 |>
left_join(planes |> select(tailnum, type, engines, seats))
2022-09-06 23:25:59 +08:00
```
When `left_join()` fails to find a match for a row in `x`, it fills in the new variables with missing values.
For example, there's no information about the plane with `N3ALAA` so the `type`, `engines`, and `seats` will be missing:
```{r}
flights2 |>
filter(tailnum == "N3ALAA") |>
left_join(planes |> select(tailnum, type, engines, seats))
```
We'll come back to this problem a few times in the rest of the chapter.
2022-09-06 23:25:59 +08:00
### Specifying join keys
By default, `left_join()` will use all variables that appear in both data frames as the join key, the so called **natural** join.
This is a useful heuristic, but it doesn't always work.
For example, what happens if we try to join `flights2` with the complete `planes`?
2022-09-06 23:25:59 +08:00
```{r}
flights2 |>
left_join(planes)
```
We get a lot of missing matches because our join is trying to use both `tailnum` and `year`.
Both `flights` and `planes` have a `year` column but they mean different things: `flights$year` is year the flight occurred and `planes$year` is the year the plane was built.
We only want to join on `tailnum` so we need to provide an explicit specification with `join_by()`:
2022-09-06 23:25:59 +08:00
```{r}
flights2 |>
left_join(planes, join_by(tailnum))
```
Note that the `year` variables are disambiguated in the output with a suffix, which you can control with the `suffix` argument.
2022-09-06 23:25:59 +08:00
`join_by(tailnum)` is short for `join_by(tailnum == tailnum)`.
2022-09-16 21:00:35 +08:00
It's important to know about this fuller form for two reasons.
Firstly, it describes the relationship between the two tables: the keys must be equal.
That's why this type of join is often called an **equi-join**.
You'll learn about non-equi-joins in @sec-non-equi-joins.
Secondly, it's how you specify different join keys in each table.
2022-09-06 23:25:59 +08:00
For example, there are two ways to join the `flight2` and `airports` table: either by `dest` or `origin:`
```{r}
flights2 |>
left_join(airports, join_by(dest == faa))
flights2 |>
left_join(airports, join_by(origin == faa))
```
In older code you might see a different way of specifying the join keys, using a character vector:
- `by = "x"` corresponds to `join_by(x)`.
2022-09-06 23:25:59 +08:00
- `by = c("a" = "x")` corresponds to `join_by(a == x)`.
2022-09-16 21:00:35 +08:00
Now that it exists, we prefer `join_by()` since it provides a clearer and more flexible specification.
2022-09-06 23:25:59 +08:00
### Filtering joins
As you might guess the primary action of a **filtering join** is to filter the rows.
There are two types: semi-joins and anti-joins.
**Semi-joins** keep all rows in `x` that have a match in `y`.
For example, we could use a semi-join to filter the `airports` dataset to show just the origin airports:
2022-09-06 23:25:59 +08:00
```{r}
airports |>
semi_join(flights2, join_by(faa == origin))
```
Or just the destinations:
```{r}
airports |>
semi_join(flights2, join_by(faa == dest))
```
**Anti-joins** are the opposite: they return all rows in `x` that don't have a match in `y`.
2022-09-19 05:23:42 +08:00
They're useful for finding missing values that are **implicit** in the data, the topic of @sec-missing-implicit.
Implicitly missing values don't show up as explicit `NA`s but instead only exist as an absence.
2022-09-16 21:00:35 +08:00
For example, we can find rows that should be in `airports` by looking for flights that don't have a matching destination:
2022-09-06 23:25:59 +08:00
```{r}
flights2 |>
2022-09-16 21:00:35 +08:00
anti_join(airports, join_by(dest == faa)) |>
distinct(dest)
2022-09-06 23:25:59 +08:00
```
2022-09-16 21:00:35 +08:00
Or we can find which `tailnum`s are missing from `planes`:
2022-09-06 23:25:59 +08:00
```{r}
flights2 |>
anti_join(planes, join_by(tailnum)) |>
distinct(tailnum)
```
### Exercises
1. Find the 48 hours (over the course of the whole year) that have the worst delays.
2022-09-06 23:25:59 +08:00
Cross-reference it with the `weather` data.
Can you see any patterns?
2. Imagine you've found the top 10 most popular destinations using this code:
2022-09-06 23:25:59 +08:00
```{r}
top_dest <- flights2 |>
count(dest, sort = TRUE) |>
head(10)
```
How can you find all flights to that destination?
3. Does every departing flight have corresponding weather data for that hour?
4. What do the tail numbers that don't have a matching record in `planes` have in common?
2022-09-06 23:25:59 +08:00
(Hint: one variable explains \~90% of the problems.)
5. Add a column to `planes` that lists every `carrier` that has flown that plane.
You might expect that there's an implicit relationship between plane and airline, because each plane is flown by a single airline.
Confirm or reject this hypothesis using the tools you've learned in previous chapters.
2022-09-06 23:25:59 +08:00
6. Add the location of the origin *and* destination (i.e. the `lat` and `lon`) to `flights`.
Is it easier to rename the columns before or after the join?
7. Compute the average delay by destination, then join on the `airports` data frame so you can show the spatial distribution of delays.
Here's an easy way to draw a map of the United States:
```{r}
#| eval: false
airports |>
semi_join(flights, join_by(faa == dest)) |>
ggplot(aes(lon, lat)) +
borders("state") +
geom_point() +
coord_quickmap()
```
You might want to use the `size` or `colour` of the points to display the average delay for each airport.
8. What happened on June 13 2013?
Display the spatial pattern of delays, and then use Google to cross-reference with the weather.
```{r}
#| eval: false
#| include: false
worst <- filter(flights, !is.na(dep_time), month == 6, day == 13)
worst |>
group_by(dest) |>
summarise(delay = mean(arr_delay), n = n()) |>
filter(n > 5) |>
inner_join(airports, by = c("dest" = "faa")) |>
ggplot(aes(lon, lat)) +
borders("state") +
geom_point(aes(size = n, colour = delay)) +
coord_quickmap()
```
## How do joins work?
Now that you've used joins a few times it's time to learn more about how they work, focusing on how each row in `x` matches zero, one, or more rows in `y`.
2022-09-08 01:30:35 +08:00
We'll begin by using @fig-join-setup to introduce a visual representation of the two simple tibbles defined below.
In these examples we'll use a single key called `key` and a single value column (`val_x` and `val_y)`, but the ideas all generalize to multiple keys and multiple values.
```{r}
2016-10-07 21:16:09 +08:00
x <- tribble(
~key, ~val_x,
1, "x1",
2, "x2",
3, "x3"
)
y <- tribble(
~key, ~val_y,
1, "y1",
2, "y2",
4, "y3"
)
```
2022-08-31 23:06:56 +08:00
```{r}
#| label: fig-join-setup
#| echo: false
#| out-width: ~
#| fig-cap: >
#| Graphical representation of two simple tables. The coloured `key`
#| columns map background colour to key value. The grey columns represent
#| the "value" columns that are carried along for the ride.
2022-08-31 23:06:56 +08:00
#| fig-alt: >
#| x and y are two data frames with 2 columns and 3 rows each. The first
#| column in each is the key and the second is the value. The contents of
2022-09-08 01:30:35 +08:00
#| these data frames are given in the previous code chunk.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/setup.png", dpi = 270)
```
@fig-join-setup2 shows all potential matches between `x` and `y` with an intersection of a pair of lines.
The rows and columns in the output are primarily determined by `x`, so the `x` table is horizontal and lines up with the output.
```{r}
2022-08-31 23:06:56 +08:00
#| label: fig-join-setup2
#| echo: false
2022-08-31 23:06:56 +08:00
#| out-width: ~
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| To understand how joins work, it's useful to think of every possible
#| match. Here we show that by drawing a grid of connecting lines.
#| fig-alt: >
2022-09-08 01:30:35 +08:00
#| x and y are placed at right-angles, with horizonal lines extending
#| from x and vertical lines extending from y. There are 3 rows in x and
#| 3 rows in y leading to 9 intersections that represent nine potential
#| matches.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/setup2.png", dpi = 270)
```
2022-08-31 23:06:56 +08:00
In an actual join, matches will be indicated with dots, as in @fig-join-inner.
2022-09-08 01:30:35 +08:00
The number of dots equals the number of matches, which in turn equals the number of rows in the output, a new data frame that contains the key, the x values, and the y values.
The join shown here is a so-called **inner join**, where rows match if the keys are equal, so that the output contains only the rows with keys that appear in both `x` and `y`.
2016-01-14 23:51:33 +08:00
```{r}
2022-08-31 23:06:56 +08:00
#| label: fig-join-inner
#| echo: false
2022-08-31 23:06:56 +08:00
#| out-width: ~
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| An inner join matches rows in `x` to rows in `y` that have the
#| same value of `key`. Each match becomes a row in the output.
#| fig-alt: >
2022-09-08 01:30:35 +08:00
#| Keys 1 and 2 appear in both x and y, so there values are equal and
#| we get a match, indicated by a dot. Each dot corresponds to a row
#| in the output, so the resulting joined data frame has two rows.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/inner.png", dpi = 270)
2016-01-14 23:51:33 +08:00
```
An **outer join** keeps observations that appear in at least one of the data frames.
2022-08-31 23:06:56 +08:00
These joins work by adding an additional "virtual" observation to each data frame.
This observation has a key that matches if no other key matches, and values filled with `NA`.
There are three types of outer joins:
2022-08-31 23:06:56 +08:00
- A **left join** keeps all observations in `x`, @fig-join-left.
2022-09-08 01:30:35 +08:00
Every row of `x` is preserved in the output because it can fall back to matching a row of `NA`s in `y`.
2022-08-31 23:06:56 +08:00
```{r}
#| label: fig-join-left
#| echo: false
#| out-width: ~
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| A visual representation of the left join where row in `x` appears
#| in the output.
2022-08-31 23:06:56 +08:00
#| fig-alt: >
2022-09-08 01:30:35 +08:00
#| Compared to the inner join, the `y` table gets a new virtual row
#| that will match any row in `x` that doesn't otherwise have a match.
#| This means that the output now has three rows. For key = 3, which
#| matches this virtual row, the value of val_y is NA.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/left.png", dpi = 270)
```
2022-08-31 23:06:56 +08:00
- A **right join** keeps all observations in `y`, @fig-join-right.
2022-09-08 01:30:35 +08:00
Every row of `y` is preserved in the output because it can fall back to matching a row of `NA`s in `x`.
Note the output will consist of all `x` rows that match a row in `y` followed by all rows of `y` that didn't match in `x`.
2016-01-13 23:06:43 +08:00
2022-08-31 23:06:56 +08:00
```{r}
#| label: fig-join-right
#| echo: false
#| out-width: ~
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| A visual representation of the right join where every row of `y`
#| appears in the output.
2022-08-31 23:06:56 +08:00
#| fig-alt: >
#| Keys 1 and 2 from x are matched to those in y, key 4 is
2022-09-08 01:30:35 +08:00
#| also carried along to the joined result since it's on the right data
#| frame, but key 3 from x is not carried along since it's on the left
#| but not on the right. The result is a data frame with 3 rows: keys
#| 1, 2, and 4, all values from val_y, and the corresponding values
#| from val_x for keys 1 and 2 with an NA for key 4, val_x.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/right.png", dpi = 270)
```
- A **full join** keeps all observations that appear in `x` or `y`, @fig-join-full.
Every row of `x` and `y` is included in the output because both `x` and `y` have a fall back row of `NA`s.
2022-09-08 01:30:35 +08:00
Note the output will consist of all `x` rows followed by the remaining `y` rows.
2022-08-31 23:06:56 +08:00
```{r}
#| label: fig-join-full
#| echo: false
#| out-width: ~
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| A visual representation of the full join where every row in `x`
#| and `y` appears in the output.
2022-08-31 23:06:56 +08:00
#| fig-alt: >
#| The result has 4 rows: keys 1, 2, 3, and 4 with all values
2022-09-08 01:30:35 +08:00
#| from val_x and val_y, however key 2, val_y and key 4, val_x are NAs
#| since those keys don't have a match in the other data frames.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/full.png", dpi = 270)
```
Another way to show how the outer joins differ is with a Venn diagram, as in @fig-join-venn.
However, this is not a great representation because while it might jog your memory about which rows are preserved, it fails to illustrate what's happening with the columns.
```{r}
2022-08-31 23:06:56 +08:00
#| label: fig-join-venn
#| echo: false
2022-08-31 23:06:56 +08:00
#| out-width: ~
#| fig-cap: >
#| Venn diagrams showing the difference between inner, left, right, and
#| full joins.
#| fig-alt: >
2022-09-08 01:30:35 +08:00
#| Venn diagrams for inner, full, left, and right joins. Each join
#| represented with two intersecting circles representing data frames x
#| and y, with x on the right and y on the left. Shading indicates the
#| result of the join.
#|
#| Inner join: Only intersection is shaded.
#| Full join: Everything is shaded.
#| Left join: All of x is shaded.
#| Right: All of y is shaded.
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/venn.png", dpi = 270)
```
2022-09-08 01:30:35 +08:00
### Row matching
2022-09-02 06:27:59 +08:00
So far we've explored what happens if a row in `x` matches zero or one rows in `y`.
What happens if it matches more than one row?
To understand what's going let's first narrow our focus to the `inner_join()` and then draw a picture, @fig-join-match-types.
2022-08-31 23:06:56 +08:00
```{r}
2022-09-06 23:25:59 +08:00
#| label: fig-join-match-types
#| echo: false
#| out-width: ~
#| fig-cap: >
#| The three ways a row in `x` can match. `x1` matches
2022-09-06 23:25:59 +08:00
#| one row in `y`, `x2` matches two rows in `y`, `x3` matches
#| zero rows in y. Note that while there are three rows in
#| `x` and three rows in the output, there isn't a direct
#| correspondence between the rows.
#| fig-alt: >
2022-09-08 01:30:35 +08:00
#| A join diagram where x has key values 1, 2, and 3, and y has
#| key values 1, 2, 2. The output has three rows because key 1 matches
#| one row, key 2 matches two rows, and key 3 matches zero rows.
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
knitr::include_graphics("diagrams/join/match-types.png", dpi = 270)
2022-09-02 06:27:59 +08:00
```
There are three possible outcomes for a row in `x`:
- If it doesn't match anything, it's dropped.
- If it matches 1 row in `y`, it's kept as is.
- If it matches more than 1 row in `y`, it's duplicated once for each match.
In principle, this means that there are no guarantees about the number of rows in the output of an `inner_join()`, compared to the number of rows in `x`.
2022-09-02 06:27:59 +08:00
2022-09-08 01:30:35 +08:00
- There might be fewer rows if some rows in `x` don't match any rows in `y`.
2022-09-06 23:25:59 +08:00
- There might be more rows if some rows in `x` match multiple rows in `y`.
2022-09-08 01:30:35 +08:00
- There might be the same number of rows if every row in `x` matches one row in `y`.
- There might be the same number of rows if some rows don't match any rows, and exactly the same number of rows match two rows in `y`!!
2022-09-02 06:27:59 +08:00
Row expansion is a fundamental property of joins, but it's dangerous because it might happen without you realizing it.
To avoid this problem, dplyr will warn whenever there are multiple matches:
2022-09-02 06:27:59 +08:00
```{r}
2022-09-06 23:25:59 +08:00
df1 <- tibble(key = c(1, 2, 3), val_x = c("x1", "x2", "x3"))
df2 <- tibble(key = c(1, 2, 2), val_y = c("y1", "y2", "y3"))
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
df1 |>
inner_join(df2, join_by(key))
2022-09-02 06:27:59 +08:00
```
This is another reason we recommend `left_join()` --- if it runs without warning, you know that every row of the output corresponds to the same row in `x`.
2022-09-02 06:27:59 +08:00
You can gain further control over row matching with two arguments:
2022-09-02 06:27:59 +08:00
- `unmatched` controls what happens when a row in `x` fails to match any rows in `y`. It defaults to `"drop"` which will silently drop any unmatched rows.
- `multiple` controls what happens when a row in `x` matches more than one row in `y`. For equi-joins, it defaults to `"warn"` which emits a warning message if any rows have multiple matches.
2022-09-02 06:27:59 +08:00
There are two common cases in which you might want to override these defaults: enforcing a one-to-one mapping or deliberately allowing the rows to increase.
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
### One-to-one mapping
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
Both `unmatched` and `multiple` can take value `"error"` which means that the join will fail unless each row in `x` matches exactly one row in `y`:
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
```{r}
#| error: true
df1 <- tibble(x = 1)
df2 <- tibble(x = c(1, 1))
df3 <- tibble(x = 3)
2022-09-06 23:25:59 +08:00
df1 |>
inner_join(df2, join_by(x), unmatched = "error", multiple = "error")
df1 |>
inner_join(df3, join_by(x), unmatched = "error", multiple = "error")
2022-09-02 06:27:59 +08:00
```
Note that `unmatched = "error"` is not useful with `left_join()` because, as described above, every row in `x` has a fallback match to a virtual row in `y`.
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
### Allow multiple rows
2022-09-02 06:27:59 +08:00
2022-09-08 01:30:35 +08:00
Sometimes it's useful to deliberately expand the number of rows in the output.
This can come about naturally if you "flip" the direction of the question you're asking.
2022-09-08 01:30:35 +08:00
For example, as we've seen above, it's natural to supplement the `flights` data with information about the plane that flew each flight:
2022-08-31 23:06:56 +08:00
```{r}
2022-09-06 23:25:59 +08:00
#| results: false
2022-09-02 06:27:59 +08:00
flights2 |>
2022-08-31 23:06:56 +08:00
left_join(planes, by = "tailnum")
```
But it's also reasonable to ask what flights did each plane fly:
2022-09-02 06:27:59 +08:00
```{r}
2022-09-06 23:25:59 +08:00
plane_flights <- planes |>
select(tailnum, type, engines, seats) |>
2022-09-06 23:25:59 +08:00
left_join(flights2, by = "tailnum")
2022-09-02 06:27:59 +08:00
```
Since this duplicates rows in `x` (the planes), we need to explicitly say we're ok with the multiple matches by setting `multiple = "all"`:
2022-08-31 23:06:56 +08:00
```{r}
2022-09-06 23:25:59 +08:00
plane_flights <- planes |>
select(tailnum, type, engines, seats) |>
2022-09-06 23:25:59 +08:00
left_join(flights2, by = "tailnum", multiple = "all")
2022-09-02 06:27:59 +08:00
2022-09-06 23:25:59 +08:00
plane_flights
2022-09-02 06:27:59 +08:00
```
2022-09-06 23:25:59 +08:00
### Filtering joins {#sec-non-equi-joins}
The number of matches also determines the behavior of the filtering joins.
2022-09-08 01:30:35 +08:00
The semi-join keeps rows in `x` that have one or more matches in `y`, as in @fig-join-semi.
The anti-join keeps rows in `x` that don't have a match in `y`, as in @fig-join-anti.
In both cases, only the existence of a match is important; it doesn't matter how many times it matches.
2022-09-06 23:25:59 +08:00
This means that filtering joins never duplicate rows like mutating joins do.
2022-08-29 23:16:18 +08:00
2022-08-31 23:06:56 +08:00
```{r}
2022-09-06 23:25:59 +08:00
#| label: fig-join-semi
2022-08-31 23:06:56 +08:00
#| echo: false
2022-09-06 23:25:59 +08:00
#| out-width: null
#| fig-cap: >
#| In a semi-join it only matters that there is a match; otherwise
#| values in `y` don't affect the output.
2022-08-31 23:06:56 +08:00
#| fig-alt: >
2022-09-06 23:25:59 +08:00
#| Diagram of a semi join. Data frame x is on the left and has two columns
#| (key and val_x) with keys 1, 2, and 3. Diagram y is on the right and also
#| has two columns (key and val_y) with keys 1, 2, and 4. Semi joining these
#| two results in a data frame with two rows and two columns (key and val_x),
#| with keys 1 and 2 (the only keys that match between the two data frames).
2022-08-31 23:06:56 +08:00
2022-09-08 01:30:35 +08:00
knitr::include_graphics("diagrams/join/semi.png", dpi = 270)
2022-08-31 23:06:56 +08:00
```
2016-01-30 21:13:46 +08:00
2022-08-31 23:06:56 +08:00
```{r}
2022-09-06 23:25:59 +08:00
#| label: fig-join-anti
#| echo: false
#| out-width: null
#| fig-cap: >
#| An anti-join is the inverse of a semi-join, dropping rows from `x`
#| that have a match in `y`.
#| fig-alt: >
#| Diagram of an anti join. Data frame x is on the left and has two columns
#| (key and val_x) with keys 1, 2, and 3. Diagram y is on the right and also
#| has two columns (key and val_y) with keys 1, 2, and 4. Anti joining these
#| two results in a data frame with one row and two columns (key and val_x),
#| with keys 3 only (the only key in x that is not in y).
2022-08-31 23:06:56 +08:00
2022-09-06 23:25:59 +08:00
knitr::include_graphics("diagrams/join/anti.png", dpi = 270)
2022-08-31 23:06:56 +08:00
```
2022-09-06 23:25:59 +08:00
## Non-equi joins
2016-01-30 21:13:46 +08:00
So far you've only seen **equi-joins**, joins where the two rows match if the `x` keys are exactly equal to the `y` keys.
2022-09-08 01:30:35 +08:00
Now we're going to relax that restriction and discuss other ways of determining if a pair of rows match.
But before we can do that, we need to revisit a simplification we made above.
In equi-joins the `x` keys and `y` are always equal, so we only need to show one in the output.
2022-09-08 01:30:35 +08:00
We can request that dplyr keep both keys with `keep = TRUE`, leading to the code below and the re-drawn `inner_join()` in @fig-inner-both.
2022-08-29 23:16:18 +08:00
```{r}
x |> left_join(y, by = "key", keep = TRUE)
```
```{r}
#| label: fig-inner-both
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| An inner join showing both `x` and `y` keys in the output.
#| fig-alt: >
#| A join diagram showing an inner join betwen x and y. The result
#| now includes four columns: key.x, val_x, key.y, and val_y. The
#| values of key.x and key.y are identical, which is why we usually
#| omit one.
2022-08-29 23:16:18 +08:00
#| echo: false
2022-08-31 23:06:56 +08:00
#| out-width: ~
2022-08-29 23:16:18 +08:00
knitr::include_graphics("diagrams/join/inner-both.png", dpi = 270)
```
When we move away from equi-joins we'll always show the keys, because the key values will often different.
For example, instead of matching only when the `x$key` and `y$key` are equal, we could match whenever the `x$key` is greater than or equal to the `y$key`, leading to @fig-join-gte.
dplyr's join functions understand this distinction so will always show both keys when you perform a non-equi-join.
2022-08-29 23:16:18 +08:00
```{r}
2022-08-31 23:06:56 +08:00
#| label: fig-join-gte
2022-08-29 23:16:18 +08:00
#| echo: false
#| fig-cap: >
2022-09-08 01:30:35 +08:00
#| A non-equi join where the `x` key must greater than or equal to
#| than the `y` key. Many rows generate multiple matches.
#| fig-alt: >
#| A join diagram illustrating join_by(key >= key). The first row
#| of x matches one row of y and the second and thirds rows each match
#| two rows. This means the output has five rows containing each of the
#| following (key.x, key.y) pairs: (1, 1), (2, 1), (2, 2), (3, 1),
#| (3, 2).
2022-08-31 23:06:56 +08:00
knitr::include_graphics("diagrams/join/gte.png", dpi = 270)
2022-08-29 23:16:18 +08:00
```
Non-equi-join isn't a particularly useful term because it only tells you what the join is not, not what it is. dplyr helps by identifying four particularly useful types of non-equi-join:
- **Cross joins** match every pair of rows.
- **Inequality joins** use `<`, `<=`, `>`, `>=` instead of `==`.
2022-09-08 01:30:35 +08:00
- **Rolling joins** are similar to inequality joins but only find the closest match.
- **Overlap joins** are a special type of inequality join designed to work with ranges.
2022-08-29 23:16:18 +08:00
2022-09-02 06:27:59 +08:00
Each of these is described in more detail in the following sections.
2022-08-29 23:16:18 +08:00
### Cross joins
2022-09-06 23:25:59 +08:00
A cross join matches everything, as in @fig-join-cross, generating the Cartesian product of rows.
2022-09-08 01:30:35 +08:00
This means the output will have `nrow(x) * nrow(y)` rows.
2022-09-06 23:25:59 +08:00
```{r}
#| label: fig-join-cross
#| echo: false
#| out-width: ~
#| fig-cap: >
#| A cross join matches each row in `x` with every row in `y`.
2022-09-08 01:30:35 +08:00
#| fig-alt: >
#| A join diagram showing a dot for every combination of x and y.
2022-09-06 23:25:59 +08:00
knitr::include_graphics("diagrams/join/cross.png", dpi = 270)
```
Cross joins are useful when generating permutations.
2022-09-08 01:30:35 +08:00
For example, the code below generates every possible pair of names.
Since we're joining `df` to itself, this is sometimes called a **self-join**.
2022-08-29 23:16:18 +08:00
```{r}
2022-09-08 01:30:35 +08:00
df <- tibble(name = c("John", "Simon", "Tracy", "Max"))
df |> left_join(df, join_by())
2022-08-29 23:16:18 +08:00
```
2022-09-08 01:30:35 +08:00
### Inequality joins
Inequality joins use `<`, `<=`, `>=`, or `>` to restrict the set of possible matches, as in @fig-join-gte and @fig-join-lt.
2022-08-29 23:16:18 +08:00
2022-09-06 23:25:59 +08:00
```{r}
#| label: fig-join-lt
2022-09-06 23:25:59 +08:00
#| echo: false
#| out-width: ~
#| fig-cap: >
#| An inequality join where `x` is joined to `y` on rows where the key
#| of `x` is less than the key of `y`.
knitr::include_graphics("diagrams/join/lt.png", dpi = 270)
2022-09-06 23:25:59 +08:00
```
2022-09-08 01:30:35 +08:00
Inequality joins are extremely general, so general that it's hard to come up with meaningful specific use cases.
One small useful technique is to use them to restrict the cross join so that instead of generating all permutations, we generate all combinations:
2022-09-08 01:30:35 +08:00
```{r}
df <- tibble(id = 1:4, name = c("John", "Simon", "Tracy", "Max"))
df |> left_join(df, join_by(id < id))
```
2022-08-29 23:16:18 +08:00
### Rolling joins
2022-09-10 21:38:12 +08:00
Rolling joins are a special type of inequality join where instead of getting *every* row that satisfies the inequality, you get just the closest row, as in @fig-join-closest.
You can turn any inequality join into a rolling join by adding `closest()`.
For example `join_by(closest(x <= y))` matches the smallest `y` that's greater than or equal to x, and `join_by(closest(x > y))` matches the biggest `y` that's less than `x`.
2022-09-06 23:25:59 +08:00
2022-08-31 23:06:56 +08:00
```{r}
2022-09-08 01:30:35 +08:00
#| label: fig-join-closest
2022-08-31 23:06:56 +08:00
#| echo: false
#| out-width: ~
#| fig-cap: >
#| A following join is similar to a greater-than-or-equal inequality join
#| but only matches the first value.
2022-09-08 01:30:35 +08:00
knitr::include_graphics("diagrams/join/closest.png", dpi = 270)
2022-08-31 23:06:56 +08:00
```
2022-09-08 01:30:35 +08:00
Rolling joins are particularly useful when you have two tables of dates that don't perfectly line up and you want to find (e.g.) the closest date in table 1 that comes before (or after) some date in table 2.
For example, imagine that you're in charge of the party planning commission for your office.
2022-09-08 01:30:35 +08:00
Your company is rather cheap so instead of having individual parties, you only have a party once each quarter.
The rules for determining when a party will be held are a little complex: parties are always on a Monday, you skip the first week of January since a lot of people are on holiday, and the first Monday of Q3 2022 is July 4, so that has to be pushed back a week.
2022-08-31 23:07:10 +08:00
That leads to the following party days:
```{r}
parties <- tibble(
q = 1:4,
party = lubridate::ymd(c("2022-01-10", "2022-04-04", "2022-07-11", "2022-10-03"))
)
```
Now imagine that you have a table of employee birthdays:
2022-08-31 23:07:10 +08:00
```{r}
employees <- tibble(
name = wakefield::name(100),
birthday = lubridate::ymd("2022-01-01") + (sample(365, 100, replace = TRUE) - 1)
)
employees
```
And for each employee we want to find the first party date that comes after (or on) their birthday.
We can express that with a rolling join:
2022-08-31 23:07:10 +08:00
```{r}
2022-09-02 06:27:59 +08:00
employees |>
left_join(parties, join_by(closest(birthday >= party)))
2022-09-06 23:25:59 +08:00
```
2022-09-02 06:27:59 +08:00
2022-08-29 23:16:18 +08:00
### Overlap joins
2022-09-08 01:30:35 +08:00
Overlap joins provide three helpers that use inequality joins to make it easier to work with intervals:
- `between(x, y_lower, y_upper)` is short for `x >= y_lower, x <= y_upper`.
- `within(x_lower, x_upper, y_lower, y_upper)` is short for `x_lower >= y_lower, x_upper <= y_upper`.
- `overlaps(x_lower, x_upper, y_lower, y_upper)` is short for `x_lower <= y_upper, x_upper >= y_lower`.
Let's continue the birthday example to see how you might use them.
There's one problem with the strategy we used above: there's no party preceding the birthdays Jan 1-9.
So it might be better to to be explicit about the date ranges that each party spans, and make a special case for those early birthdays:
2022-08-31 23:07:10 +08:00
```{r}
parties <- tibble(
q = 1:4,
party = lubridate::ymd(c("2022-01-10", "2022-04-04", "2022-07-11", "2022-10-03")),
start = lubridate::ymd(c("2022-01-01", "2022-04-04", "2022-07-11", "2022-10-03")),
end = lubridate::ymd(c("2022-04-03", "2022-07-11", "2022-10-02", "2022-12-31"))
)
parties
```
Hadley is hopelessly bad at data entry so he also wanted to check that the party periods don't overlap.
You can perform a self-join and check to see if any start-end interval overlaps with any other:
2022-08-31 23:07:10 +08:00
```{r}
2022-09-08 01:30:35 +08:00
parties |>
inner_join(parties, join_by(overlaps(start, end, start, end), q < q)) |>
select(start.x, end.x, start.y, end.y)
2022-08-31 23:07:10 +08:00
```
Ooops, there is an overlap, so let's fix that problem and continue:
2022-08-31 23:07:10 +08:00
```{r}
2022-09-08 01:30:35 +08:00
parties <- tibble(
q = 1:4,
party = lubridate::ymd(c("2022-01-10", "2022-04-04", "2022-07-11", "2022-10-03")),
start = lubridate::ymd(c("2022-01-01", "2022-04-04", "2022-07-11", "2022-10-03")),
end = lubridate::ymd(c("2022-04-03", "2022-07-10", "2022-10-02", "2022-12-31"))
)
2022-08-31 23:07:10 +08:00
```
2022-08-29 23:16:18 +08:00
2022-09-08 01:30:35 +08:00
Now we can match each employee to their party.
This is a good place to use `unmatched = "error"` because I want to quickly find out if any employees didn't get assigned a party.
2022-08-29 23:16:18 +08:00
```{r}
2022-09-08 01:30:35 +08:00
employees |>
inner_join(parties, join_by(between(birthday, start, end)), unmatched = "error")
2022-08-29 23:16:18 +08:00
```
### Exercises
1. Can you explain what's happening with the keys in this equi-join?
2022-09-06 23:25:59 +08:00
Why are they different?
2022-08-29 23:16:18 +08:00
```{r}
x |> full_join(y, by = "key")
x |> full_join(y, by = "key", keep = TRUE)
```
2022-09-08 01:30:35 +08:00
2. When finding if any party period overlapped with another party period I used `q < q` in the `join_by()`?
Why?
What happens if you remove this inequality?
## Summary
In this chapter, you've learned how to use mutating and filtering joins to combine data from a pair of data frames.
Along the way you learned how to identify keys, and the difference between primary and foreign keys.
You also understand how joins work and how to figure out how many rows the output will have.
Finally, you've gained a glimpse into the power of non-equi-joins and seen a few interesting use cases.
This chapter concludes the "Transform" part of the book where the focus was on the tools you could use with individual columns and tibbles.
You learned about dplyr and base functions for working with logical vectors, numbers, and complete tables, stringr functions for working strings, lubridate functions for working with date-times, and forcats functions for working with factors.
In the next part of the book, you'll learn more about getting various types of data into R in a tidy form.