我使用lmer（）构建了一个混合效果模型，目的是比较四个不同年龄组儿童的阅读分数增长情况。 我想在R中以置信区间绘制4个不同斜率的图形，以便可视化这种关系，但我一直陷于困境。
我已经尝试像以前的lm（）模型一样使用plot函数和ggplot的某些版本，但是到目前为止，它仍然无法正常工作。这是我尝试过的模型，我希望该模型可以查看随着时间（年龄）的阅读分数的变化如何与孩子的SESDLD分组交互（这表明孩子是否有语言问题，以及他们是高收入还是低收入）。
AgeSES.model <- lmer(ReadingMeasure ~ Age.c*SESDLD1 + (1|childid), data = reshapedomit, REML = FALSE)
ReadingMeasure是一个连续的得分，age.c是以月为单位的居中年龄。 SESDLD1是一种有4个级别的分类度量。我希望ReadingMeasure的增长有四个正斜率，截距不同，斜率也可能不同。
我非常感谢任何有关如何执行此操作的指示！
非常感谢您！
The type of plot I would like to achieve - this was done in Stata
