如何用ggplot绘制带分组X轴的qPCR数据柱状图？

核心结论 #

完全可以实现目标样式，只需针对性调整ggplot的主题、元素与统计参数即可。

关键调整方向及代码示例 #

先模拟一组qPCR常见格式的示例数据：

# 模拟qPCR示例数据（包含分组、基因、Ct值、相对表达量）
qpcr_data <- data.frame(
  Group = rep(c("Control", "Treatment1", "Treatment2"), each = 3),
  Gene = rep(c("GeneA", "GeneB", "GeneC"), 3),
  Ct_Value = c(22.1, 22.3, 22.2, 20.5, 20.7, 20.6, 18.9, 19.1, 19.0),
  Expression = c(1.0, 1.1, 0.9, 2.3, 2.5, 2.4, 5.1, 5.3, 5.2)
)

1. 复刻Excel柱状图的核心代码 #

library(ggplot2)

# 先分组统计均值与标准差（用于误差线）
summary_data <- aggregate(Expression ~ Group + Gene, qpcr_data, 
                          FUN = function(x) c(mean = mean(x), sd = sd(x)))
summary_data <- do.call(data.frame, summary_data)

# 绘制复刻Excel风格的qPCR柱状图
ggplot(summary_data, aes(x = Gene, y = Expression.mean, fill = Group)) +
  # 柱状图：调整宽度、添加黑色边框
  geom_col(position = position_dodge(width = 0.8), width = 0.7, color = "black") +
  # 误差线：与柱状图对齐，设置宽度
  geom_errorbar(aes(ymin = Expression.mean - Expression.sd, 
                    ymax = Expression.mean + Expression.sd),
                position = position_dodge(width = 0.8), width = 0.2) +
  # 数据标签：柱子上方显示保留1位小数的均值
  geom_text(aes(label = round(Expression.mean, 1)), 
            position = position_dodge(width = 0.8), vjust = -0.5) +
  # 主题调整：模拟Excel默认网格与坐标轴风格
  theme_classic() +
  theme(
    panel.grid.major.y = element_line(color = "#E0E0E0", linetype = "solid"),
    axis.text = element_text(color = "black", size = 10),
    axis.title = element_text(color = "black", size = 12),
    legend.position = "top",
    legend.title = element_text(size = 11)
  ) +
  # 坐标轴与图例标题
  labs(x = "基因名称", y = "相对表达量", fill = "处理组") +
  # 可选：使用接近Excel的配色方案
  scale_fill_brewer(palette = "Set2")

关键函数与参数说明 #

• geom_col()：绘制柱状图，position_dodge()实现分组柱子并排，width控制柱子宽度，color设置柱子边框
• geom_errorbar()：添加误差线，需提前计算均值和标准差/标准误，position_dodge()参数要与柱状图保持一致以对齐
• geom_text()：在柱子上方添加数据标签，vjust调整垂直位置避免遮挡柱子
• theme_classic() + 自定义theme：通过panel.grid.major.y添加Excel风格的Y轴网格，统一文字大小与颜色
• aggregate()：快速完成分组统计，适合qPCR数据的均值与误差值计算

细节优化扩展 #

如果需要更贴近目标Excel图的细节：

• 用scale_fill_manual()自定义填充色，完全匹配Excel的配色
• 通过element_text(family = "Arial")设置与Excel一致的字体
• 用geom_hline()添加基准线（比如Control组的1倍表达量）

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