识别lme4模型中的异常值

Question

我想识别lmer模型中的异常值（ lme4 包）。我不对删除它们感兴趣（ LMERConvenienceFunctions 包的功能） - 我只想查看列出的异常值。

我使用的模型示例：

model1<-lmer(Value~ Moisture + Planting + (day|plot), data=plants1)

Answer 1

如果您想以与LMERConvenienceFunctions相同的方式定义异常值（标准化的1级残差＆gt; 2.5SD从0开始），那么这将有效：

res1 <- resid(model1, type = "pearson") # Extract standardized residuals
plants1[which(abs(res1) > 2.5),] # Get the rows which absolute residuals > 2.5

如果plants1中缺少数据，请确保在提取残差之前在na.action = 'na.exclude'命令中设置lmer。这将确保将缺少数据的观察设置为NA而不是省略它们。