def test_holm(self): """Test function holm. Compare to the p.adjust R function. """ reject, pval_corr = holm(pvals) assert_array_equal(reject, [False, False, True, False, False]) assert_array_equal(pval_corr, [5.2e-01, 3.6e-01, 5.0e-04, 1.2e-01, 3.6e-01]) _, pval_corr = holm(pvals2) assert_array_equal(pval_corr, [0.6, 0.5, 0.5, 0.6, 0.5]) # With NaN values _, pval_corr = holm(pvals2_NA) assert_array_almost_equal(pval_corr, [0.6, np.nan, 0.4, 0.6, 0.42]) # 2D array _, pval_corr = holm(pvals_2d) pval_corr = np.round(pval_corr.ravel(), 3) assert_array_almost_equal(pval_corr, [1, 0.72, 0.001, 0.24, 0.72, 1., np.nan, 0.7, 0.9, 0.72])
def test_holm(self): """Test function holm""" holm(pvals) holm(pvals, alpha=.01) holm(pvals, alpha=.90)