def test_lladser_ci_from_r():
f = 10
t = 10
r = 4
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f)
assert_almost_equal(obs_low, 0.0806026244)
assert_almost_equal(obs_high, 0.806026244)
r = 20
t = 100
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f)
assert_almost_equal(obs_low, 0.02787923964)
assert_almost_equal(obs_high, 0.2787923964)
# make sure we test with each possible alpha
alpha = 0.99
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha)
assert_almost_equal(obs_low, 0.03184536992)
assert_almost_equal(obs_high, 0.3184536992)
alpha = 0.9
r = 3
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha)
assert_almost_equal(obs_low, 0.005635941995)
assert_almost_equal(obs_high, 0.05635941995)
# test other ci_types
ci_type = 'ULCU'
obs_low, obs_high = _lladser_ci_from_r(r=r,
t=t,
f=f,
alpha=alpha,
ci_type=ci_type)
assert_almost_equal(obs_low, 0.01095834700)
assert_almost_equal(obs_high, 0.1095834700)
alpha = 0.95
t = 10
ci_type = 'U'
obs_low, obs_high = _lladser_ci_from_r(r=r,
t=t,
f=f,
alpha=alpha,
ci_type=ci_type)
assert_almost_equal(obs_low, 0)
assert_almost_equal(obs_high, 0.6295793622)
ci_type = 'L'
obs_low, obs_high = _lladser_ci_from_r(r=r,
t=t,
f=f,
alpha=alpha,
ci_type=ci_type)
assert_almost_equal(obs_low, 0.0817691447)
assert_almost_equal(obs_high, 1)
def test_lladser_ci_from_r():
f = 10
t = 10
r = 4
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f)
assert_almost_equal(obs_low, 0.0806026244)
assert_almost_equal(obs_high, 0.806026244)
r = 20
t = 100
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f)
assert_almost_equal(obs_low, 0.02787923964)
assert_almost_equal(obs_high, 0.2787923964)
# make sure we test with each possible alpha
alpha = 0.99
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha)
assert_almost_equal(obs_low, 0.03184536992)
assert_almost_equal(obs_high, 0.3184536992)
alpha = 0.9
r = 3
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha)
assert_almost_equal(obs_low, 0.005635941995)
assert_almost_equal(obs_high, 0.05635941995)
# test other ci_types
ci_type = 'ULCU'
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha,
ci_type=ci_type)
assert_almost_equal(obs_low, 0.01095834700)
assert_almost_equal(obs_high, 0.1095834700)
alpha = 0.95
t = 10
ci_type = 'U'
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha,
ci_type=ci_type)
assert_almost_equal(obs_low, 0)
assert_almost_equal(obs_high, 0.6295793622)
ci_type = 'L'
obs_low, obs_high = _lladser_ci_from_r(r=r, t=t, f=f, alpha=alpha,
ci_type=ci_type)
assert_almost_equal(obs_low, 0.0817691447)
assert_almost_equal(obs_high, 1)
def test_lladser_ci_from_r_invalid_input():
# unsupported alpha for ci_type='U'
with assert_raises(ValueError):
_lladser_ci_from_r(r=3, t=10, f=10, alpha=0.90, ci_type='U')
# unsupported r for ci_type='U'
with assert_raises(ValueError):
_lladser_ci_from_r(r=42, t=10, f=10, alpha=0.95, ci_type='U')
# unsupported alpha for ci_type='L'
with assert_raises(ValueError):
_lladser_ci_from_r(r=3, t=10, f=10, alpha=0.90, ci_type='L')
# unsupported r for ci_type='L'
with assert_raises(ValueError):
_lladser_ci_from_r(r=50, t=10, f=10, alpha=0.95, ci_type='L')
# unknown ci_type
with assert_raises(ValueError):
_lladser_ci_from_r(r=4, t=10, f=10, alpha=0.95, ci_type='brofist')
# requesting CI for not precomputed values
with assert_raises(ValueError):
_lladser_ci_from_r(r=500, t=10, f=10)
def test_lladser_ci_from_r_invalid_input():
# unsupported alpha for ci_type='U'
with assert_raises(ValueError):
_lladser_ci_from_r(r=3, t=10, f=10, alpha=0.90, ci_type='U')
# unsupported r for ci_type='U'
with assert_raises(ValueError):
_lladser_ci_from_r(r=42, t=10, f=10, alpha=0.95, ci_type='U')
# unsupported alpha for ci_type='L'
with assert_raises(ValueError):
_lladser_ci_from_r(r=3, t=10, f=10, alpha=0.90, ci_type='L')
# unsupported r for ci_type='L'
with assert_raises(ValueError):
_lladser_ci_from_r(r=50, t=10, f=10, alpha=0.95, ci_type='L')
# unknown ci_type
with assert_raises(ValueError):
_lladser_ci_from_r(r=4, t=10, f=10, alpha=0.95, ci_type='brofist')
# requesting CI for not precomputed values
with assert_raises(ValueError):
_lladser_ci_from_r(r=500, t=10, f=10)
