R在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)

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