powered by
Huber's M-estimator
huber.M(data, q = pchisq(3, 2))
q=pchisq(k+1,k) gives maximum breakdown point