Python fold_intervalsの例

コード例 #1

ファイル: test_funcs.py プロジェクト: zkurtz/astropy

def test_histogram():
    with NumpyRNGContext(1234):
        a, b = 0.3, 3.14
        s = np.random.uniform(a, b, 10000) % 1

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

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

        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

ファイル: test_funcs.py プロジェクト: Cadair/astropy

def test_histogram():
    with NumpyRNGContext(1234):
        a, b = 0.3, 3.14
        s = np.random.uniform(a, b, 10000) % 1

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

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

        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