We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 6a629a2 commit dcb8465Copy full SHA for dcb8465
scripts/stats_wrapper.py
@@ -14,6 +14,7 @@
14
args = input_parser.parse_args()
15
16
tag = args.tag
17
+cwd = os.getcwd()
18
19
lines_per_file = 50000
20
smallfile = None
@@ -30,7 +31,7 @@
30
31
#get chunks
32
files = glob.glob('small_file_*')
33
for file in files:
- subprocess.run(f'Rscript --vanilla /home/ec2-user/workspace/regtools/scripts/compare_junctions_hist_v2.R {tag} {file}')
34
+ subprocess.run(f'Rscript --vanilla /home/ec2-user/workspace/regtools/scripts/compare_junctions_hist_v2.R {tag} {cwd}/{file}')
35
output_files = glob.glob("*_out.tsv")
36
output_files.sort() # glob lacks reliable ordering, so impose your own if output order matters
37
with open(f'junction_pvalues_{tag}.tsv', 'wb') as outfile:
0 commit comments