R函数引号使用问题：多Y变量统计分析函数报错排查

修正后的完整函数 #

library(rstatix)
library(dplyr)
library(tidyr)
library(rlang)

summarize_and_wilcoxon <- function(data, xvar, yvars) {
  # 处理分组变量xvar：同时支持裸变量名和字符串输入
  x_sym <- ensym(xvar)
  x_str <- as_name(x_sym)
  
  # 生成汇总统计表（中位数+四分位距）
  summary_tbl <- data %>%
    group_by(!!x_sym) %>%
    select(all_of(yvars), !!x_sym) %>%
    get_summary_stats(type = "median_iqr") %>%
    pivot_wider(
      names_from = !!x_sym,
      values_from = c(median, iqr),
      names_glue = "{x_str}_{.value}_{col}"
    ) %>%
    rename(variable = variable)
  
  # 批量执行Wilcoxon检验
  wilcox_tbl <- data %>%
    pivot_longer(cols = all_of(yvars), names_to = "variable", values_to = "value") %>%
    group_by(variable) %>%
    wilcox_test(value ~ !!x_sym) %>%
    add_significance() %>%
    mutate(
      comparison = paste(x_str, "分组比较"),
      p_label = case_when(
        p < 0.001 ~ "p < 0.001",
        TRUE ~ paste0("p = ", round(p, 3))
      )
    ) %>%
    select(variable, comparison, p, p_label, significance)
  
  # 合并统计结果与检验结果
  full_join(summary_tbl, wilcox_tbl, by = "variable")
}

核心问题解决说明 #

1. paste中xvar的正确引用
   之前直接在paste中使用xvar时，若传入裸变量名（如has_weight），它是一个符号对象而非字符串，无法被paste正确解析。用as_name(ensym(xvar))将符号转换为变量名字符串，即可在paste或names_glue中正常使用。

2. xvar加引号/不加引号的兼容问题
   使用ensym()替代enquo()，该函数可以自动识别输入是裸变量名还是字符串，统一转换为符号对象，再通过!!（解引用操作符）注入到dplyr的语法中，避免了分组失效或对象找不到的报错。

3. group_by失效的修复
   当传入带引号的xvar时，若直接用!!xvar会被当作字面量字符串处理，而ensym()会将字符串转换为对应的数据框列名符号，配合!!就能让group_by正确识别分组变量。

R函数参数引号的最佳实践 #

• tidyverse风格优先支持裸变量名：符合dplyr等包的语法习惯，提升代码可读性，比如调用时用xvar = has_weight而非xvar = "has_weight"。
• 用ensym/ensyms兼容字符串输入：ensym()处理单个变量，ensyms()处理多个变量，让函数同时支持两种输入方式，适配不同使用场景。
• 多变量参数用all_of()处理：对于yvars这类多变量参数，用all_of(yvars)可以同时接收字符串向量（如c("weight", "height")）和裸变量名向量（如c(weight, height)）。
• 变量名转字符串用as_name()：避免手动拼接字符串，确保变量名的准确性，尤其是动态生成输出列名时。

调用示例 #

# 构造示例数据
set.seed(123)
demo_df <- tibble(
  has_weight = sample(c("是", "否"), 100, replace = TRUE),
  weight = rnorm(100, 70, 10),
  height = rnorm(100, 170, 5),
  bmi = rnorm(100, 24, 3)
)

# 方式1：裸变量名调用
summarize_and_wilcoxon(demo_df, has_weight, c(weight, height, bmi))

# 方式2：字符串调用
summarize_and_wilcoxon(demo_df, "has_weight", c("weight", "height", "bmi"))

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