Mercurial > repos > bornea > saint_preprocessing
comparison SAINT_preprocessing.py @ 13:febb6def95cb draft
Uploaded
author | bornea |
---|---|
date | Tue, 12 Apr 2016 14:40:42 -0400 |
parents | 3e1b66d58f94 |
children | c22c1209e65f |
comparison
equal
deleted
inserted
replaced
12:3e1b66d58f94 | 13:febb6def95cb |
---|---|
55 bait_line = baits[i] + "\t" + baits[i+1] + "\t" + T_C + "\n" | 55 bait_line = baits[i] + "\t" + baits[i+1] + "\t" + T_C + "\n" |
56 read_infile = open(infile, "r") | 56 read_infile = open(infile, "r") |
57 for input_line in read_infile: | 57 for input_line in read_infile: |
58 input_line = input_line.strip() | 58 input_line = input_line.strip() |
59 temp = input_line.split('\t') | 59 temp = input_line.split('\t') |
60 if "Quantitative Variance" in str(temp): | 60 if "Accession Number" in str(temp): |
61 if baits[i] in temp: | 61 if baits[i] in temp: |
62 number_bait = temp.index(str(baits[i])) | 62 number_bait = temp.index(str(baits[i])) |
63 number_bait = number_bait - 9 | 63 number_bait = number_bait - 9 |
64 bait_cache.append((number_bait, str(bait_line))) | 64 bait_cache.append((number_bait, str(bait_line))) |
65 # Locates the Bait names in the column names and then sets the Baits in the | 65 # Locates the Bait names in the column names and then sets the Baits in the |
210 proteins.append(prot_id) | 210 proteins.append(prot_id) |
211 elif "_YEAST" in prot_id: | 211 elif "_YEAST" in prot_id: |
212 proteins.append(prot_id) | 212 proteins.append(prot_id) |
213 elif "_MOUSE" in prot_id: | 213 elif "_MOUSE" in prot_id: |
214 proteins.append(prot_id) | 214 proteins.append(prot_id) |
215 else: | 215 else: |
216 print "Accession must be uniprot ID or gene name" | 216 print "Accession must be uniprot ID or gene name" |
217 sys.exit() | |
218 return ReturnValue2(data, proteins, header) | 217 return ReturnValue2(data, proteins, header) |
219 | 218 |
220 | 219 |
221 def make_inter(Scaffold_input): | 220 def make_inter(Scaffold_input): |
222 bait = readtab(baitfile) | 221 bait = readtab(baitfile) |