【发布时间】:2012-06-08 12:45:07
【问题描述】:
我想测试矩阵的等级,有没有人可以为此推荐 R 中的包/函数?
【问题讨论】:
我想测试矩阵的等级,有没有人可以为此推荐 R 中的包/函数?
【问题讨论】:
您可以使用图书馆 pracma:实用数值数学(提供大量来自数值分析和线性代数、数值优化、微分方程、时间序列的函数,以及一些众所周知的特殊数学函数。)。
在 R 控制台中使用以下命令安装它: install.packages("pracma", repos="http://R-Forge.R-project.org") 然后您可以使用该库: 图书馆(实习) 排名(你的矩阵对象)
【讨论】:
你可以试试qr函数(“qr”,因为它执行QR decomposition):
#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]
【讨论】: