Mercurial > repos > bornea > query_crapome
diff CRAPomeQuery.py @ 16:a5444c834e72 draft
Uploaded
author | bornea |
---|---|
date | Fri, 29 Apr 2016 15:42:11 -0400 |
parents | 0c0abdfad3b7 |
children | 378fc3676c78 |
line wrap: on
line diff
--- a/CRAPomeQuery.py Mon Apr 25 11:47:05 2016 -0400 +++ b/CRAPomeQuery.py Fri Apr 29 15:42:11 2016 -0400 @@ -5,12 +5,6 @@ @author: brentkuenzi """ ################################################################################ -## Dependencies ## -import urllib2 -import sys -import numpy -import os -################################################################################ # This program will read in a SAINT formatted 'prey.txt' file or a file # containing a single column list of uniprot accessions (e.g. "P00533" or # "EGFR_HUMAN")query the CRAPome database (v1.1), and return a file specifying @@ -24,6 +18,12 @@ # A copy of the license is included in the section entitled "GNU # Free Documentation License". ################################################################################ +## Dependencies ## +import urllib2 +import sys +import numpy +import os +################################################################################ ## REQUIRED INPUT ## # 1) crappyData: Prey.txt or single column list of Uniprot accessions crappyData = sys.argv[1] # Prey file or File with single column of accessions @@ -33,9 +33,9 @@ ################################################################################ ## Global Variables ## if species == "HUMAN": - database = str(db_path) + "Human_CRAPome_v1-1.txt" + database = "Human_CRAPome_v1-1.txt" if species == "YEAST": - database = str(db_path) + "Yeast_CRAPome_v1-1.txt" + database = "Yeast_CRAPome_v1-1.txt" ################################################################################ ## CRAPomeQuery ## class ReturnValue1(object): @@ -134,7 +134,7 @@ temp.append(i[0]) # append gene name ave = [] total = len(i[3:]) # calculate total experiments - for j in i[3:]: + for j in i[3:]: if j != '0': ave.append(int(j)) # calculate Ave.SC on only experiments with ID cnt+=1 @@ -148,9 +148,6 @@ else: temp.append(i[2]) # append accession temp.append(i[3]) - temp.append("NA") - temp.append("NA") - temp.append("NA") query.append(temp) # final query results header = ["User Input","Mapped Gene Symbol","Num of Expt. (found/total)","Ave SC","Max SC"]