我如何识别不增加R中矩阵秩的行/列

Question

假设我有一个维数n x n的对称半正定矩阵A。A的秩为d ，我想找到一组长度为d的指标，使A[indices, indices]的排名为d。例如，if

A <- structure(c(3, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), .Dim = c(4L, 4L))

那么我想确定索引c(1, 2, 3)或c(1, 2, 4)中的一个集合。

我可以通过蛮力遍历所有行和列的组合来做到这一点，但我的猜测是存在更优雅和有效的方法来做到这一点。

Answer 1

您可以使用QR分解，qr

q <- qr(A)
q$pivot[seq(q$rank)]
# [1] 1 2 3
A[,q$pivot[seq(q$rank)]]
#      [,1] [,2] [,3]
# [1,]    3    2    1
# [2,]    2    2    1
# [3,]    1    1    1
# [4,]    1    1    1

这些结果列应该是线性无关的