如何用R生成n个0-1伯努利随机变量的所有排列矩阵？

嘿，这个需求在做全因子实验设计、布尔逻辑分析的时候特别常见！我给你分享几种R语言里的实现方法，从基础到进阶都有，你可以根据自己的习惯选：

方法一：用基础包的expand.grid（最直观） #

expand.grid是R基础包专门用来生成笛卡尔积的函数，简直是为这个需求量身定做的：

n <- 3  # 替换成你需要的变量个数
# 先创建一个包含n个元素的列表，每个元素是c(1, 0)
var_list <- lapply(1:n, function(x) c(1, 0))
# 生成所有组合的 data frame
full_combinations <- expand.grid(var_list, stringsAsFactors = FALSE)
# 给列起个直观的名字
colnames(full_combinations) <- paste0("var", 1:n)

print(full_combinations)

运行后会得到一个2^n行、n列的数据框，每一行就是一组0/1的排列。

方法二：用位运算生成二进制字符串（巧妙高效） #

如果你想先得到类似111、110这样的字符串形式，再拆分成变量列，位运算的方法很有意思：

n <- 3
# 生成从0到2^n - 1的所有整数
num_range <- 0:(2^n - 1)
# 把每个整数转成n位二进制字符串（补前导0）
binary_strs <- sprintf(
  paste0("%0", n, "d"),
  sapply(num_range, function(x) as.integer(intToBits(x)[1:n]))
)
# 把字符串拆分成列，转成数值型数据框
full_combinations <- as.data.frame(
  do.call(rbind, strsplit(binary_strs, "")),
  stringsAsFactors = FALSE
)
colnames(full_combinations) <- paste0("var", 1:n)
full_combinations[] <- lapply(full_combinations, as.integer)

print(full_combinations)

这个方法的核心是利用整数的二进制表示，直接生成所有可能的0/1组合，效率很高，尤其是当n比较大的时候。

方法三：tidyverse 用户专属的tidyr::crossing #

如果你平时用tidyverse生态，crossing函数比expand.grid更符合tidy风格，返回的是tibble格式：

library(tidyr)
n <- 3
# 创建命名的变量列表
var_list <- setNames(rep(list(c(1, 0)), n), paste0("var", 1:n))
# 生成所有组合
full_combinations <- crossing(!!!var_list)

print(full_combinations)

!!!是tidyeval语法，用来把列表里的元素作为参数传递给crossing，用起来很顺手。

额外调整：按你想要的顺序排列 #

上面的方法生成的顺序可能和你例子里的不完全一致，如果你需要从111到000的降序排列，可以给数据框加个排序逻辑：

# 按每行的"二进制数值"降序排列
full_combinations$binary_value <- apply(
  full_combinations, 1,
  function(row) sum(row * 2^(rev(seq_along(row)) - 1))
)
sorted_combinations <- full_combinations[order(-full_combinations$binary_value), ]
# 删掉辅助列
sorted_combinations$binary_value <- NULL

print(sorted_combinations)

这样就能得到你例子里的顺序啦！

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