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)
