r - R中expand.grid函数的修改

Question

我需要创建一个包含所有可能组合的矩阵。

L<-12 
vec <- c(0:21)
lst <- lapply(numeric(L), function(x) vec)
Mat1<-as.matrix(expand.grid(lst))

结果将是非常大的矩阵，我的计算机无法计算它。
实际上，我只需要第一列中的值大于第二列中的值，第二列中的值大于第三列中的值的组合，依此类推。我可以修改 expand.grid 函数以删除不必要的组合吗？

Answer 1

正如@AllanCameron 所确定的，您所追求的等同于采用 vec 的所有组合大小为 12 -- 您也可以使用内置函数 combn为了那个原因:

do.call(rbind, combn(vec, L, simplify = FALSE))

使用 data.table 可能会更快操作 combn 的输出转换成所需的格式:

library(data.table)
setDT(transpose(combn(vec, L, simplify=FALSE))