4

我有一个矩阵,其中的列重复字符列名。

set.seed(1)
m <- matrix(sample(1:10,12,replace=TRUE), nrow = 3, ncol = 4, byrow = TRUE,
       dimnames = list(c("s1", "s2", "s3"),c("x", "y","x","y")))

m
   x y  x  y
s1 3 4  6 10
s2 3 9 10  7
s3 7 1  3  2

我需要将具有相同列名的所有列汇总为一列,即

m <- matrix(c(9,14,13,16,10,3), nrow = 3, ncol = , byrow = TRUE,dimnames = list(c("s1", "s2", "s3"),c("x", "y")))

    x  y
s1  9 14
s2 13 16
s3 10  3

我玩过聚合函数中的简单总和,但没有任何运气。有什么建议吗?谢谢。

4

2 回答 2

8

好吧,这个解决方案不会因为代码的透明性而赢得任何奖项,但我更喜欢它:

nms <- colnames(m)
m %*% sapply(unique(nms),"==", nms)
#     x  y
# s1  9 14
# s2 13 16
# s3 10  3

它通过构造一个矩阵来工作,该矩阵形成m' 列的适当线性组合。要查看它是如何工作的,请将第二行分成两个分量矩阵,使用 相乘%*%,如下所示:

 -          -      -   -
|  3 4  6 10 |    | 1 0 |
|  3 9 10  7 |    | 0 1 |
|  7 1  3  2 |    | 1 0 |
 -          -     | 0 1 |
                   -   -
于 2012-07-16T21:26:30.303 回答
6
nms <- colnames(m)
sapply(unique(nms), function(i)rowSums(m[, nms==i]))

   x  y
s1  9 14
s2 13 16
s3 10  3
于 2012-07-16T21:24:43.580 回答