Mercurial > repos > bornea > saint_bubblebeam
comparison Bubblebeam_For_SAINT_wrapper.py @ 1:248cd769eecb draft
Uploaded
| author | bornea |
|---|---|
| date | Tue, 17 Nov 2015 10:57:19 -0500 |
| parents | |
| children | 9617d82034c6 |
comparison
equal
deleted
inserted
replaced
| 0:81bea7a47c68 | 1:248cd769eecb |
|---|---|
| 1 import os | |
| 2 import sys | |
| 3 import time | |
| 4 | |
| 5 | |
| 6 list_file = sys.argv[1] | |
| 7 prey_file = sys.argv[2] | |
| 8 crapome = sys.argv[3] | |
| 9 color = sys.argv[4] | |
| 10 label = sys.argv[5] | |
| 11 cutoff = sys.argv[6] | |
| 12 mq_sc = sys.argv[7] | |
| 13 output_file_name = sys.argv[8] | |
| 14 bub_zoom_NSAF = sys.argv[9] | |
| 15 bub_zoom_SAINT =sys.argv[10] | |
| 16 bub_SAINT = sys.argv[11] | |
| 17 bub_NSAF = sys.argv[12] | |
| 18 ins_path =sys.argv[13] | |
| 19 | |
| 20 if crapome == "None": | |
| 21 crapome = "FALSE" | |
| 22 | |
| 23 | |
| 24 if label == "false": | |
| 25 label = "FALSE" | |
| 26 elif label == "true": | |
| 27 label = "TRUE" | |
| 28 | |
| 29 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) | |
| 30 os.system(cmd) | |
| 31 time.sleep(3) | |
| 32 | |
| 33 open('./output.txt') | |
| 34 os.rename('output.txt', str(output_file_name)) | |
| 35 | |
| 36 open('./bubble_zoom_NSAF.png') | |
| 37 os.rename('bubble_zoom_NSAF.png', str(bub_zoom_NSAF)) | |
| 38 | |
| 39 open('./bubble_zoom_SAINT.png') | |
| 40 os.rename('bubble_zoom_SAINT.png', str(bub_zoom_SAINT)) | |
| 41 | |
| 42 open('./bubble_SAINT.png') | |
| 43 os.rename('bubble_SAINT.png', str(bub_SAINT)) | |
| 44 | |
| 45 open('./bubble_NSAF.png') | |
| 46 os.rename('bubble_NSAF.png', str(bub_NSAF)) | |
| 47 |
