huber_y1 = np.arange(-2.5,20,0.01)*huber.coef_[0] + 
            np.arange(-2.5,20,0.01)*huber.coef_[1] + 
            huber.intercept_