计算稀疏矩阵的零空间

Question

我找到了函数（null OR nullspace）来查找R中常规矩阵的零空间，但我找不到稀疏矩阵（sparseMatrix）的任何函数或包。

有人知道怎么做吗？

Answer 1

如果您查看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]