def test_lladser_ci_series_exact(self):
"""lladser_ci_series returns series of predictions"""
#Values are from Manuel's email of 9/11/09
#have seen RWB
urn_1 = 'RWBWWBWRRWRYWRPPZ'
results = list(lladser_ci_series(urn_1, r=4))
self.assertEqual(len(results), 3)
def test_lladser_ci_series_exact(self):
"""lladser_ci_series returns series of predictions"""
#Values are from Manuel's email of 9/11/09
#have seen RWB
urn_1 = 'RWBWWBWRRWRYWRPPZ'
results = list(lladser_ci_series(urn_1, r=4))
self.assertEqual(len(results), 3)
def test_lladser_ci_series_random(self):
"""lladser_ci_series' interval contain true prob with expected alpha."""
seq="WBWBWBWBWBWB"
observations=[]
alpha=0.95
reps = 1000
for i in range(reps):
obs = list(lladser_ci_series(seq, r=4, alpha=alpha))[0]
observations.append(obs)
tps = filter (lambda (a,b): a < 0.5 and 0.5 < b, observations)
self.assertTrue(len(tps) >= alpha*reps ) #100%-95%
def test_lladser_ci_series_random(self):
"""lladser_ci_series' interval contain true prob with expected alpha."""
seq="WBWBWBWBWBWB"
observations=[]
alpha=0.95
reps = 1000
for i in range(reps):
obs = list(lladser_ci_series(seq, r=4, alpha=alpha))[0]
observations.append(obs)
tps = filter (lambda (a,b): a < 0.5 and 0.5 < b, observations)
self.assertTrue(len(tps) >= alpha*reps ) #100%-95%
