| Previous changeset 20:b084b6a8e3ac (2018-12-14) Next changeset 22:7e570ba56b83 (2019-02-27) |
|
Commit message:
planemo upload for repository https://github.com/monikaheinzl/duplexanalysis_galaxy/tree/master/tools/hd commit b8a2f7b7615b2bcd3b602027af31f4e677da94f6 |
|
modified:
hd.py |
| b |
| diff -r b084b6a8e3ac -r 9919024d7778 hd.py --- a/hd.py Fri Dec 14 04:31:21 2018 -0500 +++ b/hd.py Fri Dec 14 05:03:24 2018 -0500 |
| [ |
| @@ -76,11 +76,15 @@ plt.ylim((0, max(numpy.bincount(numpy.concatenate(familySizeList1))) * 1.1)) plt.ylim((0, maximumY * 1.2)) - legend = "\nmax. family size: \nabsolute frequency: \nrelative frequency: " + legend = "\nfamily size: \nabsolute frequency: \nrelative frequency: " plt.text(0.15, -0.08, legend, size=12, transform=plt.gcf().transFigure) count = numpy.bincount(originalCounts) # original counts - legend1 = "{}\n{}\n{:.5f}".format(max(originalCounts), count[len(count) - 1], float(count[len(count) - 1]) / sum(count)) + if max(originalCounts) >= 20: + max_count = ">= 20" + else: + max_count = max(originalCounts) + legend1 = "{}\n{}\n{:.5f}".format(max_count, count[len(count) - 1], float(count[len(count) - 1]) / sum(count)) plt.text(0.5, -0.08, legend1, size=12, transform=plt.gcf().transFigure) legend3 = "singletons\n{:,}\n{:.5f}".format(int(counts[0][len(counts[0]) - 1][1]), float(counts[0][len(counts[0]) - 1][1]) / sum(counts[0][len(counts[0]) - 1])) plt.text(0.7, -0.08, legend3, transform=plt.gcf().transFigure, size=12) @@ -959,13 +963,13 @@ "Family size distribution separated by Hamming distance", sep, diff=False) - count = numpy.bincount(quant) # output_file.write("{}{}\n".format(sep, name1)) output_file.write("\n") - output_file.write("max. family size:{}{}\n".format(sep, max(quant))) - output_file.write("absolute frequency:{}{}\n".format(sep, count[len(count) - 1])) + max_fs = numpy.bincount(integers[result]) + output_file.write("max. family size in sample:{}{}\n".format(sep, max(integers[result]))) + output_file.write("absolute frequency:{}{}\n".format(sep, max_fs[len(max_fs) - 1])) output_file.write( - "relative frequency:{}{}\n\n".format(sep, float(count[len(count) - 1]) / sum(count))) + "relative frequency:{}{}\n\n".format(sep, float(max_fs[len(max_fs) - 1]) / sum(max_fs))) # HD within tags output_file.write( |