Mercurial > repos > lecorguille > camera_annotate

annotate lib.r @ 43:7f4ae4db39eb draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
author workflow4metabolomics
date Thu, 22 Apr 2021 10:26:15 +0000
parents 2350ab2c7e13
children db77f2b67a7f
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
1 # lib.r
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
2
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
3 #@author G. Le Corguille
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
4 # solve an issue with batch if arguments are logical TRUE/FALSE
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
5 parseCommandArgs <- function(...) {
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
6 args <- batch::parseCommandArgs(...)
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
7 for (key in names(args)) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
8 if (args[key] %in% c("TRUE", "FALSE"))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
9 args[key] <- as.logical(args[key])
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
10 }
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
11 return(args)
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
12 }
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
13
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
14 #@author G. Le Corguille
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
15 # This function will
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
16 # - load the packages
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
17 # - display the sessionInfo
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
18 loadAndDisplayPackages <- function(pkgs) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
19 for (pkg in pkgs) suppressPackageStartupMessages(stopifnot(library(pkg, quietly = TRUE, logical.return = TRUE, character.only = TRUE)))
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
20
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
21 sessioninfo <- sessionInfo()
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
22 cat(sessioninfo$R.version$version.string, "\n")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
23 cat("Main packages:\n")
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
24 for (pkg in names(sessioninfo$otherPkgs)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
25 cat(paste(pkg, packageVersion(pkg)), "\t")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
26 }
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
27 cat("\n")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
28 cat("Other loaded packages:\n")
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
29 for (pkg in names(sessioninfo$loadedOnly)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
30 cat(paste(pkg, packageVersion(pkg)), "\t")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
31 }
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
32 cat("\n")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
33 }
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
34
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
35 # This function retrieve a xset like object
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
36 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
37 getxcmsSetObject <- function(xobject) {
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
38 # XCMS 1.x
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
39 if (class(xobject) == "xcmsSet")
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
40 return(xobject)
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
41 # XCMS 3.x
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
42 if (class(xobject) == "XCMSnExp") {
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
43 # Get the legacy xcmsSet object
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
44 suppressWarnings(xset <- as(xobject, "xcmsSet"))
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
45 if (is.null(xset@phenoData$sample_group))
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
46 sampclass(xset) <- "."
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
47 else
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
48 sampclass(xset) <- xset@phenoData$sample_group
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
49 if (!is.null(xset@phenoData$sample_name))
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
50 rownames(xset@phenoData) <- xset@phenoData$sample_name
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
51 return(xset)
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
52 }
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
53 }
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
54
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
55 #@author G. Le Corguille
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
56 #The function create a pdf from the different png generated by diffreport
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
57 diffreport_png2pdf <- function(filebase) {
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
58 dir.create("pdf")
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
59
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
60 pdfEicOutput <- paste0("pdf/", filebase, "-eic_pdf.pdf")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
61 pdfBoxOutput <- paste0("pdf/", filebase, "-box_pdf.pdf")
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
62
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
63 system(paste0("gm convert ", filebase, "_eic/*.png ", pdfEicOutput))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
64 system(paste0("gm convert ", filebase, "_box/*.png ", pdfBoxOutput))
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
65
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
66 }
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
67
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
68 #@author G. Le Corguille
36
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
69 #The function create a zip archive from the different png generated by diffreport
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
70 diffreport_png2zip <- function() {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
71 zip("eic.zip", dir(pattern = "_eic"), zip = Sys.which("zip"))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
72 zip("box.zip", dir(pattern = "_box"), zip = Sys.which("zip"))
36
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
73 }
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
74
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
75 #The function create a zip archive from the different tabular generated by diffreport
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
76 diffreport_tabular2zip <- function() {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
77 zip("tabular.zip", dir(pattern = "tabular/*"), zip = Sys.which("zip"))
36
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
78 }
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
79
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
80 #@author G. Le Corguille
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
81 #This function convert if it is required the Retention Time in minutes
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
82 RTSecondToMinute <- function(variableMetadata, convertRTMinute) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
83 if (convertRTMinute) {
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
84 #converting the retention times (seconds) into minutes
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
85 print("converting the retention times into minutes in the variableMetadata")
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
86 variableMetadata[, "rt"] <- variableMetadata[, "rt"] / 60
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
87 variableMetadata[, "rtmin"] <- variableMetadata[, "rtmin"] / 60
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
88 variableMetadata[, "rtmax"] <- variableMetadata[, "rtmax"] / 60
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
89 }
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
90 return(variableMetadata)
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
91 }
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
92
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
93 #@author G. Le Corguille
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
94 #This function format ions identifiers
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
95 formatIonIdentifiers <- function(variableMetadata, numDigitsRT = 0, numDigitsMZ = 0) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
96 splitDeco <- strsplit(as.character(variableMetadata$name), "_")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
97 idsDeco <- sapply(splitDeco, function(x) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
98 deco <- unlist(x)[2]; if (is.na(deco)) return("") else return(paste0("_", deco))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
99 })
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
100 namecustom <- make.unique(paste0("M", round(variableMetadata[, "mz"], numDigitsMZ), "T", round(variableMetadata[, "rt"], numDigitsRT), idsDeco))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
101 variableMetadata <- cbind(name = variableMetadata$name, namecustom = namecustom, variableMetadata[, !(colnames(variableMetadata) %in% c("name"))])
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
102 return(variableMetadata)
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
103 }
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
104
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
105 #The function annotateDiffreport without the corr function which bugs
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
106 annotatediff <- function(xset = xset, args = args, variableMetadataOutput = "variableMetadata.tsv") {
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
107 # Resolve the bug with x11, with the function png
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
108 options(bitmapType = "cairo")
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
109
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
110 #Check if the fillpeaks step has been done previously, if it hasn't, there is an error message and the execution is stopped.
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
111 res <- try(is.null(xset@filled))
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
112
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
113 # ------ annot -------
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
114 args$calcCiS <- as.logical(args$calcCiS)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
115 args$calcIso <- as.logical(args$calcIso)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
116 args$calcCaS <- as.logical(args$calcCaS)
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
117
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
118 # common parameters
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
119 args4annotate <- list(object = xset,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
120 nSlaves = args$nSlaves, sigma = args$sigma, perfwhm = args$perfwhm,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
121 maxcharge = args$maxcharge, maxiso = args$maxiso, minfrac = args$minfrac,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
122 ppm = args$ppm, mzabs = args$mzabs, quick = args$quick,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
123 polarity = args$polarity, max_peaks = args$max_peaks, intval = args$intval)
22
65e3e1a39e16 planemo upload commit 23ad0382ec59d66220fd94684a6d158fdef2c7f5-dirty
lecorguille
parents: 21
diff changeset
124
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
125 if (args$quick == FALSE) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
126 args4annotate <- append(args4annotate,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
127 list(graphMethod = args$graphMethod, cor_eic_th = args$cor_eic_th, pval = args$pval,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
128 calcCiS = args$calcCiS, calcIso = args$calcIso, calcCaS = args$calcCaS))
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
129 # no ruleset
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
130 if (!is.null(args$multiplier)) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
131 args4annotate <- append(args4annotate,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
132 list(multiplier = args$multiplier))
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
133 }
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
134 # ruleset
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
135 else {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
136 rulset <- read.table(args$rules, h = T, sep = ";")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
137 if (ncol(rulset) < 4) rulset <- read.table(args$rules, h = T, sep = "\t")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
138 if (ncol(rulset) < 4) rulset <- read.table(args$rules, h = T, sep = ",")
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
139 if (ncol(rulset) < 4) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
140 error_message <- "Your ruleset file seems not well formatted. The column separators accepted are ; , and tabulation"
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
141 print(error_message)
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
142 stop(error_message)
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
143 }
23
2a4a811c663d planemo upload commit bf9996f1e380af6f783cc862f3e1d69f8c7b86e5
lecorguille
parents: 22
diff changeset
144
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
145 args4annotate <- append(args4annotate,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
146 list(rules = rulset))
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
147 }
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
148 }
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
149
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
150
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
151 # launch annotate
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
152 xa <- do.call("annotate", args4annotate)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
153 peakList <- getPeaklist(xa, intval = args$intval)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
154 peakList <- cbind(groupnames(xa@xcmsSet), peakList); colnames(peakList)[1] <- c("name");
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
155
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
156 # --- Multi condition : diffreport ---
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
157 diffrepOri <- NULL
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
158 if (!is.null(args$runDiffreport) & nlevels(sampclass(xset)) >= 2) {
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
159 #Check if the fillpeaks step has been done previously, if it hasn't, there is an error message and the execution is stopped.
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
160 res <- try(is.null(xset@filled))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
161 classes <- levels(sampclass(xset))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
162 x <- 1:(length(classes) - 1)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
163 for (i in seq(along = x)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
164 y <- 1:(length(classes))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
165 for (n in seq(along = y)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
166 if (i + n <= length(classes)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
167 filebase <- paste(classes[i], class2 = classes[i + n], sep = "-vs-")
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
168
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
169 diffrep <- diffreport(
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
170 object = xset, class1 = classes[i], class2 = classes[i + n],
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
171 filebase = filebase, eicmax = args$eicmax, eicwidth = args$eicwidth,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
172 sortpval = TRUE, value = args$value, h = args$h, w = args$w, mzdec = args$mzdec, missing = 0)
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
173
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
174 diffrepOri <- diffrep
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
175
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
176 # renamming of the column rtmed to rt to fit with camera peaklist function output
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
177 colnames(diffrep)[colnames(diffrep) == "rtmed"] <- "rt"
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
178 colnames(diffrep)[colnames(diffrep) == "mzmed"] <- "mz"
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
179
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
180 # combines results and reorder columns
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
181 diffrep <- merge(peakList, diffrep[, c("name", "fold", "tstat", "pvalue")], by.x = "name", by.y = "name", sort = F)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
182 diffrep <- cbind(diffrep[, !(colnames(diffrep) %in% c(sampnames(xa@xcmsSet)))], diffrep[, (colnames(diffrep) %in% c(sampnames(xa@xcmsSet)))])
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
183
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
184 diffrep <- RTSecondToMinute(diffrep, args$convertRTMinute)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
185 diffrep <- formatIonIdentifiers(diffrep, numDigitsRT = args$numDigitsRT, numDigitsMZ = args$numDigitsMZ)
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
186
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
187 if (args$sortpval) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
188 diffrep <- diffrep[order(diffrep$pvalue), ]
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
189 }
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
190
37
636e36a64d31 planemo upload commit a5f400649acb64e1c66a92379027ac5abdb4f546
lecorguille
parents: 36
diff changeset
191 dir.create("tabular", showWarnings = FALSE)
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
192 write.table(diffrep, sep = "\t", quote = FALSE, row.names = FALSE, file = paste("tabular/", filebase, "_tsv.tabular", sep = ""))
35
b02797fbead1 planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents: 33
diff changeset
193
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
194 if (args$eicmax != 0) {
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
195 if (args$png2 == "pdf")
36
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
196 diffreport_png2pdf(filebase)
41
2350ab2c7e13 "planemo upload commit f978cc0a7aa81c244591e560f44970f0d61bf74e"
workflow4metabolomics
parents: 40
diff changeset
197 if (args$png2 == "zip")
2350ab2c7e13 "planemo upload commit f978cc0a7aa81c244591e560f44970f0d61bf74e"
workflow4metabolomics
parents: 40
diff changeset
198 diffreport_png2zip()
33
3b2b642542f6 planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents: 30
diff changeset
199 }
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
200 }
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
201 }
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
202 }
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
203 if (args$tabular2 == "zip")
36
b6a8143c397c planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents: 35
diff changeset
204 diffreport_tabular2zip()
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
205 }
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
206
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
207 # --- variableMetadata ---
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
208 variableMetadata <- peakList[, !(make.names(colnames(peakList)) %in% c(make.names(sampnames(xa@xcmsSet))))]
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
209 variableMetadata <- RTSecondToMinute(variableMetadata, args$convertRTMinute)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
210 variableMetadata <- formatIonIdentifiers(variableMetadata, numDigitsRT = args$numDigitsRT, numDigitsMZ = args$numDigitsMZ)
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
211 # if we have 2 conditions, we keep stat of diffrep
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
212 if (!is.null(args$runDiffreport) & nlevels(sampclass(xset)) == 2) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
213 variableMetadata <- merge(variableMetadata, diffrep[, c("name", "fold", "tstat", "pvalue")], by.x = "name", by.y = "name", sort = F)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
214 if (exists("args[[\"sortpval\"]]")) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
215 variableMetadata <- variableMetadata[order(variableMetadata$pvalue), ]
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
216 }
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
217 }
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
218
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
219 variableMetadataOri <- variableMetadata
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
220 write.table(variableMetadata, sep = "\t", quote = FALSE, row.names = FALSE, file = variableMetadataOutput)
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
221
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
222 return(list("xa" = xa, "diffrep" = diffrepOri, "variableMetadata" = variableMetadataOri));
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
223
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
224 }
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
225
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
226
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
227 combinexsAnnos_function <- function(xaP, xaN, diffrepP = NULL, diffrepN = NULL,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
228 pos = TRUE, tol = 2, ruleset = NULL, keep_meta = TRUE, convertRTMinute = F, numDigitsMZ = 0,
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
229 numDigitsRT = 0, variableMetadataOutput = "variableMetadata.tsv") {
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
230
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
231 #Load the two Rdata to extract the xset objects from positive and negative mode
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
232 cat("\tObject xset from positive mode\n")
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
233 print(xaP)
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
234 cat("\n")
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
235
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
236 cat("\tObject xset from negative mode\n")
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
237 print(xaN)
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
238 cat("\n")
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
239
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
240 cat("\n")
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
241 cat("\tCombining...\n")
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
242 #Convert the string to numeric for creating matrix
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
243 row <- as.numeric(strsplit(ruleset, ",")[[1]][1])
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
244 column <- as.numeric(strsplit(ruleset, ",")[[1]][2])
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
245 ruleset <- cbind(row, column)
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
246 #Test if the file comes from an older version tool
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
247 if ((!is.null(xaP)) & (!is.null(xaN))) {
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
248 #Launch the combinexsannos function from CAMERA
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
249 cAnnot <- combinexsAnnos(xaP, xaN, pos = pos, tol = tol, ruleset = ruleset)
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
250 } else {
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
251 stop("You must relauch the CAMERA.annotate step with the lastest version.")
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
252 }
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
253
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
254 if (pos) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
255 xa <- xaP
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
256 mode <- "neg. Mode"
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
257 } else {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
258 xa <- xaN
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
259 mode <- "pos. Mode"
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
260 }
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
261
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
262 peakList <- getPeaklist(xa)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
263 peakList <- cbind(groupnames(xa@xcmsSet), peakList); colnames(peakList)[1] <- c("name");
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
264 variableMetadata <- cbind(peakList, cAnnot[, c("isotopes", "adduct", "pcgroup", mode)]);
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
265 variableMetadata <- variableMetadata[, !(colnames(variableMetadata) %in% c(sampnames(xa@xcmsSet)))]
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
266
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
267 #Test if there are more than two classes (conditions)
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
268 if (nlevels(sampclass(xaP@xcmsSet)) == 2 & (!is.null(diffrepN)) & (!is.null(diffrepP))) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
269 diffrepP <- diffrepP[, c("name", "fold", "tstat", "pvalue")]; colnames(diffrepP) <- paste("P.", colnames(diffrepP), sep = "")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
270 diffrepN <- diffrepN[, c("name", "fold", "tstat", "pvalue")]; colnames(diffrepN) <- paste("N.", colnames(diffrepN), sep = "")
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
271
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
272 variableMetadata <- merge(variableMetadata, diffrepP, by.x = "name", by.y = "P.name")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
273 variableMetadata <- merge(variableMetadata, diffrepN, by.x = "name", by.y = "N.name")
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
274 }
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
275
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
276 rownames(variableMetadata) <- NULL
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
277 #TODO: checker colnames(variableMetadata)[1:2] = c("name", "mz/rt");
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
278
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
279 variableMetadata <- RTSecondToMinute(variableMetadata, convertRTMinute)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
280 variableMetadata <- formatIonIdentifiers(variableMetadata, numDigitsRT = numDigitsRT, numDigitsMZ = numDigitsMZ)
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
281
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
282 #If the user want to keep only the metabolites which match a difference
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
283 if (keep_meta) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
284 variableMetadata <- variableMetadata[variableMetadata[, c(mode)] != "", ]
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
285 }
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
286
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
287 #Write the output into a tsv file
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
288 write.table(variableMetadata, sep = "\t", quote = FALSE, row.names = FALSE, file = variableMetadataOutput)
25
230f0bc9e792 planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents: 23
diff changeset
289 return(variableMetadata);
9
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
290
c02d80efba80 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
291 }
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
292
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
293 # This function get the raw file path from the arguments
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
294 getRawfilePathFromArguments <- function(singlefile, zipfile, args) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
295 if (!is.null(args$zipfile)) zipfile <- args$zipfile
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
296 if (!is.null(args$zipfilePositive)) zipfile <- args$zipfilePositive
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
297 if (!is.null(args$zipfileNegative)) zipfile <- args$zipfileNegative
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
298
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
299 if (!is.null(args$singlefile_galaxyPath)) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
300 singlefile_galaxyPaths <- args$singlefile_galaxyPath;
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
301 singlefile_sampleNames <- args$singlefile_sampleName
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
302 }
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
303 if (!is.null(args$singlefile_galaxyPathPositive)) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
304 singlefile_galaxyPaths <- args$singlefile_galaxyPathPositive;
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
305 singlefile_sampleNames <- args$singlefile_sampleNamePositive
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
306 }
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
307 if (!is.null(args$singlefile_galaxyPathNegative)) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
308 singlefile_galaxyPaths <- args$singlefile_galaxyPathNegative;
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
309 singlefile_sampleNames <- args$singlefile_sampleNameNegative
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
310 }
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
311 if (exists("singlefile_galaxyPaths")) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
312 singlefile_galaxyPaths <- unlist(strsplit(singlefile_galaxyPaths, ","))
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
313 singlefile_sampleNames <- unlist(strsplit(singlefile_sampleNames, ","))
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
314
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
315 singlefile <- NULL
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
316 for (singlefile_galaxyPath_i in seq_len(length(singlefile_galaxyPaths))) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
317 singlefile_galaxyPath <- singlefile_galaxyPaths[singlefile_galaxyPath_i]
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
318 singlefile_sampleName <- singlefile_sampleNames[singlefile_galaxyPath_i]
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
319 singlefile[[singlefile_sampleName]] <- singlefile_galaxyPath
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
320 }
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
321 }
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
322 for (argument in c("zipfile", "zipfilePositive", "zipfileNegative",
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
323 "singlefile_galaxyPath", "singlefile_sampleName",
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents: 37
diff changeset
324 "singlefile_galaxyPathPositive", "singlefile_sampleNamePositive",
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
325 "singlefile_galaxyPathNegative", "singlefile_sampleNameNegative")) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
326 args[[argument]] <- NULL
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
327 }
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
328 return(list(zipfile = zipfile, singlefile = singlefile, args = args))
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
329 }
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
330
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
331
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
332 # This function retrieve the raw file in the working directory
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
333 # - if zipfile: unzip the file with its directory tree
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
334 # - if singlefiles: set symlink with the good filename
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
335 retrieveRawfileInTheWorkingDir <- function(singlefile, zipfile) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
336 if (!is.null(singlefile) && (length("singlefile") > 0)) {
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
337 for (singlefile_sampleName in names(singlefile)) {
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
338 singlefile_galaxyPath <- singlefile[[singlefile_sampleName]]
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
339 if (!file.exists(singlefile_galaxyPath)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
340 error_message <- paste("Cannot access the sample:", singlefile_sampleName, "located:", singlefile_galaxyPath, ". Please, contact your administrator ... if you have one!")
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
341 print(error_message); stop(error_message)
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
342 }
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
343
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
344 file.symlink(singlefile_galaxyPath, singlefile_sampleName)
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
345 }
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
346 directory <- "."
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
347
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
348 }
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
349 if (!is.null(zipfile) && (zipfile != "")) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
350 if (!file.exists(zipfile)) {
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
351 error_message <- paste("Cannot access the Zip file:", zipfile, ". Please, contact your administrator ... if you have one!")
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
352 print(error_message)
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
353 stop(error_message)
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
354 }
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
355
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
356 #unzip
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
357 suppressWarnings(unzip(zipfile, unzip = "unzip"))
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
358
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
359 #get the directory name
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
360 filesInZip <- unzip(zipfile, list = T);
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
361 directories <- unique(unlist(lapply(strsplit(filesInZip$Name, "/"), function(x) x[1])));
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
362 directories <- directories[!(directories %in% c("__MACOSX")) & file.info(directories)$isdir]
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
363 directory <- "."
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
364 if (length(directories) == 1) directory <- directories
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
365
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
366 cat("files_root_directory\t", directory, "\n")
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
367
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
368 }
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 41
diff changeset
369 return(directory)
28
591d26b9027e planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents: 26
diff changeset
370 }