R中编写自定义函数执行多变量配对t检验的报错排查

错误原因 #

你写的函数里直接用{{parameter}} ~ {{grouping_variable}}构造公式，rstatix的t_test无法正确解析这种方式，导致内部pull()找不到对应列，抛出错误。

正确的自定义函数写法 #

方法1：用rlang构造公式 #

利用rlang::ensym()把参数转成符号，再用new_formula()生成合法的公式对象：

library(rstatix)
library(rlang)

my_t_test <- function(dataset, parameter, grouping_variable) {
  # 将传入的裸变量名转为符号
  param_sym <- ensym(parameter)
  group_sym <- ensym(grouping_variable)
  # 构造公式
  test_formula <- new_formula(param_sym, group_sym)
  
  dataset %>% 
    t_test(formula = test_formula, paired = TRUE, detailed = TRUE) %>% 
    add_significance()
}

# 测试（设种子保证结果可复现）
set.seed(123)
df <- data.frame(grouping = c(rep("left", 50), rep("right", 50)), 
                 var1 = rnorm(100, mean=21, sd=3))

my_t_test(df, var1, grouping)

方法2：字符串拼接构造公式 #

把参数转成字符串后拼接成公式：

library(rstatix)
library(rlang)

my_t_test <- function(dataset, parameter, grouping_variable) {
  param_str <- as.character(ensym(parameter))
  group_str <- as.character(ensym(grouping_variable))
  # 字符串转公式
  test_formula <- as.formula(paste0(param_str, " ~ ", group_str))
  
  dataset %>% 
    t_test(formula = test_formula, paired = TRUE, detailed = TRUE) %>% 
    add_significance()
}

dplyr风格函数的符号与作用域规则 #

• {{ }}（大括号插值）：仅在dplyr核心动词（如mutate、select）中完美生效，它会自动处理非标准评估，把函数参数映射到数据框的列。但不是所有tidyverse工具都支持这种直接拼接公式的方式，比如rstatix的t_test对公式的解析逻辑不同。
• ensym()/enquo()：用于把函数接收的裸变量名（比如调用时写的var1）转为符号/引用对象，方便手动构造公式或适配其他非标准评估函数。
• 数据掩码作用域：tidyverse函数默认使用数据掩码，优先从传入的数据框环境中查找变量，而非全局环境。自定义函数时，需要明确告诉工具如何在数据掩码环境中解析参数——要么用{{ }}适配支持的函数，要么手动构造符号/公式来匹配目标函数的参数要求。

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