L在lme中绘制边际效应

Question

我有以下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))

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

Answer 1

尝试ggeffects-package。这应该起作用：

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

详细信息可以在上述链接网站的文章中找到。