Projection

If $$A$$ is skinny and full rank, then the projection operation is $$\mathcal P_{\mathcal R(A)}(y) = A(A^TA)^{-1}A^T y$$

Properties of projection matrix $$P$$

 * If $$P$$ is a projection matrix, then so is $$I-P$$.
 * TBC

TBC