@@ -57,21 +57,21 @@ def dataminer(indfile_name, fmt, popfile=None):
5757 qvalues = np .array ([])
5858
5959 # Start file parsing
60- r = False
60+ parse = False
6161 with open (indfile_name ) as file_handle :
6262 for line in file_handle :
6363 if line .strip ().lower ().startswith ("inferred ancestry of "
6464 "individuals:" ):
6565 # Enter parse mode ON
66- r = True
66+ parse = True
6767 # Skip subheader
6868 next (file_handle )
6969 elif line .strip ().lower ().startswith ("estimated allele "
7070 "frequencies in each "
7171 "cluster" ):
7272 # parse mode OFF
73- r = False
74- elif r :
73+ parse = False
74+ elif parse :
7575 if line .strip () != "" :
7676 fields = line .strip ().split ()
7777 # Get cluster values
@@ -130,27 +130,27 @@ def plotter(qvalues, poplist, outfile):
130130 if i == 0 :
131131 ax .bar (range (numinds ), qvalues [:, i ], facecolor = clr ,
132132 edgecolor = "none" , width = 1 )
133- formerQ = qvalues [:, i ]
133+ former_q = qvalues [:, i ]
134134 else :
135- ax .bar (range (numinds ), qvalues [:, i ], bottom = formerQ ,
135+ ax .bar (range (numinds ), qvalues [:, i ], bottom = former_q ,
136136 facecolor = clr , edgecolor = "none" , width = 1 )
137- formerQ = formerQ + qvalues [:, i ]
137+ former_q = former_q + qvalues [:, i ]
138138
139139 # Annotate population info
140140 if poplist :
141- c = 1
142- for p , vals in enumerate (poplist ):
141+ count = 1
142+ for ppl , vals in enumerate (poplist ):
143143 # Add population delimiting lines
144144 plt .axvline (x = vals [0 ], linewidth = 1.5 , color = 'black' )
145145 # Add population labels
146146 # Determine x pos
147- xpos = vals [0 ] - ((vals [0 ] - poplist [p - 1 ][0 ]) / 2 ) if p > 0 \
147+ xpos = vals [0 ] - ((vals [0 ] - poplist [ppl - 1 ][0 ]) / 2 ) if ppl > 0 \
148148 else vals [0 ] / 2
149149 # Draw text
150- ax .text (xpos , - 0.05 , vals [1 ] if vals [1 ] else "Pop{}" .format (c ),
150+ ax .text (xpos , - 0.05 , vals [1 ] if vals [1 ] else "Pop{}" .format (count ),
151151 rotation = 45 , va = "top" , ha = "right" , fontsize = 14 ,
152152 weight = "bold" )
153- c += 1
153+ count += 1
154154
155155 for axis in ["top" , "bottom" , "left" , "right" ]:
156156 ax .spines [axis ].set_linewidth (2 )
@@ -168,9 +168,9 @@ def main(result_files, fmt, outdir, popfile=None):
168168 :return:
169169 """
170170
171- for f in result_files :
172- data , pops = dataminer (f , fmt , popfile )
171+ for files in result_files :
172+ data , pops = dataminer (files , fmt , popfile )
173173 # Get output file name from input file name
174- outfile = os .path .join (outdir , f .split (os .sep )[- 1 ])
174+ outfile = os .path .join (outdir , files .split (os .sep )[- 1 ])
175175 # Create plots
176176 plotter (data , pops , outfile )
0 commit comments