示例#1
0
def test_skewness_basic():
    """ circ.skewness: test basic call... """
    data = np.array([1.80044838, 2.02938314, 1.03534016, 4.84225057,
                     1.54256458, 5.19290675, 2.18474784,
                     4.77054777, 1.51736933, 0.72727580])
    mp = pycircstat.skewness(data)
    assert_allclose(mp, 0.40660, rtol=1e-5)
def test_skewness_basic():
    """ circ.skewness: test basic call... """
    data = np.array([1.80044838, 2.02938314, 1.03534016, 4.84225057,
                     1.54256458, 5.19290675, 2.18474784,
                     4.77054777, 1.51736933, 0.72727580])
    mp = pycircstat.skewness(data)
    assert_allclose(mp, 0.40660, rtol=1e-5)
def circvVals(rad, w, d):
    vStr = pcirc.vector_strength(rad, w=w, d=d/2, ci=None)
    vDir = pcirc.mean(rad, w=w, d=d/2, ci=None)
    astd = pcirc.astd(rad, w=w, d=d/2, ci=None)
    skew = pcirc.skewness(rad, w=w, ci=None)
    kurt = pcirc.kurtosis(rad, w=w, ci=None)
    return vStr, vDir, astd, skew, kurt
示例#4
0
def test_skewness_2d_data_axis0_fisher():
    data = np.array([
                    [0.58429, 0.88333],
                    [1.14892, 2.22854],
                    [2.87128, 3.06369],
                    [1.07677, 1.49836],
                    [2.96969, 1.51748],
                    ])
    mp = pycircstat.skewness(data, axis=0, mode='fisher')
    assert_allclose(mp, [-0.999794, -1.157187], rtol=1e-4)
def test_skewness_2d_data_axis0_fisher():
    data = np.array([
        [0.58429, 0.88333],
        [1.14892, 2.22854],
        [2.87128, 3.06369],
        [1.07677, 1.49836],
        [2.96969, 1.51748],
    ])
    mp = pycircstat.skewness(data, axis=0, mode='fisher')
    assert_allclose(mp, [-0.999794, -1.157187], rtol=1e-4)
def test_skewness_2d_data_axis1():
    data = np.array([
                    [0.58429, 0.88333],
                    [1.14892, 2.22854],
                    [2.87128, 3.06369],
                    [1.07677, 1.49836],
                    [2.96969, 1.51748],
                    ])
    mp = pycircstat.skewness(data, axis=1)
    assert_allclose(mp,  [-8.3267e-17, 0.0000e+00, -1.3878e-17, 1.3878e-16, 0.0000e+00], atol=1e-4)
示例#7
0
def test_skewness_2d_data_axisNone_bootstrap():
    "basically only test whether boostrapping does not throw an error"
    data = np.array([
                    [0.58429, 0.88333],
                    [1.14892, 2.22854],
                    [2.87128, 3.06369],
                    [1.07677, 1.49836],
                    [2.96969, 1.51748],
                    ])
    mp, (low, high) = pycircstat.skewness(data, ci=0.95)
    assert_allclose(mp, -0.18157, rtol=1e-4)
示例#8
0
def test_skewness_2d_data_axisNone():
    """circ.skewness: test 2D data (axis=0)..."""
    data = np.array([
                    [0.58429, 0.88333],
                    [1.14892, 2.22854],
                    [2.87128, 3.06369],
                    [1.07677, 1.49836],
                    [2.96969, 1.51748],
                    ])
    mp = pycircstat.skewness(data)
    assert_allclose(mp, -0.18157, rtol=1e-4)
def test_skewness_2d_data_axisNone_bootstrap():
    "basically only test whether boostrapping does not throw an error"
    data = np.array([
        [0.58429, 0.88333],
        [1.14892, 2.22854],
        [2.87128, 3.06369],
        [1.07677, 1.49836],
        [2.96969, 1.51748],
    ])
    mp, (low, high) = pycircstat.skewness(data, ci=0.95)
    assert_allclose(mp, -0.18157, rtol=1e-4)
def test_skewness_2d_data_axisNone():
    """circ.skewness: test 2D data (axis=0)..."""
    data = np.array([
        [0.58429, 0.88333],
        [1.14892, 2.22854],
        [2.87128, 3.06369],
        [1.07677, 1.49836],
        [2.96969, 1.51748],
    ])
    mp = pycircstat.skewness(data)
    assert_allclose(mp, -0.18157, rtol=1e-4)
示例#11
0
def test_skewness_2d_data_axis1():
    data = np.array([
                    [0.58429, 0.88333],
                    [1.14892, 2.22854],
                    [2.87128, 3.06369],
                    [1.07677, 1.49836],
                    [2.96969, 1.51748],
                    ])
    mp = pycircstat.skewness(data, axis=1)
    assert_allclose(
        mp, [-8.3267e-17, 0.0000e+00, -1.3878e-17, 1.3878e-16, 0.0000e+00], atol=1e-4)