Mercurial > repos > drosofff > mi_rna_parser

changeset 5:59dfb33ca70e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author drosofff
date Mon, 23 Jun 2014 05:14:35 -0400
parents ac8fd41b8c98
children 11d042a0fbdf
files smRtools.py
diffstat 1 files changed, 3 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/smRtools.py	Sun Jun 22 18:34:29 2014 -0400
+++ b/smRtools.py	Mon Jun 23 05:14:35 2014 -0400
@@ -304,6 +304,8 @@
     for offset in self.readDict:
       for size in self.readDict[offset]:
         dicsize[size] = dicsize.get(size, 0) + 1
+    for offset in range (min(dicsize.keys()), max(dicsize.keys())+1):
+      dicsize[size] = dicsize.get(size, 0) # to fill offsets with null values
     return dicsize
     
   def statsizes (self, upstream_coord=None, downstream_coord=None):
@@ -393,7 +395,7 @@
     return mylist
 
   def readcoverage (self, upstream_coord=None, downstream_coord=None, windowName=None):
-    '''This method has not been tested yet 15-11-2013'''
+    '''Use by MirParser tool'''
     upstream_coord = upstream_coord or 1
     downstream_coord = downstream_coord or self.size
     windowName = windowName or "%s_%s_%s" % (self.gene, upstream_coord, downstream_coord)