File tree Expand file tree Collapse file tree 1 file changed +6
-1
lines changed Expand file tree Collapse file tree 1 file changed +6
-1
lines changed Original file line number Diff line number Diff line change 10
10
'tag' ,
11
11
help = "Variant tag parameter used to run RegTools." ,
12
12
)
13
+ input_parser .add_argument (
14
+ 'cohort_directory' ,
15
+ help = "Path to directory containing cohort files." ,
16
+ )
13
17
14
18
args = input_parser .parse_args ()
15
19
16
20
tag = args .tag
21
+ cohort_dir = args .cohort_directory
17
22
18
23
lines_per_file = 50000
19
24
smallfile = None
30
35
#get chunks
31
36
files = glob .glob ('small_file_*' )
32
37
for file in files :
33
- subprocess .run (f'Rscript --vanilla /home/ec2-user/workspace/regtools/scripts/compare_junctions_hist_v2.R { tag } { file } ' )
38
+ subprocess .run (f'Rscript --vanilla /home/ec2-user/workspace/regtools/scripts/compare_junctions_hist_v2.R { tag } { cohort_dir } / { file } ' )
34
39
output_files = glob .glob ("*_out.tsv" )
35
40
output_files .sort () # glob lacks reliable ordering, so impose your own if output order matters
36
41
with open (f'junction_pvalues_{ tag } .tsv' , 'wb' ) as outfile :
You can’t perform that action at this time.
0 commit comments