提取行和列的最大值

时间:2017-05-23 09:29:03

标签: r

我需要查找行max和列max对于相同位置的值。测试数据(真实数据不需要是方阵):

scores<-structure(c(0.4, 0.6, 0.222222222222222, 0.4, 0.4, 0, 0.25, 0.5, 
0.285714285714286), .Dim = c(3L, 3L), .Dimnames = list(c("a", 
"b", "c"), c("d", "e", "f")))

我已经找到了具有该行/列的最大值的列/行。

rows<-structure(list(a = c("d", "e"), b = "d", c = "f"), .Names = c("a", 
"b", "c"))
cols<-structure(list(d = "b", e = c("a", "b"), f = "b"), .Names = c("d", 
"e", "f"))

但我没有设法从矩阵中获取值。问题是当相同(最大)值出现两次或更多时。我不知道如何检查该案件的指数。我尝试使用mapply:

mapply(function(x, y) {
    cols[x] == rows[y] 
    }, rows, cols)

但是当行或列有多个元素时,这会停止。

预期输出:c(0.6, 0.4)
第一个是第1列和第2行的最大值,第二个值是第1行和第2列的最大值。

            d   e         f   | Max
a   0.4000000 0.4 0.2500000   0.4
b   0.6000000 0.4 0.5000000   0.6
c   0.2222222 0.0 0.2857143   0.2857
Max:  0.6     0.4 0.5

正如您在第2行和第1列中看到的那样,最大值是相同的,对于第1行和第1列,它们是相同的值,但对于第3行和第3列,它不是&#39; t

2 个答案:

答案 0 :(得分:2)

我想,我明白你要做什么。虽然不是最佳解决方案。

我们在行和列中找出最大值的索引,然后找出intersect的索引并显示数据帧中的相应值。

a1 <- which(apply(scores, 1, function(x) x == max(x)))
a2 <- which(apply(scores, 2, function(x) x == max(x)))
scores[intersect(a1, a2)]

#[1] 0.6 0.4

并在一行

scores[intersect(which(apply(scores, 1, function(x) x == max(x))), 
                 which(apply(scores, 2, function(x) x == max(x))))]

答案 1 :(得分:2)

这就是你想要的:

# Compute rows and columns max and rows max positions
row_max<-apply(scores, 1, max)
row_max_pos<-apply(scores, 1, which.max)
col_max<-apply(scores, 2, max)

# For each row, check if max is equal to corresponding column max
res <- sapply(1:length(row_max), 
              function(i) ifelse(row_max[i] == col_max[row_max_pos[i]], T, F))
row_max[res]

它也可以在多行/列上使用相同的最大值,例如使用此数据:

scores <- structure(c(0.4, 0.6, 0.222222222222222, 0.4, 0.4, 0, 0.25, 0.5, 
                      0.285714285714286, 0.13, 0.2, 0.6), .Dim = c(4L, 3L), 
                      .Dimnames = list(c("a", "b", "c", "d"), c("e", "f", "g")))