Mercurial > repos > bornea > saint_preprocessing

changeset 29:bd71998aec8d draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bornea
date Tue, 26 Apr 2016 15:40:35 -0400
parents dbd1af88f060
children e6e456d3ac14
files SAINT_preprocessing_mq_pep.py
diffstat 1 files changed, 11 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/SAINT_preprocessing_mq_pep.py	Tue Apr 26 14:42:16 2016 -0400
+++ b/SAINT_preprocessing_mq_pep.py	Tue Apr 26 15:40:35 2016 -0400
@@ -31,15 +31,22 @@
 mq_file = sys.argv[1]
 ins_path = sys.argv[8]
 names_path = str(ins_path) + r"uniprot_names.txt"
-cmd = (r"Rscript "+ str(ins_path) +"pre_process_protein_name_set.R " + str(mq_file) +
-       " " + str(names_path))
-os.system(cmd)
+fasta_db = sys.argv[3]
+
+# Uses faster names list for filtering when default db used. 
+if fasta_db == "None":
+    cmd = (r"Rscript "+ str(ins_path) +"pre_process_protein_name_set.R " + str(mq_file) +
+        " " + str(names_path))
+    os.system(cmd)
+else:
+    cmd = (r"Rscript "+ str(ins_path) +"pre_process_protein_name_set.R " + str(mq_file) +
+        " " + str(fasta_db)
+    os.system(cmd)
 
 infile = "./tukeys_output.txt" 
 # The MaxQuant "Samples Report" output.
 prey = sys.argv[2] 
 # Y or N boolean from Galaxy.
-fasta_db = sys.argv[3]
 if fasta_db == "None":
     fasta_db = str(ins_path)  + "SwissProt_HUMAN_2014_08.fasta"
 make_bait = sys.argv[6]