R语言for循环遍历数据框列运行prop_test无法识别变量i如何解决

报错原因 #

报错的核心原因是循环中的i是字符串格式的列名，直接写入公式i~sex时，infer包的函数会直接把符号i当做要查找的变量名，而非取i对应的字符串值作为响应变量名。

可行实现方案 #

方案1：for循环+动态公式构建 #

通过reformulate()函数动态生成每轮检验对应的公式，即可正确识别列名：

cols <- gss %>% select(-sex) %>% names()

for (i in cols){
  # 构建公式：响应变量为i对应的列，自变量为sex
  test_formula <- reformulate(termlabels = "sex", response = i)
  test_res <- infer::prop_test(
    gss,
    test_formula,
    order = c("female", "male")
  )
  # 按需打印或存储结果
  print(paste("当前检验变量：", i))
  print(test_res)
}

方案2：purrr迭代批量输出结果（更推荐） #

用purrr包的迭代函数替代显式for循环，可直接把所有检验结果整理成结构化数据框，方便后续分析：

library(purrr)
library(dplyr)

cols <- gss %>% select(-sex) %>% names()

# 批量运行检验，结果合并为数据框
all_test_res <- map_dfr(cols, function(col_name) {
  test_formula <- reformulate(termlabels = "sex", response = col_name)
  infer::prop_test(
    gss,
    test_formula,
    order = c("female", "male")
  ) %>%
    mutate(test_variable = col_name, .before = 1) # 新增列标注检验对应的变量
})

方案3：符号注入写法 #

也可以通过rlang的符号注入语法把字符串列名转为公式可识别的符号：

library(rlang)

cols <- gss %>% select(-sex) %>% names()

for (i in cols){
  col_sym <- sym(i)
  test_res <- infer::prop_test(
    gss,
    !!col_sym ~ sex,
    order = c("female", "male")
  )
  print(test_res)
}

注意事项 #

• 上述所有方案都保留了单变量检验时的order参数，确保检验逻辑和你之前的单变量运行逻辑完全一致
• 如果不需要打印结果，可自行定义列表存储每轮循环的检验结果即可

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