Mercurial > repos > bornea > saint_bubblebeam

comparison Bubblebeam_For_SAINT_wrapper.py @ 13:e63b61738581 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bornea
date Fri, 20 Nov 2015 11:00:42 -0500
parents 9617d82034c6
children b11a6884496e
comparison
equal deleted inserted replaced
12:7b797749e8c4 13:e63b61738581
16 bub_zoom_NSAF = sys.argv[11] 16 bub_zoom_NSAF = sys.argv[11]
17 bub_zoom_SAINT =sys.argv[12] 17 bub_zoom_SAINT =sys.argv[12]
18 bub_SAINT = sys.argv[13] 18 bub_SAINT = sys.argv[13]
19 bub_NSAF = sys.argv[14] 19 bub_NSAF = sys.argv[14]
20 ins_path =sys.argv[15] 20 ins_path =sys.argv[15]
21 r_path = sys.argv[16]
22 r_path = str(r_path) + r"bin/"
21 23
22 if crapome == "None": 24 if crapome == "None":
23 crapome = "FALSE" 25 crapome = "FALSE"
24 26
25 27
26 if label == "false": 28 if label == "false":
27 label = "FALSE" 29 label = "FALSE"
28 elif label == "true": 30 elif label == "true":
29 label = "TRUE" 31 label = "TRUE"
30 32
31 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) + r" " + str(inc_file) + r" " + str(exc_file) 33 cmd = str(r_path) + 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) + r" " + str(inc_file) + r" " + str(exc_file)
32 os.system(cmd) 34 os.system(cmd)
33 time.sleep(3) 35 time.sleep(3)
34 36
35 open('./output.txt') 37 open('./output.txt')
36 os.rename('output.txt', str(output_file_name)) 38 os.rename('output.txt', str(output_file_name))