从数据框中提取重复行

4

我有一个大型数据框正在处理,前几行如下:

      Assay   Genotype   Sample    Result
1     001        G         1         0
2     001        A         2         1
3     001        G         3         0 
4     001        NA        1         NA
5     002        T         1         0
6     002        G         2         1
7     002        T         2         0 
8     002        T         4         0
9     003        NA        1         NA

总共我将处理2000个样本,每个样本有168个测定。

我想提取出其中同一测定和样本都有多个条目的行。我希望结果数据以数据框的形式呈现,包含所有重复的条目,并按照重复项相邻排序。从上面的示例中,结果将如下所示:

      Assay   Genotype   Sample    Result
1     001        G         1         0
4     001        NA        1         NA
6     002        G         2         1
7     002        T         2         0 
1个回答

6

用于简便加载的演示数据:

df <- structure(list(Assay = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L), Genotype = structure(c(2L, 1L, 2L, NA, 3L, 2L, 3L, 3L, NA), .Label = c("A", "G", "T"), class = "factor"), Sample = c(1L, 2L, 3L, 1L, 1L, 2L, 2L, 4L, 1L), Result = c(0L, 1L, 0L, NA, 0L, 1L, 0L, 0L, NA)), .Names = c("Assay", "Genotype", "Sample", "Result"), class = "data.frame", row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9"))

您可以使用duplicated轻松获取重复的测定/样品对:
vars <- c('Assay', 'Sample')
dup <- df[duplicated(x[, vars]), vars]

导致:
> dup
  Assay Sample
4     1      1
7     2      2

需要简单合并以获取所需结果:merge
> merge(dup, df)
  Assay Sample Genotype Result
1     1      1     <NA>     NA
2     1      1        G      0
3     2      2        G      1
4     2      2        T      0

网页内容由stack overflow 提供, 点击上面的
可以查看英文原文,
原文链接