R找到所有可能的唯一组合

Question

我试图在 R 中找到所有可能的独特组合。似乎有很多类似的问题被问到，但我找不到相同的问题。

我的问题是从向量 x 中找到 m 个元素的组合，但 m 可能大于 x。例如，从字母 [1:2] 中选取 3 个元素，希望可以返回：

combn(letters[1:2],3)
       [,1]  [,2]  [,3]  [,4] 
[1,]   "a"   "a"   "a"    "b"
[2,]   "a"   "a"   "b"    "b"
[3,]   "a"   "b"   "b"    "b"

但是组合函数 n < m 中的错误。有类似的功能，包括 gtools:permutations、expand.grid。

如果之前有人问过同样的问题，但我没有听懂，再次道歉。谢谢。

Answer 1

有一些专门为此构建的软件包。基本前提是我们需要长度重复的组合，m其中m可能大于输入向量。我们从经典开始gtools：

library(gtools)
combinations(2, 3, letters[1:2], repeats.allowed = TRUE)
    [,1] [,2] [,3]
[1,] "a"  "a"  "a" 
[2,] "a"  "a"  "b" 
[3,] "a"  "b"  "b" 
[4,] "b"  "b"  "b"

然后是arrangements替代品iterpc（@Gregor 在上面的评论中链接的包）：

library(arrangements)
arrangements::combinations(2, 3, letters[1:2], replace = TRUE)
     [,1] [,2] [,3]
[1,] "a"  "a"  "a" 
[2,] "a"  "a"  "b" 
[3,] "a"  "b"  "b" 
[4,] "b"  "b"  "b"

最后是RcppAlgos，我写的：

library(RcppAlgos)
comboGeneral(letters[1:2], 3, TRUE)
     [,1] [,2] [,3]
[1,] "a"  "a"  "a" 
[2,] "a"  "a"  "b" 
[3,] "a"  "b"  "b" 
[4,] "b"  "b"  "b"

combn是一个很棒的函数，它作为带有 的基本包之一提供R，但是它的一个缺点是它不允许重复（这是这里要求的）。我为与此完全相同的问题编写了一份非常全面的概述，可在此处找到：A Walk Through a Slice of Combinatorics in R。