Esempio n. 1
0
def create_prob2a_full():
    
    ubg = UndirectedBinaryGraph()
    
    #construct graph
    ubg.add_node(1)
    ubg.add_node(2)
    ubg.add_node(3)
    ubg.add_node(4)    
    
    #construct single node potentials    
    for n in ubg.g.nodes():
        config = {}
        rnum = np.random.rand()
        config[(0,)] = rnum
        config[(1,)] = 1 - rnum
        ftable = FunctionTable(config)    
        ubg.add_clique([n], ftable)    
    
    #construct and set potential functions
    #cliques = [ (1, 2, 3, 4)]
    cliques = [ (1, 2), (1, 3), (1, 4), (2, 3), (2, 4), (3, 4)]
    for c in cliques:
        bcombs = get_all_binary_combos(len(c))
        rnums = np.random.rand(len(bcombs))
        probs = rnums / rnums.sum()
        config = {}            
        for k,bc in enumerate(bcombs):            
            config[tuple(bc)] = probs[k]
        ftable = FunctionTable(config)
        ubg.add_clique(c, ftable)
    
    data = read_data_file('../IPF.dat')
        
    return (ubg, data)
Esempio n. 2
0
def create_simple_graph():
    
    ubg = UndirectedBinaryGraph()
    
    #construct graph
    ubg.add_node(1)
    
    #construct and set potential functions
    config = {}
    config[(0,)] = 0.8
    config[(1,)] = 0.2
    ftable = FunctionTable(config)    
    ubg.add_clique([1], ftable)
    
    #create data to go with it
    N = 100
    data = []    
    rnums = np.random.rand(N)
    for r in rnums:
        if r <= 0.30:
            data.append( [0] )
        else:
            data.append( [1] )
        
    return (ubg, data)