Instead of just minimizing the vector u = y-Xβ, we