如何修改R代码为ggplot2柱状图添加第三个emotionality子图

核心思路 #

原有代码基于宽格式数据生成两个facet子图，新增变量后需要将数据重塑为长格式，让ggplot自动识别第三个特质（emotionality），再调整facet布局实现垂直堆叠。

具体修改步骤 #

假设你的数据框名为dt，包含新增的rG_emotionality、rE_emotionality及对应误差列（如se_rG_emotionality、se_rE_emotionality），按以下步骤修改：

1. 将宽格式数据转为长格式
   用data.table的melt函数整合三类特质变量：

library(data.table)
library(ggplot2)

# 重塑数据：提取特质、分组、数值及误差
dt_long <- melt(dt, 
                id.vars = "age",
                measure.vars = patterns(rG = "^rG_", rE = "^rE_", se_rG = "^se_rG_", se_rE = "^se_rE_"),
                variable.name = "trait",
                value.name = c("rG_val", "rE_val", "se_rG", "se_rE"))

# 拆分出rG/rG分组，统一数值与误差列
dt_long <- dt_long[, .(age, trait,
                       group = rep(c("rG", "rE"), each = .N),
                       value = c(rG_val, rE_val),
                       se = c(se_rG, se_rE))]

2. 修改ggplot代码生成垂直堆叠的三个子图
   替换原有硬编码逻辑，基于trait列实现facet映射：

ggplot(dt_long, aes(x = age, y = value, fill = group)) +
  geom_col(position = position_dodge(0.8), width = 0.7) +
  # 添加误差棒
  geom_errorbar(aes(ymin = value - se, ymax = value + se),
                position = position_dodge(0.8), width = 0.2) +
  # 添加数值标签
  geom_text(aes(label = round(value, 2)),
            position = position_dodge(0.8), vjust = -0.5, size = 3) +
  # 垂直堆叠三个特质子图
  facet_grid(trait ~ ., scales = "free_y") +
  labs(x = "年龄", y = "数值", fill = "分组") +
  theme_bw()

关键说明 #

• 若误差列命名规则不同，只需调整melt函数中se_rG、se_rE的匹配正则即可。
• scales="free_y"让每个子图y轴自适应数值范围，若需统一y轴，删除该参数即可。

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