使用metafor包循环调用rma.glmm函数报错求助

错误根源 #

你遇到的Error in ni - xi : non-numeric argument to binary operator，本质是循环中直接用列名字符串赋值给xi，R把它当成了普通字符向量，而非数据框里的数值列，导致后续运算时类型不匹配。

修正方案 #

用数据框的[[ ]]索引提取目标列的数值向量，替代直接赋值列名字符串。以下是完整的可运行示例：

假设你的数据框结构 #

假设数据框名为meta_df，包含：

• 多个待分析的比例/成功数列（比如prop_a、prop_b）
• 总样本量列total_n

修正后的循环代码 #

library(metafor)

# 定义要分析的目标列名
target_cols <- c("prop_a", "prop_b", "prop_c")

# 创建空列表存储每个列的分析结果
analysis_results <- list()

# 循环遍历目标列
for (col_name in target_cols) {
  # 提取当前列的数值：如果是比例，先转成整数成功数；如果已是成功数，直接取列
  xi <- round(meta_df[[col_name]] * meta_df$total_n)  # 比例转成功数
  # xi <- meta_df[[col_name]]  # 若列本身是成功数，用这行替代上面
  
  ni <- meta_df$total_n
  
  # 调用rma.glmm做比例元分析
  res <- rma.glmm(measure = "PLO", xi = xi, ni = ni, data = meta_df)
  
  # 将结果存入列表，用列名作为键
  analysis_results[[col_name]] <- res
  
  # 可选：打印当前列的分析摘要
  cat("--- 列", col_name, "分析结果 ---\n")
  print(summary(res))
}

关键注意事项 #

• 确认target_cols中的列名在数据框中存在，且列类型为数值型（用str(meta_df)检查，避免字符/因子类型的列）
• 如果目标列是比例值，必须转换为整数成功数（xi要求是事件发生的整数次数）
• 若总样本量列名不是total_n，记得替换成你实际的列名

内容的提问来源于stack exchange，提问作者Mariano C Giglio