Exemplo n.º 1
0
def rejection_sample_demo():
    """Show estimation of probabilities by rejection sampling"""
    from Samplers import BNSampler

    sampler = BNSampler(asia)
    sampler.condition({'XRay':'abnormal'})
    print "Samples with 'None's are rejected" 
    for i in xrange(100):
        print sampler.rejection_sample()