Mercurial > repos > artbio > small_rna_maps
comparison small_rna_maps.py @ 27:fe1a9cfaf5c3 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_maps commit d234ef45deb84f70c85b3372b8e0137df0be4e29
author | artbio |
---|---|
date | Wed, 24 Apr 2019 11:18:24 -0400 |
parents | 376fae7c9f32 |
children | 14adf24603b6 |
comparison
equal
deleted
inserted
replaced
26:376fae7c9f32 | 27:fe1a9cfaf5c3 |
---|---|
372 def main(inputs, samples, methods, outputs, minsize, maxsize, cluster, | 372 def main(inputs, samples, methods, outputs, minsize, maxsize, cluster, |
373 nostrand, bedfile=None, bed_skipsize=0): | 373 nostrand, bedfile=None, bed_skipsize=0): |
374 for method, output in zip(methods, outputs): | 374 for method, output in zip(methods, outputs): |
375 out = open(output, 'w') | 375 out = open(output, 'w') |
376 if method == 'Size': | 376 if method == 'Size': |
377 header = ["Dataset", "Chromosome", "Polarity", method, "Counts", | 377 header = ["# Dataset", "Chromosome", "Polarity", method, "Counts", |
378 "Strandness", "z-score"] | 378 "Strandness", "z-score"] |
379 elif cluster: | 379 elif cluster: |
380 header = ["Dataset", "Chromosome", "Chrom_length", "Coordinate", | 380 header = ["# Dataset", "Chromosome", "Chrom_length", "Coordinate", |
381 "Polarity", method, "Start-End", "Cluster Size", | 381 "Polarity", method, "Start-End", "Cluster Size", |
382 "density"] | 382 "density"] |
383 else: | 383 else: |
384 header = ["Dataset", "Chromosome", "Chrom_length", "Coordinate", | 384 header = ["# Dataset", "Chromosome", "Chrom_length", "Coordinate", |
385 "Polarity", method] | 385 "Polarity", method] |
386 out.write('\t'.join(header) + '\n') | 386 out.write('\t'.join(header) + '\n') |
387 for input, sample in zip(inputs, samples): | 387 for input, sample in zip(inputs, samples): |
388 mapobj = Map(input, sample, minsize, maxsize, cluster, nostrand) | 388 mapobj = Map(input, sample, minsize, maxsize, cluster, nostrand) |
389 token = {"Counts": mapobj.compute_readcount, | 389 token = {"Counts": mapobj.compute_readcount, |