where ψ (·) is partial derivative of ρ (·) and xij is jth entry in the ith row of matrix X with xi0 = 1.
Solution to these k equations is obtained by iterative reweighted least squares method (see Draper
and Smith [4]). At convergence, M-estimator may be given as