annotate viz_graphs.py @ 17:a09d238416ba draft

Uploaded
author glogobyte
date Wed, 28 Oct 2020 07:42:50 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
17
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
1 import sys
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
2
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
3 # Read a file and return it as a list
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
4 def read(path, flag):
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
5 if flag == 0:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
6 with open(path) as fp:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
7 file=fp.readlines()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
8 fp.close()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
9 return file
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
10
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
11 if flag == 1:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
12 with open(path) as fp:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
13 file = fp.read().splitlines()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
14 fp.close()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
15 return file
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
16
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
17
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
18 #################################################################################################
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
19 def pdf_before_DE(analysis):
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
20
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
21 # Import FPDF class
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
22 from fpdf import FPDF, fpdf
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
23
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
24 # Import glob module to find all the files matching a pattern
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
25 import glob
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
26
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
27 # Image extensions
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
28 if analysis=="2":
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
29 image_extensions = ("c_hist_red.png","t_hist_red.png","pie_non.png","spider_red.png","spider_non_red.png","c_logo.png","t_logo.png","c_bar.png","t_bar.png")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
30 else:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
31 image_extensions = ("c_hist_red.png","t_hist_red.png","pie_tem.png","spider_red.png","spider_non_red.png")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
32 # This list will hold the images file names
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
33 images = []
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
34
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
35 # Build the image list by merging the glob results (a list of files)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
36 # for each extension. We are taking images from current folder.
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
37 for extension in image_extensions:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
38 images.extend(glob.glob(extension))
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
39 #sys.exit(images)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
40 # Create instance of FPDF class
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
41 pdf = FPDF('P', 'in', 'A4')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
42 # Add new page. Without this you cannot create the document.
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
43 pdf.add_page()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
44 # Set font to Arial, 'B'old, 16 pts
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
45 pdf.set_font('Arial', 'B', 20.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
46
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
47 # Page header
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
48 pdf.cell(pdf.w-0.5, 0.5, 'IsomiR Profile Report',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
49 pdf.ln(0.7)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
50 pdf.set_font('Arial','', 16.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
51 pdf.cell(pdf.w-0.5, 0.5, 'sRNA Length Distribution',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
52
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
53 # Smaller font for image captions
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
54 pdf.set_font('Arial', '', 11.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
55
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
56 # Image caption
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
57 pdf.ln(0.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
58
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
59 yh=FPDF.get_y(pdf)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
60 pdf.image(images[0],x=0.3,w=4, h=3)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
61 pdf.image(images[1],x=4,y=yh, w=4, h=3)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
62 pdf.ln(0.3)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
63
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
64 # Image caption
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
65 pdf.cell(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
66 pdf.cell(3.0, 0.0, " Mapped and unmapped reads to custom precussor arm reference DB (5p and 3p arms) in Control (left)")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
67 pdf.ln(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
68 pdf.cell(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
69 pdf.cell(3.0, 0.0, " and Treated (right) groups")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
70
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
71
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
72 pdf.ln(0.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
73 h1=FPDF.get_y(pdf)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
74 pdf.image(images[2],x=1, w=6.5, h=5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
75 h2=FPDF.get_y(pdf)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
76 FPDF.set_y(pdf,h1+0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
77 pdf.set_font('Arial','', 14.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
78 pdf.cell(pdf.w-0.5, 0.5, 'Template and non-template IsomiRs',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
79 pdf.set_font('Arial', '', 11.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
80 FPDF.set_y(pdf,h2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
81 FPDF.set_y(pdf,9.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
82 # Image caption
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
83 pdf.cell(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
84 if analysis=="2":
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
85 pdf.cell(3.0, 0.0, " Template, non-template, miRNA reference and unmapped sequences as percentage of total sRNA")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
86 else:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
87 pdf.cell(3.0, 0.0, " Template, miRNA reference and unmapped sequences as percentage of total sRNA")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
88 pdf.ln(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
89 pdf.cell(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
90 pdf.cell(3.0, 0.0, " reads in Control (left) and treated (right) groups")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
91
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
92
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
93
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
94 pdf.add_page()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
95 pdf.set_font('Arial', 'B', 16.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
96 pdf.cell(pdf.w-0.5, 0.5, "Reference form and isomiR among total miRNA reads",align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
97 pdf.ln(0.7)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
98 pdf.set_font('Arial', 'B', 12.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
99 pdf.cell(pdf.w-0.5, 0.5, "Template isomiR profile (redundant)",align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
100 pdf.ln(0.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
101 pdf.image(images[3],x=1.5, w=5.5, h=4)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
102 pdf.ln(0.6)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
103 pdf.cell(pdf.w-0.5, 0.0, "Template isomiR profile (non-redundant)",align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
104 pdf.set_font('Arial', '', 12.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
105 pdf.ln(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
106 pdf.image(images[4],x=1.5, w=5.5, h=4)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
107 pdf.ln(0.3)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
108 pdf.set_font('Arial', '', 11.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
109 pdf.cell(0.2)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
110 pdf.cell(3.0, 0.0, " * IsomiRs potentialy initiated from multiple loci")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
111
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
112
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
113 if analysis=="2":
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
114 pdf.add_page('L')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
115
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
116 pdf.set_font('Arial', 'B', 16.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
117 pdf.cell(pdf.w-0.5, 0.5, "Non-template IsomiRs",align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
118 pdf.ln(0.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
119 pdf.set_font('Arial', 'B', 12.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
120 pdf.cell(pdf.w-0.5, 0.5, "3' Additions of reference of isomiR sequence",align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
121 pdf.ln(0.7)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
122
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
123 yh=FPDF.get_y(pdf)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
124 pdf.image(images[5],x=1.5,w=3.65, h=2.65)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
125 pdf.image(images[7],x=6.5,y=yh, w=3.65, h=2.65)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
126 pdf.ln(0.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
127 yh=FPDF.get_y(pdf)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
128 pdf.image(images[6],x=1.5,w=3.65, h=2.65)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
129 pdf.image(images[8],x=6.5,y=yh, w=3.65, h=2.65)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
130
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
131 pdf.close()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
132 pdf.output('report1.pdf','F')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
133
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
134
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
135
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
136
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
137 #############################################################################################################################################################3
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
138
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
139 def pdf_after_DE(analysis,top):
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
140
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
141 # Import FPDF class
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
142 from fpdf import FPDF
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
143
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
144 # Import glob module to find all the files matching a pattern
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
145 import glob
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
146
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
147 # Image extensions
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
148 if analysis=="2":
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
149 image_extensions = ("tem.png","a2.png","non.png")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
150 else:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
151 image_extensions = ("tem.png","a2.png")
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
152
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
153 # This list will hold the images file names
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
154 images = []
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
155
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
156 # Build the image list by merging the glob results (a list of files)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
157 # for each extension. We are taking images from current folder.
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
158 for extension in image_extensions:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
159 images.extend(glob.glob(extension))
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
160 #sys.exit(images)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
161 # Create instance of FPDF class
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
162 pdf = FPDF('P', 'in', 'letter')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
163 # Add new page. Without this you cannot create the document.
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
164 pdf.add_page()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
165 # Set font to Arial, 'B'old, 16 pts
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
166 pdf.set_font('Arial', 'B', 16.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
167
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
168 # Page header
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
169 pdf.cell(pdf.w-0.5, 0.5, 'Differential expression of miRNAs and Isoforms',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
170 #pdf.ln(0.25)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
171
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
172 pdf.ln(0.7)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
173 pdf.set_font('Arial','B', 12.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
174 if "tem.png" in images:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
175 pdf.cell(pdf.w-0.5, 0.5, 'Top '+top+' most differentially expressed miRNA and template isoforms',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
176 # Smaller font for image captions
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
177 pdf.set_font('Arial', '', 10.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
178 # Image caption
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
179 pdf.ln(0.4)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
180 pdf.image(images[images.index("tem.png")],x=0.8, w=7, h=8)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
181 pdf.ln(0.3)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
182 pdf.set_font('Arial','B', 12.0)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
183 else:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
184 print("WARNING: There aren't miRNAs which fullfiled these criteria" )
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
185
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
186 if "non.png" in images and analysis=="2":
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
187 if "tem.png" in images: pdf.add_page()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
188 pdf.ln(0.7)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
189 pdf.cell(pdf.w-0.5, 0.5, 'Top '+top+' most differentially expressed non-template isomiRs',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
190 pdf.ln(0.4)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
191 pdf.image(images[images.index("non.png")],x=0.5, w=7.5, h=6.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
192 else:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
193 print("WARNING: There aren't non-template miRNAs which fullfiled these criteria" )
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
194
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
195
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
196 if "a2.png" in images:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
197 if len(images)>=2: pdf.add_page()
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
198 pdf.ln(0.5)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
199 pdf.cell(pdf.w-0.5, 0.5, 'Top '+top+' most differentially expressed miRNAs and isomiRs grouped by arm',align='C')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
200 pdf.ln(0.4)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
201 pdf.image(images[images.index("a2.png")],x=0.8, w=7, h=8)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
202 pdf.ln(0.3)
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
203 else:
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
204 print("WARNING: There aren't non-template miRNAs which fullfiled these criteria" )
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
205
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
206
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
207 pdf.output('report2.pdf', 'F')
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
208
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
209
a09d238416ba Uploaded
glogobyte
parents:
diff changeset
210