如何在自定义函数中使用purrr::map与aov并传递参数至aov？

核心问题 #

你的函数存在三个关键问题：

• 未将metals_list作为参数传入，函数内部无法识别该变量，导致无法正确提取要分析的Y变量列
• ANOVA公式里直接写x.var * z.var会被当作字面量字符串，而非你传入的实际变量名，引发公式解析错误
• 函数没有返回计算结果，即使逻辑正确，调用后也无法获取输出

修正后的函数（两种写法） #

写法一：字符串拼接构造公式 #

anova_fun <- function(df, y_vars, x_var, z_var) {
  # 遍历每个Y变量，动态构造ANOVA公式
  map(df[y_vars], function(y_col) {
    tryCatch({
      # 将变量名拼接成公式字符串，再转换为公式对象
      formula <- as.formula(paste(deparse(substitute(y_col)), "~", x_var, "*", z_var))
      aov(formula, data = df)
    }, error = function(e) NULL)
  })
}

# 调用示例：x和z需传入字符串形式的列名
aov_res <- anova_fun(my_data, metals_list, "my_x", "my_z")

写法二：用reformulate构造公式（更安全） #

anova_fun <- function(df, y_vars, x_var, z_var) {
  # 遍历Y变量列表，逐个构造公式并执行ANOVA
  map(y_vars, function(y_var) {
    tryCatch({
      # 直接生成公式：响应变量为y_var，预测变量包含x、z及二者交互项
      form <- reformulate(c(x_var, z_var, paste(x_var, z_var, sep = ":")), response = y_var)
      aov(form, data = df)
    }, error = function(e) NULL)
  }) %>% set_names(y_vars)  # 给结果列表添加变量名，方便后续查看
}

# 调用示例
aov_res <- anova_fun(my_data, metals_list, "my_x", "my_z")

关键修改说明 #

1. 新增y_vars参数：将原metals_list作为参数传入，明确函数要分析的Y变量集合
2. 动态构造公式：根据传入的变量名生成合法的ANOVA公式，避免硬编码导致的解析错误
3. 保留变量名：第二种写法通过set_names给结果列表添加原Y变量名称，提升结果可读性
4. 自动返回结果：R函数默认返回最后一行计算结果，无需额外编写return()语句

关于警告的说明 #

你遇到的「restarting interrupted promise evaluation」警告，本质是函数内部变量解析失败导致延迟求值中断。修正变量传递和公式构造逻辑后，该警告会自动消失。

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