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
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)
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)
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)
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)