|
0
|
1 ---
|
|
1
|
2 title: 'InterProScan analysis'
|
|
|
3 output: html_document
|
|
0
|
4 ---
|
|
|
5
|
|
|
6 ```{r setup, include=FALSE, warning=FALSE, message=FALSE}
|
|
|
7 knitr::opts_chunk$set(
|
|
|
8 echo = TRUE,
|
|
|
9 error = TRUE
|
|
|
10 )
|
|
|
11 ```
|
|
|
12
|
|
|
13
|
|
1
|
14 # InterProScan job script
|
|
0
|
15
|
|
|
16 ```{bash echo=FALSE}
|
|
|
17 #---- build job script ----
|
|
|
18 cat >temp.sh <<EOL
|
|
|
19 /opt/main/InterProScan/5.28-67.0/interproscan.sh \\
|
|
|
20 -i ${X_A} \\
|
|
|
21 -f ${X_B} \\
|
|
|
22 --disable-precalc ${X_C} \\
|
|
|
23 --iprlookup ${X_D} \\
|
|
|
24 --goterms ${X_E} \\
|
|
|
25 --pathways ${X_F} \\
|
|
|
26 --tempdir ${X_D}/temp_dir \\
|
|
|
27 --output-dir ${X_D}/output_dir > /dev/null 2>&1
|
|
1
|
28
|
|
0
|
29 EOL
|
|
|
30
|
|
|
31 # remove empty input lines
|
|
1
|
32 grep -v 'None' temp.sh |\
|
|
|
33 grep -v 'NO_ARGUMENT_NO' |\
|
|
|
34 sed 's/NO_ARGUMENT_YES//g' > job-script.sh
|
|
0
|
35
|
|
|
36 # copy job-script.sh to OUTPUT_DIR
|
|
1
|
37 cp job-script.sh ${X_d}/job-script.sh
|
|
0
|
38 ```
|
|
|
39
|
|
|
40 ```{r echo=FALSE}
|
|
|
41 # display job script
|
|
|
42 job_script = paste0(opt$X_d, '/job-script.sh')
|
|
1
|
43 tags$code(tags$pre(readChar(job_script, file.info(job_script)$size )))
|
|
0
|
44 ```
|
|
|
45
|
|
|
46 ```{bash echo=FALSE}
|
|
|
47 #run job-script.sh
|
|
|
48 sh job-script.sh
|
|
|
49 ```
|
|
|
50
|
|
|
51
|
|
|
52 ```{r, 'display output directory contents', results='asis', echo=FALSE}
|
|
|
53 ## after the job is done, we list all files from the output directory.
|
|
|
54 ## full relative path to the output directory needs to be displayed.
|
|
|
55
|
|
|
56 cat('##All output files')
|
|
|
57 cat('\n\n')
|
|
1
|
58 all_files = list.files(path = OUTPUT_DIR,
|
|
0
|
59 full.names = TRUE,
|
|
|
60 recursive = TRUE)
|
|
|
61
|
|
1
|
62 for (f in sub(OUTPUT_DIR, '.', all_files) ) {
|
|
0
|
63 cat('* [', f, '](', f, ')\n')
|
|
|
64 }
|
|
|
65 cat('\n')
|
|
|
66 ``` |
