如何绘制混合效应线性回归？

Question

我正在尝试制作一个箱线图来可视化这个回归模型

library(lme4)
lmer(dv1 ~ intervention + (1|id/area),
            data=data,
            REML=T)

在这个实验中，控制和治疗干预都应用于离散区域内的受试者。

这是我正在使用的数据

data <- data.frame("id" = 1:2, 
           "intervention" = c(rep("a",27),rep("b", 27)),
           "area" = 1:3,
           "dv1" = rnorm(54),
           "dv2" = rnorm(54),
           "dv3" = rnorm(54))
data$area <- as.factor(data$area)
data$id <- as.factor(data$id)

这是我尝试过的

library(ggplot2)
ggplot(data,aes(x=area,y=dv1,col=intervention)) +
  geom_point() +
  geom_boxplot(alpha=0.2) + 
  facet_wrap(~id) + 
  ggtitle("DV1") +
  xlab("Intervention") +
  ylab("DV1")

不是覆盖在红色箱线图上的红点，而是到处都是。我该如何解决这个问题？

编辑：我使用了 u/eipi10 建议的抖动选项，这就是我现在所拥有的。

ggplot(data,aes(x=area,y=dv1,col=intervention)) +
  geom_point(position=position_jitterdodge(dodge.width=0.75, jitter.height=0, jitter.width=0.25), alpha=0.6) +
  geom_boxplot(alpha=0.2, size=0.3) + 
  facet_wrap(~id) + 
  ggtitle("DV1") +
  xlab("Area") +
  ylab("DV1")