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

有谁知道如何计算（或提取）mer类对象（通过lme4包获得）的杠杆作用和库克距离？我想将它们绘制出来以进行残差分析。

您应该看看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")