¿Alguien sabe cómo calcular (o extraer) el apalancamiento y las distancias de Cook para un mer
objeto de clase (obtenido a través del lme4
paquete)? Me gustaría trazar estos para un análisis de residuos.
¿Alguien sabe cómo calcular (o extraer) el apalancamiento y las distancias de Cook para un mer
objeto de clase (obtenido a través del lme4
paquete)? Me gustaría trazar estos para un análisis de residuos.
Respuestas:
Usted debe tener una mirada en el paquete R influence.ME
. Le permite calcular medidas de datos influyentes para modelos de efectos mixtos generados por lme4
.
Un modelo de ejemplo:
library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)
La función influence
es la base de todos los pasos adicionales:
library(influence.ME)
infl <- influence(model, obs = TRUE)
Calcule la distancia de Cook:
cooks.distance(infl)
Trazar la distancia del cocinero:
plot(infl, which = "cook")
influence.ME
paquete. Desafortunadamente, no tengo una solución para esta tarea.
infl <- influence(model, group = "cyl")
porque especificaste el efecto aleatorio como (1|cyl)
? No sé, no entiendo en absoluto, que acaba de instalar influencia ... pero yo realmente no sé cuándo usar obs = TRUE
y cuándo usar group
...
cooksD_data<-as.data.frame(cooks.distance(ft1)) cooksD_data_select<-cooksd[cooksD_data>0.1,drop=FALSE,] cooksD_oultiers<-as.numeric(rownames(cooksD_data_select))]
hatvalues()
función recomendada aquí ?