The psi function for Huber’s t estimator
The analytic derivative of rho
1d array
ndarray
psi(z) = z for |z| <= t
psi(z) = sign(z)*t for |z| > t
HuberT.psi