（R）生成半正交随机矩阵

Question

我需要生成随机NxK矩阵（其中N> K），其中列K是正交随机向量。

我尝试的一个选项是生成一个大小为NxN的平方正交矩阵，然后选择前K列，但我想知道是否有更有效的方法。

当前代码（在R中）：

library(pracma)
Z <- rortho(N)[,1:K]

Answer 1

您可以使用rstiefel软件包来做到这一点：

library(rstiefel)
rustiefel(N, K)