Mercurial > repos > bornea > saint_preprocessing

comparison mzID_process2.py @ 71:f6d916d1d304 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bornea
date Sat, 27 Aug 2016 23:22:20 -0400
parents 71e47a3e1bf5
children 5ec0b997fb13
comparison
equal deleted inserted replaced
70:71e47a3e1bf5 71:f6d916d1d304
182 for i in file_list: 182 for i in file_list:
183 cmd = (r"Rscript "+ str(ins_path) +"flatten_mzIdentML.R " + i) 183 cmd = (r"Rscript "+ str(ins_path) +"flatten_mzIdentML.R " + i)
184 os.system(cmd) 184 os.system(cmd)
185 mzIdentML = read_tab("flat_mzIdentML.txt") 185 mzIdentML = read_tab("flat_mzIdentML.txt")
186 inter = inter + make_inter(mzIdentML,bait[cnt][0],bait[cnt][1]).inter 186 inter = inter + make_inter(mzIdentML,bait[cnt][0],bait[cnt][1]).inter
187 print inter
187 accessions.append(make_inter(mzIdentML,bait[cnt][0],bait[cnt][1]).accessions) 188 accessions.append(make_inter(mzIdentML,bait[cnt][0],bait[cnt][1]).accessions)
189 print accessions
188 cnt+=1 190 cnt+=1
189 191
190 with open("inter.txt","w") as x: 192 with open("inter.txt","w") as x:
191 x.write(inter) 193 x.write(inter)
192 if make_prey == "Y": 194 if make_prey == "Y":