使用metafor::escalc()封装函数优化Meta分析工作流遇数值错误

问题根源 #

原函数存在两个核心问题：

1. dplyr::select()中直接使用varname1:varname4无法识别字符串类型的变量名，导致后续数据框列选择错误；
2. escalc()参数中直接传入字符串变量名，函数无法将其解析为数据框中的数值列，从而触发"非数值"错误。

修复后的完整代码 #

library(tidyverse)
library(metafor)

# 可复现数据集
dat <- tibble(
  study = c("study1", "study2", "study3"),
  int1 = c("GA", "GA", "GA"),
  int2 = c("Spinal", "Epidural", "Spinal"),
  n_int1 = c(20, 25, 30),
  n_int2 = c(20, 25, 30),
  dynamic_pain_mean_24h_int1 = c(4.5, 4.7, 4.2),
  dynamic_pain_sd_24h_int1 = c(0.5, 0.7, 1.0),
  dynamic_pain_mean_24h_int2 = c(2.2, 2.1, 1.5),
  dynamic_pain_sd_24h_int2 = c(0.4, 0.6, 0.7),
  ome_mean_24h_int1 = c(54, 65, 70),
  ome_sd_24h_int1 = c(15, 20, 18),
  ome_mean_24h_int2 = c(40, 35, 50),
  ome_sd_24h_int2 = c(12, 17, 15),
)

# 修复后的变量选择与分析函数
var_select <- function(df, x, y, z){
  # 统一变量名拼接逻辑，避免重复判断
  build_var <- function(stat, int){
    parts <- c(if(!missing(x)) x, y, stat, z, int)
    paste(parts, sep = "_")
  }
  
  varname1 <- build_var("mean", "int1")
  varname2 <- build_var("sd", "int1")
  varname3 <- build_var("mean", "int2")
  varname4 <- build_var("sd", "int2")
  
  # 使用all_of()引用字符串列名，正确选择变量
  dat_processed <- df %>% 
    dplyr::select(study, int1, int2, n_int1, n_int2, all_of(c(varname1, varname2, varname3, varname4))) %>% 
    drop_na()
  
  # 使用!!sym()将字符串转换为列名符号，让escalc正确识别数值列
  dat_escaled <- metafor::escalc(measure = "MD", 
                                 n1i = n_int2, n2i = n_int1, 
                                 m1i = !!sym(varname1), sd1i = !!sym(varname2), 
                                 m2i = !!sym(varname3), sd2i = !!sym(varname4), 
                                 data = dat_processed)
  return(dat_escaled)
}

关键修复点说明 #

• 简化变量名拼接：通过内部函数build_var()统一处理有无x参数的情况，减少重复代码，提升可读性；
• 正确选择列名：使用dplyr::all_of()包裹字符串向量，让select()正确识别动态生成的变量名，替代原函数中无效的varname1:varname4写法；
• 解析列名给escalc：通过!!sym()将字符串变量名转换为R语言的符号对象，让escalc()直接调用数据框中的数值列，解决"非数值"参数错误。

测试调用示例 #

# 带x参数的场景（如dynamic_pain变量）
result1 <- var_select(dat, x = "dynamic", y = "pain", z = "24h")
print(result1)

# 不带x参数的场景（如ome变量）
result2 <- var_select(dat, y = "ome", z = "24h")
print(result2)

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