def test_hampel_sigma():
x = np.random.random(100)
# append a value slightly higher than the others
x_new = np.append(x, 10)
# v. high sigma value
out = stats.hampel(x_new, sigma=100)
ans = np.zeros(101).astype(np.int)
assert all(out == ans)
def test_hampel_negative():
x = np.random.random(100)
x_new = np.append(x, -100)
out = stats.hampel(x_new)
ans = np.zeros(100).astype(np.int)
ans = np.append(ans, -1.0)
assert len(out) == 101
assert sum(out) == -1.0
assert all(out == ans)