如何在R中通过循环批量绘制20个同格式的OR森林图？

第一步：批量获取目标文件路径 #

先一次性获取所有待处理文件的完整路径，假设你的20个文件都存在/data/logistic_reg/目录下，文件名都是log_xxx.txt格式（比如log_diab.txt），用list.files就能快速抓取：

# 指定数据文件所在目录
data_dir <- "/data/logistic_reg/"
# 匹配所有符合命名规则的txt文件，full.names=TRUE返回完整路径
file_list <- list.files(data_dir, pattern = "^log_.*\\.txt$", full.names = TRUE)

如果你的文件名不是这个规律，修改pattern参数即可，比如所有文件都是*_reg.txt，就改成pattern = ".*_reg\\.txt$"。

第二步：封装绘图函数 #

把你原有的绘图逻辑打包成函数，方便循环调用，还能自动生成对应文件名的图片：

library(ggplot2)

# 定义绘图函数：输入数据和图片名称前缀，输出并保存森林图
plot_forest <- function(data, plot_name) {
  p <- ggplot(data, aes(x = OR, y = Group)) + 
    geom_vline(xintercept = 1, linewidth = 0.25, linetype = "dashed") + 
    geom_errorbarh(aes(xmax = CI_high, xmin = CI_low), size = .5, height = .2, color = "gray50") +
    geom_point(size = 2.5, color = "orange") +
    scale_x_continuous(breaks = seq(0.2, 4.0, 0.2), labels = seq(0.2, 4.0, 0.2), limits = c(0.5,4.0)) +
    theme_bw() +
    theme(panel.grid.minor = element_blank()) +
    ylab("") +
    xlab("Adjusted odds ratio")
  
  # 保存图片，这里用png格式，也可替换为pdf
  ggsave(paste0(plot_name, "_forest.png"), plot = p, width = 8, height = 6, dpi = 300)
  return(p)
}

如果不同数据集的OR范围差异较大，不想固定x轴，可将scale_x_continuous的limits改为自适应：

scale_x_continuous(breaks = seq(0.2, max(data$CI_high)*1.1, 0.2), 
                   limits = c(min(data$CI_low)*0.9, max(data$CI_high)*1.1))

第三步：循环处理所有文件 #

用for循环遍历每个文件，完成读数据、绘图、保存的流程：

# 逐个处理列表中的文件
for (file in file_list) {
  # 读取当前文件数据
  lr_data <- read.table(file, sep = "\t", header = TRUE)
  # 提取纯文件名（去掉路径和后缀）作为图片名称前缀
  plot_prefix <- tools::file_path_sans_ext(basename(file))
  # 调用绘图函数生成并保存图片
  plot_forest(lr_data, plot_prefix)
}

运行后，每个数据对应的森林图会以xxx_forest.png的格式保存在当前工作目录，也可在ggsave中指定保存路径，比如paste0(data_dir, plot_prefix, "_forest.png")。

小提示 #

• 确保所有文件表头完全一致，避免读取或绘图报错；
• 可调整ggsave的width、height参数控制图片长宽比例；
• 若想在循环中实时查看图片，可在函数内添加print(p)。

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