在R中按研究分组循环执行Meta分析时遇报错的技术问询

示例数据 #

创建包含研究、阶段、效应量(b)和标准误(se)的数据集：

dat <- data.frame(
  Study = rep(c("Study 1", "Study 2", "Study 3", "Study 4", "Study 5"), each = 2),
  Stage = rep(c("Discovery", "Replication"), 5),
  b = c(0.1, 0.2, 0.15, 0.3, 0.05, 0.1, 0.25, 0.4, 0.08, 0.15),
  se = c(0.05, 0.06, 0.07, 0.08, 0.09, 0.06, 0.07, 0.08, 0.1, 0.12)
)

打印数据集：

print(dat)

需求 #

对所有研究的Discovery阶段和Replication阶段分别执行Meta分析。

尝试代码 #

# 创建空列表存储Meta分析结果
meta_results <- list()

# 按阶段执行Meta分析
for (stage in c("Discovery", "Replication")) {
  # 筛选当前阶段的数据
  stage_data <- dat[dat$Stage == stage, ]
  
  # 使用meta包的函数执行Meta分析
  meta_result <- meta::metagen(
    estimate = stage_data$b,
    se = stage_data$se,
    data = stage_data
  )
  
  # 将结果存入列表
  meta_results[[stage]] <- meta_result
}

# 打印Meta分析结果
print(meta_results)

报错信息 #

Error in meta::metagen(estimate = stage_data$b, se = stage_data$se, data = stage_data) : 
  argument 2 matches multiple formal arguments

问题原因与解决方法 #

原因 #

meta::metagen函数的参数传递逻辑冲突：

• 同时指定data参数（传入数据集）和直接传递向量（如stage_data$b）时，函数无法正确匹配参数；
• 此外，metagen的标准参数为TE（效应量）和seTE（标准误），使用别名estimate/se搭配data参数易引发匹配歧义。

解决方法 #

两种修正方式任选其一即可：

方式1：移除data参数，直接传递向量

# 创建空列表存储Meta分析结果
meta_results <- list()

# 按阶段执行Meta分析
for (stage in c("Discovery", "Replication")) {
  # 筛选当前阶段的数据
  stage_data <- dat[dat$Stage == stage, ]
  
  # 直接传递向量，不指定data参数
  meta_result <- meta::metagen(
    estimate = stage_data$b,
    se = stage_data$se
  )
  
  meta_results[[stage]] <- meta_result
}

# 查看结果
print(meta_results)

方式2：使用data参数，传入列名（推荐更规范的TE/seTE参数）

# 创建空列表存储Meta分析结果
meta_results <- list()

# 按阶段执行Meta分析
for (stage in c("Discovery", "Replication")) {
  # 筛选当前阶段的数据
  stage_data <- dat[dat$Stage == stage, ]
  
  # 指定data参数，传入列名（用TE和seTE更符合函数规范）
  meta_result <- meta::metagen(
    TE = b,
    seTE = se,
    data = stage_data,
    studlab = Study  # 可选：添加研究标签，便于结果解读
  )
  
  meta_results[[stage]] <- meta_result
}

# 查看结果
print(meta_results)

说明 #

两种方法都能正确完成分阶段Meta分析：

• 筛选stage_data时，每个阶段会得到5个研究的对应效应量和标准误，符合Meta分析的输入要求；
• 若需要在结果中显示研究名称，可添加studlab = Study参数（方式2中已包含）。

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