Esempio n. 1
0
def check_prior_recovery(net_fn, priors_fn):
    (net, genes, tfs) = load_network(net_fn)
    (p, sign) = load_priors(priors_fn, '1')
    pcon = fl.priors_to_constraints(['1'], [genes],[tfs],p,0.5)
    prior_inter = []
    for con in pcon:
        prior_inter.append(net[con.c1.r, con.c1.c])
    
    return (np.array(prior_inter), np.array(sign))
Esempio n. 2
0
def check_prior_corr(expr_fn, tfs_fn, priors_fn):
    (e, t, genes, tfs) = load_bacteria(expr_fn, tfs_fn)
    (p, s) = load_priors(priors_fn,'1')
    pcon = fl.priors_to_constraints(['1'], [genes],[tfs],p,0.5)

    prior_corr = []
    for con in pcon:
        tfi = con.c1.r
        gi = con.c1.c
        tf_expr = t[:, tfi]
        g_expr = e[:, gi]
        corr  = np.corrcoef(tf_expr, g_expr)[0,1]
        prior_corr.append(corr)
    return (prior_corr, s)