def test_binom_rejection_interval():
# consistency check with binom_test
# some code duplication but limit checks are different
alpha = 0.05
nobs = 200
prop = 12. / 20
alternative = 'smaller'
ci_low, ci_upp = smprop.binom_test_reject_interval(prop,
nobs,
alpha=alpha,
alternative=alternative)
assert_equal(ci_upp, nobs)
pval = smprop.binom_test(ci_low, nobs, prop=prop, alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_low + 1,
nobs,
prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
alternative = 'larger'
ci_low, ci_upp = smprop.binom_test_reject_interval(prop,
nobs,
alpha=alpha,
alternative=alternative)
assert_equal(ci_low, 0)
pval = smprop.binom_test(ci_upp, nobs, prop=prop, alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_upp - 1,
nobs,
prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
alternative = 'two-sided'
ci_low, ci_upp = smprop.binom_test_reject_interval(prop,
nobs,
alpha=alpha,
alternative=alternative)
pval = smprop.binom_test(ci_upp, nobs, prop=prop, alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_upp - 1,
nobs,
prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
pval = smprop.binom_test(ci_upp, nobs, prop=prop, alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_upp - 1,
nobs,
prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
def test_binom_rejection_interval():
# consistency check with binom_test
# some code duplication but limit checks are different
alpha = 0.05
nobs = 200
prop = 12./20
alternative='smaller'
ci_low, ci_upp = smprop.binom_test_reject_interval(prop, nobs, alpha=alpha,
alternative=alternative)
assert_equal(ci_upp, nobs)
pval = smprop.binom_test(ci_low, nobs, prop=prop,
alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_low + 1, nobs, prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
alternative='larger'
ci_low, ci_upp = smprop.binom_test_reject_interval(prop, nobs, alpha=alpha,
alternative=alternative)
assert_equal(ci_low, 0)
pval = smprop.binom_test(ci_upp, nobs, prop=prop,
alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_upp - 1, nobs, prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
alternative='two-sided'
ci_low, ci_upp = smprop.binom_test_reject_interval(prop, nobs, alpha=alpha,
alternative=alternative)
pval = smprop.binom_test(ci_upp, nobs, prop=prop,
alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_upp - 1, nobs, prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)
pval = smprop.binom_test(ci_upp, nobs, prop=prop,
alternative=alternative)
assert_array_less(pval, alpha)
pval = smprop.binom_test(ci_upp - 1, nobs, prop=prop,
alternative=alternative)
assert_array_less(alpha, pval)