在ggplot中按多水平多因子分组绘制均值与置信区间图

我尝试为数据集绘制均值与置信区间图，但需要按两个不同因子分组。我设置了4种处理（Control、Virgin、Low EE2、High EE2），并在两个不同日期（D14、D21）测量了这4种处理下的总移动距离。同时希望各因子按上述指定顺序排列。我是R语言新手，恳请帮助！

当前使用的代码 #

f <- factor(c("Control", "Virgin", "Low EE2", "High EE2"), levels=c("D14", "D21"))

ggplot(Working_Exp_1_Data_summary, aes(x=Juvenile.Treatment, y=Distance.moved, color=Juvenile.Treatment, fct_relevel(f))) +
  geom_errorbar(aes(ymin=Distance.moved-ci, ymax=Distance.moved+ci), width=0.2) +
  geom_line() +
  geom_point()

当前生成的图 #

修正方案 #

原代码的问题 #

• 因子定义错误：把处理组名称和日期搞混了，f的levels应该对应处理组顺序，日期是独立的分组因子
• 没有利用日期因子做分组，导致折线无法正确连接不同日期的同处理组数据点
• aes里的fct_relevel(f)位置错误，且未起到调整顺序的作用

修正步骤与代码 #

1. 先正确设置两个因子的顺序（确保图表按指定顺序展示）

# 设置处理组的显示顺序
Working_Exp_1_Data_summary$Juvenile.Treatment <- factor(
  Working_Exp_1_Data_summary$Juvenile.Treatment,
  levels = c("Control", "Virgin", "Low EE2", "High EE2")
)

# 设置日期的显示顺序
Working_Exp_1_Data_summary$Date <- factor(
  Working_Exp_1_Data_summary$Date,
  levels = c("D14", "D21")
)

2. 修正ggplot代码，实现双因子分组可视化

library(ggplot2)

ggplot(Working_Exp_1_Data_summary, aes(x = Juvenile.Treatment, y = Distance.moved)) +
  # 绘制置信区间误差棒，按日期分组并错开位置避免重叠
  geom_errorbar(
    aes(ymin = Distance.moved - ci, ymax = Distance.moved + ci, group = Date, color = Date),
    width = 0.2,
    position = position_dodge(0.5)
  ) +
  # 用折线连接同一处理组不同日期的均值点
  geom_line(
    aes(group = Juvenile.Treatment, color = Date),
    position = position_dodge(0.5)
  ) +
  # 绘制均值点，按日期错开位置
  geom_point(
    aes(color = Date),
    position = position_dodge(0.5),
    size = 2
  ) +
  # 强制x轴按指定处理组顺序展示
  scale_x_discrete(limits = c("Control", "Virgin", "Low EE2", "High EE2")) +
  # 添加图表标签
  labs(
    title = "不同处理组在D14、D21的总移动距离",
    x = "处理组",
    y = "总移动距离",
    color = "测量日期"
  ) +
  # 用简洁的主题优化图表样式
  theme_bw()

代码说明 #

• position_dodge(0.5)：让同一处理组下的两个日期数据点、误差棒横向错开，避免重叠
• 用Date作为color和group变量：通过颜色区分日期，同时让折线正确连接同一处理组的不同日期数据
• 提前设置因子levels：确保x轴处理组、日期图例都按你指定的顺序展示

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