7

遺伝子発現データ(40000遺伝子(変数)×30観測)から40000×40000の共分散行列を作成したい。これは間違いなく私のRAMよりも大きいです。パッケージ「ff」を使用して、相関のために40000x40000の空のマトリックスを事前に割り当てることができました。ただし、「cov」または「cor」関数は、私のシステムでは 5000x5000 の共分散行列のみを管理するため、ブロックごとに 1:5000、5001:10000 などの共分散計算を行い、事前に割り当てられた行列を対角線に沿って埋める必要があります。行列の「欠落しているパッチ」、つまり共分散 (または相関) 1 と 22000 を埋めるアルゴリズムを知っている人はいますか。 ' は非常に高速です... では、既に計算された共分散を使用して 1/22000 の cov (または cor) を計算する方法はありますか?

前もって感謝します!

4

1 に答える 1

1

cov2 つの引数を使用して、非対角ブロックを計算できます。

cov( x[,1:5000], x[,5001:10000] )
于 2013-02-16T15:39:41.940 に答える