lme4モデルの外れ値の特定

Question

lmerモデル（lme4パッケージ）の外れ値を特定したいと考えています。私はではないそれらを削除することに興味があります（LMERConvenienceFunctionsパッケージは何ですか？） - 私は単に外れ値が表示されていることを確認したいと思います。

私が使用したモデルの例：

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

Answer 1

あなたはLMERConvenienceFunctionsは（標準化されたレベル1の残差> 0から2.5SD）と同じ方法で異常値を定義する場合、これは動作します：

を

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

あなたはplants1で欠落データがある場合には、残差を抽出する前にlmerコマンドでna.action = 'na.exclude'を設定してください。これにより、データが欠落している観測値を省略するのではなく、NAに設定することができます。