R语言按Periodo分组过滤循环生成直方图的异常排查

你的代码会出现生成图片数量远超预期的问题，核心有3处逻辑错误：

• 循环遍历范围设置错误
  你写的for (i in 1:length(data_tot$Periodo)) 是沿着数据集的每一行触发循环，不是遍历去重后的Periodo值。如果data_tot有上千条记录，循环就会执行上千次，完全符合你观察到的“按每行数据生成一张图”的现象，也是图片数量超标最直接的原因。
• 管道传递的过滤结果未被使用
  前面通过filter(Periodo == 对应值)得到的分组数据，在写ggplot时被你硬编码指定读取全量data_tot，等于过滤操作完全失效，不管过滤条件是什么，画出来的永远是全量数据的直方图，没有分组统计效果。
• 统计值取数范围错误
  绘制参考线、标注均值时，你直接取了全量data_tot$b1的均值，不是当前过滤分组内的b1均值，哪怕前面的过滤逻辑生效，最终画出来的参考线位置、标注数值也不是对应分组的结果。

修正后的可运行代码 #

# 先提取去重后的Periodo值作为循环遍历对象，从根源上避免按行循环的问题
target_periods <- unique(data_tot$Periodo)

for (year in target_periods) {
  # 单独过滤出当前分组的数据存为对象，避免管道传参混乱
  sub_group <- data_tot %>%
    filter(Periodo == year)
  
  # 明确传入过滤后的分组数据作为ggplot的输入源
  p <- ggplot(sub_group, aes(x = b1)) +
    geom_histogram(fill = "steelblue", aes(y = ..density.., alpha = ..count..), bins = 60) +
    labs(title = "Beneficios", subtitle = paste0("Año ", year)) +
    ylab("Densidad") +
    xlab("Beneficios ($millones)") +
    # 用当前分组的均值画参考线
    geom_vline(aes(xintercept = mean(b1)), color = "red4", linetype = "dashed") +
    theme(legend.position = "none") +
    annotate("text",
             x = mean(sub_group$b1),
             y = 0,
             label = round(mean(sub_group$b1), 2),
             color = "red4",
             size = 3.5,
             vjust = -1.5,
             hjust = -0.5)
  
  # 用实际年份做文件名，比循环序号更方便后续识别
  ggsave(p, filename = paste0(year, "_histograma.png"))
}

运行前可以打印一行print(target_periods) 确认提取到的周期值刚好覆盖2020-2029、没有缺失值或异常值，避免生成无意义的图片。

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