1
A
回答
2
如果你看看ggm::null
的代码,你会发现它是基于输入矩阵的QR分解。
另一方面,Matrix
程序包提供了自己的方法来计算稀疏矩阵的QR
分解。
例如:
require(Matrix)
A <- matrix(rep(0:1, 3), 3, 2)
As <- Matrix(A, sparse = TRUE)
qr.Q(qr(A), complete=TRUE)[, 2:3]
qr.Q(qr(As), complete=TRUE)[, 2:3]
相关问题
- 1. 加快稀疏矩阵计算
- 2. 如何在Matlab中计算大型稀疏矩阵的零点?
- 3. 稀疏矩阵
- 4. 稀疏矩阵和矩阵
- 5. 计算矩阵的零空间
- 6. Python稀疏矩阵非零值交点
- 7. 本征稀疏矩阵乘法似乎计算全矩阵
- 8. 以稀疏矩阵
- 9. CUDA上的稀疏矩阵运算
- 10. 将列放入空的稀疏矩阵
- 11. 特征3.3中的稀疏矩阵非零计数
- 12. 98%稀疏矩阵的矩阵完成
- 13. python稀疏矩阵的矩阵功率
- 14. 在Python中连接零矩阵和稀疏矩阵
- 15. 基于另一个矩阵(矩阵包)的稀疏矩阵的清零元素
- 16. R矩阵包:Demean稀疏矩阵
- 17. 稀疏矩阵乘法执行时间
- 18. 用稀疏矩阵乘二次形式矩阵的算法
- 19. 稀疏矩阵的特征向量的快速计算
- 20. 反相稀疏矩阵
- 21. R:在稀疏矩阵
- 22. 点产品稀疏矩阵
- 23. 切片稀疏(scipy)矩阵
- 24. 多稀疏矩阵与SQL
- 25. matlab:稀疏矩阵分解
- 26. scipy稀疏矩阵分裂
- 27. 稀疏矩阵内存
- 28. 朱莉娅稀疏矩阵
- 29. 转身稀疏矩阵
- 30. 稀疏矩阵对齐
嗨文森特。我用你的评论,这很好。虽然矩阵和它的转置是分解的,但在稀疏情况下,当行数小于列数时,它不分解。作为矩阵(A,sparse = TRUE)和qr.Q(qr(A),complete = TRUE),可以使用以下公式来计算矩阵A和矩阵B的相关系数矩阵:A <矩阵(1:6,nc = 2,nr = 3)和 。和qr.Q(qr(A),qr.Q(qr(As),complete = TRUE)和qr.Q(qr(t(A)),complete = TRUE) ),complete = TRUE)不起作用 .local(x,...)中的错误:必须有#{rows}> =#{columns}。非常感谢。最好的 – Mahin
我得到相同的错误信息。我不知道为什么它在密集矩阵上工作,但对于稀疏矩阵失败... –
非常感谢您的关注。 – Mahin