r - 将混合模型的预测绘制为 R 中的一条线

Question

我试图在predict()下面绘制我的混合模型的预测（

我试图绘制我的模型的预测，但我没有达到我想要的情节。有没有更好的方法来定义predict()，这样我就可以实现我想要的情节？

library(lme4)

dat3 <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/dat3.csv')

m4 <- lmer(math~pc1+pc2+discon+(pc1+pc2+discon|id), data=dat3)

newdata <- with(dat3, expand.grid(pc1=unique(pc1), pc2=unique(pc2), discon=unique(discon)))

y <- predict(m4, newdata=newdata, re.form=NA)

plot(newdata$pc1+newdata$pc2, y)

标签： rplotregressionlme4mixed-models

Answer 1

更多sjPlot。请参阅参数grid以将多个预测变量包装在一个图中。

library(lme4)
library(sjPlot)
library(patchwork)

dat3 <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/dat3.csv')

m4 <- lmer(math~pc1+pc2+discon+(pc1+pc2+discon|id), data=dat3) # Does not converge

m4 <- lmer(math~pc1+pc2+discon+(1|id), data=dat3) # Converges

# To remove discon
a <- plot_model(m4,type = 'pred')[[1]]
b <- plot_model(m4,type = 'pred',title = '')[[2]]

a + b 

编辑 1：我在框架中删除该dropcon术语时遇到了一些麻烦。sjPlot我放弃了，又重新开始了patchwork。我相信丹尼尔可以知道正确的方法。