列方向「外積」

Question

は私が20 x 10行列Xを持っている、と私は1:10にi,jためA[i,j,] = X[,i] * X[,j]と10 x 10 x 20配列Aを作成したいです。ここでforループを使用して、簡単な解決策は以下のとおりです。

X <- matrix(seq(10*20),20,10) 
A <- array(dim=c(10,10,20)) 
for (i in 1:10) 
    for (j in 1:10) 
    A[i,j,] <- X[,i] * X[,j] 

どのように私はベクトル化でこれを実装するのですか？

ありがとうございました。

Answer 1

outerを使用してください、これは何ですか？

A2 = array(apply(X, 1L, function (x) outer(x,x)), dim=c(10,10,20)) 

# test 
X <- matrix(seq(10*20),20,10) 
A <- array(dim=c(10,10,20)) 
for (i in 1:10){ 
    for (j in 1:10){ 
    A[i,j,] <- X[,i] * X[,j] 
    } 
} 

all.equal(A, A2) # TRUE 

EDIT：問題の質問の例の後更新答えが変更されました。