fsd: fsd.py annotate

annotate fsd.py @ 6:c4b8222dce29 draft

planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit 1667d50b521325eea1eada1b16029b1ba1379c5f

author	mheinzl
date	Wed, 09 May 2018 09:20:06 -0400
parents	69f47e0b804e
children	43b16247c35a

rev	line source
0 9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	1 #!/usr/bin/env python
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	2
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	3 # Family size distribution of SSCSs
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	4 #
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	5 # Author: Monika Heinzl, Johannes-Kepler University Linz (Austria)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	6 # Contact: monika.heinzl@edumail.at
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	7 #
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	8 # Takes at least one TABULAR file with tags before the alignment to the SSCS, but up to 4 files can be provided, as input.
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	9 # The program produces a plot which shows the distribution of family sizes of the all SSCSs from the input files and
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	10 # a CSV file with the data of the plot, as well as a TXT file with all tags of the DCS and their family sizes.
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	11 # If only one file is provided, then a family size distribution, which is separated after SSCSs without a partner and DCSs, is produced.
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	12 # Whereas a family size distribution with multiple data in one plot is produced, when more than one file (up to 4) is given.
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	13
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	14 # USAGE: python FSD_Galaxy_1.4_commandLine_FINAL.py filename --inputFile2 filename2 --inputFile3 filename3 --inputFile4 filename4 /
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	15 # --title_file outputFileName --sep "characterWhichSeparatesCSVFile"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	16
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	17 import numpy
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	18 import matplotlib.pyplot as plt
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	19 from matplotlib.backends.backend_pdf import PdfPages
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	20 import argparse
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	21 import sys
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	22 import os
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	23 import re
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	24
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	25 def readFileReferenceFree(file):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	26 with open(file, 'r') as dest_f:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	27 data_array = numpy.genfromtxt(dest_f, skip_header=0, delimiter='\t', comments='#', dtype='string')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	28 return(data_array)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	29
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	30 def make_argparser():
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	31 parser = argparse.ArgumentParser(description='Family Size Distribution of duplex sequencing data')
6 c4b8222dce29 planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit 1667d50b521325eea1eada1b16029b1ba1379c5f mheinzl parents: 5 diff changeset	32 parser.add_argument('--inputFile1',
c4b8222dce29 planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit 1667d50b521325eea1eada1b16029b1ba1379c5f mheinzl parents: 5 diff changeset	33 help='Tabular File with three columns: ab or ba, tag and family size.')
0 9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	34 parser.add_argument('--inputName1')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	35 parser.add_argument('--inputFile2',default=None,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	36 help='Tabular File with three columns: ab or ba, tag and family size.')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	37 parser.add_argument('--inputName2')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	38 parser.add_argument('--inputFile3',default=None,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	39 help='Tabular File with three columns: ab or ba, tag and family size.')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	40 parser.add_argument('--inputName3')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	41 parser.add_argument('--inputFile4',default=None,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	42 help='Tabular File with three columns: ab or ba, tag and family size.')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	43 parser.add_argument('--inputName4')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	44 parser.add_argument('--sep', default=",",
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	45 help='Separator in the csv file.')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	46 parser.add_argument('--output_csv', default="data.csv",type=str,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	47 help='Name of the pdf and csv file.')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	48 parser.add_argument('--output_pdf', default="data.pdf",type=str,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	49 help='Name of the pdf and csv file.')
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	50 return parser
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	51
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	52 def compare_read_families(argv):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	53 parser = make_argparser()
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	54 args=parser.parse_args(argv[1:])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	55
6 c4b8222dce29 planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit 1667d50b521325eea1eada1b16029b1ba1379c5f mheinzl parents: 5 diff changeset	56 firstFile = args.inputFile1
0 9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	57 name1 = args.inputName1
6 c4b8222dce29 planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit 1667d50b521325eea1eada1b16029b1ba1379c5f mheinzl parents: 5 diff changeset	58 #firstFile = args.inputName1
5 69f47e0b804e planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit 9829d657c07703b117bed2f259ddba432b244312 mheinzl parents: 3 diff changeset	59
0 9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	60 secondFile = args.inputFile2
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	61 name2 = args.inputName2
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	62 thirdFile = args.inputFile3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	63 name3 = args.inputName3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	64 fourthFile = args.inputFile4
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	65 name4 = args.inputName4
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	66
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	67 title_file = args.output_csv
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	68 title_file2 = args.output_pdf
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	69 sep = args.sep
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	70
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	71 if type(sep) is not str or len(sep)>1:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	72 print("Error: --sep must be a single character.")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	73 exit(4)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	74
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	75 plt.rc('figure', figsize=(11.69, 8.27)) # A4 format
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	76 plt.rcParams['patch.edgecolor'] = "black"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	77 plt.rcParams['axes.facecolor'] = "E0E0E0" # grey background color
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	78 plt.rcParams['xtick.labelsize'] = 12
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	79 plt.rcParams['ytick.labelsize'] = 12
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	80
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	81 list_to_plot = []
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	82 label = []
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	83 data_array_list = []
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	84
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	85 with open(title_file, "w") as output_file, PdfPages(title_file2) as pdf:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	86 fig = plt.figure()
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	87 plt.subplots_adjust(bottom=0.25)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	88 if firstFile != str(None):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	89 file1 = readFileReferenceFree(firstFile)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	90 integers = numpy.array(file1[:, 0]).astype(int) ## keep original family sizes
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	91
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	92 # for plot: replace all big family sizes by 22
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	93 data1 = numpy.array(file1[:, 0]).astype(int)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	94 bigFamilies = numpy.where(data1 > 20)[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	95 data1[bigFamilies] = 22
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	96
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	97 name1 = name1.split(".tabular")[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	98 list_to_plot.append(data1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	99 label.append(name1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	100 data_array_list.append(file1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	101
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	102 legend = "\n\n\n{}".format(name1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	103 plt.text(0.1, 0.11, legend, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	104 legend1 = "singletons:\nabsolute nr.\n{:,}".format(numpy.bincount(data1)[1])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	105 plt.text(0.4, 0.11, legend1, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	106
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	107 legend3 = "rel. freq\n{:.3f}".format(float(numpy.bincount(data1)[1]) / len(data1))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	108 plt.text(0.5, 0.11, legend3, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	109
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	110 legend4 = "family size > 20:\nabsolute nr.\n{:,}".format(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	111 numpy.bincount(data1)[len(numpy.bincount(data1)) - 1].astype(int))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	112 plt.text(0.6, 0.11, legend4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	113
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	114 legend5 = "rel. freq\n{:.3f}".format(float(numpy.bincount(data1)[len(numpy.bincount(data1)) - 1]) / len(data1))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	115 plt.text(0.7, 0.11, legend5, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	116
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	117 legend6 = "total length\n{:,}".format(len(data1))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	118 plt.text(0.8, 0.11, legend6, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	119
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	120 if secondFile != str(None):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	121 file2 = readFileReferenceFree(secondFile)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	122 data2 = numpy.asarray(file2[:, 0]).astype(int)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	123 bigFamilies2 = numpy.where(data2 > 20)[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	124 data2[bigFamilies2] = 22
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	125
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	126 list_to_plot.append(data2)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	127 name2 = name2.split(".tabular")[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	128 label.append(name2)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	129 data_array_list.append(file2)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	130
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	131 plt.text(0.1, 0.09, name2, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	132
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	133 legend1 = "{:,}".format(numpy.bincount(data2)[1])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	134 plt.text(0.4, 0.09, legend1, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	135
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	136 legend3 = "{:.3f}".format(float(numpy.bincount(data2)[1]) / len(data2))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	137 plt.text(0.5, 0.09, legend3, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	138
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	139 legend4 = "{:,}".format(numpy.bincount(data2)[len(numpy.bincount(data2)) - 1].astype(int))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	140 plt.text(0.6, 0.09, legend4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	141
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	142 legend5 = "{:.3f}".format(float(numpy.bincount(data2)[len(numpy.bincount(data2)) - 1]) / len(data2))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	143 plt.text(0.7, 0.09, legend5, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	144
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	145 legend6 = "{:,}".format(len(data2))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	146 plt.text(0.8, 0.09, legend6, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	147
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	148 if thirdFile != str(None):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	149 file3 = readFileReferenceFree(thirdFile)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	150
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	151 data3 = numpy.asarray(file3[:, 0]).astype(int)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	152 bigFamilies3 = numpy.where(data3 > 20)[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	153 data3[bigFamilies3] = 22
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	154
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	155 list_to_plot.append(data3)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	156 name3 = name3.split(".tabular")[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	157 label.append(name3)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	158 data_array_list.append(file3)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	159
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	160 plt.text(0.1, 0.07, name3, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	161
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	162 legend1 = "{:,}".format(numpy.bincount(data3)[1])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	163 plt.text(0.4, 0.07, legend1, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	164
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	165 legend3 = "{:.3f}".format(float(numpy.bincount(data3)[1]) / len(data3))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	166 plt.text(0.5, 0.07, legend3, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	167
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	168 legend4 = "{:,}".format(numpy.bincount(data3)[len(numpy.bincount(data3)) - 1].astype(int))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	169 plt.text(0.6, 0.07, legend4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	170
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	171 legend5 = "{:.3f}".format(float(numpy.bincount(data3)[len(numpy.bincount(data3)) - 1]) / len(data3))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	172 plt.text(0.7, 0.07, legend5, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	173
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	174 legend6 = "{:,}".format(len(data3))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	175 plt.text(0.8, 0.07, legend6, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	176
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	177 if fourthFile != str(None):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	178 file4 = readFileReferenceFree(fourthFile)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	179
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	180 data4 = numpy.asarray(file4[:, 0]).astype(int)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	181 bigFamilies4 = numpy.where(data4 > 20)[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	182 data4[bigFamilies4] = 22
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	183
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	184 list_to_plot.append(data4)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	185 name4 = name4.split(".tabular")[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	186 label.append(name4)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	187 data_array_list.append(file4)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	188
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	189 plt.text(0.1, 0.05, name4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	190
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	191 legend1 = "{:,}".format(numpy.bincount(data4)[1])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	192 plt.text(0.4, 0.05, legend1, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	193
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	194 legend4 = "{:.3f}".format(float(numpy.bincount(data4)[1]) / len(data4))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	195 plt.text(0.5, 0.05, legend4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	196
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	197 legend4 = "{:,}".format(numpy.bincount(data4)[len(numpy.bincount(data4)) - 1].astype(int))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	198 plt.text(0.6, 0.05, legend4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	199
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	200 legend5 = "{:.3f}".format(float(numpy.bincount(data4)[len(numpy.bincount(data4)) - 1]) / len(data4))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	201 plt.text(0.7, 0.05, legend5, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	202
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	203 legend6 = "{:,}".format(len(data4))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	204 plt.text(0.8, 0.05, legend6, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	205
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	206 maximumX = numpy.amax(numpy.concatenate(list_to_plot))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	207 minimumX = numpy.amin(numpy.concatenate(list_to_plot))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	208
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	209 counts = plt.hist(list_to_plot, bins=range(minimumX, maximumX + 1), stacked=False, edgecolor="black",
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	210 linewidth=1, label=label, align="left", alpha=0.7, rwidth=0.8)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	211
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	212 ticks = numpy.arange(minimumX - 1, maximumX, 1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	213 ticks1 = map(str, ticks)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	214 ticks1[len(ticks1) - 1] = ">20"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	215 plt.xticks(numpy.array(ticks), ticks1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	216
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	217 plt.legend(loc='upper right', fontsize=14, frameon=True, bbox_to_anchor=(0.9, 1))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	218 plt.title("Family Size Distribution", fontsize=14)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	219 plt.xlabel("No. of Family Members", fontsize=14)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	220 plt.ylabel("Absolute Frequency", fontsize=14)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	221 plt.margins(0.01, None)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	222 plt.grid(b=True, which="major", color="#424242", linestyle=":")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	223 pdf.savefig(fig)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	224 plt.close()
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	225
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	226 # write data to CSV file
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	227 output_file.write("Values from family size distribution with all datasets\n")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	228 output_file.write("\nFamily size")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	229 for i in label:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	230 output_file.write("{}{}".format(sep, i))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	231 output_file.write("{}sum".format(sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	232 output_file.write("\n")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	233 j = 0
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	234 for fs in counts[1][0:len(counts[1]) - 1]:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	235 if fs == 21:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	236 fs = ">20"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	237 else:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	238 fs = "={}".format(fs)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	239 output_file.write("FS{}{}".format(fs, sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	240 values_of_fs = []
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	241 if len(label) == 1:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	242 output_file.write("{}{}".format(int(counts[0][j]), sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	243 values_of_fs.append(int(counts[0][j]))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	244 else:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	245 for n in range(len(label)):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	246 output_file.write("{}{}".format(int(counts[0][n][j]), sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	247 values_of_fs.append(int(counts[0][n][j]))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	248 output_file.write("{}\n".format(sum(values_of_fs)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	249 j += 1
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	250 output_file.write("sum{}".format(sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	251 values_for_sum = []
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	252 if len(label) == 1:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	253 output_file.write("{}{}".format(int(sum(counts[0])), sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	254 values_for_sum.append(int(sum(counts[0])))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	255 else:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	256 for i in counts[0]:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	257 output_file.write("{}{}".format(int(sum(i)), sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	258 values_for_sum.append(int(sum(i)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	259
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	260 output_file.write("{}\n".format(sum(values_for_sum)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	261
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	262 ### Family size distribution after DCS and SSCS
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	263 for dataset, data, name_file in zip(list_to_plot, data_array_list, label):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	264 maximumX = numpy.amax(dataset)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	265 minimumX = numpy.amin(dataset)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	266
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	267 tags = numpy.array(data[:, 2])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	268 seq = numpy.array(data[:, 1])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	269 data = numpy.array(dataset)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	270
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	271 # find all unique tags and get the indices for ALL tags, but only once
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	272 u, index_unique, c = numpy.unique(numpy.array(seq), return_counts=True, return_index=True)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	273 d = u[c > 1]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	274
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	275 # get family sizes, tag for duplicates
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	276 duplTags_double = data[numpy.in1d(seq, d)]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	277 duplTags = duplTags_double[0::2] # ab of DCS
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	278 duplTagsBA = duplTags_double[1::2] # ba of DCS
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	279
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	280 duplTags_double_tag = tags[numpy.in1d(seq, d)]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	281 duplTags_double_seq = seq[numpy.in1d(seq, d)]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	282
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	283 # get family sizes for SSCS with no partner
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	284 ab = numpy.where(tags == "ab")[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	285 abSeq = seq[ab]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	286 ab = data[ab]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	287 ba = numpy.where(tags == "ba")[0]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	288 baSeq = seq[ba]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	289 ba = data[ba]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	290
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	291 dataAB = ab[numpy.in1d(abSeq, d, invert=True)]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	292 dataBA = ba[numpy.in1d(baSeq, d, invert=True)]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	293
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	294 # write DCS tags to file
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	295 # with open("DCS information_{}.txt".format(firstFile), "w") as file:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	296 # for t, s, f in zip(duplTags_double_tag, duplTags_double_seq, duplTags_double):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	297 # file.write("{}\t{}\t{}\n".format(t, s, f))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	298
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	299 list1 = [duplTags_double, dataAB, dataBA] # list for plotting
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	300
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	301 ## information for family size >= 3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	302 dataAB_FS3 = dataAB[dataAB >= 3]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	303 dataBA_FS3 = dataBA[dataBA >= 3]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	304 ab_FS3 = ab[ab >= 3]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	305 ba_FS3 = ba[ba >= 3]
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	306
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	307 duplTags_FS3 = duplTags[(duplTags >= 3) & (duplTagsBA >= 3)] # ab+ba with FS>=3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	308 duplTags_FS3_BA = duplTagsBA[(duplTags >= 3) & (duplTagsBA >= 3)] # ba+ab with FS>=3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	309 duplTags_double_FS3 = len(duplTags_FS3)+len(duplTags_FS3_BA) # both ab and ba strands with FS>=3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	310
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	311 fig = plt.figure()
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	312
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	313 plt.subplots_adjust(bottom=0.3)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	314 counts = plt.hist(list1, bins=range(minimumX, maximumX + 1), stacked=True,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	315 label=["duplex", "ab", "ba"], edgecolor="black", linewidth=1,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	316 align="left", color=["#FF0000", "#5FB404", "#FFBF00"])
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	317 # tick labels of x axis
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	318 ticks = numpy.arange(minimumX - 1, maximumX, 1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	319 ticks1 = map(str, ticks)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	320 ticks1[len(ticks1) - 1] = ">20"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	321 plt.xticks(numpy.array(ticks), ticks1)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	322 singl = counts[0][2][0] # singletons
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	323 last = counts[0][2][len(counts[0][0]) - 1] # large families
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	324
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	325 plt.legend(loc='upper right', fontsize=14, bbox_to_anchor=(0.9, 1), frameon=True)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	326 plt.title(name1, fontsize=14)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	327 plt.xlabel("No. of Family Members", fontsize=14)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	328 plt.ylabel("Absolute Frequency", fontsize=14)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	329 plt.margins(0.01, None)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	330 plt.grid(b=True, which="major", color="#424242", linestyle=":")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	331
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	332 ## extra information beneath the plot
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	333 legend = "SSCS ab= \nSSCS ba= \nDCS (total)= \nlength of dataset="
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	334 plt.text(0.1, 0.09, legend, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	335
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	336 legend = "absolute numbers\n\n{:,}\n{:,}\n{:,} ({:,})\n{:,}" \
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	337 .format(len(dataAB), len(dataBA), len(duplTags), len(duplTags_double),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	338 (len(dataAB) + len(dataBA) + len(duplTags)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	339 plt.text(0.35, 0.09, legend, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	340
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	341 legend = "relative frequencies\nunique\n{:.3f}\n{:.3f}\n{:.3f}\n{:,}" \
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	342 .format(float(len(dataAB)) / (len(dataAB) + len(dataBA) + len(duplTags)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	343 float(len(dataBA)) / (len(dataAB) + len(dataBA) + len(duplTags)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	344 float(len(duplTags)) / (len(dataAB) + len(dataBA) + len(duplTags)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	345 (len(dataAB) + len(dataBA) + len(duplTags)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	346 plt.text(0.54, 0.09, legend, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	347
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	348 legend = "total\n{:.3f}\n{:.3f}\n{:.3f} ({:.3f})\n{:,}" \
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	349 .format(float(len(dataAB)) / (len(ab) + len(ba)), float(len(dataBA)) / (len(ab) + len(ba)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	350 float(len(duplTags)) / (len(ab) + len(ba)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	351 float(len(duplTags_double)) / (len(ab) + len(ba)), (len(ab) + len(ba)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	352 plt.text(0.64, 0.09, legend, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	353
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	354 legend1 = "\nsingletons:\nfamily size > 20:"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	355 plt.text(0.1, 0.03, legend1, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	356
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	357 legend4 = "{:,}\n{:,}".format(singl.astype(int), last.astype(int))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	358 plt.text(0.35, 0.03, legend4, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	359
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	360 legend3 = "{:.3f}\n{:.3f}".format(singl / len(data),last / len(data))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	361 plt.text(0.54, 0.03, legend3, size=12, transform=plt.gcf().transFigure)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	362
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	363 pdf.savefig(fig)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	364 plt.close()
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	365
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	366 # write same information to a csv file
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	367 count = numpy.bincount(integers) # original counts of family sizes
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	368 output_file.write("\nDataset:{}{}\n".format(sep, name_file))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	369 output_file.write("max. family size:{}{}\n".format(sep, max(integers)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	370 output_file.write("absolute frequency:{}{}\n".format(sep, count[len(count) - 1]))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	371 output_file.write("relative frequency:{}{:.3f}\n\n".format(sep, float(count[len(count) - 1]) / sum(count)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	372
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	373 output_file.write("{}singletons:{}{}family size > 20:\n".format(sep, sep, sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	374 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	375 "{}absolute nr.{}rel. freq{}absolute nr.{}rel. freq{}total length\n".format(sep, sep, sep, sep, sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	376 output_file.write("{}{}{}{}{:.3f}{}{}{}{:.3f}{}{}\n\n".format(name_file, sep, singl.astype(int), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	377 singl / len(data), sep,last.astype(int), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	378 last / len(data), sep, len(data)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	379
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	380 ## information for FS >= 1
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	381 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	382 "The unique frequencies were calculated from the dataset where the tags occured only once (=ab without DCS, ba without DCS)\n" \
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	383 "Whereas the total frequencies were calculated from the whole dataset (=including the DCS).\n\n")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	384 output_file.write("FS >= 1{}{}unique:{}total:\n".format(sep, sep, sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	385 output_file.write("nr./rel. freq of ab={}{}{}{:.3f}{}{:.3f}\n".format(sep, len(dataAB), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	386 float(len(dataAB)) / (len(dataAB) + len(dataBA) + len( duplTags)), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	387 float(len(dataAB)) / (len(ab) + len(ba))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	388 output_file.write("nr./rel. freq of ba={}{}{}{:.3f}{}{:.3f}\n".format(sep, len(dataBA), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	389 float(len(dataBA)) / (len(dataBA) + len(dataBA) + len(duplTags)), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	390 float(len(dataBA)) / (len(ba) + len(ba))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	391 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	392 "nr./rel. freq of DCS (total)={}{} ({}){}{:.3f}{}{:.3f} ({:.3f})\n".format(sep, len(duplTags), len(duplTags_double), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	393 float(len(duplTags)) / ( len(dataAB) + len( dataBA) + len(duplTags)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	394 sep, float(len(duplTags)) / ( len(ab) + len(ba)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	395 float( len(duplTags_double)) / (len(ab) + len(ba))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	396 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	397 "length of dataset={}{}{}{}{}{}\n".format(sep, (len(dataAB) + len(dataBA) + len(duplTags)), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	398 (len(dataAB) + len(dataBA) + len(duplTags)), sep,(len(ab) + len(ba))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	399 ## information for FS >= 3
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	400 output_file.write("FS >= 3{}{}unique:{}total:\n".format(sep, sep, sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	401 output_file.write("nr./rel. freq of ab={}{}{}{:.3f}{}{:.3f}\n".format(sep, len(dataAB_FS3), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	402 float(len(dataAB_FS3)) / (len(dataAB_FS3) + len(dataBA_FS3) + len(duplTags_FS3)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	403 sep, float(len(dataAB_FS3)) / ( len(ab_FS3) + len(ba_FS3))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	404 output_file.write("nr./rel. freq of ba={}{}{}{:.3f}{}{:.3f}\n".format(sep, len(dataBA_FS3), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	405 float(len(dataBA_FS3)) / ( len(dataBA_FS3) + len(dataBA_FS3) + len(duplTags_FS3)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	406 sep,float(len(dataBA_FS3)) / (len(ba_FS3) + len(ba_FS3))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	407 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	408 "nr./rel. freq of DCS (total)={}{} ({}){}{:.3f}{}{:.3f} ({:.3f})\n".format(sep, len(duplTags_FS3),duplTags_double_FS3,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	409 sep, float(len( duplTags_FS3)) / (len(dataBA_FS3) + len(duplTags_FS3)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	410 sep, float(len(duplTags_FS3)) / (len(ab_FS3) + len(ba_FS3)),
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	411 float(duplTags_double_FS3) / (len(ab_FS3) + len(ba_FS3))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	412 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	413 "length of dataset={}{}{}{}{}{}\n".format(sep, (len(dataAB_FS3) + len(dataBA_FS3) + len(duplTags_FS3)), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	414 (len(dataAB_FS3) + len(dataBA_FS3) + len(duplTags_FS3)), sep,
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	415 (len(ab_FS3) + len(ba_FS3))))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	416
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	417 output_file.write("\nValues from family size distribution\n")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	418 output_file.write("{}duplex{}ab{}ba{}sum\n".format(sep, sep, sep, sep))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	419 for dx, ab, ba, fs in zip(counts[0][0], counts[0][1], counts[0][2], counts[1]):
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	420 if fs == 21:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	421 fs = ">20"
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	422 else:
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	423 fs = "={}".format(fs)
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	424 ab1 = ab - dx
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	425 ba1 = ba - ab
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	426 output_file.write(
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	427 "FS{}{}{}{}{}{}{}{}{}\n".format(fs, sep, int(dx), sep, int(ab1), sep, int(ba1), sep, int(ba)))
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	428
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	429 print("Files successfully created!")
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	430
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	431 if __name__ == '__main__':
9736b9d04a0b planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/fsd commit f674213e798956531c935e7b9eb7f444286d0a5e-dirty mheinzl parents: diff changeset	432 sys.exit(compare_read_families(sys.argv))

Mercurial > repos > mheinzl > fsd

annotate fsd.py @ 6:c4b8222dce29 draft