Prerequisites: Understanding partitioned least squares and its solution.
Setup
Given the linear model:
yy=Xbb+ee
We can partition the columns of X into X1 and X2 to get:
yy=X1bb1+X2bb2+ee
Partitioned Least Squares Result
Recall that bb2=(X~2′X~2)−1X~2′yy~ and that each variable with a tilde is the result of being partialled out by M1=I−X1(X1′X1)−1X1′.
Applying result to inclusion of a column of 1's
Consider the case where X1=ιι, then each column of X~2 is the form:
xx~2=(M1)xx2=(I−ιι(ιι′ιι)−1ιι′)xx2=M0xx=xx¨
We also have M0yy so yy~=yy¨ so all variables in the regression are effectively converted to deviations from their means in the regression. (See here for more on M0)
Which brings us to the result including a coefficient has on the solution to regression:
bb2=(X¨2′X¨2)−1X¨2′yy¨
Deriving the coefficient on ιι
Recall that bb1=(X1′X1)−1X1′(y−X2bb2).
Also, notice that bb1 is just a a scalar b1 (because X1 is just one column).
We have:
b1=(ιι′ιι)−1ιι′(yy−X2bb2)=(ιι′ιι)−1(ιι′yy−ιι′X2bb2)=(1/n)(nyˉ−nxxˉ′bb2)=yˉ−xxˉ′bb2