Extrahieren / Berechnen der Hebelwirkung und der Cook-Abstände für lineare Modelle mit gemischten Effekten

10

Weiß jemand, wie man Hebelwirkung und Cooks Abstände für ein merKlassenobjekt berechnet (oder extrahiert) (erhalten durch lme4Paket)? Ich möchte diese für eine Residuenanalyse zeichnen.

r mixed-model linear-model residuals leverage Roey Angel
quelle

Antworten:

13

Sie sollten sich das R-Paket ansehen influence.ME. Sie können damit einflussreiche Daten für Modelle mit gemischten Effekten berechnen, die von generiert wurden lme4.

Ein Beispielmodell:

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

Die Funktion influenceist die Basis für alle weiteren Schritte:

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

Berechnen Sie die Entfernung von Cook:

cooks.distance(infl)

Plot Cooks Entfernung:

plot(infl, which = "cook")

Geben Sie hier die Bildbeschreibung ein

Sven Hohenstein
quelle
Vielen Dank! Das hilft sicherlich. Wie wäre es mit der Berechnung der Hebelwirkung für den Abstand zwischen Cook und Distanz?
Roey Angel
@RoeyAngel Ich nehme an, dass dies mit dem influence.MEPaket nicht möglich ist . Leider habe ich keine Lösung für diese Aufgabe.
Sven Hohenstein
Sollte es nicht sein infl <- influence(model, group = "cyl"), weil Sie zufälligen Effekt als angegeben haben (1|cyl)? Ich weiß nicht, ich verstehe das überhaupt nicht, ich habe nur Einfluss installiert ... aber ich weiß nicht wirklich, wann obs = TRUEund wann ich es verwenden soll group...
Neugierig
Ich möchte Folgendes hinzufügen: Wenn Sie die Zeilennummer erhalten möchten, in der Cooks D-Abstände auftreten - dieselbe Zahl, die im Plot ohne Darstellung vorkommt, können Sie die folgende r-Formel für Cooks 'D-Abstandsnummern mit einem Schnitt verwenden Aus Wert von zB 0,1cooksD_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
Ist das besser als die hier empfohlenehatvalues() Funktion ?
Neugierig