Mercurial > repos > eschen42 > w4mclassfilter
annotate w4mclassfilter_wrapper.R @ 15:08d4ca8bc6dd draft
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
author | eschen42 |
---|---|
date | Thu, 11 Mar 2021 20:46:26 +0000 |
parents | |
children |
rev | line source |
---|---|
15
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
1 #!/usr/bin/env Rscript |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
2 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
3 library(batch) ## parseCommandArgs |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
4 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
5 ######## |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
6 # MAIN # |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
7 ######## |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
8 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
9 argVc <- unlist(parseCommandArgs(evaluate=FALSE)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
10 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
11 ##------------------------------ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
12 ## Initializing |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
13 ##------------------------------ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
14 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
15 ## options |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
16 ##-------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
17 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
18 strAsFacL <- options()$stringsAsFactors |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
19 options(stringsAsFactors = FALSE) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
20 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
21 ## libraries |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
22 ##---------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
23 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
24 suppressMessages(library(w4mclassfilter)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
25 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
26 expected_version <- "0.98.18" |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
27 actual_version <- packageVersion("w4mclassfilter") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
28 if(packageVersion("w4mclassfilter") < expected_version) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
29 stop( |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
30 sprintf( |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
31 "Unrecoverable error: Version %s of the 'w4mclassfilter' R package was loaded instead of expected version %s", |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
32 actual_version, expected_version |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
33 ) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
34 ) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
35 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
36 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
37 ## constants |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
38 ##---------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
39 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
40 modNamC <- "w4mclassfilter" ## module name |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
41 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
42 topEnvC <- environment() |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
43 flgC <- "\n" |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
44 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
45 ## functions |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
46 ##---------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
47 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
48 flgF <- function(tesC, |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
49 envC = topEnvC, |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
50 txtC = NA) { ## management of warning and error messages |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
51 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
52 tesL <- eval(parse(text = tesC), envir = envC) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
53 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
54 if(!tesL) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
55 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
56 #sink(NULL) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
57 stpTxtC <- ifelse(is.na(txtC), |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
58 paste0(tesC, " is FALSE"), |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
59 txtC) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
60 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
61 stop(stpTxtC, |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
62 call. = FALSE) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
63 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
64 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
65 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
66 } ## flgF |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
67 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
68 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
69 ## log file |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
70 ##--------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
71 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
72 my_print <- function(x, ...) { cat(c(x, ...))} |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
73 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
74 my_print("\nStart of the '", modNamC, "' Galaxy module call: ", |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
75 format(Sys.time(), "%a %d %b %Y %X"), "\n", sep="") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
76 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
77 ## arguments |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
78 ##---------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
79 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
80 # files |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
81 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
82 dataMatrix_in <- as.character(argVc["dataMatrix_in"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
83 dataMatrix_out <- as.character(argVc["dataMatrix_out"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
84 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
85 sampleMetadata_in <- as.character(argVc["sampleMetadata_in"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
86 sampleMetadata_out <- as.character(argVc["sampleMetadata_out"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
87 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
88 variableMetadata_in <- as.character(argVc["variableMetadata_in"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
89 variableMetadata_out <- as.character(argVc["variableMetadata_out"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
90 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
91 # other parameters |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
92 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
93 transformation <- as.character(argVc["transformation"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
94 my_imputation_label <- as.character(argVc["imputation"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
95 my_imputation_function <- if (my_imputation_label == "zero") { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
96 w4m_filter_zero_imputation |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
97 } else if (my_imputation_label == "center") { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
98 w4m_filter_median_imputation |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
99 } else if (my_imputation_label == "none") { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
100 w4m_filter_no_imputation |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
101 } else { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
102 stop(sprintf("Unknown value %s supplied for 'imputation' parameter. Expected one of {zero,center,none}.")) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
103 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
104 wildcards <- as.logical(argVc["wildcards"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
105 sampleclassNames <- as.character(argVc["sampleclassNames"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
106 sampleclassNames <- strsplit(x = sampleclassNames, split = ",", fixed = TRUE)[[1]] |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
107 if (wildcards) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
108 sampleclassNames <- gsub("[.]", "[.]", sampleclassNames) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
109 sampleclassNames <- utils::glob2rx(sampleclassNames, trim.tail = FALSE) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
110 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
111 inclusive <- as.logical(argVc["inclusive"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
112 classnameColumn <- as.character(argVc["classnameColumn"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
113 samplenameColumn <- as.character(argVc["samplenameColumn"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
114 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
115 order_vrbl <- as.character(argVc["order_vrbl"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
116 centering <- as.character(argVc["centering"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
117 order_smpl <- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
118 if (centering == 'centroid' || centering == 'median') { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
119 "sampleMetadata" |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
120 } else { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
121 as.character(argVc["order_smpl"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
122 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
123 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
124 variable_range_filter <- as.character(argVc["variable_range_filter"]) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
125 variable_range_filter <- strsplit(x = variable_range_filter, split = ",", fixed = TRUE)[[1]] |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
126 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
127 ## ----------------------------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
128 ## Transformation and imputation |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
129 ## ----------------------------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
130 my_transformation_and_imputation <- if (transformation == "log10") { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
131 function(m) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
132 # convert negative intensities to missing values |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
133 m[m < 0] <- NA |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
134 if (!is.matrix(m)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
135 stop("Cannot transform and impute data - the supplied data is not in matrix form") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
136 if (nrow(m) == 0) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
137 stop("Cannot transform and impute data - data matrix has no rows") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
138 if (ncol(m) == 0) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
139 stop("Cannot transform and impute data - data matrix has no columns") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
140 suppressWarnings({ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
141 # suppress warnings here since non-positive values will produce NaN's that will be fixed in the next step |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
142 m <- log10(m) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
143 m[is.na(m)] <- NA |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
144 }) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
145 return ( my_imputation_function(m) ) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
146 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
147 } else if (transformation == "log2") { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
148 function(m) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
149 # convert negative intensities to missing values |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
150 m[m < 0] <- NA |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
151 if (!is.matrix(m)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
152 stop("Cannot transform and impute data - the supplied data is not in matrix form") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
153 if (nrow(m) == 0) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
154 stop("Cannot transform and impute data - data matrix has no rows") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
155 if (ncol(m) == 0) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
156 stop("Cannot transform and impute data - data matrix has no columns") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
157 suppressWarnings({ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
158 # suppress warnings here since non-positive values will produce NaN's that will be fixed in the next step |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
159 m <- log2(m) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
160 m[is.na(m)] <- NA |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
161 }) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
162 return ( my_imputation_function(m) ) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
163 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
164 } else { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
165 function(m) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
166 # convert negative intensities to missing values |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
167 m[m < 0] <- NA |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
168 if (!is.matrix(m)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
169 stop("Cannot transform and impute data - the supplied data is not in matrix form") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
170 if (nrow(m) == 0) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
171 stop("Cannot transform and impute data - data matrix has no rows") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
172 if (ncol(m) == 0) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
173 stop("Cannot transform and impute data - data matrix has no columns") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
174 suppressWarnings({ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
175 # suppress warnings here since non-positive values will produce NaN's that will be fixed in the next step |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
176 m[is.na(m)] <- NA |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
177 }) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
178 return ( my_imputation_function(m) ) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
179 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
180 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
181 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
182 ##------------------------------ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
183 ## Computation |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
184 ##------------------------------ |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
185 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
186 result <- w4m_filter_by_sample_class( |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
187 dataMatrix_in = dataMatrix_in |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
188 , sampleMetadata_in = sampleMetadata_in |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
189 , variableMetadata_in = variableMetadata_in |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
190 , dataMatrix_out = dataMatrix_out |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
191 , sampleMetadata_out = sampleMetadata_out |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
192 , variableMetadata_out = variableMetadata_out |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
193 , classes = sampleclassNames |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
194 , include = inclusive |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
195 , class_column = classnameColumn |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
196 , samplename_column = samplenameColumn |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
197 , order_vrbl = order_vrbl |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
198 , order_smpl = order_smpl |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
199 , centering = centering |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
200 , variable_range_filter = variable_range_filter |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
201 , failure_action = my_print |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
202 , data_imputation = my_transformation_and_imputation |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
203 ) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
204 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
205 my_print("\nResult of '", modNamC, "' Galaxy module call to 'w4mclassfilter::w4m_filter_by_sample_class' R function: ", |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
206 as.character(result), "\n", sep = "") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
207 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
208 ##-------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
209 ## Closing |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
210 ##-------- |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
211 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
212 my_print("\nEnd of '", modNamC, "' Galaxy module call: ", |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
213 as.character(Sys.time()), "\n", sep = "") |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
214 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
215 #sink() |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
216 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
217 if (!file.exists(dataMatrix_out)) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
218 print(sprintf("ERROR %s::w4m_filter_by_sample_class - file '%s' was not created", modNamC, dataMatrix_out)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
219 }# else { print(sprintf("INFO %s::w4m_filter_by_sample_class - file '%s' was exists", modNamC, dataMatrix_out)) } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
220 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
221 if (!file.exists(variableMetadata_out)) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
222 print(sprintf("ERROR %s::w4m_filter_by_sample_class - file '%s' was not created", modNamC, variableMetadata_out)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
223 } # else { print(sprintf("INFO %s::w4m_filter_by_sample_class - file '%s' was exists", modNamC, variableMetadata_out)) } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
224 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
225 if (!file.exists(sampleMetadata_out)) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
226 print(sprintf("ERROR %s::w4m_filter_by_sample_class - file '%s' was not created", modNamC, sampleMetadata_out)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
227 } # else { print(sprintf("INFO %s::w4m_filter_by_sample_class - file '%s' was exists", modNamC, sampleMetadata_out)) } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
228 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
229 if( !result ) { |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
230 stop(sprintf("ERROR %s::w4m_filter_by_sample_class - method failed", modNamC)) |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
231 } |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
232 |
08d4ca8bc6dd
"planemo upload for repository https://github.com/HegemanLab/w4mclassfilter_galaxy_wrapper/tree/master commit 9639dde5737c9aa2330bb603c2299345939407cf"
eschen42
parents:
diff
changeset
|
233 rm(list = ls()) |