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 )