我有一个大数据框,我正在使用,前几行如下:
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和Sample的多个条目的行。我希望结果数据位于包含所有重复条目的数据框中,进行排序以使重复项彼此相邻。从上面的例子可以看出结果如下:
Assay Genotype Sample Result
1 001 G 1 0
4 001 NA 1 NA
6 002 G 2 1
7 002 T 2 0
答案 0 :(得分:5)
演示数据以便于加载:
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