R语言for循环中向tabyl管道传入列名批量执行卡方检验的方法

报错的核心是tabyl()函数不支持直接传入字符串格式的列名，会将字符串识别为不存在的列名，需要先将字符串转为R语言的符号对象才能被正常识别。

首先确保已加载依赖包：

library(dplyr)
library(janitor)

修正后的循环代码如下，可直接适配你60+症状列的数据集：

# 初始化存储p值的列表
p_result <- list()

# 遍历所有症状列（第2列到最后一列，自动适配列数）
for (i in 2:ncol(df)) {
  # 将字符串格式的症状列名转为符号对象
  symptom_col <- sym(colnames(df)[i])
  # 生成列联表、执行卡方检验、提取p值
  chi_p <- chisq.test(df %>% tabyl(Race, !!symptom_col, show_na = FALSE))$p.value
  # 存入列表，用列名作为键
  p_result[[colnames(df)[i]]] <- chi_p
}

# 转为要求格式的结果数据框
final_result <- data.frame(p_result, row.names = "p.value")

# 查看结果
print(final_result)

• 用sym()函数将字符串列名转为符号，搭配!!运算符注入到tabyl()的参数中，解决列名识别错误问题
• 仅提取p值存入列表，避免存入完整卡方检验对象导致后续转数据框出错
• 遍历范围设为2:ncol(df)，不需要手动调整列数范围，任意数量的症状列都可直接运行

如果习惯tidyverse语法，也可以用purrr::map实现更简洁的批量计算，效果和for循环完全一致：

library(purrr)
final_result <- df %>%
  select(-Race) %>%
  map_dbl(~ chisq.test(tabyl(df, Race, .x, show_na = FALSE))$p.value) %>%
  t() %>%
  as.data.frame() %>%
  `rownames<-`("p.value")

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