def test_median_ci():
alpha = np.ones((2, 10))
m1 = np.ones(2)
m0 = np.ones(10)
mout1, ci_1 = pycircstat.median(alpha, axis=1, ci=.8)
mout0, ci_0 = pycircstat.median(alpha, axis=0, ci=.8)
moutall, ci_all = pycircstat.median(alpha, axis=0, ci=.8)
assert_allclose(mout1, m1)
assert_allclose(mout0, m0)
assert_allclose(moutall, 1.)
assert_allclose(ci_0.lower, m0)
assert_allclose(ci_0.upper, m0)
assert_allclose(ci_1.lower, m1)
assert_allclose(ci_1.upper, m1)
assert_allclose(ci_all.lower, 1.)
assert_allclose(ci_all.upper, 1.)
def test_median_ci():
alpha = np.ones((2, 10))
m1 = np.ones(2)
m0 = np.ones(10)
mout1, ci_1 = pycircstat.median(alpha, axis=1, ci=.8)
mout0, ci_0 = pycircstat.median(alpha, axis=0, ci=.8)
moutall, ci_all = pycircstat.median(alpha, axis=0, ci=.8)
assert_allclose(mout1, m1)
assert_allclose(mout0, m0)
assert_allclose(moutall, 1.)
assert_allclose(ci_0.lower, m0)
assert_allclose(ci_0.upper, m0)
assert_allclose(ci_1.lower, m1)
assert_allclose(ci_1.upper, m1)
assert_allclose(ci_all.lower, 1.)
assert_allclose(ci_all.upper, 1.)
def test_median():
alpha = np.array([
[3.73153000, 1.63904879, 4.03175622, 3.90422402, 4.61029613,
4.04117818, 5.79313473, 5.50863002, 5.81530225, 2.44973903],
[2.12868554, 0.09073566, 0.05581025, 5.10673712, 1.68712454,
3.72915575, 4.45439608, 4.70694685, 3.58470730, 2.49742028]
])
m0 = np.array([2.93010777, 0.86489223, -1.09780942, -1.77770474,
-3.13447497, -2.39801834, -1.15941990, -1.17539688,
-1.58318053, 2.47357966]) % (2 * np.pi)
m1 = np.array([-2.24671810, -1.24910966]) % (2 * np.pi)
m11 = np.array([-2.24200713, -1.82878923]) % (2 * np.pi)
mall = -2.2467 % (2 * np.pi)
assert_allclose(pycircstat.median(alpha, axis=1), m1)
assert_allclose(pycircstat.median(alpha[:, :-1], axis=1), m11)
assert_allclose(pycircstat.median(alpha, axis=0), m0)
assert_allclose(pycircstat.median(alpha), mall, atol=1e-4)
def test_median():
alpha = np.array([
[3.73153000, 1.63904879, 4.03175622, 3.90422402, 4.61029613,
4.04117818, 5.79313473, 5.50863002, 5.81530225, 2.44973903],
[2.12868554, 0.09073566, 0.05581025, 5.10673712, 1.68712454,
3.72915575, 4.45439608, 4.70694685, 3.58470730, 2.49742028]
])
m0 = np.array([2.93010777, 0.86489223, -1.09780942, -1.77770474,
-3.13447497, -2.39801834, -1.15941990, -1.17539688,
-1.58318053, 2.47357966]) % (2 * np.pi)
m1 = np.array([-2.24671810, -1.24910966]) % (2 * np.pi)
m11 = np.array([-2.24200713, -1.82878923]) % (2 * np.pi)
mall = -2.2467 % (2 * np.pi)
assert_allclose(pycircstat.median(alpha, axis=1), m1)
assert_allclose(pycircstat.median(alpha[:, :-1], axis=1), m11)
assert_allclose(pycircstat.median(alpha, axis=0), m0)
assert_allclose(pycircstat.median(alpha), mall, atol=1e-4)