如何使用for循环对多变量执行配对Wilcoxon检验

云原生阿岚

2026-9-24

问题核心原因

报错本质是wilcox.test()的公式接口在循环调用时，字符串形式的变量名无法被公式环境直接识别，以及循环索引取值时未正确匹配数据框列、配对分组的样本长度对应逻辑出错。

修正方案

前置说明

默认你的数据符合配对检验的结构要求：要么每一行对应1组配对样本，1列标识分组（仅2个水平，每组样本量完全相等），其余列为待检验的指标；要么同一配对的两组数值分别存在两个列中，行位置一一对应。

可用代码实现

方案1：公式接口适配（适合批量多变量场景）

# 1. 先定义待检验的变量名列表，替换为你实际的变量名
test_vars <- c("var1", "var2", "var3")
result_list <- list()

for (i in seq_along(test_vars)) {
  # 用as.formula将拼接的字符串转为公式对象，指定数据环境避免识别错误
  f <- as.formula(paste0(test_vars[i], " ~ group")) # group替换为你实际的分组列名
  res <- wilcox.test(f, data = df, paired = TRUE) # df替换为你的数据框名
  result_list[[test_vars[i]]] <- res
}

方案2：直接传向量（规避公式识别问题，报错率更低）

如果分组为两个独立的列，或担心公式识别出错，直接提取两组向量传入：

test_vars <- c("var1", "var2", "var3")
result_list <- list()

for (v in test_vars) {
  # 提取两组对应数值，确保一一对应、无缺失
  g1_val <- df[df$group == "level1", v] # level1替换为你第一组的实际分组值
  g2_val <- df[df$group == "level2", v] # level2替换为你第二组的实际分组值
  res <- wilcox.test(x = g1_val, y = g2_val, paired = TRUE)
  result_list[[v]] <- res
}