annotate tests/run-tests.py @ 13:df1fb577db0d

Version 1.2: Add strand bias column, rename minor allele frequency column.
author me <nmapsy@gmail.com>
date Tue, 03 Dec 2013 13:45:16 -0500
parents 97b772e3a0f1
children 44c3abd1b767
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
1 #!/usr/bin/env python
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
2 import os
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
3 import sys
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
4 import subprocess
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
5
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
6 DATASETS = [
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
7 'artificial',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
8 'artificial-samples',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
9 'artificial-nofilt',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
10 'real',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
11 'real-mit',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
12 'real-mit-s',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
13 'real-nofilt',
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
14 ]
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
15 IN_EXT = '.vcf.in'
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
16 OUT_EXT = '.csv.out'
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
17 ARGS_KEY = '##comment="ARGS='
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
18
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
19 def main():
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
20
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
21 test_dir = os.path.dirname(os.path.relpath(sys.argv[0]))
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
22 if test_dir:
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
23 test_dir += os.sep
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
24
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
25 for dataset in DATASETS:
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
26 infile = test_dir+dataset+IN_EXT
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
27 outfile = test_dir+dataset+OUT_EXT
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
28
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
29 if not os.path.exists(infile):
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
30 sys.stderr.write("Error: file not found: "+infile+"\n")
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
31 continue
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
32 if not os.path.exists(outfile):
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
33 sys.stderr.write("Error: file not found: "+outfile+"\n")
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
34 continue
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
35
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
36 options = read_options(infile)
12
97b772e3a0f1 Uploaded tarball.
nick
parents: 11
diff changeset
37 script_cmd = 'allele-counts.py '+options+' -i '+infile
97b772e3a0f1 Uploaded tarball.
nick
parents: 11
diff changeset
38 bash_cmd = 'diff '+outfile+' <('+script_cmd+')'
97b772e3a0f1 Uploaded tarball.
nick
parents: 11
diff changeset
39 # print infile+":"
97b772e3a0f1 Uploaded tarball.
nick
parents: 11
diff changeset
40 print script_cmd
97b772e3a0f1 Uploaded tarball.
nick
parents: 11
diff changeset
41 subprocess.call(['bash', '-c', bash_cmd])
10
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
42
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
43
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
44 def read_options(infile):
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
45 with open(infile, 'r') as infilehandle:
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
46 for line in infilehandle:
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
47 line.strip()
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
48 if ARGS_KEY == line[:len(ARGS_KEY)]:
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
49 return line[len(ARGS_KEY):-2]
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
50 return ''
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
51
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
52
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
53 if __name__ == '__main__':
db6f217dc45a Uploaded tarball.
nick
parents:
diff changeset
54 main()