Mercurial > repos > mheinzl > fsd_regions

comparison fsd_regions.py @ 13:63432e6f6a61 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/monikaheinzl/duplexanalysis_galaxy/tree/master/tools/fsd_regions commit edbe0509294e5c23c26cdae659ce5cf3ff69c363
author mheinzl
date Fri, 07 Dec 2018 07:23:07 -0500
parents 37db9decb5d0
children 6879295d3f11
comparison
equal deleted inserted replaced
12:fa3ff1c5280d 13:63432e6f6a61
199 199
200 legend = "BA\n{}\n{}\n{:.5f}" \ 200 legend = "BA\n{}\n{}\n{:.5f}" \
201 .format(max(map(int, quant_ba)), count2[len(count2) - 1], float(count2[len(count2) - 1]) / sum(count2)) 201 .format(max(map(int, quant_ba)), count2[len(count2) - 1], float(count2[len(count2) - 1]) / sum(count2))
202 plt.text(0.45, 0.1475, legend, size=11, transform=plt.gcf().transFigure) 202 plt.text(0.45, 0.1475, legend, size=11, transform=plt.gcf().transFigure)
203 203
204 plt.text(0.55, 0.2125, "total nr. of tags:", size=11, transform=plt.gcf().transFigure) 204 plt.text(0.53, 0.2125, "total nr. of tags:", size=11, transform=plt.gcf().transFigure)
205 plt.text(0.8, 0.2125, "{:,} ({:,})".format(length_regions, length_regions / 2), size=11, 205 plt.text(0.85, 0.2125, "{:,} ({:,})".format(length_regions, length_regions / 2), size=11,
206 transform=plt.gcf().transFigure) 206 transform=plt.gcf().transFigure)
207 207
208 legend4 = "* In the plot, both family sizes of the ab and ba strands were used.\nWhereas the total numbers indicate only the single count of the tags per region.\n" 208 legend4 = "* In the plot, both family sizes of the ab and ba strands were used.\nWhereas the total numbers indicate only the single count of the tags per region.\n"
209 plt.text(0.1, 0.01, legend4, size=11, transform=plt.gcf().transFigure) 209 plt.text(0.1, 0.01, legend4, size=11, transform=plt.gcf().transFigure)
210 210
211 space = 0 211 space = 0
212 for i, count in zip(group, quantAfterRegion): 212 for i, count in zip(group, quantAfterRegion):
213 plt.text(0.55, 0.15 - space, "{}:\n".format(i), size=11, transform=plt.gcf().transFigure) 213 plt.text(0.53, 0.15 - space, "{}:\n".format(i), size=11, transform=plt.gcf().transFigure)
214 plt.text(0.8, 0.15 - space, "{:,}\n".format(len(count) / 2), size=11, transform=plt.gcf().transFigure) 214 plt.text(0.85, 0.15 - space, "{:,}\n".format(len(count) / 2), size=11, transform=plt.gcf().transFigure)
215 space = space + 0.02 215 space = space + 0.02
216 216
217 plt.legend(loc='upper right', fontsize=14, bbox_to_anchor=(0.9, 1), frameon=True) 217 plt.legend(loc='upper right', fontsize=14, bbox_to_anchor=(0.9, 1), frameon=True)
218 plt.xlabel("Family size", fontsize=14) 218 plt.xlabel("Family size", fontsize=14)
219 plt.ylabel("Absolute Frequency", fontsize=14) 219 plt.ylabel("Absolute Frequency", fontsize=14)