如何提取/计算线性混合效应模型的杠杆和库克距离

机器算法验证 r 混合模式 线性模型 残差 杠杆作用
2022-03-10 03:05:01

有谁知道如何计算(或提取)mer类对象(通过lme4包获得)的杠杆作用和库克距离?我想绘制这些以进行残差分析。

1个回答

你应该看看 R 包influence.ME它允许您计算由lme4.

示例模型:

library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)

该函数influence是所有进一步步骤的基础:

library(influence.ME)
infl <- influence(model, obs = TRUE)

计算库克距离:

cooks.distance(infl)

绘制库克距离:

plot(infl, which = "cook")

在此处输入图像描述