对于我的一个项目,我想创建几个具有完全排名的随机矩阵。有人知道在R中快速做到这一点还是知道如何继续?
答案 0 :(得分:5)
如果使用iid元素生成矩阵,并且没有其他约束,则极有可能获得满秩矩阵:
library(Matrix)
set.seed(101)
r <- replicate(1000,rankMatrix(matrix(rnorm(10000),100)))
table(r) ## all values are equal to 100
(花费更多时间在数学上的人可能能够证明这个矩阵空间内的降秩矩阵集实际上有0 ...)