|
6
|
1 import os
|
|
|
2 import sys
|
|
|
3 import subprocess
|
|
|
4 import shutil
|
|
|
5
|
|
|
6 def main():
|
|
|
7
|
|
|
8 input_file=sys.argv[1]
|
|
|
9 tmp_dir=sys.argv[4]
|
|
|
10 script_dir=os.path.dirname(os.path.abspath(__file__))
|
|
|
11 plot=sys.argv[11]
|
|
|
12 pdffigures=sys.argv[13]
|
|
|
13 outputlog=sys.argv[14]
|
|
|
14 log=sys.argv[15]
|
|
|
15 user=sys.argv[16]
|
|
|
16 package=sys.argv[17]
|
|
|
17
|
|
|
18 iFile=open(input_file,'r')
|
|
|
19 dataSetLine=iFile.readline()
|
|
|
20 dataset=dataSetLine.split("\t")[1]
|
|
|
21 iFile.close()
|
|
|
22
|
|
|
23 if (outputlog=="TRUE"):
|
|
|
24 errfile=open(log,'w')
|
|
|
25 else:
|
|
|
26 errfile=open(os.path.join(tmp_dir,"errfile.log"),'w')
|
|
|
27
|
|
|
28 retcode=subprocess.call(["Rscript", os.path.join(script_dir,"selection.R"), dataset, sys.argv[2], sys.argv[3], sys.argv[4], sys.argv[5], sys.argv[6], sys.argv[7], sys.argv[8], sys.argv[9], sys.argv[10], sys.argv[11], sys.argv[12],sys.argv[16],package], stdout = errfile, stderr = errfile)
|
|
|
29
|
|
|
30 if (plot=="TRUE"):
|
|
|
31 shutil.copy(os.path.join(tmp_dir,"mpagenomics",user,"Rplots.pdf"), pdffigures)
|
|
|
32
|
|
|
33 errfile.close()
|
|
|
34
|
|
|
35 sys.exit(retcode)
|
|
|
36
|
|
|
37 if __name__ == "__main__":
|
|
|
38 main()
|
