ggplot2按州批量生成分区指标图并保存PDF报错求助

问题背景 #

有一个包含state（州）、district（分区）及数值变量的数据集，共32个州、712个分布不均的分区。需求是按州生成包含下属分区指标的图表，自动添加州名标题，并将每个图表保存为PDF到指定文件夹。原模拟代码运行报错，同时不清楚如何实现自动保存。

报错原因分析 #

原代码中使用purrr::map时，直接传入ggplot(., aes(...))的写法有误——map的参数需要是一个函数，而不是直接调用ggplot。group_split返回的是数据框列表，直接把列表元素传给ggplot会导致数据格式不匹配，触发fortify()的报错。

修正后的完整代码 #

library(tidyverse)

# 生成模拟数据
set.seed(321)
statenames <- c("state1","state2","state1","state2","state1","state2")
distnames <- c("dist1","dist2","dist3","dist4","dist5","dist6")
hcrv <- rnorm(n=6,mean = 10, sd=2)
intv <- rnorm(n=6,mean = 100, sd=20)
rdiv <- rnorm(n=6,mean = 150, sd=2)

plotdata <- cbind.data.frame(statenames, distnames, hcrv,intv,rdiv)

# 创建输出文件夹（如果不存在）
output_dir <- "state_plots"
if (!dir.exists(output_dir)) {
  dir.create(output_dir)
}

# 批量生成并保存图表
plotdata %>%
  pivot_longer(cols = -c(statenames, distnames), 
               names_to = "indicatorName", 
               values_to = "indicatorValue") %>%
  group_split(statenames) %>%
  purrr::walk(function(state_data) {
    # 获取当前州的名称
    state_name <- unique(state_data$statenames)
    
    # 绘制图表
    p <- ggplot(state_data, aes(x=indicatorName, y=indicatorValue)) +
      geom_col() +
      facet_wrap(~distnames) +
      labs(title = state_name) +
      theme_bw() # 可选：添加主题美化
    
    # 保存为PDF
    ggsave(
      filename = file.path(output_dir, paste0(state_name, ".pdf")),
      plot = p,
      device = "pdf",
      width = 8,
      height = 6
    )
  })

代码说明 #

• group_split(statenames)：将数据按州拆分为数据框列表，每个元素对应一个州的所有数据
• purrr::walk：遍历列表执行绘图和保存操作，适合无返回值的场景，比map更高效
• 动态获取州名：通过unique(state_data$statenames)自动提取当前州名称，用于图表标题和文件名
• ggsave：指定输出路径、文件名，设置图表尺寸；file.path确保路径在不同操作系统下正常工作
• 可选：添加theme_bw()等主题函数美化图表

关键修正点 #

• 将map改为传入匿名函数，确保每个州的数据框被正确传给ggplot
• 使用walk替代map，避免不必要的对象返回
• 新增文件夹创建逻辑，解决保存时路径不存在的错误
• 动态生成州名标题和文件名，实现完全自动化

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