Mercurial > repos > greg > cwpair2
changeset 5:2e0ddcc726f9 draft
Uploaded
author | greg |
---|---|
date | Tue, 24 Nov 2015 08:15:40 -0500 |
parents | 9ed566138ecb |
children | 001551406848 |
files | cwpair2_util.py |
diffstat | 1 files changed, 3 insertions(+), 19 deletions(-) [+] |
line wrap: on
line diff
--- a/cwpair2_util.py Sat Nov 21 09:01:46 2015 -0500 +++ b/cwpair2_util.py Tue Nov 24 08:15:40 2015 -0500 @@ -225,7 +225,7 @@ def process_file(dataset_path, galaxy_hid, method, threshold, up_distance, - down_distance, binsize, output_files, sort_chromosome, sort_score): + down_distance, binsize, output_files, sort_score): if method == 'all': match_methods = METHODS.keys() else: @@ -240,7 +240,6 @@ down_distance, binsize, output_files, - sort_chromosome, sort_score) statistics.append(stats) if output_files == 'all' and method == 'all': @@ -251,7 +250,7 @@ def perform_process(dataset_path, galaxy_hid, method, threshold, up_distance, - down_distance, binsize, output_files, sort_chromosome, sort_score): + down_distance, binsize, output_files, sort_score): output_details = output_files in ["all", "simple_orphan_detail"] output_plots = output_files in ["all"] output_orphans = output_files in ["all", "simple_orphan", "simple_orphan_detail"] @@ -367,27 +366,12 @@ orphan_output.writerow((cname, cpeak[0], cpeak[1], cpeak[2], cpeak[3])) # Keep track of orphans for statistics. orphans += len(crick) - # Sort output by chromosome if specified. - if sort_chromosome == "asc": - try: - x.sort(key=lambda data: int(data[3])) - x.sort(key=lambda data: int(data[0])) - except: - # Cannot sort because chromosome number is not a numeric. - pass - elif sort_chromosome == "desc": - try: - x.sort(key=lambda data: int(data[0]), reverse=True) - x.sort(key=lambda data: int(data[3]), reverse=True) - except: - # Cannot sort because chromosome number is not a numeric. - pass # Sort output by score if specified. if sort_score == "desc": x.sort(key=lambda data: float(data[5]), reverse=True) elif sort_score == "asc": x.sort(key=lambda data: float(data[5])) - # Writing a summary to txt or gff format file + # Writing a summary to gff format file for row in x: row_tmp = list(row) # Dataset in tuple cannot be modified in Python, so row will