在R中矩阵的秩

Question

我想测试矩阵的秩，是否有人能推荐一个在R中实现这个功能的包或函数？

- user1274212

Answer 1

你可以尝试使用函数qr（"qr"因为它执行了QR分解）：

#define a matrix for this example
M <- matrix(data = rnorm(12), ncol = 3)

#run the function qr() 
qr(M)$rank

#Alternative: load the Matrix package...
require(Matrix)

#...and run the function rankMatrix()
rankMatrix(M)[1]

- Qaswed

Answer 2

http://cran.r-project.org/web/packages/Matrix/Matrix.pdf，第101页

http://cran.r-project.org/web/packages/matrixcalc/matrixcalc.pdf，第12页

- Andrius Naruševičius

Answer 3

你可以使用pracma库：实用数值数学（提供大量的数值分析和线性代数、数值优化、微分方程、时间序列以及一些著名的特殊数学函数）。

在R控制台中使用以下命令进行安装： install.packages("pracma", repos="http://R-Forge.R-project.org")

然后你就可以使用该库了： library(pracma) Rank(你的矩阵对象)

- Abdul Razzak