0
|
1 #######################################################################################
|
|
2 # Python-code: NSAF scoring wrapper
|
|
3 # Author: Adam L Borne
|
|
4 # Contributers: Paul A Stewart, Brent Kuenzi
|
|
5 #######################################################################################
|
|
6 # This program runs the R script that generates a nsaf scoring table. Python script simply
|
|
7 # handles arguments and interacts with Galaxy.
|
|
8 #######################################################################################
|
|
9 # Copyright (C) Adam Borne.
|
|
10 # Permission is granted to copy, distribute and/or modify this document
|
|
11 # under the terms of the GNU Free Documentation License, Version 1.3
|
|
12 # or any later version published by the Free Software Foundation;
|
|
13 # with no Invariant Sections, no Front-Cover Texts, and no Back-Cover Texts.
|
|
14 # A copy of the license is included in the section entitled "GNU
|
|
15 # Free Documentation License".
|
|
16 #######################################################################################
|
|
17 ## REQUIRED INPUT ##
|
|
18
|
|
19 # 1) listfile: SaintExpress output file.
|
|
20 # 2) preyfile: Prey file listing gene name, sequence legnth, and gene id.
|
|
21 # 3) crapfile: Crapome file can be created at http://crapome.org. (default = "None")
|
|
22 # 4) interfile: inter file from preprocessing script.
|
|
23 #######################################################################################
|
|
24 import os
|
|
25 import sys
|
|
26 import time
|
|
27
|
|
28
|
|
29 listfile = sys.argv[1]
|
|
30 preyfile = sys.argv[2]
|
|
31 crapfile = sys.argv[3]
|
|
32 interfile = sys.argv[4]
|
|
33 ins_path = sys.argv[6]
|
|
34
|
|
35
|
|
36 cmd = (r"Rscript "+ str(ins_path) + r"/nsaf_scoring.R " + str(listfile) + r" "
|
|
37 + str(preyfile) + r" " + str(crapfile) + r" " + str(inter_file))
|
|
38 os.system(cmd)
|
|
39
|
|
40 os.rename("SaintTable.txt", sys.argv[5])
|