Hi.

I know that R computes sums of squares based on the diagonal of

t(Q) %*% y %*% t(y) %*% Q,

where Q comes from the QR-decomposition of the model matrix.

Does anyone know where I can find a proof for this result?

All Best and Happy New Year,
Ethan

	[[alternative HTML version deleted]]

