#!/usr/bin/python import os, csv, sys from TF_Module import * #Example of using the class TF_Module #Data of TF-genes binding relaton TF_gene_relation = 'protein-dna_interactions_Gid.txt' #Create an object for the calss TF_Module ObjectTF = TF_Module(TF_gene_relation) #A list of co-regulated genes Co_reg_genes = ['AGA1', 'FAR1', 'FUS1', 'GPA1', 'KAR4', 'SST2', 'STE12', 'STE2', 'STE6', 'TEC1'] #Find the TF module that regulates the co-regulated genes TF_module = ObjectTF. TF_Module_Greedy(Co_reg_genes) #print result: print 'The TF module is:', TF_module