2017-12-09 6 views
6

考慮し、一般的な行列のn×n - 例えば:和 - R

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

がどのように「右下」の三角形の行の合計を計算するとの情報を表示することができベクター?

example

答えて

1

ここに解決策があります。キーは、lower.triまたはupper.tri関数を使用し、applyを使用して各行または列を並べ替えることです。最後に、各行の合計を計算します。

# Create example data frame 
m <- matrix(1:16, 4) 
m 
#  [,1] [,2] [,3] [,4] 
# [1,] 1 5 9 13 
# [2,] 2 6 10 14 
# [3,] 3 7 11 15 
# [4,] 4 8 12 16 

# Calculate the sum 
rowSums(m * apply(lower.tri(m), 1, sort)) 
# [1] 0 14 26 36 

rowSums(m * apply(upper.tri(m), 2, sort)) 
# [1] 0 14 26 36 
2
A = matrix(1:16, 4) 
A 
#  [,1] [,2] [,3] [,4] 
# [1,] 1 5 9 13 
# [2,] 2 6 10 14 
# [3,] 3 7 11 15 
# [4,] 4 8 12 16 

sapply(1:NROW(A), function(i) sum(tail(A[i,], i - 1))) 
#[1] 0 14 26 36 
6

データが必要なカプセル化のマスク検索:

> 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 

乗にこのマスクをし、合計を計算する:ちょうど `rowSums(A * lower.triについて

> 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 
+5

方法を(A)[、ncol(A):1]) ' –