如何修复循环生成ggplot高级统计柱状图时的报错问题

原代码存在嵌套循环冗余、变量引用错误、参数过时等问题，以下是修复后的完整代码，可实现批量生成包含均值、标准差、样本量、显著性字母的柱状图：

library(ggplot2)
library(agricolae)
library(dplyr)
library(tibble)

# 数据预处理，避免使用attach引发变量冲突
dat <- mtcars
dat$cyl <- factor(dat$cyl)

selected_names <- c("mpg", "disp", "hp", "drat", "wt", "qsec")
ylab_names <- c(expression("MPG (unit "[A]* ")"), 
                expression("disp (unit "[B]* ")"), 
                expression("hp (unit "[C]* ")"),  
                expression("drat (unit "[D]* ")"), 
                expression("wt (unit "[E]* ")"), 
                expression("qsec (unit "[F]* ")"))

# 按索引循环，让响应变量与y轴标签一一对应
for (idx in seq_along(selected_names)) {
  i <- selected_names[idx]
  j <- ylab_names[idx]
  
  # 拟合线性模型并执行多重比较
  model_1 <- lm(dat[[i]] ~ cyl, data = dat)
  out <- HSD.test(model_1, "cyl", group = TRUE)
  out_names <- out$groups %>%
    rownames_to_column("cyl") %>%
    arrange(match(cyl, levels(dat$cyl))) %>%
    mutate(groups = as.character(groups))
  
  statistics_letter <- out_names$groups
  
  # 计算y轴上限：均值+2倍标准差，适配误差棒范围
  ylimit_max <- mean(dat[[i]], na.rm = TRUE) + 2 * sd(dat[[i]], na.rm = TRUE)
  
  # 汇总分组统计量
  df <- dat %>%
    group_by(cyl) %>%
    summarise(
      n = n(),
      mean = round(mean(.data[[i]], na.rm = TRUE), 2),
      sd = sd(.data[[i]], na.rm = TRUE),
      .groups = "drop"
    )
  
  # 绘制图表
  rplot <- ggplot(df, aes(x = cyl, y = mean, fill = cyl)) +
    geom_bar(stat = "identity", lwd = 1.25, width = 0.6) +
    geom_errorbar(aes(ymin = mean - sd, ymax = mean + sd),
                  width = .2, lwd = 0.5, col = "black", na.rm = TRUE) +
    # 替换过时的fun.y参数为当前支持的fun，用after_stat获取统计值
    stat_summary(aes(label = after_stat(y)), fun = mean, 
                 geom = "text", colour = "darkblue", size = 5, 
                 position = position_dodge(0.6), vjust = -1) +
    geom_text(aes(label = statistics_letter, y = mean + sd), 
              position = position_dodge(0.6), color = "black", size = 6, vjust = -0.5) +
    geom_text(aes(y = max(mean)/10, label = paste0("n = ", n)), 
              position = position_dodge(0.6), size = 5, colour = "darkred") +
    labs(x = "cyl", y = j) +
    scale_y_continuous(limits = c(0, ylimit_max)) +
    theme_bw() # 添加主题优化图表美观度
  
  # 保存生成的图表
  ggsave(rplot, file = paste0(i, ".png"), width = 10, height = 8.5, units = "in")
}

关键修改说明 #

• 移除冗余嵌套循环：selected_names与ylab_names为一一对应关系，改用索引循环避免重复生成图表。
• 修正列引用逻辑：用dat[[i]]和.data[[i]]（dplyr环境中）正确调用响应变量列，替代原代码中错误的paste0(i)字符串引用。
• 修复y轴标签调用：直接使用j调用ylab_names中的表达式对象，而非错误的paste0("j")字符串。
• 调整y轴上限计算：基于真实数据列计算均值与标准差，设置为均值+2倍标准差适配误差棒显示范围。
• 更新过时参数：将stat_summary中的fun.y替换为当前支持的fun，用after_stat(y)获取统计值。
• 统一位置对齐：position_dodge的数值与geom_bar的width保持一致，避免文本与柱子错位。
• 移除attach：避免全局变量命名冲突，直接通过dat引用数据更安全。

内容的提问来源于stack exchange，提问作者Baltazár Tivadar