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 )
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 )