求大神啊,在R中可以根据qr分解来求得某一矩阵的秩和pivot,如矩阵:
[,1] [,2] [,3] [,4] [,5]
[1,] 1 1 0 0.000000 50
[2,] 1 1 0 0.000000 51
[3,] 1 1 0 0.000000 52
[4,] 1 1 0 0.000000 54
[5,] 1 1 0 0.000000 53
[6,] 1 0 2 0.000000 60
[7,] 1 0 2 0.000000 59
[8,] 1 0 2 0.000000 65
[9,] 1 0 2 0.000000 67
[10,] 1 0 2 0.000000 70
[11,] 1 0 0 3.000000 70
[12,] 1 0 0 3.000000 73
[13,] 1 0 0 3.000000 74
[14,] 1 0 0 3.000000 78
[15,] 1 0 0 3.000000 82
[16,] 1 0 0 3.000000 80
[17,] 1 0 0 3.000000 87
[18,] 1 0 0 3.000000 84
[19,] 1 0 0 3.000000 88
[20,] 1 0 0 3.000001 92
qr(X)$rank=4,
qr(X)$pivot: 1 2 3 5 4
我想在C++中求得某一矩阵值包含X矩阵的1 2 3 5列,不包含第4列。在C++中可以利用Eigen中的ColPivHouseholderQR求矩阵的秩,如何求没有多重共线性列的矩阵,求高人指点啊.