如何在R语言中基于结果数据绘制带误差棒的分组柱状图

问题原因 #

你的代码出错来自三个核心问题：

1. 原始数据中的ph.mean、ph.sd、hc.mean、hc.sd列均为字符类型，无法直接用于数值计算和绘图
2. 调用gather转换长表时，指定的phm、hcm列不存在于原数据框，无法正确提取均值数据
3. 仅转换了均值列，没有将标准差和对应分组、分类绑定，误差棒缺少匹配的数据源

解决代码 #

首先确保你已经加载tidyverse包（包含dplyr、tidyr、ggplot2核心功能）：

library(tidyverse)

第一步：修正数据类型+转换长表 #

# 原始数据
df<-data.frame(
  chq=c("Eb","Pb","Sb","En","Bn"),
  ph.mean=c("9.8","8.0","7.9","13.9","10.0"),
  ph.sd=c("3.3","4.2","3.3","6.3","4.2"),
  hc.mean=c("6.6","7.4","6.1","11.1","8.3"),
  hc.sd=c("1.3","2.6","2.6","3.4","1.5"),
  p.val=c("0.001","0.120","0.098","0.021","0.054")
)

# 转换数值列+长表重构
df.long <- df %>%
  # 将所有均值、标准差列从字符转为数值类型
  mutate(across(c(ph.mean:hc.sd), as.numeric)) %>%
  # 保留分类列chq、显著性列p.val，其余列转长格式
  pivot_longer(cols = -c(chq, p.val),
               # 拆分原列名，提取分组名（ph/hc）和指标类型（mean/sd）
               names_to = c("group", ".value"),
               names_sep = "\\.")

转换后的df.long每一行对应1个chq分类+1个分组（ph/hc）的组合，同时包含对应的均值mean和标准差sd，可直接用于绘图。

第二步：绘制带误差棒的分组柱状图 #

ggplot(df.long, aes(x = chq, y = mean, fill = group)) +
  # 分组并列柱状图，加黑色边框区分
  geom_col(position = position_dodge(0.9), color = "black") +
  # 误差棒，范围为均值±标准差，和柱状图位置对齐
  geom_errorbar(aes(ymin = mean - sd, ymax = mean + sd),
                position = position_dodge(0.9), width = 0.25) +
  # 基础主题调整
  theme_bw() +
  labs(x = "chq分类", y = "测量值", fill = "分组")

关键逻辑说明 #

• 使用pivot_longer的.value参数，可以直接将同分组的均值、标准差自动匹配到同一行，无需额外关联操作
• position_dodge(0.9)用于保证误差棒和对应柱状图的中心对齐，避免错位
• 优先使用geom_col代替geom_bar(stat = "identity")，语法更简洁

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