motifs_tfbs = MotifTools.load(file_tfbs) match_dict = {} for unknown in motifs_unknown: tf_list = [] for tfbs in motifs_tfbs: #print #print "Comparing motifs:" #print " %s vs %s" % (unknown.source, tfbs.source) #print " Unknown motif ( %s ) vs TFBS ( %s ) " % (unknown, tfbs) #print joined_motifs = [] joined_motifs.append(unknown) joined_motifs.append(tfbs) print joined_motifs Dmat = MotifCompare.computeDmat(joined_motifs) cluster = Kmedoids.bestaveKMedoids_cluster(Dmat,kmax=2,min_dist=0.2) for i in cluster[1].keys(): cluster_list = cluster[1][i] if len(cluster_list) > 1: dist = Dmat[0][1] tf_list.append( (dist, tfbs.source, tfbs.oneletter) ) #print #print "*** Motif match found! ***" #tfbs.giflogo("%s" % tfbs.oneletter) #for i in cluster_list: #print joined_motifs[i] novel_id = unknown.source tf_list.sort() match_dict[novel_id] = {'query_motif': unknown.oneletter, 'subject': tf_list }