r - 连续一些位置的总和 - R

Question

考虑一个通用矩阵 nxn - 例如:

A <- matrix(1:16, nrow = 4, ncol = 4)

如何计算“右下”三角形中行的总和并在向量中显示信息？

Answer 1

找到封装所需数据的掩码:

> mask <- apply(lower.tri(A, diag = FALSE), 1, rev)
> mask
      [,1]  [,2]  [,3]  [,4]
[1,] FALSE FALSE FALSE FALSE
[2,] FALSE FALSE FALSE  TRUE
[3,] FALSE FALSE  TRUE  TRUE
[4,] FALSE  TRUE  TRUE  TRUE

乘以这个掩码并计算总和:

> A * mask
     [,1] [,2] [,3] [,4]
[1,]    0    0    0    0
[2,]    0    0    0   14
[3,]    0    0   11   15
[4,]    0    8   12   16

>  rowSums(A * mask)
[1]  0 14 26 36