我有这样的数据:
x <- matrix(c(2,2,3,3,3,4,4,20,33,2,3,45,6,9,45,454,7,4,6,7,5), nrow = 7, ncol = 3)
在真实数据集中,我有一个包含大量列的巨大矩阵。 我想提取相对于第一列(Id)和第三列最小值的唯一行。例如,对于这个矩阵,我希望
y <- matrix(c(2,3,4,20,3,9,45,4,5), nrow = 3, ncol = 3)
我尝试了很多东西,但我无法理解。 任何帮助表示赞赏。
提前致谢, Zeray
答案 0 :(得分:5)
这是一个更复杂的版本,但比Chase的ddply解决方案更快 - 大约快200倍: - )
uniqueMin <- function(m, idCol = 1L, minCol = ncol(m)) {
t(vapply(split(1:nrow(m), m[,idCol]), function(i, x, minCol) x[i, , drop=FALSE][which.min(x[i,minCol]),], m[1,], x=m, minCol=minCol))
}
以下测试代码:
nRows <- 10000
nCols <- 100
ids <- nRows/5
m <- cbind(sample(ids, nRows, T), matrix(runif(nRows*nCols), nRows))
system.time( a<-uniqueMin(m, minCol=3L) ) # 0.07
system.time(ddply(as.data.frame(m), "V1", function(x) x[which.min(x$V3) ,])) # 15.72
答案 1 :(得分:2)
您可以使用package plyr。转换为data.frame,以便您可以对第一列进行分组,然后使用which.min按组提取最小行:
library(plyr)
ddply(as.data.frame(x), "V1", function(x) x[which.min(x$V3) ,])
V1 V2 V3
1 2 20 45
2 3 3 4
3 4 9 5