ggplot2：带fill、alpha与分面的小提琴图中，匹配dodge连接散点

现有示例数据集（实际规模更大）：

df <- data.frame(ID   = rep(c(1:8),each = 8 ),
                 var1 = factor(rep(rep(c('1','2'), each = 4),8)),
                 var2 = factor(rep(c("EXP", "CTR"), 32)), 
                 var3 = factor(rep(c('X','Z','Z','X'), 16)),
                 var4 = factor(rep(rep(c('G1','G1','G2','G2'),each=8),2)),
                 var5 = factor(rep(rep(c('GA','GB','GB','GA'),each=8),2)),
                 value = sample(0:100, 64, rep = TRUE)) %>%
  arrange(ID,var1,var2,var3,var4,var5) %>% 
  ungroup()

需求是连接图中抖动散点的单个数据点，尤其是var3的X与Z组之间的点。当前绘图代码生成的连线位置随机，推测是fill和alpha的设置导致的问题。优先寻求无需手动在数据中设置dodge的解决方案；若不可行，需明确dodge的设置方法并提供通用公式。

当前使用的绘图代码：

#set position of the jitter
pos_jit=position_jitterdodge(jitter.width=0.3, dodge.width = 0.9, seed=1)

ggplot(df, aes(x=var1, y=value, fill=var2, alpha=var3)) +theme_classic() +
  facet_wrap(vars(var4,var5), nrow=4) +
  scale_fill_brewer(palette='Set1')+ # scale_fill_manual(values=c("#49111c","#bddbd0"))+
  scale_alpha_discrete(range = c(1, 0.3),guide = guide_legend(override.aes = list(fill = "black"))) +
  geom_violin(position = 'dodge') + 
  geom_point(               size=2, color='grey',  position=pos_jit) +
  #geom_path(aes(group=interaction(ID,var1)), position=pos_jit) +
  geom_line(aes(group=ID), position=pos_jit) +
  stat_summary(fun = "mean", geom = "crossbar", width = .5, size=0.4, color = "black", position = position_dodge(.9), show.legend=FALSE) +
  theme(legend.position= "top")+theme(legend.title=element_blank()) + 
  ylab('RT in ms') + xlab('') 

优先方案：统一分组与位置映射 #

问题核心是geom_line的group仅设为ID，但position_jitterdodge会根据fill（var2）和alpha（var3）进行分组偏移，导致同一ID的点被分到不同偏移组，连线错位。

修改方法：

1. 将geom_line的group设置为interaction(ID, var2, var3)，确保同一ID、同一var2+var3组合的点归为同一组，跟随对应偏移位置。
2. 保留position_jitterdodge的固定种子，保证抖动位置一致。

修改后的代码：

pos_jit=position_jitterdodge(jitter.width=0.3, dodge.width = 0.9, seed=1)

ggplot(df, aes(x=var1, y=value, fill=var2, alpha=var3)) +theme_classic() +
  facet_wrap(vars(var4,var5), nrow=4) +
  scale_fill_brewer(palette='Set1')+
  scale_alpha_discrete(range = c(1, 0.3),guide = guide_legend(override.aes = list(fill = "black"))) +
  geom_violin(position = 'dodge') + 
  geom_point(size=2, color='grey', position=pos_jit) +
  # 修正group为ID+var2+var3的交互项
  geom_line(aes(group=interaction(ID, var2, var3)), position=pos_jit) +
  stat_summary(fun = "mean", geom = "crossbar", width = .5, size=0.4, color = "black", position = position_dodge(.9), show.legend=FALSE) +
  theme(legend.position= "top")+theme(legend.title=element_blank()) + 
  ylab('RT in ms') + xlab('') 

此方法无需修改原始数据，直接通过调整group参数让连线跟随散点的偏移逻辑，解决错位问题。

备选方案：手动计算dodge位置 #

如果上述方法不适用，可手动给每个分组计算偏移后的x坐标，替代position_jitterdodge：

步骤1：计算分组偏移量 #

1. 将var1转换为数值型x_num
2. 生成var2+var3的组合分组，计算每个分组的偏移索引
3. 添加固定种子的随机抖动，保证位置一致性

代码实现 #

set.seed(1)
df <- df %>%
  mutate(
    x_num = as.numeric(var1),
    # 生成var2和var3的组合分组
    dodge_group = interaction(var2, var3),
    # 获取分组数量，计算每个分组的偏移量
    n_dodge = n_distinct(dodge_group),
    dodge_pos = (as.numeric(dodge_group) - (n_dodge + 1)/2) * 0.9 / n_dodge,
    # 添加抖动（jitter.width=0.3对应±0.15）
    x_jitter = x_num + dodge_pos + runif(n(), -0.15, 0.15)
  )

# 绘图时用手动计算的x_jitter作为x轴
ggplot(df, aes(x=x_jitter, y=value, fill=var2, alpha=var3)) +theme_classic() +
  facet_wrap(vars(var4,var5), nrow=4) +
  scale_fill_brewer(palette='Set1')+
  scale_alpha_discrete(range = c(1, 0.3),guide = guide_legend(override.aes = list(fill = "black"))) +
  # 小提琴图基于原始var1，单独设置x和position
  geom_violin(aes(x=x_num), position=position_dodge(0.9)) + 
  geom_point(size=2, color='grey') +
  geom_line(aes(group=ID)) +
  stat_summary(aes(x=x_num), fun = "mean", geom = "crossbar", width = .5, size=0.4, color = "black", position = position_dodge(.9), show.legend=FALSE) +
  # 还原x轴刻度为原始var1标签
  scale_x_continuous(breaks=unique(df$x_num), labels=unique(df$var1)) +
  theme(legend.position= "top")+theme(legend.title=element_blank()) + 
  ylab('RT in ms') + xlab('') 

通用公式说明 #

手动dodge偏移量公式：

dodge_pos = (分组索引 - (分组总数 + 1)/2) * dodge_width / 分组总数

参数说明：

• 分组索引：每个dodge_group的数值型编号（从1开始）
• 分组总数：需要dodge的分组数量（即var2×var3的组合数）
• dodge_width：与position_jitterdodge中参数一致（通常设为0.9）

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