changeset 0:46a582eb1ea9 draft

Uploaded
author bornea
date Wed, 11 May 2016 10:50:52 -0400
parents
children e3625c710dd3
files nsaf_scoring.py
diffstat 1 files changed, 40 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/nsaf_scoring.py	Wed May 11 10:50:52 2016 -0400
@@ -0,0 +1,40 @@
+#######################################################################################
+# Python-code: NSAF scoring wrapper
+# Author: Adam L Borne
+# Contributers: Paul A Stewart, Brent Kuenzi
+#######################################################################################
+# This program runs the R script that generates a nsaf scoring table. Python script simply 
+# handles arguments and interacts with Galaxy.
+#######################################################################################
+# Copyright (C)  Adam Borne.
+# Permission is granted to copy, distribute and/or modify this document
+# under the terms of the GNU Free Documentation License, Version 1.3
+# or any later version published by the Free Software Foundation;
+# with no Invariant Sections, no Front-Cover Texts, and no Back-Cover Texts.
+# A copy of the license is included in the section entitled "GNU
+# Free Documentation License".
+#######################################################################################
+## REQUIRED INPUT ##
+
+# 1) listfile: SaintExpress output file.
+# 2) preyfile: Prey file listing gene name, sequence legnth, and gene id.
+# 3) crapfile: Crapome file can be created at http://crapome.org. (default = "None")
+# 4) interfile: inter file from preprocessing script.
+#######################################################################################
+import os 
+import sys
+import time
+
+
+listfile = sys.argv[1]
+preyfile = sys.argv[2]
+crapfile = sys.argv[3]
+interfile = sys.argv[4]
+ins_path = sys.argv[6]
+
+
+cmd = (r"Rscript "+ str(ins_path) + r"/nsaf_scoring.R " + str(listfile) + r" "
+    + str(preyfile) + r" " + str(crapfile) + r" " + str(inter_file))
+os.system(cmd)
+
+os.rename("SaintTable.txt", sys.argv[5])