huber.py
residuals = np.abs(model.predict(x) - y)