¿Qué hacen los "efectos" de función en R?

17

No entiendo la explicación en Rel archivo de ayuda de efectos () :

Para un modelo lineal ajustado por lmo aov, los efectos son los valores no correlacionados de un solo grado de libertad obtenidos al proyectar los datos en los sucesivos subespacios ortogonales generados por la descomposición QR durante el proceso de ajuste.

¿Alguien puede explicar lo que esto significa?

¿Se aluden los subespacios ortogonales a los subespacios unidimensionales abarcados por las columnas de la parte Q de la descomposición QR (y, por lo tanto, ortogonales entre sí)? ¿O se supone que son ortogonales a otra cosa?

Vincent
fuente

Respuestas:

3

Dado el vector de respuesta , la matriz variable explicativa X y su descomposición QR X = Q R , los efectos devueltos por es el vector Q T y .yXX=QRRQTy

Aquí está el ejemplo numérico que confirma lo anterior:

> 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
fuente