def test_hdr_plot(close_figures):
fig = plt.figure()
ax = fig.add_subplot(111)
hdrboxplot(data, labels=labels.tolist(), ax=ax, threshold=1, seed=12345)
ax.set_xlabel("Month of the year")
ax.set_ylabel("Sea surface temperature (C)")
ax.set_xticks(np.arange(13, step=3) - 1)
ax.set_xticklabels(["", "Mar", "Jun", "Sep", "Dec"])
ax.set_xlim([-0.2, 11.2])
def test_hdr_plot(close_figures):
fig = plt.figure()
ax = fig.add_subplot(111)
try:
hdrboxplot(data, labels=labels.tolist(), ax=ax, threshold=1,
seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
ax.set_xlabel("Month of the year")
ax.set_ylabel("Sea surface temperature (C)")
ax.set_xticks(np.arange(13, step=3) - 1)
ax.set_xticklabels(["", "Mar", "Jun", "Sep", "Dec"])
ax.set_xlim([-0.2, 11.2])
def test_hdr_plot(close_figures):
fig = plt.figure()
ax = fig.add_subplot(111)
try:
hdrboxplot(data, labels=labels.tolist(), ax=ax, threshold=1,
seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
ax.set_xlabel("Month of the year")
ax.set_ylabel("Sea surface temperature (C)")
ax.set_xticks(np.arange(13, step=3) - 1)
ax.set_xticklabels(["", "Mar", "Jun", "Sep", "Dec"])
ax.set_xlim([-0.2, 11.2])
def test_hdr_bw(close_figures):
fig, hdr = hdrboxplot(data, bw='cv_ml')
median_t = [
24.25, 25.64, 25.99, 25.04, 23.71, 22.38, 21.31, 20.44, 20.24, 20.51,
21.19, 22.38
]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_alpha(close_figures):
_, hdr = hdrboxplot(data, alpha=[0.7], seed=12345)
extra_quant_t = np.vstack([[25.1, 26.5, 27.0, 26.4, 25.4, 24.1,
23.0, 22.0, 21.7, 22.1, 22.7, 23.8],
[23.4, 24.8, 25.0, 23.9, 22.4, 21.1,
20.0, 19.3, 19.2, 19.4, 20.1, 21.3]])
assert_almost_equal(hdr.extra_quantiles, extra_quant_t, decimal=0)
def test_hdr_basic(close_figures):
_, hdr = hdrboxplot(data, labels=labels, seed=12345)
assert len(hdr.extra_quantiles) == 0
median_t = [24.247, 25.625, 25.964, 24.999, 23.648, 22.302,
21.231, 20.366, 20.168, 20.434, 21.111, 22.299]
assert_almost_equal(hdr.median, median_t, decimal=2)
quant = np.vstack([hdr.outliers, hdr.hdr_90, hdr.hdr_50])
quant_t = np.vstack([[24.36, 25.42, 25.40, 24.96, 24.21, 23.35,
22.50, 21.89, 22.04, 22.88, 24.57, 25.89],
[27.25, 28.23, 28.85, 28.82, 28.37, 27.43,
25.73, 23.88, 22.26, 22.22, 22.21, 23.19],
[23.70, 26.08, 27.17, 26.74, 26.77, 26.15,
25.59, 24.95, 24.69, 24.64, 25.85, 27.08],
[28.12, 28.82, 29.24, 28.45, 27.36, 25.19,
23.61, 22.27, 21.31, 21.37, 21.60, 22.81],
[25.48, 26.99, 27.51, 27.04, 26.23, 24.94,
23.69, 22.72, 22.26, 22.64, 23.33, 24.44],
[23.11, 24.50, 24.66, 23.44, 21.74, 20.58,
19.68, 18.84, 18.76, 18.99, 19.66, 20.86],
[24.84, 26.23, 26.67, 25.93, 24.87, 23.57,
22.46, 21.45, 21.26, 21.57, 22.14, 23.41],
[23.62, 25.10, 25.34, 24.22, 22.74, 21.52,
20.40, 19.56, 19.63, 19.67, 20.37, 21.76]])
assert_almost_equal(quant, quant_t, decimal=0)
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape),
axis=1)
outliers = labels[labels_pos]
assert_equal([1982, 1983, 1997, 1998], outliers)
assert_equal(labels[hdr.outliers_idx], outliers)
def test_hdr_ncomp(close_figures):
fig, hdr = hdrboxplot(data, ncomp=3)
median_t = [
24.33, 25.71, 26.04, 25.08, 23.74, 22.40, 21.32, 20.45, 20.25, 20.53,
21.2, 22.39
]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_multiple_alpha(close_figures):
try:
_, hdr = hdrboxplot(data, alpha=[0.4, 0.92], seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
extra_quant_t = [[
25.712, 27.052, 27.711, 27.200, 26.162, 24.833, 23.639, 22.378, 22.250,
22.640, 23.472, 24.649
],
[
22.973, 24.526, 24.608, 23.343, 21.908, 20.655,
19.750, 19.046, 18.812, 18.989, 19.520, 20.685
],
[
24.667, 26.033, 26.416, 25.584, 24.308, 22.849,
21.684, 20.948, 20.483, 21.019, 21.751, 22.890
],
[
23.873, 25.371, 25.667, 24.644, 23.177, 21.923,
20.791, 20.015, 19.697, 19.951, 20.622, 21.858
]]
assert_almost_equal(hdr.extra_quantiles,
np.vstack(extra_quant_t),
decimal=0)
def test_hdr_threshold():
fig, hdr = hdrboxplot(data, alpha=[0.8], threshold=0.93)
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape),
axis=1)
outliers = labels[labels_pos]
assert_equal([1968, 1982, 1983, 1997, 1998], outliers)
plt.close(fig)
def test_hdr_alpha():
fig, hdr = hdrboxplot(data, alpha=[0.7])
extra_quant_t = np.vstack([[25.1, 26.5, 27.0, 26.4, 25.4, 24.1,
23.0, 22.0, 21.7, 22.1, 22.7, 23.8],
[23.4, 24.8, 25.0, 23.9, 22.4, 21.1,
20.0, 19.3, 19.2, 19.4, 20.1, 21.3]])
assert_almost_equal(hdr.extra_quantiles, extra_quant_t, decimal=0)
plt.close(fig)
def test_hdr_ncomp(close_figures):
try:
_, hdr = hdrboxplot(data, ncomp=3, seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
median_t = [24.33, 25.71, 26.04, 25.08, 23.74, 22.40,
21.32, 20.45, 20.25, 20.53, 21.20, 22.39]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_bw(close_figures):
try:
_, hdr = hdrboxplot(data, bw='cv_ml', seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
median_t = [24.25, 25.64, 25.99, 25.04, 23.71, 22.38,
21.31, 20.44, 20.24, 20.51, 21.19, 22.38]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_basic_brute(close_figures):
_, hdr = hdrboxplot(data, labels=labels, use_brute=True)
assert len(hdr.extra_quantiles) == 0
median_t = [24.247, 25.625, 25.964, 24.999, 23.648, 22.302,
21.231, 20.366, 20.168, 20.434, 21.111, 22.299]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_bw(close_figures):
try:
_, hdr = hdrboxplot(data, bw='cv_ml', seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
median_t = [24.25, 25.64, 25.99, 25.04, 23.71, 22.38,
21.31, 20.44, 20.24, 20.51, 21.19, 22.38]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_ncomp(close_figures):
try:
_, hdr = hdrboxplot(data, ncomp=3, seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
median_t = [24.33, 25.71, 26.04, 25.08, 23.74, 22.40,
21.32, 20.45, 20.25, 20.53, 21.20, 22.39]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_threshold(close_figures):
try:
_, hdr = hdrboxplot(data, alpha=[0.8], threshold=0.93, seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape),
axis=1)
outliers = labels[labels_pos]
assert_equal([1968, 1982, 1983, 1997, 1998], outliers)
def test_hdr_alpha(close_figures):
try:
_, hdr = hdrboxplot(data, alpha=[0.7], seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
extra_quant_t = np.vstack([[25.1, 26.5, 27.0, 26.4, 25.4, 24.1,
23.0, 22.0, 21.7, 22.1, 22.7, 23.8],
[23.4, 24.8, 25.0, 23.9, 22.4, 21.1,
20.0, 19.3, 19.2, 19.4, 20.1, 21.3]])
assert_almost_equal(hdr.extra_quantiles, extra_quant_t, decimal=0)
def test_hdr_threshold(close_figures):
try:
_, hdr = hdrboxplot(data, alpha=[0.8], threshold=0.93,
seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape),
axis=1)
outliers = labels[labels_pos]
assert_equal([1968, 1982, 1983, 1997, 1998], outliers)
def test_hdr_alpha(close_figures):
try:
_, hdr = hdrboxplot(data, alpha=[0.7], seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
extra_quant_t = np.vstack([[25.1, 26.5, 27.0, 26.4, 25.4, 24.1,
23.0, 22.0, 21.7, 22.1, 22.7, 23.8],
[23.4, 24.8, 25.0, 23.9, 22.4, 21.1,
20.0, 19.3, 19.2, 19.4, 20.1, 21.3]])
assert_almost_equal(hdr.extra_quantiles, extra_quant_t, decimal=0)
def test_hdr_basic_brute(close_figures, reset_randomstate):
try:
_, hdr = hdrboxplot(data, ncomp=2, labels=labels, use_brute=True)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
assert len(hdr.extra_quantiles) == 0
median_t = [24.247, 25.625, 25.964, 24.999, 23.648, 22.302,
21.231, 20.366, 20.168, 20.434, 21.111, 22.299]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_plot():
fig = plt.figure()
ax = fig.add_subplot(111)
fig, res = hdrboxplot(data, labels=labels.tolist(), ax=ax, threshold=1)
ax.set_xlabel("Month of the year")
ax.set_ylabel("Sea surface temperature (C)")
ax.set_xticks(np.arange(13, step=3) - 1)
ax.set_xticklabels(["", "Mar", "Jun", "Sep", "Dec"])
ax.set_xlim([-0.2, 11.2])
plt.close(fig)
def test_hdr_basic_brute(close_figures, reset_randomstate):
try:
_, hdr = hdrboxplot(data, ncomp=2, labels=labels, use_brute=True)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
assert len(hdr.extra_quantiles) == 0
median_t = [24.247, 25.625, 25.964, 24.999, 23.648, 22.302,
21.231, 20.366, 20.168, 20.434, 21.111, 22.299]
assert_almost_equal(hdr.median, median_t, decimal=2)
def test_hdr_multiple_alpha():
fig, hdr = hdrboxplot(data, alpha=[0.4, 0.92])
extra_quant_t = [[25.712, 27.052, 27.711, 27.200,
26.162, 24.833, 23.639, 22.378,
22.250, 22.640, 23.472, 24.649],
[22.973, 24.526, 24.608, 23.343,
21.908, 20.655, 19.750, 19.046,
18.812, 18.989, 19.520, 20.685],
[24.667, 26.033, 26.416, 25.584,
24.308, 22.849, 21.684, 20.948,
20.483, 21.019, 21.751, 22.890],
[23.873, 25.371, 25.667, 24.644,
23.177, 21.923, 20.791, 20.015,
19.697, 19.951, 20.622, 21.858]]
assert_almost_equal(hdr.extra_quantiles, np.vstack(extra_quant_t), decimal=0)
plt.close(fig)
def test_hdr_multiple_alpha():
fig, hdr = hdrboxplot(data, alpha=[0.4, 0.92])
extra_quant_t = [[25.712, 27.052, 27.711, 27.200,
26.162, 24.833, 23.639, 22.378,
22.250, 22.640, 23.472, 24.649],
[22.973, 24.526, 24.608, 23.343,
21.908, 20.655, 19.750, 19.046,
18.812, 18.989, 19.520, 20.685],
[24.667, 26.033, 26.416, 25.584,
24.308, 22.849, 21.684, 20.948,
20.483, 21.019, 21.751, 22.890],
[23.873, 25.371, 25.667, 24.644,
23.177, 21.923, 20.791, 20.015,
19.697, 19.951, 20.622, 21.858]]
assert_almost_equal(hdr.extra_quantiles, np.vstack(extra_quant_t), decimal=0)
plt.close(fig)
def test_hdr_multiple_alpha(close_figures):
try:
_, hdr = hdrboxplot(data, alpha=[0.4, 0.92], seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
extra_quant_t = [[25.712, 27.052, 27.711, 27.200,
26.162, 24.833, 23.639, 22.378,
22.250, 22.640, 23.472, 24.649],
[22.973, 24.526, 24.608, 23.343,
21.908, 20.655, 19.750, 19.046,
18.812, 18.989, 19.520, 20.685],
[24.667, 26.033, 26.416, 25.584,
24.308, 22.849, 21.684, 20.948,
20.483, 21.019, 21.751, 22.890],
[23.873, 25.371, 25.667, 24.644,
23.177, 21.923, 20.791, 20.015,
19.697, 19.951, 20.622, 21.858]]
assert_almost_equal(hdr.extra_quantiles, np.vstack(extra_quant_t),
decimal=0)
def test_hdr_basic(close_figures):
try:
_, hdr = hdrboxplot(data, labels=labels, seed=12345)
except WindowsError:
pytest.xfail('Multiprocess randomly crashes in Windows testing')
assert len(hdr.extra_quantiles) == 0
median_t = [24.247, 25.625, 25.964, 24.999, 23.648, 22.302,
21.231, 20.366, 20.168, 20.434, 21.111, 22.299]
assert_almost_equal(hdr.median, median_t, decimal=2)
quant = np.vstack([hdr.outliers, hdr.hdr_90, hdr.hdr_50])
quant_t = np.vstack([[24.36, 25.42, 25.40, 24.96, 24.21, 23.35,
22.50, 21.89, 22.04, 22.88, 24.57, 25.89],
[27.25, 28.23, 28.85, 28.82, 28.37, 27.43,
25.73, 23.88, 22.26, 22.22, 22.21, 23.19],
[23.70, 26.08, 27.17, 26.74, 26.77, 26.15,
25.59, 24.95, 24.69, 24.64, 25.85, 27.08],
[28.12, 28.82, 29.24, 28.45, 27.36, 25.19,
23.61, 22.27, 21.31, 21.37, 21.60, 22.81],
[25.48, 26.99, 27.51, 27.04, 26.23, 24.94,
23.69, 22.72, 22.26, 22.64, 23.33, 24.44],
[23.11, 24.50, 24.66, 23.44, 21.74, 20.58,
19.68, 18.84, 18.76, 18.99, 19.66, 20.86],
[24.84, 26.23, 26.67, 25.93, 24.87, 23.57,
22.46, 21.45, 21.26, 21.57, 22.14, 23.41],
[23.62, 25.10, 25.34, 24.22, 22.74, 21.52,
20.40, 19.56, 19.63, 19.67, 20.37, 21.76]])
assert_almost_equal(quant, quant_t, decimal=0)
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape),
axis=1)
outliers = labels[labels_pos]
assert_equal([1982, 1983, 1997, 1998], outliers)
assert_equal(labels[hdr.outliers_idx], outliers)
def test_hdr_threshold():
fig, hdr = hdrboxplot(data, alpha=[0.8], threshold=0.93)
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape), axis=1)
outliers = labels[labels_pos]
assert_equal([1968, 1982, 1983, 1997, 1998], outliers)
plt.close(fig)
def test_hdr_threshold(close_figures):
_, hdr = hdrboxplot(data, alpha=[0.8], threshold=0.93, seed=12345)
labels_pos = np.all(np.in1d(data, hdr.outliers).reshape(data.shape),
axis=1)
outliers = labels[labels_pos]
assert_equal([1968, 1982, 1983, 1997, 1998], outliers)
def test_hdr_ncomp():
fig, hdr = hdrboxplot(data, ncomp=3)
median_t = [24.33, 25.71, 26.04, 25.08, 23.74, 22.40,
21.32, 20.45, 20.25, 20.53, 21.2 , 22.39]
assert_almost_equal(hdr.median, median_t, decimal=2)
plt.close(fig)
def test_hdr_bw():
fig, hdr = hdrboxplot(data, bw='cv_ml')
median_t = [24.25, 25.64, 25.99, 25.04, 23.71, 22.38,
21.31, 20.44, 20.24, 20.51, 21.19, 22.38]
assert_almost_equal(hdr.median, median_t, decimal=2)
plt.close(fig)