Gradient descent involves significant computational effort, whereas the method of least squares enables direct and accurate calculation. Does gradient descent offer any advantages over least squares estimation?
The only potential advantage of gradient descent that I can identify is that it allows us to use other loss functions. However, I'm not sure when and why we would need to use a different loss function, since by Gauss–Markov theorem the least squares estimates should be optimal.