Пример #1
0
def test_histogram():
    a, b = 0.3, 3.14
    s = np.random.uniform(a, b, 10000) % 1

    b, w = kuiper.fold_intervals([(a, b, 1. / (b - a))])

    h = kuiper.histogram_intervals(16, b, w)
    nn, bb = np.histogram(s, bins=len(h), range=(0, 1), new=True)

    uu = np.sqrt(nn)
    nn, uu = len(h) * nn / h / len(s), len(h) * uu / h / len(s)

    c2 = np.sum(((nn - 1) / uu)**2)

    assert scipy.stats.chi2(len(h)).cdf(c2) > 0.01
    assert scipy.stats.chi2(len(h)).sf(c2) > 0.01
Пример #2
0
def test_histogram():
    a, b = 0.3, 3.14
    s = np.random.uniform(a,b,10000) % 1
    
    b, w = kuiper.fold_intervals([(a,b,1./(b-a))])

    h = kuiper.histogram_intervals(16,b,w)
    nn, bb = np.histogram(s, bins=len(h), range=(0,1), new=True)

    uu = np.sqrt(nn)
    nn, uu = len(h)*nn/h/len(s), len(h)*uu/h/len(s)

    c2 = np.sum(((nn-1)/uu)**2)
    
    assert scipy.stats.chi2(len(h)).cdf(c2)>0.01
    assert scipy.stats.chi2(len(h)).sf(c2)>0.01
Пример #3
0
def check_histogram_intervals_known(ii, rr):
    assert_array_almost_equal(kuiper.histogram_intervals(*ii), rr)
Пример #4
0
def check_histogram_intervals_known(ii, rr):
    assert_array_almost_equal(kuiper.histogram_intervals(*ii),rr)