Example #1
0
def look_for_factor_pair( name, all_pairs, pssms1, pssms2 ):
    both_pairs = biopsy.get_pairs_from( all_pairs, pssms1, pssms2 )
    biopsy.sort_pairs_by_log_odds( both_pairs )
    print 'Found %d pairs' % len( both_pairs )
    for p in both_pairs:
        print p.binder_pair, p.log_odds_ratio
    biopsy.write_pair_separation_histograms( both_pairs, name )
Example #2
0
def test_get_pairs_for_pssm():
    pssm = 'M00736'
    pairs = biopsy.load(
            'c:/data/remos/100/100.top_pairs',
            display = 'pairs log odds'
    )
    pssm_pairs = biopsy.get_pairs_for_pssm(
            pairs,
            pssm
    )
    print 'Got %d pairs for %s from total of %d' % \
            ( len(pssm_pairs), pssm, len(pairs) )
    for p in pssm_pairs[:20]:
        print p.binder_pair, p.log_odds_ratio
    biopsy.write_pair_separation_histograms( pssm_pairs[:20], pssm )