def test_huber_huberT(self): n = scale.norms.HuberT() n.t = 1.5 h = scale.Huber(norm=n) assert_almost_equal(scale.huber(self.chem)[0], h(self.chem)[0], DECIMAL) assert_almost_equal(scale.huber(self.chem)[1], h(self.chem)[1], DECIMAL)
def test_huber_Hampel(self): hh = scale.Huber(norm=scale.norms.Hampel()) assert_almost_equal(hh(self.chem)[0], 3.17434, DECIMAL) assert_almost_equal(hh(self.chem)[1], 0.66782, DECIMAL)
def setup_class(cls): np.random.seed(54321) cls.X = standard_normal((40, 10, 30)) cls.h = scale.Huber(maxiter=1000, tol=1.0e-05)
def test_huber_result_shape(self): h = scale.Huber(maxiter=100) m, s = h(self.X) assert_equal(m.shape, (10, ))
def basic_functionality(self): h = scale.Huber(maxiter=100) m, s = h(self.X) assert_equal(m.shape, (10,))
def __init__(self): np.random.seed(54321) self.X = standard_normal((40, 10, 30)) self.h = scale.Huber(maxiter=1000, tol=1.0e-05)