sys.exit(1)
else:
print "Input file is from MACS %s" % macs_version
# Don't try to convert output from MACS2
if macs_version.startswith("2."):
logging.error("input XLS comes from MACS %s, this version only handles 1.4" %
macs_version)
sys.exit(1)
# Sort into order by fold_enrichment and then by -10*log10(pvalue) column
data.sort(lambda line: line['fold_enrichment'],reverse=True)
data.sort(lambda line: line['-10*log10(pvalue)'],reverse=True)
# Restore first line
data.insert(0,tabdata=header_line)
# Insert "order" column
data.appendColumn("order")
# Perhaps confusingly must also insert initial value "#order"
data[0]['order'] = "#order"
for i in range(1,len(data)):
data[i]['order'] = i
# Reorder columns to put it at the start
data = data.reorderColumns(['order','chr','start','end','length','summit','tags',
'-10*log10(pvalue)','fold_enrichment','FDR(%)'])
# Legnds text
legends_text = """order\tSorting order Pvalue and FE
chr\tChromosome location of binding region
start\tStart coordinate of binding region
else:
print "Input file is from MACS %s" % macs_version
# Don't try to convert output from MACS2
if macs_version.startswith("2."):
logging.error(
"input XLS comes from MACS %s, this version only handles 1.4" %
macs_version)
sys.exit(1)
# Sort into order by fold_enrichment and then by -10*log10(pvalue) column
data.sort(lambda line: line['fold_enrichment'], reverse=True)
data.sort(lambda line: line['-10*log10(pvalue)'], reverse=True)
# Restore first line
data.insert(0, tabdata=header_line)
# Insert "order" column
data.appendColumn("order")
# Perhaps confusingly must also insert initial value "#order"
data[0]['order'] = "#order"
for i in range(1, len(data)):
data[i]['order'] = i
# Reorder columns to put it at the start
data = data.reorderColumns([
'order', 'chr', 'start', 'end', 'length', 'summit', 'tags',
'-10*log10(pvalue)', 'fold_enrichment', 'FDR(%)'
])
# Legnds text
legends_text = """order\tSorting order Pvalue and FE
