Je souhaite identifier les valeurs aberrantes dans les modèles lmer (lme4). Je suis pas intéressé par les enlever (ce que fait le paquet LMERConvenienceFunctions) - Je veux simplement voir les valeurs aberrantes répertoriées.Identification des valeurs aberrantes dans les modèles lme4
Un exemple d'un modèle que j'utilise:
model1<-lmer(Value~ Moisture + Planting + (day|plot), data=plants1)
Si vous souhaitez définir des valeurs aberrantes de la même manière que LMERConvenienceFunctions fait (résidus de niveau normalisé 1> 2.5SD de 0), alors cela fonctionnera:
res1 <- resid(model1, type = "pearson") # Extract standardized residuals
plants1[which(abs(res1) > 2.5),] # Get the rows which absolute residuals > 2.5
si vous avez des données manquantes dans plants1 assurez-vous de mettre na.action = 'na.exclude' dans la commande lmer avant d'extraire les résidus. Cela garantira que les observations avec des données manquantes sont définies sur NA au lieu de les omettre.