我在 R 中找到了用于查找正则矩阵的零空间的函数 (null OR nullspace),但我找不到用于稀疏矩阵 (sparseMatrix) 的任何函数或包。
有人知道怎么做吗?
最佳答案
如果你看一下 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]
关于r - 计算稀疏矩阵的零空间,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/46036316/