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)
