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)