r - (R) 生成半正交随机矩阵

Question

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

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

当前代码（在 R 中）：

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

Answer 1

你可以用这个rstiefel包做到这一点：

library(rstiefel)
rustiefel(N, K)