annotate BC/batch_correction_docker_wrapper.R @ 3:2e3a23dd6c24 draft default tip

Uploaded
author melpetera
date Thu, 28 Feb 2019 05:12:34 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
1 #!/usr/bin/env Rscript
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
2
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
3 ################################################################################################
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
4 # batch_correction_main_wrapper #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
5 # #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
6 # Author: Nils Paulhe #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
7 # User: Galaxy #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
8 # Original data: -- #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
9 # Starting date: 2017-12-11 #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
10 # Version 1: 2017-12-11 #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
11 # #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
12 # #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
13 # #
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
14 ################################################################################################
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
15
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
16 library(batch) #necessary for parseCommandArgs function
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
17
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
18 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
19 ## init. prog. constants
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
20 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
21
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
22 argv.wrapper <- commandArgs(trailingOnly = FALSE)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
23 script.path <- sub("--file=", "", argv.wrapper[grep("--file=", argv.wrapper)])
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
24 prog.name <- basename(script.path)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
25
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
26 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
27 ## init. functions
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
28 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
29
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
30 script_bypass <- function(other.script.name) {
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
31 initial.options <- commandArgs(trailingOnly = FALSE)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
32 file.arg.name <- "--file="
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
33 script.name <- sub(file.arg.name, "", initial.options[grep(file.arg.name, initial.options)])
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
34 script.basename <- dirname(script.name)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
35 other.script.fullpath <- paste(sep="/", script.basename, other.script.name)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
36 other.script.cmd <- paste(sep=" ", "Rscript", other.script.fullpath, "-h")
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
37 system(other.script.cmd, wait=TRUE)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
38 }
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
39
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
40 source_wrapper <- function(other.script.name){
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
41 initial.options <- commandArgs(trailingOnly = FALSE)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
42 file.arg.name <- "--file="
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
43 script.name <- sub(file.arg.name, "", initial.options[grep(file.arg.name, initial.options)])
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
44 script.basename <- dirname(script.name)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
45 other.script.fullpath <- paste(sep="/", script.basename, other.script.name)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
46 source(other.script.fullpath)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
47 }
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
48
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
49 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
50 ## Test Help
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
51 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
52
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
53 if (length(grep('-h', argv.wrapper)) > 0) {
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
54 cat("Usage: Rscript ",
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
55 prog.name,
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
56 "{args} \n",
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
57 "parameters: \n",
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
58 "\t-h: display this help message, call all scripts with the same option and exit (optional) \n",
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
59 "\t--loess \"TRUE\": call the script as \"batch_correction_all_loess_wrapper.R\"; otherwise call it as \"batch_correction_wrapper.R\" one (optional) \n",
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
60 "for other parameters, please refer to each script specific options and parameters. \n",
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
61 "\n")
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
62 script_bypass("batch_correction_all_loess_wrapper.R")
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
63 script_bypass("batch_correction_wrapper.R")
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
64 quit(status = 0)
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
65 }
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
66
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
67 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
68 ## check if loess or normal
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
69 ##------------------------------
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
70
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
71 if (length(grep('--loess', argv.wrapper)) > 0) {
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
72 source_wrapper("batch_correction_all_loess_wrapper.R")
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
73 } else {
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
74 source_wrapper("batch_correction_wrapper.R")
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
75 }
2e3a23dd6c24 Uploaded
melpetera
parents:
diff changeset
76