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)

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