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,