2

私はベクトルからこのような行列を作成しようとしています：

    vec= c(2, 5, 9)
    > A
            [,1] [,2] [,3] [,4]
    [1,]     2    0    0    0
    [2,]     5    3    0    0
    [3,]     9    7    4    0

実際には、常に最初の列はベクトル要素であり、2番目の列は0で始まり、次に（5-2 = 3）、次に2番目の列の3番目の要素は（9-2 = 7）です。次に、3番目の列は0で始まり、次に0と（9-5 = 4）で始まり、最後の列は常にゼロです。vecの変更の長さは、たとえば4、5、...などの任意の数になります。この行列を作成するための効率的な関数またはコードを作成するにはどうすればよいですか。

4

2 に答える 2

7

効率についてはわかりませんが、forループを使用しない2つの解決策を次に示します。

n <- length(vec)    
A <- replicate(n+1, vec) - cbind(0, t(replicate(n, vec)))
A[upper.tri(A)] <- 0

これは長いですが、1つのマトリックスしか作成しません

n <- length(vec)
A <- replicate(n, vec)
A <- A - t(A)
A <- cbind(vec, A)
A[upper.tri(A)] <- 0
于 2012-06-01T05:25:00.573 に答える
5

私はこれがあなたが望むことをするだろうと思います：

f = function(vec)
{
   n = length(vec)  
   M = matrix(0,n,n+1)
   M[,1] = vec
   for(i in 1:n) M[,i+1] = c(rep(0,i),vec[-c(1:i)]-vec[i])
   return(M)
}

vec = c(2,5,9)
f(vec)
     [,1] [,2] [,3] [,4]
[1,]    2    0    0    0
[2,]    5    3    0    0
[3,]    9    7    4    0
于 2012-06-01T02:45:04.050 に答える