Mercurial > repos > bornea > saint_bubblebeam
changeset 1:248cd769eecb draft
Uploaded
author | bornea |
---|---|
date | Tue, 17 Nov 2015 10:57:19 -0500 |
parents | 81bea7a47c68 |
children | d3be2b91e8d4 |
files | Bubblebeam_For_SAINT_wrapper.py |
diffstat | 1 files changed, 47 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/Bubblebeam_For_SAINT_wrapper.py Tue Nov 17 10:57:19 2015 -0500 @@ -0,0 +1,47 @@ +import os +import sys +import time + + +list_file = sys.argv[1] +prey_file = sys.argv[2] +crapome = sys.argv[3] +color = sys.argv[4] +label = sys.argv[5] +cutoff = sys.argv[6] +mq_sc = sys.argv[7] +output_file_name = sys.argv[8] +bub_zoom_NSAF = sys.argv[9] +bub_zoom_SAINT =sys.argv[10] +bub_SAINT = sys.argv[11] +bub_NSAF = sys.argv[12] +ins_path =sys.argv[13] + +if crapome == "None": + crapome = "FALSE" + + +if label == "false": + label = "FALSE" +elif label == "true": + label = "TRUE" + +cmd = r"Rscript "+ str(ins_path) + r"/bubbles_v9_NSAF_natural_log.R " + str(list_file) + r" " + str(prey_file) + r" " + str(crapome) + r" " + str(color) + r" " + str(label) + r" " + str(cutoff) + r" " + str(mq_sc) +os.system(cmd) +time.sleep(3) + +open('./output.txt') +os.rename('output.txt', str(output_file_name)) + +open('./bubble_zoom_NSAF.png') +os.rename('bubble_zoom_NSAF.png', str(bub_zoom_NSAF)) + +open('./bubble_zoom_SAINT.png') +os.rename('bubble_zoom_SAINT.png', str(bub_zoom_SAINT)) + +open('./bubble_SAINT.png') +os.rename('bubble_SAINT.png', str(bub_SAINT)) + +open('./bubble_NSAF.png') +os.rename('bubble_NSAF.png', str(bub_NSAF)) +