識別lme4模型中的異常值

Question

我想識別lmer模型中的異常值（lme4包）。我是而不是有興趣去除它們（什麼LMERConvenienceFunctions包） - 我只是想看到列出的異常值。

一個模式，我用的例子：

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

Answer 1

如果你想以同樣的方式是LMERConvenienceFunctions做（1標準化水平殘差> 0 2.5SD）來定義異常值，那麼這將工作：

如果你有 plants1丟失的數據

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

確保提取的殘留物之前設置na.action = 'na.exclude'在lmer命令。這將確保將缺失數據的觀測值設置爲NA而不是忽略它們。