# HG changeset patch # User bornea # Date 1447173558 18000 # Node ID 333b4c487fc356517d221ae599a641927edfa3e7 Uploaded diff -r 000000000000 -r 333b4c487fc3 saint_wrapper.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/saint_wrapper.py Tue Nov 10 11:39:18 2015 -0500 @@ -0,0 +1,47 @@ +import os +import sys + +inter_file = sys.argv[1] +prey_file = sys.argv[2] +bait_file = sys.argv[3] +num_of_rep = sys.argv[4] +vc_bool = sys.argv[5] +vc_num = sys.argv[6] +go_bool = sys.argv[7] +go_file = sys.argv[8] +output_file = sys.argv[9] + +def default_run(inter_file1,prey_file1,bait_file1,output_file1,num_of_rep1): + cmd = r"SAINTexpress-spc " + r"-R" + str(num_of_rep1) + " " + str(inter_file1) + " " + str(prey_file1) + " " + str(bait_file1) + os.system(cmd) + open('list.txt') + os.rename('list.txt', str(output_file1)) + +def with_L(inter_file1,prey_file1,bait_file1,output_file1,vc_num1,num_of_rep1): + cmd = r"SAINTexpress-spc "+ r"-R" + str(num_of_rep1) + " " + r"-L" + str(vc_num1) + " " + str(inter_file1) + " " + str(prey_file1) + " " + str(bait_file1) + os.system(cmd) + open('list.txt') + os.rename('list.txt', str(output_file1)) + +def external_data_no_L(inter_file1,prey_file1,bait_file1,output_file1,go_file1,num_of_rep1): + cmd = r"SAINTexpress-spc "+ r"-R" + str(num_of_rep1) + " " + str(inter_file1) + " " + str(prey_file1) + " " + str(bait_file1) + " " + str(go_file1) + os.system(cmd) + open('list.txt') + os.rename('list.txt', str(output_file1)) + +def external_data_with_L(inter_file1,prey_file1,bait_file1,output_file1,go_file1,num_of_rep1,vc_num1): + cmd = r"SAINTexpress-spc "+ r"-R" + str(num_of_rep1) + " " + r"-L" + str(vc_num1) + " " + str(inter_file1) + " " + str(prey_file1) + " " + str(bait_file1) + " " + str(go_file1) + os.system(cmd) + open('list.txt') + os.rename('list.txt', str(output_file1)) + +if (vc_bool == "true"): + if (go_bool == "false"): + with_L(inter_file, prey_file, bait_file, output_file, vc_num, num_of_rep) + elif (go_bool == "true"): + external_data_with_L(inter_file, prey_file, bait_file, output_file, go_file, num_of_rep, vc_num) +elif (vc_bool == "false"): + if (go_bool == "false"): + default_run(inter_file, prey_file, bait_file, output_file, num_of_rep) + elif (go_bool == "true"): + external_data_no_L(inter_file, prey_file, bait_file, output_file, go_file, num_of_rep)