annotate htseq_count_render.R @ 0:10a460f61ac2 draft

planemo upload
author mingchen0919
date Tue, 06 Mar 2018 15:16:40 -0500
parents
children cc33706a1d4b
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
1 ##============ Sink warnings and errors to a file ==============
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
2 ## use the sink() function to wrap all code within it.
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
3 ##==============================================================
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
4 zz = file('warnings_and_errors.txt')
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
5 sink(zz)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
6 sink(zz, type = 'message')
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
7
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
8 #------------import libraries--------------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
9 options(stringsAsFactors = FALSE)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
10
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
11 library(getopt)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
12 library(rmarkdown)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
13 #------------------------------------------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
14
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
15
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
16 #------------get arguments into R--------------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
17 # library(dplyr)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
18 # getopt_specification_matrix(extract_short_flags('')) %>%
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
19 # write.table(file = 'spec.txt', sep = ',', row.names = FALSE, col.names = TRUE, quote = FALSE)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
20
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
21
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
22 spec_matrix = as.matrix()
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
23 opt = getopt(spec_matrix)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
24 #----------------------------------------------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
25
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
26
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
27 #-----------using passed arguments in R
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
28 # to define system environment variables---
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
29 do.call(Sys.setenv, opt[-1])
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
30 #----------------------------------------------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
31
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
32 #---------- often used variables ----------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
33 # OUTPUT_DIR: path to the output associated directory, which stores all outputs
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
34 # TOOL_DIR: path to the tool installation directory
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
35 OUTPUT_DIR = ''
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
36 TOOL_DIR = ''
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
37 RMD_NAME = ''
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
38 OUTPUT_REPORT = opt$X_o
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
39
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
40 # create the output associated directory to store all outputs
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
41 dir.create(OUTPUT_DIR, recursive = TRUE)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
42
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
43 #-----------------render Rmd--------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
44 render(paste0(TOOL_DIR, RMD_NAME, sep = '/'), output_file = OUTPUT_REPORT)
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
45 #------------------------------------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
46
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
47 #==============the end==============
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
48
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
49
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
50 ##--------end of code rendering .Rmd templates----------------
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
51 sink()
10a460f61ac2 planemo upload
mingchen0919
parents:
diff changeset
52 ##=========== End of sinking output=============================