R语言中为分组柱状图添加误差棒时遇因子运算错误

你遇到的警告核心原因是数据源不匹配+数据逻辑错误：

• 你用原始数据集R做geom_bar(stat="identity")，实际是把所有原始int值叠加展示，而非你计算的均值；
• geom_errorbar里直接用原始列int和整个汇总表sd_int做运算，不仅单位和数据结构不兼容，还触发了因子运算的警告。

直接用你已经汇总好的combined_data画图，确保柱状图展示均值、误差棒对应标准差：

# 你已完成的合并步骤保留
combined_data <- merge(mean_int, sd_int, by = c("PSm", "Day"))

ggplot(combined_data, aes(x = factor(PSm), y = mean/60, fill = factor(PSm))) +
  geom_bar(stat = "identity") +
  # 误差棒和y轴单位对齐（转成分钟）
  geom_errorbar(aes(ymin = (mean - sd)/60, ymax = (mean + sd)/60), 
                width = 0.15, size = 0.8) +
  facet_wrap(~Day) +
  theme_gray() +
  guides(fill = FALSE) +
  scale_fill_gradient(low = "#CD6600", high = "#EEAD0E") +
  ylab("Interaction duration in minutes") +
  # 直接设置x轴标签，替代annotate更规范
  scale_x_discrete(NULL, labels = c("no problemsolver", "problemsolver")) +
  ggtitle("Interactions between males and females by \nmale problem solving ability per day") +
  theme(plot.title = element_text(hjust = 0.5, size = 19),
        axis.text = element_text(size = 10.3),
        axis.title.y = element_text(size = 18),
        strip.text = element_text(size = 13.5))

1. 数据源替换：将ggplot的数据源从原始数据R换成汇总后的combined_data，确保柱状图展示你计算的均值，而非原始值叠加。
2. 误差棒单位对齐：把标准差sd也除以60，和y轴的分钟单位保持一致，避免数值偏差。
3. x轴标签优化：用scale_x_discrete直接设置x轴类别标签，替代原来的annotate，避免位置偏移问题。
4. 因子类型处理：将PSm转为因子，确保x轴两个类别匹配正确，消除因子运算警告。

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