R 用于说明线性混合效应模型结果的图-续

我正在尝试创建一个类似于此处的绘图：

（为便于参考，以下是该链接的代码：

library(nlme)
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1|Subject)

newdat <- expand.grid(Sex=unique(Orthodont$Sex),
                  age=c(min(Orthodont$age),
                        max(Orthodont$age)))

library(ggplot2)
p <- ggplot(Orthodont, aes(x=age, y=distance, colour=Sex)) +
  geom_point(size=3) +
  geom_line(aes(y=predict(fm2), group=Subject, size="Subjects")) +
  geom_line(data=newdat, aes(y=predict(fm2, level=0, newdata=newdat), size="Population")) +
  scale_size_manual(name="Predictions", values=c("Subjects"=0.5, "Population"=3)) +
  theme_bw(base_size=22) 
print(p)

我在x轴上绘制刺激性曲线，在y轴上绘制反应曲线。我想要每个参与者的单独线条（按参与者性别进行颜色编码）和两条更粗的线条，显示男性和女性参与者的总体平均值。以下是我的代码：

lmer1<- lmer(Response~StimSex*Sex + (Sex|StimID) + (StimSex|Subject), data = dat)


p1<- ggplot(dat, aes(x=StimSex, y = Response, colour = Sex))+
  geom_point(size = 1) +
  geom_line(aes(y=predict(lmer1), group = Subject, size = "Subjects"))+
  geom_line(aes(y=predict(lmer1), group = Sex, size = "Population Average"))+
  scale_size_manual(name="Predictions", values=c("Subjects"=0.5, "Population Average"=2))

p1

---

lmer1'How to ____; in R'问题通常与主题无关。它们可以在主题上（例如，这一个应该是）。如果您可以等待，我们可以为您迁移它。阅读
帮助（“predict.merMod”）
。您需要
re.form
参数。感谢您将此问题迁移到stackoverflow-抱歉造成混淆。@Roland谢谢！这很有效！如何在R中进行问题通常是离题的。它们可以在主题上（例如，这一个应该是）。如果您可以等待，我们可以为您迁移此问题。请阅读
帮助（“predict.merMod”）
。您需要
re.form
参数。感谢您将此问题迁移到stackoverflow-很抱歉造成混淆。@Roland谢谢！这很有效！
lmer1<- lmer(Response~StimSex*Sex + (Sex|StimID) + (StimSex|Subject), data = dat)


p1<- ggplot(dat, aes(x=StimSex, y = Response, colour = Sex))+
  geom_point(size = 1) +
  geom_line(aes(y=predict(lmer1), group = Subject, size = "Subjects"))+
  geom_line(aes(y=predict(lmer1), group = Sex, size = "Population Average"))+
  scale_size_manual(name="Predictions", values=c("Subjects"=0.5, "Population Average"=2))

p1