R语言如何通过循环/函数批量绘制多数据框的ggplot2图表

你之前用数据框列表实现的逻辑是目前批量绘图最易维护的方案，以下是可直接运行的完整实现，两种形式可选：

方案1：列表统一管理（推荐，即你之前用过的实现逻辑） #

该方案不会在全局环境生成大量离散变量，后续批量调整、导出图表效率更高。

• 第一步：加载依赖、整合散列的染色体数据框为命名列表

library(ggplot2)

# 读取全局环境中chr1~chr18数据框，整合为命名列表
chr_list <- lapply(1:18, function(i) get(paste0("chr", i), envir = .GlobalEnv))
names(chr_list) <- paste0("chr", 1:18)

• 第二步：编写通用绘图函数，自动匹配数据源和标题

# 自定义绘图函数：传入染色体编号、对应数据集，返回配置好的ggplot对象
generate_chr_plot <- function(chr_index, dataset) {
  ggplot(dataset, aes(V4, colour = factor(V6))) +
    geom_freqpoly(binwidth = 1, size = 0.8) +
    labs(
      x = "score",
      y = "Number of Mutations",
      title = paste0("Chromosome ", chr_index) # 自动生成对应染色体的标题
    ) +
    theme_bw()
}

• 第三步：批量生成所有绘图对象，存入命名列表

# 批量生成18张图，按p1~p18命名存在列表中
p_list <- mapply(
  FUN = generate_chr_plot,
  chr_index = 1:18,
  dataset = chr_list,
  SIMPLIFY = FALSE
)
names(p_list) <- paste0("p", 1:18)

调用时直接按索引取即可，比如要获取1号染色体的图，直接运行p_list$p1就能输出。如果需要批量导出所有图片，直接遍历列表即可：

# 批量导出所有图到单个PDF示例
pdf("chromosome_mut_freq.pdf", width = 8, height = 6)
invisible(lapply(p_list, print))
dev.off()

方案2：直接生成p1~p18独立全局对象 #

如果不需要列表管理，就是要在全局环境生成独立命名的p1到p18变量，可以直接用循环赋值实现：

library(ggplot2)

generate_chr_plot <- function(chr_index, dataset) {
  ggplot(dataset, aes(V4, colour = factor(V6))) +
    geom_freqpoly(binwidth = 1, size = 0.8) +
    labs(
      x = "score",
      y = "Number of Mutations",
      title = paste0("Chromosome ", chr_index)
    ) +
    theme_bw()
}

# 循环给全局环境赋值p1~p18
for (i in 1:18) {
  assign(
    x = paste0("p", i),
    value = generate_chr_plot(
      chr_index = i,
      dataset = get(paste0("chr", i), envir = .GlobalEnv)
    ),
    envir = .GlobalEnv
  )
}

运行完成后，直接在环境里就能看到p1到p18共18个独立的ggplot对象，和你手动逐张写代码生成的对象完全一致。

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