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.
Extrahieren / Berechnen der Hebelwirkung und der Cook-Abstände für lineare Modelle mit gemischten Effekten
Antworten:
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")
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
—
Sven Hohenstein
influence.MEPaket nicht möglich ist . Leider habe ich keine Lösung für diese Aufgabe.
Sollte es nicht sein
—
Neugierig
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...
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,1
—
Elias Estatistics
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))]