如何用ggplot对同一数据框循环生成多张绘图并导出为PDF

前置准备 #

先补全原代码依赖（示例中用到的randomStrings函数来自random包，需提前安装加载）：

# 加载所有依赖
library(random)
library(dplyr)
library(ggplot2)
library(ggrepel)

# 生成示例数据（和你提供的逻辑一致）
set.seed(15)
res_df <- data.frame(exp = randomStrings(N = 148, string_size = 4))
res_df <- data.frame(res_df[rep(seq_len(nrow(res_df)), each = 148), ])
colnames(res_df)[1] <- "exp"
res_df <- mutate(res_df, y = randomStrings(N = 148, string_size = 5),
                 logp = abs(rnorm(n = 148, mean = 5, sd = 6)),
                 r = rnorm(n = 148, mean = 0.5, sd = 0.1))

步骤1：封装绘图函数 #

把单张图的生成逻辑封装为可复用函数，动态传入exposure名称：

plot_single_exp <- function(current_exp, df) {
  # 筛选当前exposure的所有数据
  plot_df <- filter(df, exp == current_exp)
  
  ggplot(plot_df, aes(x = r, y = logp, label = y)) +
    geom_point(data = filter(plot_df, logp < 10), color = "grey50") +
    geom_text_repel(data = filter(plot_df, logp > 10), box.padding = 0.5, max.overlaps = Inf) +
    geom_point(data = filter(plot_df, logp > 10), color = "red") +
    xlab("Variance explained (%)") + 
    ylab("-log10(pvalue)") +
    # 动态设置标题，显示当前exposure名称
    ggtitle(paste0("Outcome ~ ", current_exp))
}

步骤2：批量导出方案 #

方案A：导出为单份多页PDF（所有图在同一个PDF文件，每页1张） #

直接调用PDF绘图设备循环打印即可，代码最简单，内存占用低：

# 打开PDF设备，每页尺寸可按需调整
pdf("所有exposure分析结果.pdf", width = 8, height = 6)

# 循环遍历所有exposure生成并打印图
all_exp <- unique(res_df$exp)
for (exp_name in all_exp) {
  p <- plot_single_exp(exp_name, res_df)
  print(p)
}

# 关闭设备，完成导出
dev.off()

方案B：导出为独立PDF文件（每个exposure对应1个单独PDF） #

命名时带上exposure名称避免重名，可指定导出文件夹：

# 可选：新建文件夹存所有独立PDF，避免文件散乱
dir.create("单独exposure结果", showWarnings = FALSE)

all_exp <- unique(res_df$exp)
for (exp_name in all_exp) {
  p <- plot_single_exp(exp_name, res_df)
  # 自定义文件名，替换特殊字符避免导出报错
  file_name <- paste0("单独exposure结果/", gsub("[^[:alnum:]]", "_", exp_name), "_分析结果.pdf")
  ggsave(file_name, p, width = 8, height = 6, dpi = 300)
}

优化提示 #

如果数据量很大，不需要后续修改绘图对象，可以直接在循环中生成并导出，不需要把所有148个ggplot对象都存在内存中，能大幅降低内存占用。如果需要后续统一调整所有图，也可以用purrr::map(all_exp, ~plot_single_exp(.x, res_df))把所有图存在列表中批量处理。

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