2

我想计算上/下三角矩阵(对角线)的行/列和。例子:

m <- matrix(1:9, nrow=3)
#result:
upperRowSums(m)
## [1] 12 13 9
lowerRowSums(m)
## [1] 1 7 18

我知道这可以通过简单的for循环来完成,但我想避免这种情况。我正在寻找纯粹的R解决方案。

4

1 回答 1

4

离开:

rowSums(m * upper.tri(m, diag=TRUE))
#[1] 12 13  9
于 2014-06-29T21:03:05.557 回答