def loadZslv(options, zslv_type, tohash=None): zslv_file = getBasename(os.path.basename(options.dbname)) + getFileSuffix(options, zslv_type, tohash) + ".zslv" abs_zslv_file = "%s/%s" % (options.output_path, zslv_file) if os.path.exists(abs_zslv_file): return zshelve.open("%s/%s" % (options.output_path, zslv_file)) else: print "%s not exists, run Stats.py first" % abs_zslv_file sys.exit(1)
help="The range of distance of motif pairs", metavar="INT") parser.add_option("-o", "--output", dest="output_path", help="Result output path", default="", metavar="DIR") return parser if __name__ == '__main__': parser = getOptions() (options, args) = parser.parse_args() if not options.output_path: abs_path = os.path.abspath(options.list_file) options.output_path = getBasename(os.path.basename(abs_path)) # abs_path = os.path.dirname(abs_path) # options.output_path = "%s/%s" % (abs_path, basename) if not os.path.exists(options.output_path): os.makedirs(options.output_path) united = UnitiedDB(options.dbname) REF_SeqNames = [united.SeqName2REF_SeqName[seqname] for seqname in loadList(options.list_file)] if options.motif_list_file: REF_uids = [united.Motif2uid[motif] for motif in loadList(options.motif_list_file)] else: REF_uids = [] # Motifs = zshelve.open(options.motif_zslv) # PairedMotifs = zshelve.open(options.paired_zslv) # Result = zshelve.open(options.result_file)