R：因子をバイナリ行列に変換する効率的な方法

Question

サイズfactorをfactorのi番目とj番目の要素が1の場合は、（i、j）要素が1である×Nのバイナリに変換したいそれ以外は0です。

以下は私がやりたいことを実装する素朴な方法ですが、このコードはかなり遅いです。同じことをするもっと効率的なやり方はありますか？

size <- 100 
id <- factor(sample(3, size, replace=TRUE)) 
mat <- matrix(0, nrow=size, ncol=size) 

for(i in 1:size){ 
    for(j in 1:size){ 
    if(id[i] == id[j]){ 
     mat[i, j] <- 1 
    } 
    } 
} 

Answer 1

outerを使用できます。

mat <- outer(id, id, "==") 

出力はバイナリ行列であり、O（N^2）オブジェクトは一種の大であるので、これはスパース行列のための良好なユースケースである：

library(Matrix) 
mat <- Matrix(nrow=100, ncol=100) 
mat[] <- outer(id, id, "==") # [] means to assign into the existing 'mat' matrix 

Answer 2

別の代替、これべきあなたはまた、スパース行列

library(Matrix) 
tcrossprod(sparse.model.matrix(~ id + 0)) 

012を使用することができ香港大井の答えに

tcrossprod(model.matrix(~ id + 0)) 

同様に、比較的高速であります