The robust criterion function for Tukey’s biweight estimator
z (array_like) – 1d array
rho – rho(z) = -(1 - (z/c)**2)**3 * c**2/6. for |z| <= R
rho(z) = 0 for |z| > R
ndarray
TukeyBiweight.rho