7

从基因表达数据(40000 个基因(变量)x 30 个观察值)我想创建一个 40000 x 40000 协方差矩阵。这绝对比我的 RAM 大。使用包“ff”,我设法为相关性预分配了一个 40000x40000 的空矩阵。然而,'cov' 或 'cor' 函数在我的系统上只能管理一个 5000x5000 的协方差矩阵,所以我必须按块进行 1:5000、5001:10000 等协方差计算并沿对角线填充预分配的矩阵。有谁知道填充矩阵中“缺失补丁”的算法,即 1 和 22000 之间的协方差(或相关性)。我知道我可以做所有成对组合并一个接一个地填充矩阵,但是 'cor ' 非常快......那么,有没有办法通过使用已经计算的协方差来计算 1/22000 的 cov(或 cor)?

提前致谢!

4

1 回答 1

1

您可以使用cov2 个参数来计算非对角块。

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