Come estrarre / calcolare la leva e le distanze di Cook per modelli lineari di effetti misti

Qualcuno sa come calcolare (o estrarre) la leva e le distanze di Cook per un meroggetto di classe (ottenuto attraverso il lme4pacchetto)? Vorrei tracciare questi per un'analisi dei residui.

— Roey Angel
fonte

Si dovrebbe dare un'occhiata al pacchetto R influence.ME. Consente di calcolare misure di dati influenti per modelli di effetti misti generati da lme4.

Un modello di esempio:

library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)

La funzione influenceè la base per tutti gli ulteriori passaggi:

library(influence.ME)
infl <- influence(model, obs = TRUE)

Calcola la distanza di Cook:

cooks.distance(infl)

Traccia la distanza di Cook:

plot(infl, which = "cook")

inserisci qui la descrizione dell'immagine

— Sven Hohenstein
fonte

Grazie! Questo sicuramente aiuta. Che ne dici di calcolare la leva per la distanza di un cuoco rispetto alla trama della leva?

— Roey Angel,

@RoeyAngel Suppongo che questo non sia possibile con il influence.MEpacchetto. Sfortunatamente, non ho una soluzione per questo compito.

— Sven Hohenstein,

Non dovrebbe essere infl <- influence(model, group = "cyl"), perché hai specificato l'effetto casuale come (1|cyl)? Non lo so, non lo capisco affatto, ho appena installato influenza ... ma non so davvero quando usare obs = TRUEe quando usare group...

— Curioso,

Vorrei aggiungere quanto segue: Se si desidera ottenere il numero di riga in cui si verificano le distanze D di Cook - lo stesso numero presente nella trama senza tracciare, è possibile utilizzare la seguente formula r sui numeri delle distanze D di Cook con un taglio valore fuori di esempio 0.1

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))]

— Elias Estatistics

È meglio della hatvalues()funzione consigliata qui ?

— Curioso