r - R在lme中绘制边际效应

Question

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

• Y=依赖V.，
• V1-V3=个体水平变量（1级），
• V4=国家级变量（第 3 级）

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

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)

详情可在上述链接网站的文章中找到。