[R] residuals from lm
CharlieTheBrown77 at gmail.com
Mon Jul 2 20:21:56 CEST 2012
I forgot the column of 1's for the intercept term
> I was playing around with something else and I noticed this matrix code
> for residuals in a linear model doesn't say what lm() says. Please tell
> me if I am completely misguided here.
> Y <- as.matrix(mtcars[,1])
> X <- as.matrix(mtcars[,c(2:11)])
> # shouldnt this:
> H <- X %*% solve(t(X) %*% X) %*% t(X)
> (diag(dim(H)) - H) %*% Y
> # be equal to this:
> # ???
> # thanks
View this message in context: http://r.789695.n4.nabble.com/residuals-from-lm-tp4635155p4635183.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help