如何在R中高效实现多变量配对的Fisher精确检验？

针对你需要批量处理39个名义变量两两Fisher精确检验的需求，用R可以通过以下步骤高效实现，避免手动重复编码：

1. 准备数据 #

假设你的问卷数据存储在数据框df中，所有39个名义变量都是因子类型（如果不是，先转成因子：df <- lapply(df, factor) %>% as.data.frame()）。

2. 生成所有不重复的变量配对 #

用combn()函数生成所有两两不重复的变量组合，避免重复计算（比如变量A&变量B与变量B&变量A是同一个列联表）：

# 获取所有变量名的两两组合
var_pairs <- combn(names(df), 2)

3. 定义检验函数 #

写一个自定义函数，输入一对变量名，生成列联表并执行Fisher精确检验，返回关键结果：

fisher_test_pair <- function(var1, var2) {
  # 生成列联表
  cont_table <- table(df[[var1]], df[[var2]])
  # 执行Fisher精确检验
  test_result <- fisher.test(cont_table)
  # 返回需要的结果，可按需调整字段
  return(data.frame(
    Var1 = var1,
    Var2 = var2,
    P_value = test_result$p.value,
    Odds_Ratio = test_result$estimate,
    Conf_Low = test_result$conf.int[1],
    Conf_High = test_result$conf.int[2]
  ))
}

4. 批量执行检验并整理结果 #

用apply()遍历所有变量对，应用上面的函数，最后把结果合并成一个数据框：

# 批量执行检验
all_results <- apply(var_pairs, 2, function(x) fisher_test_pair(x[1], x[2]))
# 合并成单个数据框
final_results <- do.call(rbind, all_results)

5. 查看结果 #

直接打印final_results，或者用View(final_results)在RStudio中打开查看，方便筛选显著的变量对。

内容的提问来源于stack exchange，提问作者Jorge Guerrero De León