我有一个很大的相关矩阵,62589x62589。我已经将矩阵二值化到某个阈值之上,我没有遇到任何问题,但我对基本计算时间的显着差异感到有些困惑。
我第一次这样做.... 1 的数量:425,491 ... 0 的数量:3,916,957,430 这两个数字的总和 == 62589^2,这意味着矩阵是真正二值化的。我将其保存为 Rdata 对象 (31Mb)。执行矩阵的基本计算大约需要 3.5 分钟。
fooB <- foo %*% foo
第二次,阈值较低...... 1 的数量:30,384,683 ... 0 的数量:3,886,998,238。这些总和又是 62589^2,因此是真正的二值化。Rdata 对象为 84Mb。执行与上述相同的乘法步骤目前仍在计算一个小时后。
最新矩阵中增加的 1 数量是否应该如此大幅度地增加文件大小和处理时间?
谢谢阅读
编辑:对第二个矩阵进行相同计算的最后时间是 65 分钟
Edit2:执行 is() 结果:矩阵数组结构向量