在R中对稀疏向量进行对角化的有效方法

时间:2017-07-30 14:40:21

标签: r sparse-matrix diagonal

我正在处理各种值的矢量(~​​14000x1),我想把它放在稀疏矩阵的对角线上,我在那里使用矩阵库。我想这样做,同时避免创建一个完整的矩阵,然后转换回稀疏矩阵。

到目前为止,我可以使用for循环执行此操作,但这需要很长时间。你能想到一种更有效,记忆力更强的方式吗?

这是一个简单的可重复示例:

library(Matrix)
x = Matrix(matrix(1,14000,1),sparse=TRUE)

X = Diagonal(14000)

for(i in 1:13383){
  X[i,i]=aa[i]
  print(i)

}

0 个答案:

没有答案