lme中的R图边际效应

我有以下lme模型（在3个级别上随机截取）：

• Y=从属V
• V1-V3=单个级别变量（级别1）
• V4=国家级变量（3级）

V1的随机斜率 V1和V4的跨级别交互

modelCI1= lme(Y~ V1+V2+V3+V1*V4+V4,
            data=total,
            method="ML",
            na.action="na.omit",
            random=list(~1|level2,~V1|level3),
            control=lmeControl(msMaxIter = 200))

我想为V4的每个级别绘制V1的边际效应

试试这个。这应该起作用：

library(ggeffects)
p <- ggpredict(nodelCI1, c("V1", "V4"))
plot(p)

库（ggeffects）
P