如何在ggplot中批量绘制多列小提琴图并批量保存PDF

核心实现步骤 #

• 加载绘图依赖包
• 创建输出文件夹（避免路径不存在报错）
• 遍历所有Q列生成并保存图表

完整代码示例 #

假设你的数据框名为df，首列是SYSTEM，其余为Q1至Q100的0/1数值列：

# 加载必要包
library(ggplot2)

# 设置输出文件夹路径
output_dir <- "小提琴图结果"

# 若文件夹不存在则创建，recursive参数支持创建多级目录
if (!dir.exists(output_dir)) {
  dir.create(output_dir, recursive = TRUE)
}

# 获取所有Q列的列名（排除首列SYSTEM）
q_columns <- colnames(df)[-1]

# 循环处理每个Q列
for (col_name in q_columns) {
  # 提取当前Q列和SYSTEM列的临时数据
  temp_data <- df[, c("SYSTEM", col_name)]
  # 重命名列，简化ggplot语法（避免动态变量名报错）
  colnames(temp_data) <- c("SYSTEM", "value")
  
  # 绘制小提琴图，搭配箱线图增强可读性（0/1数据的小提琴图形态较特殊）
  plot <- ggplot(temp_data, aes(x = SYSTEM, y = value)) +
    geom_violin(fill = "#619CFF", alpha = 0.6) +
    geom_boxplot(width = 0.12, fill = "white", outlier.size = 1.5) +
    labs(
      title = col_name,
      x = "系统分类",
      y = "0/1数值"
    ) +
    theme_minimal() +
    theme(
      plot.title = element_text(hjust = 0.5, size = 14, face = "bold"),
      axis.title = element_text(size = 12)
    )
  
  # 保存为PDF文件，文件名对应Q列名
  save_path <- file.path(output_dir, paste0(col_name, ".pdf"))
  ggsave(
    filename = save_path,
    plot = plot,
    width = 7,
    height = 5,
    dpi = 300
  )
}

关键细节说明 #

• 数据处理：每次循环提取当前Q列和SYSTEM列并重新命名，避免ggplot中动态变量名引用的报错（这是新手代码失效的常见原因）
• 图表优化：针对0/1二元数据，添加箱线图可直观展示每组的中位数、四分位数，弥补小提琴图形态扁平的不足
• 路径处理：使用file.path()生成跨系统兼容的文件路径，避免Windows和Linux路径格式冲突
• 保存参数：通过ggsave()的width/height/dpi控制输出PDF的尺寸和清晰度

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