# HG changeset patch # User mheinzl # Date 1565768972 14400 # Node ID 8beced3064e38ccc87c1eedf8ac89b9aabd02dad # Parent 46bfbec0f9e6f2821def84a1f2e07809f7e81c82 planemo upload for repository https://github.com/monikaheinzl/duplexanalysis_galaxy/tree/master/tools/hd commit 033dd7b750f68e8aa68f327d7d72bd311ddbee4e-dirty diff -r 46bfbec0f9e6 -r 8beced3064e3 hd.py --- 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)) diff -r 46bfbec0f9e6 -r 8beced3064e3 hd.xml --- a/hd.xml Wed Aug 07 04:01:32 2019 -0400 +++ b/hd.xml Wed Aug 14 03:49:32 2019 -0400 @@ -1,5 +1,5 @@ - + hamming distance analysis of duplex tags python diff -r 46bfbec0f9e6 -r 8beced3064e3 test-data/hd_output.pdf Binary file test-data/hd_output.pdf has changed