To increase the speed of convergence.

To handle cases of complete separation where MLE does not exist.

To account for non-normality of the residuals.

To adjust for heteroscedasticity.