如何用R为300个同结构数据框生成正确的分组条形图

问题背景 #

拥有300个结构完全一致的数据框，需为每个数据框生成独立分组条形图：

• X轴：Mz_Round列
• Y轴：HW Intensity和Reg Intensity两列数值

数据框示例：

precursorMz Mz_Round HW Intensity Reg Intensity diff1 diff2
1    256.6814   141.10         4216          3994  0.96  1.00
2    256.6814   142.10         7184          5988  1.00  1.02
3    256.6814   143.12        44510         30020  1.02  1.00
4    256.6814   144.12         1858          1312  1.00  0.00
5    256.6814   260.20        43010         23230  4.52  1.00
6    256.6814   261.20         9452          6388  1.00  0.99

用户尝试的错误代码：

intensities <- table(split1$`HW Intensity`, split1$`Reg Intensity`)
barplot(intensities,
        main = "Intensity Compared",
        xlab = "M/z", ylab = "Intensity",
        col = c("darkgrey", "blue"),
        rownames(split1$Mz_Round),
        beside = TRUE)

单数据框的正确绘图代码 #

错误原因是误用table()函数（该函数用于生成频数表，不适合直接绘制数值对比条形图），正确做法是将需要对比的两列转换为矩阵并转置，指定X轴标签为Mz_Round：

# 以单个数据框split1为例
# 提取数值列并转置为矩阵（行是强度类型，列是Mz_Round）
plot_matrix <- t(as.matrix(split1[, c("HW Intensity", "Reg Intensity")]))
# 设置列名为Mz_Round的值，作为X轴标签
colnames(plot_matrix) <- split1$Mz_Round

# 绘制分组条形图
barplot(plot_matrix,
        main = "Intensity Comparison",
        xlab = "M/z",
        ylab = "Intensity",
        col = c("darkgrey", "blue"),
        beside = TRUE,  # 分组并排显示条形
        legend.text = rownames(plot_matrix),  # 添加图例
        args.legend = list(x = "topright"))  # 设置图例位置

批量处理300个数据框的方法 #

假设所有数据框已存入列表df_list，可通过循环批量生成并保存图表：

# 遍历列表中的每个数据框
for (idx in seq_along(df_list)) {
  current_df <- df_list[[idx]]
  # 构建绘图矩阵
  plot_matrix <- t(as.matrix(current_df[, c("HW Intensity", "Reg Intensity")]))
  colnames(plot_matrix) <- current_df$Mz_Round
  
  # 保存为PNG文件，文件名带序号区分
  png(filename = paste0("intensity_plot_", idx, ".png"), width = 800, height = 600)
  barplot(plot_matrix,
          main = paste("Intensity Comparison - No.", idx),
          xlab = "M/z",
          ylab = "Intensity",
          col = c("darkgrey", "blue"),
          beside = TRUE,
          legend.text = rownames(plot_matrix),
          args.legend = list(x = "topright"))
  dev.off()  # 关闭绘图设备
}

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