def run_merge_rnaseq_metrics(input_files, summary_file):
"""Merge the outputs of collectrnaseqmetrics into one tab-separated file.
"""
metrics = []
for fn in input_files:
metrics.extend(picard_helpers.parse_picard_rnaseq_metrics(fn))
fieldnames = metrics[0].keys()
with open(summary_file, 'w') as fou:
dw = csv.DictWriter(fou, delimiter='\t', fieldnames=fieldnames)
dw.writeheader()
dw.writerows(metrics)
def run_merge_rnaseq_metrics(input_files, summary_file):
"""Merge the outputs of collectrnaseqmetrics into one tab-separated file.
"""
metrics = []
for fn in input_files:
## Only take the first non-comment line; there is some histogram data after it!
metrics.extend(picard_helpers.parse_picard_rnaseq_metrics(fn)[0:1])
fieldnames = metrics[0].keys()
with open(summary_file, 'w') as fou:
dw = csv.DictWriter(fou, delimiter='\t', fieldnames=fieldnames)
dw.writeheader()
dw.writerows(metrics)
