Was bewirkt die Funktion „Effekte“ in R?

17

Ich verstehe die Erklärung in Rder Hilfedatei für effects () nicht :

Für ein lineares Modell, das mit lmoder angepasst wurde aov, sind die Effekte die nicht korrelierten Einzelfreiheitsgradwerte, die durch Projektion der Daten auf die aufeinanderfolgenden orthogonalen Teilräume erhalten wurden, die durch die QR-Zerlegung während des Anpassungsprozesses generiert wurden.

Kann jemand erklären, was das bedeutet?

Sind die orthogonalen Teilräume auf die eindimensionalen Teilräume bezogen, die von den Spalten des Q-Teils der QR-Zerlegung aufgespannt werden (und daher orthogonal zueinander sind)? Oder sollen sie orthogonal zu etwas anderem sein?

Vincent
quelle

Antworten:

3

yXX=Q.RRQ.Ty

Hier ist das numerische Beispiel, das das oben Gesagte bestätigt:

> set.seed(1001)
> x<-rnorm(100)
> y<-1+2*x+rnorm(100)
> mod<-lm(y~x)
> xqr<-qr(cbind(1,x))
> sum(abs(qr.qty(xqr,y)-effects(mod)))
[1] 0
mpiktas
quelle