使用bigalgebra構造對角矩陣

我想從矢量構造更大的對角矩陣。我安裝了bigalgebra包，但它沒有diag功能。另外，如何製作逆矩陣（solve）和轉置矩陣（t）。使用bigalgebra構造對角矩陣

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

如果稀疏矩陣是一個選項，您可以使用Matrix包（隨R提供）。

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

2017-07-13 01:34:19

謝謝你！但是，如何解決這個問題：'../Core/cholmod_memory.c，line 147'文件中的'Cholmod error'內存不足'。它發生在我認識產品'sqrt（V）％*％X％*％V％*％t（X）％*％sqrt（V）'時，其中'X'也是一個稀疏矩陣。 –

回答