gtools的combinations/permutations处理含重复元素向量失效原因及解决方案

gtools包的combinations()和permutations()函数原生逻辑没有考虑源向量内元素重复的频数限制：

• 当未开启repeats.allowed时，函数默认要求源向量唯一值数量≥抽取数量，你的源向量只有H、M、A3个唯一值，抽取数量为5，直接触发“too few different elements”报错
• 就算开启repeats.allowed，函数逻辑是允许元素无限制重复抽取，不会匹配你原始向量中H仅4个、M仅2个、A仅2个的实际库存，哪怕绕过报错生成的结果也不符合你的需求。

方法1：base R 最简实现 #

不需要手动做循环替换，只需要基于索引生成组合后映射回原向量即可：

# 原始向量
vec <- c("H", "H", "H", "H", "M", "M", "A", "A")
# 生成所有5个元素的索引组合
idx_comb <- combn(seq_along(vec), 5)
# 索引映射为原向量值，转置为和示例一致的矩阵格式
result <- t(apply(idx_comb, 2, function(x) vec[x]))
# 查看前6行
head(result)

输出结果和你给出的示例完全一致。

方法2：第三方包更高效实现（支持频数控制、排列/组合切换） #

用arrangements包可以直接指定每个元素的可用数量，不需要处理索引逻辑，如果需要去重的组合/排列会更方便：

# 安装包（首次使用执行）
# install.packages("arrangements")
library(arrangements)

# 统计每个元素的可用频数
freq <- table(vec)
# 生成符合频数限制的不重复组合（相同元素的不同位置不重复输出）
comb_res <- combinations(
  x = names(freq),
  k = 5,
  freq = as.integer(freq)
)

# 如果需要生成排列，替换为permutations函数即可
perm_res <- permutations(
  x = names(freq),
  k = 5,
  freq = as.integer(freq)
)

内容的提问来源于stack exchange，提问作者mmw