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