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.

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


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
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.