R中的巨大对角矩阵

Question

以下代码导致内存错误：

diag(1:100000)

有没有替代diag允许产生巨大的对角矩阵？

Answer 1

更长的答案：我建议不要创建对角矩阵，因为在大多数情况下你可以不用它。为了清楚起见，请考虑最典型的矩阵运算：

1. 将对角矩阵D乘以向量v以产生Dv。不要保持矩阵，而是将“矩阵”保持为对角线元素的向量d，然后将d元素乘以v。结果相同。
2. 反转矩阵。同样，简单：反转每个元素（当然，仅对于对角矩阵，这通常是正确的反转）。
3. 各种分解/特征值/决定因素/跟踪。同样，这些都可以在矢量d上完成。

简而言之，虽然它需要在代码中引起一些注意，但您始终可以将对角矩阵表示为向量，这应该可以解决您的内存问题。

更短的回答：现在，说了这么多，当然人们已经使用稀疏矩阵隐式地实现了上述步骤，这些步骤完成了上述步骤。在R中，Matrix包非常适合稀疏矩阵：https://cran.r-project.org/web/packages/Matrix/Matrix.pdf