def test_kurtosis_basic(): """ circ.kurtosis: 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.kurtosis(data) assert_allclose(mp, 0.55260, 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_kurtosis_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.kurtosis(data, axis=0, mode='fisher') assert_allclose(mp, [-1.76462, -0.97441], rtol=1e-4)
def test_kurtosis_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.kurtosis(data, axis=1) assert_allclose(mp, [0.95562, 0.47166, 0.98155, 0.91244, 0.11831], rtol=1e-4)
def test_kurtosis_2d_data_axisNone_bootstrap_kw(): "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.kurtosis(alpha=data, ci=0.95) assert_allclose(mp, -0.013320, rtol=1e-4)
def test_kurtosis_2d_data_axisNone(): """circ.kurtosis: 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.kurtosis(data) assert_allclose(mp, -0.013320, rtol=1e-4)