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%