使用bigalgebra构造对角矩阵

Question

我想从矢量构造一个更大的对角矩阵。我安装了bigalgebra包，但它没有diag功能。此外，如何将逆（solve）和转置（t）转换为大矩阵。

v <- runif(42109)
V <- diag(v)
Error: cannot allocate vector of size 13.2 Gb

Answer 1

如果稀疏矩阵是一个选项，您可以使用Matrix包（随R提供）。

library(Matrix)
V <- Matrix(0, nrow=42109, ncol=42109)
diag(V) <- v