Mercurial > repos > eschen42 > w4mclassfilter

annotate w4mclassfilter_wrapper.R @ 13:b24ca78a425b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
author eschen42
date Mon, 03 Sep 2018 22:28:03 -0400
parents 518cc205f289
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
1 #!/usr/bin/env Rscript
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
2
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
3 library(batch) ## parseCommandArgs
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
4
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
5 ########
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
6 # MAIN #
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
7 ########
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
8
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
9 argVc <- unlist(parseCommandArgs(evaluate=FALSE))
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
10
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
11 ##------------------------------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
12 ## Initializing
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
13 ##------------------------------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
14
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
15 ## options
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
16 ##--------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
17
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
18 strAsFacL <- options()$stringsAsFactors
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
19 options(stringsAsFactors = FALSE)
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
20
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
21 ## libraries
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
22 ##----------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
23
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
24 suppressMessages(library(w4mclassfilter))
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
25
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
26 if(packageVersion("w4mclassfilter") < "0.98.0")
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
27 stop("Please use 'w4mclassfilter' versions of 0.98.0 and above")
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
28
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
29 ## constants
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
30 ##----------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
31
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
32 modNamC <- "w4mclassfilter" ## module name
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
33
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
34 topEnvC <- environment()
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
35 flgC <- "\n"
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
36
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
37 ## functions
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
38 ##----------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
39
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
40 flgF <- function(tesC,
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
41 envC = topEnvC,
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
42 txtC = NA) { ## management of warning and error messages
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
43
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
44 tesL <- eval(parse(text = tesC), envir = envC)
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
45
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
46 if(!tesL) {
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
47
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
48 #sink(NULL)
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
49 stpTxtC <- ifelse(is.na(txtC),
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
50 paste0(tesC, " is FALSE"),
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
51 txtC)
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
52
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
53 stop(stpTxtC,
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
54 call. = FALSE)
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
55
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
56 }
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
57
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
58 } ## flgF
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
59
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
60
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
61 ## log file
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
62 ##---------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
63
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
64 my_print <- function(x, ...) { cat(c(x, ...))}
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
65
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
66 my_print("\nStart of the '", modNamC, "' Galaxy module call: ",
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
67 format(Sys.time(), "%a %d %b %Y %X"), "\n", sep="")
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
68
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
69 ## arguments
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
70 ##----------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
71
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
72 # files
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
73
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
74 dataMatrix_in <- as.character(argVc["dataMatrix_in"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
75 dataMatrix_out <- as.character(argVc["dataMatrix_out"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
76
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
77 sampleMetadata_in <- as.character(argVc["sampleMetadata_in"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
78 sampleMetadata_out <- as.character(argVc["sampleMetadata_out"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
79
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
80 variableMetadata_in <- as.character(argVc["variableMetadata_in"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
81 variableMetadata_out <- as.character(argVc["variableMetadata_out"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
82
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
83 # other parameters
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
84
12
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
85 transformation <- as.character(argVc["transformation"])
6
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
86 wildcards <- as.logical(argVc["wildcards"])
0
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
87 sampleclassNames <- as.character(argVc["sampleclassNames"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
88 sampleclassNames <- strsplit(x = sampleclassNames, split = ",", fixed = TRUE)[[1]]
5
ae791fe4fbe9 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit a2d17eac4a1343a34cf2908d9ab31b3202a21f64
eschen42
parents: 1
diff changeset
89 if (wildcards) {
ae791fe4fbe9 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit a2d17eac4a1343a34cf2908d9ab31b3202a21f64
eschen42
parents: 1
diff changeset
90 sampleclassNames <- gsub("[.]", "[.]", sampleclassNames)
ae791fe4fbe9 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit a2d17eac4a1343a34cf2908d9ab31b3202a21f64
eschen42
parents: 1
diff changeset
91 sampleclassNames <- utils::glob2rx(sampleclassNames, trim.tail = FALSE)
ae791fe4fbe9 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit a2d17eac4a1343a34cf2908d9ab31b3202a21f64
eschen42
parents: 1
diff changeset
92 }
0
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
93 inclusive <- as.logical(argVc["inclusive"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
94 classnameColumn <- as.character(argVc["classnameColumn"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
95 samplenameColumn <- as.character(argVc["samplenameColumn"])
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
96
6
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
97 variable_range_filter <- as.character(argVc["variable_range_filter"])
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
98 variable_range_filter <- strsplit(x = variable_range_filter, split = ",", fixed = TRUE)[[1]]
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
99
12
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
100 ## -----------------------------
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
101 ## Transformation and imputation
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
102 ## -----------------------------
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
103 my_w4m_filter_imputation <- if (transformation == "log10") {
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
104 function(m) {
13
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
105 if (!is.matrix(m))
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
106 stop("Cannot impute and transform data - the supplied data is not in matrix form")
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
107 if (nrow(m) == 0)
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
108 stop("Cannot impute and transform data - data matrix has no rows")
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
109 if (ncol(m) == 0)
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
110 stop("Cannot impute and transform data - data matrix has no columns")
12
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
111 suppressWarnings(
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
112 # suppress warnings here since non-positive values will produce NaN's that will be fixed in the next step
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
113 m <- log10(m)
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
114 )
13
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
115 return ( w4m_filter_imputation(m) )
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
116 }
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
117 } else if (transformation == "log2") {
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
118 function(m) {
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
119 if (!is.matrix(m))
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
120 stop("Cannot impute and transform data - the supplied data is not in matrix form")
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
121 if (nrow(m) == 0)
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
122 stop("Cannot impute and transform data - data matrix has no rows")
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
123 if (ncol(m) == 0)
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
124 stop("Cannot impute and transform data - data matrix has no columns")
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
125 suppressWarnings(
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
126 # suppress warnings here since non-positive values will produce NaN's that will be fixed in the next step
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
127 m <- log2(m)
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
128 )
b24ca78a425b planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 32cf96ed6c39a93eb9c6cf00494b128bd3c3febe
eschen42
parents: 12
diff changeset
129 return ( w4m_filter_imputation(m) )
12
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
130 }
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
131 } else {
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
132 # use the method from the w4mclassfilter class
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
133 w4m_filter_imputation
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
134 }
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
135
0
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
136 ##------------------------------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
137 ## Computation
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
138 ##------------------------------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
139
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
140 result <- w4m_filter_by_sample_class(
6
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
141 dataMatrix_in = dataMatrix_in
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
142 , sampleMetadata_in = sampleMetadata_in
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
143 , variableMetadata_in = variableMetadata_in
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
144 , dataMatrix_out = dataMatrix_out
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
145 , sampleMetadata_out = sampleMetadata_out
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
146 , variableMetadata_out = variableMetadata_out
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
147 , classes = sampleclassNames
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
148 , include = inclusive
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
149 , class_column = classnameColumn
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
150 , samplename_column = samplenameColumn
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
151 , variable_range_filter = variable_range_filter
f1eabb5973b1 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 08a878ac05ac0c05eafac27b9ebaadb88ae21eb7
eschen42
parents: 5
diff changeset
152 , failure_action = my_print
12
518cc205f289 planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit aaa20ca94614124d11723bb906dee91636144d05
eschen42
parents: 6
diff changeset
153 , data_imputation = my_w4m_filter_imputation
0
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
154 )
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
155
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
156 my_print("\nResult of '", modNamC, "' Galaxy module call to 'w4mclassfilter::w4m_filter_by_sample_class' R function: ",
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
157 as.character(result), "\n", sep = "")
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
158
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
159 ##--------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
160 ## Closing
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
161 ##--------
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
162
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
163 my_print("\nEnd of '", modNamC, "' Galaxy module call: ",
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
164 as.character(Sys.time()), "\n", sep = "")
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
165
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
166 #sink()
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
167
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
168 if (!file.exists(dataMatrix_out)) {
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
169 print(sprintf("ERROR %s::w4m_filter_by_sample_class - file '%s' was not created", modNamC, dataMatrix_out))
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
170 }# else { print(sprintf("INFO %s::w4m_filter_by_sample_class - file '%s' was exists", modNamC, dataMatrix_out)) }
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
171
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
172 if (!file.exists(variableMetadata_out)) {
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
173 print(sprintf("ERROR %s::w4m_filter_by_sample_class - file '%s' was not created", modNamC, variableMetadata_out))
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
174 } # else { print(sprintf("INFO %s::w4m_filter_by_sample_class - file '%s' was exists", modNamC, variableMetadata_out)) }
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
175
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
176 if (!file.exists(sampleMetadata_out)) {
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
177 print(sprintf("ERROR %s::w4m_filter_by_sample_class - file '%s' was not created", modNamC, sampleMetadata_out))
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
178 } # else { print(sprintf("INFO %s::w4m_filter_by_sample_class - file '%s' was exists", modNamC, sampleMetadata_out)) }
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
179
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
180 if( !result ) {
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
181 stop(sprintf("ERROR %s::w4m_filter_by_sample_class - method failed", modNamC))
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
182 }
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
183
9549cc9e8ec6 planemo upload
eschen42
parents:
diff changeset
184 rm(list = ls())