Rのゼロと1の大きな行列を素早く効率的に計算する方法

Question

私はしばらくRを使用していましたが、パフォーマンスのために自分のコードを最適化するのは得策ではありません。そのi番目の行は、すべてj = 1, 2, ..., tため$ 1_ {Y_I < Y_j} $に対応するように長さtのベクトルY、与えられ

、Iは、順tの正方行列を生成しようとしています。これは0/1の2進行列です。

しかし、私はこれを実行する最も効率的な方法ではないと主に考えています。主にオブジェクトのサイズを見ています。これは、密な行列と同じサイズです（バイト単位）。私の目標はこれを大きなtに拡張することですが、計算速度とそのようなマトリックスのストレージスペースは長すぎると考えています。この行列を保存して計算するためのより迅速かつ効率的な方法がありますか？

ありがとうございます！

Answer 1

私が正しくあなたを理解している場合は、我々は使用することができます。

outer(t, t, "<=") + 0L 
## operator `"<="` returns TRUE/FALSE logical, 
## add an integer 0 to make a 0/1 binary integer matrix 

はのは、あなたの参照forループと比較して、簡単なテストを持ってみましょう。

set.seed(0); t <- runif(5) ## reproducible data 

outer(t, t, "<=") + 0L 
#  [,1] [,2] [,3] [,4] [,5] 
#[1,] 1 0 0 0 1 
#[2,] 1 1 1 1 1 
#[3,] 1 0 1 1 1 
#[4,] 1 0 0 1 1 
#[5,] 0 0 0 0 1 

Y <- matrix(0, nrow=length(t), ncol=length(t)) 
for(i in 1:length(t)) { 
    Y[i, ] <- t[i] <= t 
} 

Y 
#  [,1] [,2] [,3] [,4] [,5] 
#[1,] 1 0 0 0 1 
#[2,] 1 1 1 1 1 
#[3,] 1 0 1 1 1 
#[4,] 1 0 0 1 1 
#[5,] 0 0 0 0 1 

同じ結果が得られます。ビンゴ！！

---

注：

人々はこの回答の編集の歴史を読みました。 sparseMatrixの方法が示されている。私は当初、疎な行列を使うと考えていました。そして、カラムピボットは効率を改善するでしょうが、そうではありません。ベンチマークでは、outerは6倍以上高速です。したがって、私はかなり悲しいと感じていますが、私はそれらの大きなチャンクを私の答えの中に入れなければなりません。