diff hd.py @ 31:8beced3064e3 draft

planemo upload for repository https://github.com/monikaheinzl/duplexanalysis_galaxy/tree/master/tools/hd commit 033dd7b750f68e8aa68f327d7d72bd311ddbee4e-dirty
author mheinzl
date Wed, 14 Aug 2019 03:49:32 -0400
parents 46bfbec0f9e6
children b432cfa895ee
line wrap: on
line diff
--- a/hd.py	Wed Aug 07 04:01:32 2019 -0400
+++ b/hd.py	Wed Aug 14 03:49:32 2019 -0400
@@ -240,7 +240,10 @@
             lenTags, len_sample, len(numpy.concatenate(list1)))
     plt.text(0.14, -0.07, legend, size=12, transform=plt.gcf().transFigure)
 
-    legend2 = "SSCS ab = {:,}\nSSCS ba = {:,}\nDCS = {:,}".format(len(list1[1]), len(list1[2]), len(list1[0]))
+    legend2 = "SSCS ab = {:,} ({:.5f})\nSSCS ba = {:,} ({:.5f})\nDCS = {:,} ({:.5f})".format(
+        len(list1[1]), len(list1[1]) / float(nr_unique_chimeras),
+        len(list1[2]), len(list1[2]) / float(nr_unique_chimeras),
+        len(list1[0]), len(list1[0]) / float(nr_unique_chimeras))
     plt.text(0.6, -0.047, legend2, size=12, transform=plt.gcf().transFigure)
 
     pdf.savefig(fig, bbox_inches="tight")
@@ -1137,7 +1140,7 @@
                 sample_half_b = tag1[len(tag1) / 2:len(tag1)]
 
                 max_tags = data_chimeraAnalysis[i, 1]
-                if len(max_tags) > 1 and type(max_tags) is not numpy.ndarray:
+                if len(max_tags) > 1 and len(max_tags) !=  len(data_chimeraAnalysis[0, 1]) and type(max_tags) is not numpy.ndarray :
                     max_tags = numpy.concatenate(max_tags)
                 max_tags = numpy.unique(max_tags)
                 stat_maxTags.append(len(max_tags))