我正在处理各种值的矢量(~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)
}