Mutating joins add columns from y to x, matching observations based on
the keys. There are four mutating joins: the inner join, and the three outer
joins.
Inner join
An inner_join() only keeps observations from x that have a matching key
in y.
The most important property of an inner join is that unmatched rows in either input are not included in the result. This means that generally inner joins are not appropriate in most analyses, because it is too easy to lose observations.
Outer joins
The three outer joins keep observations that appear in at least one of the data frames:
A
left_join()keeps all observations inx.A
right_join()keeps all observations iny.A
full_join()keeps all observations inxandy.
Usage
# S3 method for class 'Seurat'
left_join(x, y, by = NULL, copy = FALSE, suffix = c(".x", ".y"), ...)Arguments
- x, y
A pair of data frames, data frame extensions (e.g. a tibble), or lazy data frames (e.g. from dbplyr or dtplyr). See Methods, below, for more details.
- by
A join specification created with
join_by(), or a character vector of variables to join by.If
NULL, the default,*_join()will perform a natural join, using all variables in common acrossxandy. A message lists the variables so that you can check they're correct; suppress the message by supplyingbyexplicitly.To join on different variables between
xandy, use ajoin_by()specification. For example,join_by(a == b)will matchx$atoy$b.To join by multiple variables, use a
join_by()specification with multiple expressions. For example,join_by(a == b, c == d)will matchx$atoy$bandx$ctoy$d. If the column names are the same betweenxandy, you can shorten this by listing only the variable names, likejoin_by(a, c).join_by()can also be used to perform inequality, rolling, and overlap joins. See the documentation at ?join_by for details on these types of joins.For simple equality joins, you can alternatively specify a character vector of variable names to join by. For example,
by = c("a", "b")joinsx$atoy$aandx$btoy$b. If variable names differ betweenxandy, use a named character vector likeby = c("x_a" = "y_a", "x_b" = "y_b").To perform a cross-join, generating all combinations of
xandy, seecross_join().- copy
If
xandyare not from the same data source, andcopyisTRUE, thenywill be copied into the same src asx. This allows you to join tables across srcs, but it is a potentially expensive operation so you must opt into it.- suffix
If there are non-joined duplicate variables in
xandy, these suffixes will be added to the output to disambiguate them. Should be a character vector of length 2.- ...
Other parameters passed onto methods.
Value
An object of the same type as x (including the same groups). The order of
the rows and columns of x is preserved as much as possible. The output has
the following properties:
The rows are affect by the join type.
inner_join()returns matchedxrows.left_join()returns allxrows.right_join()returns matched ofxrows, followed by unmatchedyrows.full_join()returns allxrows, followed by unmatchedyrows.
Output columns include all columns from
xand all non-key columns fromy. Ifkeep = TRUE, the key columns fromyare included as well.If non-key columns in
xandyhave the same name,suffixes are added to disambiguate. Ifkeep = TRUEand key columns inxandyhave the same name,suffixes are added to disambiguate these as well.If
keep = FALSE, output columns included inbyare coerced to their common type betweenxandy.
Many-to-many relationships
By default, dplyr guards against many-to-many relationships in equality joins by throwing a warning. These occur when both of the following are true:
A row in
xmatches multiple rows iny.A row in
ymatches multiple rows inx.
This is typically surprising, as most joins involve a relationship of one-to-one, one-to-many, or many-to-one, and is often the result of an improperly specified join. Many-to-many relationships are particularly problematic because they can result in a Cartesian explosion of the number of rows returned from the join.
If a many-to-many relationship is expected, silence this warning by
explicitly setting relationship = "many-to-many".
In production code, it is best to preemptively set relationship to whatever
relationship you expect to exist between the keys of x and y, as this
forces an error to occur immediately if the data doesn't align with your
expectations.
Inequality joins typically result in many-to-many relationships by nature, so they don't warn on them by default, but you should still take extra care when specifying an inequality join, because they also have the capability to return a large number of rows.
Rolling joins don't warn on many-to-many relationships either, but many
rolling joins follow a many-to-one relationship, so it is often useful to
set relationship = "many-to-one" to enforce this.
Note that in SQL, most database providers won't let you specify a many-to-many relationship between two tables, instead requiring that you create a third junction table that results in two one-to-many relationships instead.
Methods
These functions are generics, which means that packages can provide implementations (methods) for other classes. See the documentation of individual methods for extra arguments and differences in behaviour.
Methods available in currently loaded packages:
inner_join(): dplyr (data.frame), tidyseurat (Seurat) .left_join(): dplyr (data.frame), tidyseurat (Seurat) .right_join(): dplyr (data.frame), tidyseurat (Seurat) .full_join(): dplyr (data.frame), tidyseurat (Seurat) .
See also
Other joins:
cross_join(),
filter-joins,
nest_join()
Examples
data(pbmc_small)
tt <- pbmc_small
tt |> left_join(tt |>
distinct(groups) |>
mutate(new_column=1:2))
#> tidyseurat says: A data frame is returned for independent data analysis.
#> Joining with `by = join_by(groups)`
#> # A Seurat-tibble abstraction: 80 × 16
#> # Features=230 | Cells=80 | Active assay=RNA | Assays=RNA
#> .cell orig.ident nCount_RNA nFeature_RNA RNA_snn_res.0.8 letter.idents groups
#> <chr> <fct> <dbl> <int> <fct> <fct> <chr>
#> 1 ATGC… SeuratPro… 70 47 0 A g2
#> 2 CATG… SeuratPro… 85 52 0 A g1
#> 3 GAAC… SeuratPro… 87 50 1 B g2
#> 4 TGAC… SeuratPro… 127 56 0 A g2
#> 5 AGTC… SeuratPro… 173 53 0 A g2
#> 6 TCTG… SeuratPro… 70 48 0 A g1
#> 7 TGGT… SeuratPro… 64 36 0 A g1
#> 8 GCAG… SeuratPro… 72 45 0 A g1
#> 9 GATA… SeuratPro… 52 36 0 A g1
#> 10 AATG… SeuratPro… 100 41 0 A g1
#> # ℹ 70 more rows
#> # ℹ 9 more variables: RNA_snn_res.1 <fct>, new_column <int>, PC_1 <dbl>,
#> # PC_2 <dbl>, PC_3 <dbl>, PC_4 <dbl>, PC_5 <dbl>, tSNE_1 <dbl>, tSNE_2 <dbl>