# HG changeset patch # User artbio # Date 1500795820 14400 # Node ID f924a33e1eef62b55e526418668061c07fb90bd5 # Parent d65045e976e634d0c5a344384ea3de92a7642d73 planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/small_rna_map commit 19a447cd21c746c852cf7434b5df423504baf383 diff -r d65045e976e6 -r f924a33e1eef small_rna_map.py --- a/small_rna_map.py Sat Jul 22 11:45:52 2017 -0400 +++ b/small_rna_map.py Sun Jul 23 03:43:40 2017 -0400 @@ -121,13 +121,15 @@ Takes a map_dictionary and returns a dictionary of sizes: {chrom: {size: {polarity: nbre of reads}}} ''' - size_dictionary = defaultdict(lambda: defaultdict(lambda: defaultdict( int ))) + size_dictionary = defaultdict(lambda: defaultdict( + lambda: defaultdict( int ))) for key in map_dictionary: + if len(map_dictionary) == 0: + # to track empty chromosomes + size_dictionary[key[0]][key[2]][size] = 0 + continue for size in map_dictionary[key]: - try: - size_dictionary[key[0]][key[2]][size] += 1 - except KeyError: - size_dictionary[key[0]][key[2]][size] = 1 + size_dictionary[key[0]][key[2]][size] += 1 return size_dictionary def write_size_table(self, out): diff -r d65045e976e6 -r f924a33e1eef small_rna_map.xml --- a/small_rna_map.xml Sat Jul 22 11:45:52 2017 -0400 +++ b/small_rna_map.xml Sun Jul 23 03:43:40 2017 -0400 @@ -1,4 +1,4 @@ - + numpy