def test_compare_zm2():
data = np.random.random(1000)
H, M, fpp = htest.h_test(data)
Z, zm2fpp = zm2.zm2(data, M)
assert np.abs(Z - (H+4*M-4))<1e-8
assert fpp>zm2fpp
def test_null():
for N in 100, 1000, 10000:
for m in [1,2,5]:
yield check_uniform, lambda x: zm2.zm2(x,m)[1], N
def test_non_null():
Z, fpp = zm2.zm2(np.random.random(1000)/2,5)
assert fpp<0.01
def test_fpp():
for N in 10, 100, 500:
for m in [1,2,5]:
yield check_fpp, lambda x: zm2.zm2(x,m)[1], N, 1000, 0.05
def test_compare_zm2():
data = np.random.random(1000)
H, M, fpp = htest.h_test(data)
Z, zm2fpp = zm2.zm2(data, M)
assert np.abs(Z - (H + 4 * M - 4)) < 1e-8
assert fpp > zm2fpp
