I'm trying to show that $\hat{e}$ and $\hat{\beta}$, the residuals and coefficient vector for the OLS problem $y=X\beta+\epsilon$, are uncorrelated.
The proof I see online doesn't make sense to me because it doesn't subtract $E[\hat{e}]E[\hat{\beta}]$ when computing the covariance.
From page 21 of http://homepage.ntu.edu.tw/~ckuan/pdf/et01/et_Ch3.pdf it says $$ cov(\hat{e},\hat{\beta}) = E[(I-P)yy'X(X'X)^{-1}] \\ = (I-P)E[yy']X(X'X)^{-1} = \sigma_0^2(I-P)X(X'X)^{-1} = 0. $$ where $P=X(X'X)^{-1}X'$ is the orthogonal projection onto $span(X)$. This seems wrong to me because $cov(A,B)=E[AB^{T}]-E[A]E[B^T]$, but above is just $E[\hat{e}\hat{\beta}^T].$