Mercurial > repos > lecorguille > xcms_retcor
annotate lib.r @ 25:aa252eec9229 draft
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
author | workflow4metabolomics |
---|---|
date | Wed, 07 Apr 2021 12:08:25 +0000 |
parents | c8011370e90f |
children | 36480435e92b |
rev | line source |
---|---|
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
1 #@authors ABiMS TEAM, Y. Guitton |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
2 # lib.r for Galaxy Workflow4Metabolomics xcms tools |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
3 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
4 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
5 # solve an issue with batch if arguments are logical TRUE/FALSE |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
6 parseCommandArgs <- function(...) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
7 args <- batch::parseCommandArgs(...) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
8 for (key in names(args)) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
9 if (args[key] %in% c("TRUE", "FALSE")) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
10 args[key] <- as.logical(args[key]) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
11 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
12 return(args) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
13 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
14 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
15 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
16 # This function will |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
17 # - load the packages |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
18 # - display the sessionInfo |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
19 loadAndDisplayPackages <- function(pkgs) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
20 for (pkg in pkgs) suppressPackageStartupMessages(stopifnot(library(pkg, quietly = TRUE, logical.return = TRUE, character.only = TRUE))) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
21 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
22 sessioninfo <- sessionInfo() |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
23 cat(sessioninfo$R.version$version.string, "\n") |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
24 cat("Main packages:\n") |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
25 for (pkg in names(sessioninfo$otherPkgs)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
26 cat(paste(pkg, packageVersion(pkg)), "\t") |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
27 } |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
28 cat("\n") |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
29 cat("Other loaded packages:\n") |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
30 for (pkg in names(sessioninfo$loadedOnly)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
31 cat(paste(pkg, packageVersion(pkg)), "\t") |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
32 } |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
33 cat("\n") |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
34 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
35 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
36 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
37 # This function merge several chromBPI or chromTIC into one. |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
38 mergeChrom <- function(chrom_merged, chrom) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
39 if (is.null(chrom_merged)) return(NULL) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
40 chrom_merged@.Data <- cbind(chrom_merged@.Data, chrom@.Data) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
41 return(chrom_merged) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
42 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
43 |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
44 #@author G. Le Corguille |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
45 # This function merge several xdata into one. |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
46 mergeXData <- function(args) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
47 chromTIC <- NULL |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
48 chromBPI <- NULL |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
49 chromTIC_adjusted <- NULL |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
50 chromBPI_adjusted <- NULL |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
51 md5sumList <- NULL |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
52 for (image in args$images) { |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
53 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
54 load(image) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
55 # Handle infiles |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
56 if (!exists("singlefile")) singlefile <- NULL |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
57 if (!exists("zipfile")) zipfile <- NULL |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
58 rawFilePath <- retrieveRawfileInTheWorkingDir(singlefile, zipfile, args) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
59 zipfile <- rawFilePath$zipfile |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
60 singlefile <- rawFilePath$singlefile |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
61 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
62 if (exists("raw_data")) xdata <- raw_data |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
63 if (!exists("xdata")) stop("\n\nERROR: The RData doesn't contain any object called 'xdata'. This RData should have been created by an old version of XMCS 2.*") |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
64 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
65 cat(sampleNamesList$sampleNamesOrigin, "\n") |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
66 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
67 if (!exists("xdata_merged")) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
68 xdata_merged <- xdata |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
69 singlefile_merged <- singlefile |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
70 md5sumList_merged <- md5sumList |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
71 sampleNamesList_merged <- sampleNamesList |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
72 chromTIC_merged <- chromTIC |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
73 chromBPI_merged <- chromBPI |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
74 chromTIC_adjusted_merged <- chromTIC_adjusted |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
75 chromBPI_adjusted_merged <- chromBPI_adjusted |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
76 } else { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
77 if (is(xdata, "XCMSnExp")) xdata_merged <- c(xdata_merged, xdata) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
78 else if (is(xdata, "OnDiskMSnExp")) xdata_merged <- xcms:::.concatenate_OnDiskMSnExp(xdata_merged, xdata) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
79 else stop("\n\nERROR: The RData either a OnDiskMSnExp object called raw_data or a XCMSnExp object called xdata") |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
80 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
81 singlefile_merged <- c(singlefile_merged, singlefile) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
82 md5sumList_merged$origin <- rbind(md5sumList_merged$origin, md5sumList$origin) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
83 sampleNamesList_merged$sampleNamesOrigin <- c(sampleNamesList_merged$sampleNamesOrigin, sampleNamesList$sampleNamesOrigin) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
84 sampleNamesList_merged$sampleNamesMakeNames <- c(sampleNamesList_merged$sampleNamesMakeNames, sampleNamesList$sampleNamesMakeNames) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
85 chromTIC_merged <- mergeChrom(chromTIC_merged, chromTIC) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
86 chromBPI_merged <- mergeChrom(chromBPI_merged, chromBPI) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
87 chromTIC_adjusted_merged <- mergeChrom(chromTIC_adjusted_merged, chromTIC_adjusted) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
88 chromBPI_adjusted_merged <- mergeChrom(chromBPI_adjusted_merged, chromBPI_adjusted) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
89 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
90 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
91 rm(image) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
92 xdata <- xdata_merged; rm(xdata_merged) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
93 singlefile <- singlefile_merged; rm(singlefile_merged) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
94 md5sumList <- md5sumList_merged; rm(md5sumList_merged) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
95 sampleNamesList <- sampleNamesList_merged; rm(sampleNamesList_merged) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
96 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
97 if (!is.null(args$sampleMetadata)) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
98 cat("\tXSET PHENODATA SETTING...\n") |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
99 sampleMetadataFile <- args$sampleMetadata |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
100 sampleMetadata <- getDataFrameFromFile(sampleMetadataFile, header = F) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
101 xdata@phenoData@data$sample_group <- sampleMetadata$V2[match(xdata@phenoData@data$sample_name, sampleMetadata$V1)] |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
102 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
103 if (any(is.na(pData(xdata)$sample_group))) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
104 sample_missing <- pData(xdata)$sample_name[is.na(pData(xdata)$sample_group)] |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
105 error_message <- paste("Those samples are missing in your sampleMetadata:", paste(sample_missing, collapse = " ")) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
106 print(error_message) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
107 stop(error_message) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
108 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
109 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
110 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
111 if (!is.null(chromTIC_merged)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
112 chromTIC <- chromTIC_merged; chromTIC@phenoData <- xdata@phenoData |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
113 } |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
114 if (!is.null(chromBPI_merged)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
115 chromBPI <- chromBPI_merged; chromBPI@phenoData <- xdata@phenoData |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
116 } |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
117 if (!is.null(chromTIC_adjusted_merged)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
118 chromTIC_adjusted <- chromTIC_adjusted_merged; chromTIC_adjusted@phenoData <- xdata@phenoData |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
119 } |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
120 if (!is.null(chromBPI_adjusted_merged)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
121 chromBPI_adjusted <- chromBPI_adjusted_merged; chromBPI_adjusted@phenoData <- xdata@phenoData |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
122 } |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
123 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
124 return(list("xdata" = xdata, "singlefile" = singlefile, "md5sumList" = md5sumList, "sampleNamesList" = sampleNamesList, "chromTIC" = chromTIC, "chromBPI" = chromBPI, "chromTIC_adjusted" = chromTIC_adjusted, "chromBPI_adjusted" = chromBPI_adjusted)) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
125 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
126 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
127 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
128 # This function convert if it is required the Retention Time in minutes |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
129 RTSecondToMinute <- function(variableMetadata, convertRTMinute) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
130 if (convertRTMinute) { |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
131 #converting the retention times (seconds) into minutes |
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
132 print("converting the retention times into minutes in the variableMetadata") |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
133 variableMetadata[, "rt"] <- variableMetadata[, "rt"] / 60 |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
134 variableMetadata[, "rtmin"] <- variableMetadata[, "rtmin"] / 60 |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
135 variableMetadata[, "rtmax"] <- variableMetadata[, "rtmax"] / 60 |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
136 } |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
137 return(variableMetadata) |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
138 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
139 |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
140 #@author G. Le Corguille |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
141 # This function format ions identifiers |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
142 formatIonIdentifiers <- function(variableMetadata, numDigitsRT = 0, numDigitsMZ = 0) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
143 splitDeco <- strsplit(as.character(variableMetadata$name), "_") |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
144 idsDeco <- sapply(splitDeco, |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
145 function(x) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
146 deco <- unlist(x)[2]; if (is.na(deco)) return("") else return(paste0("_", deco)) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
147 } |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
148 ) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
149 namecustom <- make.unique(paste0("M", round(variableMetadata[, "mz"], numDigitsMZ), "T", round(variableMetadata[, "rt"], numDigitsRT), idsDeco)) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
150 variableMetadata <- cbind(name = variableMetadata$name, namecustom = namecustom, variableMetadata[, !(colnames(variableMetadata) %in% c("name"))]) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
151 return(variableMetadata) |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
152 } |
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
153 |
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
154 #@author G. Le Corguille |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
155 # This function convert the remain NA to 0 in the dataMatrix |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
156 naTOzeroDataMatrix <- function(dataMatrix, naTOzero) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
157 if (naTOzero) { |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
158 dataMatrix[is.na(dataMatrix)] <- 0 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
159 } |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
160 return(dataMatrix) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
161 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
162 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
163 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
164 # Draw the plotChromPeakDensity 3 per page in a pdf file |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
165 getPlotChromPeakDensity <- function(xdata, param = NULL, mzdigit = 4) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
166 pdf(file = "plotChromPeakDensity.pdf", width = 16, height = 12) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
167 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
168 par(mfrow = c(3, 1), mar = c(4, 4, 1, 0.5)) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
169 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
170 if (length(unique(xdata$sample_group)) < 10) { |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
171 group_colors <- brewer.pal(length(unique(xdata$sample_group)), "Set1") |
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
172 }else{ |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
173 group_colors <- hcl.colors(length(unique(xdata$sample_group)), palette = "Dark 3") |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
174 } |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
175 names(group_colors) <- unique(xdata$sample_group) |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
176 col_per_samp <- as.character(xdata$sample_group) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
177 for (i in seq_len(length(group_colors))) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
178 col_per_samp[col_per_samp == (names(group_colors)[i])] <- group_colors[i] |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
179 } |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
180 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
181 xlim <- c(min(featureDefinitions(xdata)$rtmin), max(featureDefinitions(xdata)$rtmax)) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
182 for (i in seq_len(nrow(featureDefinitions(xdata)))) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
183 mzmin <- featureDefinitions(xdata)[i, ]$mzmin |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
184 mzmax <- featureDefinitions(xdata)[i, ]$mzmax |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
185 plotChromPeakDensity(xdata, param = param, mz = c(mzmin, mzmax), col = col_per_samp, pch = 16, xlim = xlim, main = paste(round(mzmin, mzdigit), round(mzmax, mzdigit))) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
186 legend("topright", legend = names(group_colors), col = group_colors, cex = 0.8, lty = 1) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
187 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
188 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
189 dev.off() |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
190 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
191 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
192 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
193 # Draw the plotChromPeakDensity 3 per page in a pdf file |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
194 getPlotAdjustedRtime <- function(xdata) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
195 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
196 pdf(file = "raw_vs_adjusted_rt.pdf", width = 16, height = 12) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
197 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
198 # Color by group |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
199 if (length(unique(xdata$sample_group)) < 10) { |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
200 group_colors <- brewer.pal(length(unique(xdata$sample_group)), "Set1") |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
201 } else { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
202 group_colors <- hcl.colors(length(unique(xdata$sample_group)), palette = "Dark 3") |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
203 } |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
204 if (length(group_colors) > 1) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
205 names(group_colors) <- unique(xdata$sample_group) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
206 plotAdjustedRtime(xdata, col = group_colors[xdata$sample_group]) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
207 legend("topright", legend = names(group_colors), col = group_colors, cex = 0.8, lty = 1) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
208 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
209 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
210 # Color by sample |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
211 plotAdjustedRtime(xdata, col = rainbow(length(xdata@phenoData@data$sample_name))) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
212 legend("topright", legend = xdata@phenoData@data$sample_name, col = rainbow(length(xdata@phenoData@data$sample_name)), cex = 0.8, lty = 1) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
213 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
214 dev.off() |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
215 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
216 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
217 #@author G. Le Corguille |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
218 # value: intensity values to be used into, maxo or intb |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
219 getPeaklistW4M <- function(xdata, intval = "into", convertRTMinute = F, numDigitsMZ = 4, numDigitsRT = 0, naTOzero = T, variableMetadataOutput, dataMatrixOutput, sampleNamesList) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
220 dataMatrix <- featureValues(xdata, method = "medret", value = intval) |
12
f85ebbed9ad0
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit e3956bf3393be5bff67b4d90f1007da00b0ec9b2-dirty
lecorguille
parents:
10
diff
changeset
|
221 colnames(dataMatrix) <- make.names(tools::file_path_sans_ext(colnames(dataMatrix))) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
222 dataMatrix <- cbind(name = groupnames(xdata), dataMatrix) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
223 variableMetadata <- featureDefinitions(xdata) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
224 colnames(variableMetadata)[1] <- "mz"; colnames(variableMetadata)[4] <- "rt" |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
225 variableMetadata <- data.frame(name = groupnames(xdata), variableMetadata) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
226 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
227 variableMetadata <- RTSecondToMinute(variableMetadata, convertRTMinute) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
228 variableMetadata <- formatIonIdentifiers(variableMetadata, numDigitsRT = numDigitsRT, numDigitsMZ = numDigitsMZ) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
229 dataMatrix <- naTOzeroDataMatrix(dataMatrix, naTOzero) |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
230 |
14
e95f0dcdea3b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 37b0a6a7686f701e4bf00db97ae2c1b82cd6e989
lecorguille
parents:
12
diff
changeset
|
231 # FIX: issue when the vector at peakidx is too long and is written in a new line during the export |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
232 variableMetadata[, "peakidx"] <- vapply(variableMetadata[, "peakidx"], FUN = paste, FUN.VALUE = character(1), collapse = ",") |
14
e95f0dcdea3b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 37b0a6a7686f701e4bf00db97ae2c1b82cd6e989
lecorguille
parents:
12
diff
changeset
|
233 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
234 write.table(variableMetadata, file = variableMetadataOutput, sep = "\t", quote = F, row.names = F) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
235 write.table(dataMatrix, file = dataMatrixOutput, sep = "\t", quote = F, row.names = F) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
236 |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
237 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
238 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
239 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
240 # It allow different of field separators |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
241 getDataFrameFromFile <- function(filename, header = T) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
242 myDataFrame <- read.table(filename, header = header, sep = ";", stringsAsFactors = F) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
243 if (ncol(myDataFrame) < 2) myDataFrame <- read.table(filename, header = header, sep = "\t", stringsAsFactors = F) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
244 if (ncol(myDataFrame) < 2) myDataFrame <- read.table(filename, header = header, sep = ",", stringsAsFactors = F) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
245 if (ncol(myDataFrame) < 2) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
246 error_message <- "Your tabular file seems not well formatted. The column separators accepted are ; , and tabulation" |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
247 print(error_message) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
248 stop(error_message) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
249 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
250 return(myDataFrame) |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
251 } |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
252 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
253 #@author G. Le Corguille |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
254 # Draw the BPI and TIC graphics |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
255 # colored by sample names or class names |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
256 getPlotChromatogram <- function(chrom, xdata, pdfname = "Chromatogram.pdf", aggregationFun = "max") { |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
257 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
258 if (aggregationFun == "sum") |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
259 type <- "Total Ion Chromatograms" |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
260 else |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
261 type <- "Base Peak Intensity Chromatograms" |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
262 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
263 adjusted <- "Raw" |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
264 if (hasAdjustedRtime(xdata)) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
265 adjusted <- "Adjusted" |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
266 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
267 main <- paste(type, ":", adjusted, "data") |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
268 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
269 pdf(pdfname, width = 16, height = 10) |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
270 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
271 # Color by group |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
272 if (length(unique(xdata$sample_group)) < 10) { |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
273 group_colors <- brewer.pal(length(unique(xdata$sample_group)), "Set1") |
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
274 }else{ |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
275 group_colors <- hcl.colors(length(unique(xdata$sample_group)), palette = "Dark 3") |
24
c8011370e90f
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 0dcc8bcb1e645574c7f81ec1a43f86be38acd065"
lecorguille
parents:
22
diff
changeset
|
276 } |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
277 if (length(group_colors) > 1) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
278 names(group_colors) <- unique(xdata$sample_group) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
279 plot(chrom, col = group_colors[chrom$sample_group], main = main, peakType = "none") |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
280 legend("topright", legend = names(group_colors), col = group_colors, cex = 0.8, lty = 1) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
281 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
282 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
283 # Color by sample |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
284 plot(chrom, col = rainbow(length(xdata@phenoData@data$sample_name)), main = main, peakType = "none") |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
285 legend("topright", legend = xdata@phenoData@data$sample_name, col = rainbow(length(xdata@phenoData@data$sample_name)), cex = 0.8, lty = 1) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
286 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
287 dev.off() |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
288 } |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
289 |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
290 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
291 # Get the polarities from all the samples of a condition |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
292 #@author Misharl Monsoor misharl.monsoor@sb-roscoff.fr ABiMS TEAM |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
293 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr ABiMS TEAM |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
294 getSampleMetadata <- function(xdata = NULL, sampleMetadataOutput = "sampleMetadata.tsv") { |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
295 cat("Creating the sampleMetadata file...\n") |
7
bb602a5b8819
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 08e7f269a5c59687a7768be8db5fcb4e4d736093
lecorguille
parents:
2
diff
changeset
|
296 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
297 #Create the sampleMetada dataframe |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
298 sampleMetadata <- xdata@phenoData@data |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
299 rownames(sampleMetadata) <- NULL |
14
e95f0dcdea3b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 37b0a6a7686f701e4bf00db97ae2c1b82cd6e989
lecorguille
parents:
12
diff
changeset
|
300 colnames(sampleMetadata) <- c("sample_name", "class") |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
301 |
14
e95f0dcdea3b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 37b0a6a7686f701e4bf00db97ae2c1b82cd6e989
lecorguille
parents:
12
diff
changeset
|
302 sampleNamesOrigin <- sampleMetadata$sample_name |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
303 sampleNamesMakeNames <- make.names(sampleNamesOrigin) |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
304 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
305 if (any(duplicated(sampleNamesMakeNames))) { |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
306 write("\n\nERROR: Usually, R has trouble to deal with special characters in its column names, so it rename them using make.names().\nIn your case, at least two columns after the renaming obtain the same name, thus XCMS will collapse those columns per name.", stderr()) |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
307 for (sampleName in sampleNamesOrigin) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
308 write(paste(sampleName, "\t->\t", make.names(sampleName)), stderr()) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
309 } |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
310 stop("\n\nERROR: One or more of your files will not be import by xcmsSet. It may due to bad characters in their filenames.") |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
311 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
312 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
313 if (!all(sampleNamesOrigin == sampleNamesMakeNames)) { |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
314 cat("\n\nWARNING: Usually, R has trouble to deal with special characters in its column names, so it rename them using make.names()\nIn your case, one or more sample names will be renamed in the sampleMetadata and dataMatrix files:\n") |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
315 for (sampleName in sampleNamesOrigin) { |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
316 cat(paste(sampleName, "\t->\t", make.names(sampleName), "\n")) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
317 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
318 } |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
319 |
14
e95f0dcdea3b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 37b0a6a7686f701e4bf00db97ae2c1b82cd6e989
lecorguille
parents:
12
diff
changeset
|
320 sampleMetadata$sample_name <- sampleNamesMakeNames |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
321 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
322 |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
323 #For each sample file, the following actions are done |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
324 for (fileIdx in seq_len(length(fileNames(xdata)))) { |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
325 #Check if the file is in the CDF format |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
326 if (!mzR:::netCDFIsFile(fileNames(xdata))) { |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
327 |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
328 # If the column isn't exist, with add one filled with NA |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
329 if (is.null(sampleMetadata$polarity)) sampleMetadata$polarity <- NA |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
330 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
331 #Extract the polarity (a list of polarities) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
332 polarity <- fData(xdata)[fData(xdata)$fileIdx == fileIdx, "polarity"] |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
333 #Verify if all the scans have the same polarity |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
334 uniq_list <- unique(polarity) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
335 if (length(uniq_list) > 1) { |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
336 polarity <- "mixed" |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
337 } else { |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
338 polarity <- as.character(uniq_list) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
339 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
340 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
341 #Set the polarity attribute |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
342 sampleMetadata$polarity[fileIdx] <- polarity |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
343 } |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
344 |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
345 } |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
346 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
347 write.table(sampleMetadata, sep = "\t", quote = FALSE, row.names = FALSE, file = sampleMetadataOutput) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
348 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
349 return(list("sampleNamesOrigin" = sampleNamesOrigin, "sampleNamesMakeNames" = sampleNamesMakeNames)) |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
350 |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
351 } |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
352 |
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
353 |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
354 # This function will compute MD5 checksum to check the data integrity |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
355 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
356 getMd5sum <- function(files) { |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
357 cat("Compute md5 checksum...\n") |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
358 library(tools) |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
359 return(as.matrix(md5sum(files))) |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
360 } |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
361 |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
362 # This function retrieve the raw file in the working directory |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
363 # - if zipfile: unzip the file with its directory tree |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
364 # - if singlefiles: set symlink with the good filename |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
365 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
366 retrieveRawfileInTheWorkingDir <- function(singlefile, zipfile, args, prefix = "") { |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
367 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
368 if (!(prefix %in% c("", "Positive", "Negative", "MS1", "MS2"))) stop("prefix must be either '', 'Positive', 'Negative', 'MS1' or 'MS2'") |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
369 |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
370 # single - if the file are passed in the command arguments -> refresh singlefile |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
371 if (!is.null(args[[paste0("singlefile_galaxyPath", prefix)]])) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
372 singlefile_galaxyPaths <- unlist(strsplit(args[[paste0("singlefile_galaxyPath", prefix)]], "\\|")) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
373 singlefile_sampleNames <- unlist(strsplit(args[[paste0("singlefile_sampleName", prefix)]], "\\|")) |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
374 |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
375 singlefile <- NULL |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
376 for (singlefile_galaxyPath_i in seq_len(length(singlefile_galaxyPaths))) { |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
377 singlefile_galaxyPath <- singlefile_galaxyPaths[singlefile_galaxyPath_i] |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
378 singlefile_sampleName <- singlefile_sampleNames[singlefile_galaxyPath_i] |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
379 # In case, an url is used to import data within Galaxy |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
380 singlefile_sampleName <- tail(unlist(strsplit(singlefile_sampleName, "/")), n = 1) |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
381 singlefile[[singlefile_sampleName]] <- singlefile_galaxyPath |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
382 } |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
383 } |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
384 # zipfile - if the file are passed in the command arguments -> refresh zipfile |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
385 if (!is.null(args[[paste0("zipfile", prefix)]])) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
386 zipfile <- args[[paste0("zipfile", prefix)]] |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
387 |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
388 # single |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
389 if (!is.null(singlefile) && (length("singlefile") > 0)) { |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
390 files <- vector() |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
391 for (singlefile_sampleName in names(singlefile)) { |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
392 singlefile_galaxyPath <- singlefile[[singlefile_sampleName]] |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
393 if (!file.exists(singlefile_galaxyPath)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
394 error_message <- paste("Cannot access the sample:", singlefile_sampleName, "located:", singlefile_galaxyPath, ". Please, contact your administrator ... if you have one!") |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
395 print(error_message); stop(error_message) |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
396 } |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
397 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
398 if (!suppressWarnings(try(file.link(singlefile_galaxyPath, singlefile_sampleName), silent = T))) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
399 file.copy(singlefile_galaxyPath, singlefile_sampleName) |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
400 files <- c(files, singlefile_sampleName) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
401 } |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
402 } |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
403 # zipfile |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
404 if (!is.null(zipfile) && (zipfile != "")) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
405 if (!file.exists(zipfile)) { |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
406 error_message <- paste("Cannot access the Zip file:", zipfile, ". Please, contact your administrator ... if you have one!") |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
407 print(error_message) |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
408 stop(error_message) |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
409 } |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
410 suppressWarnings(unzip(zipfile, unzip = "unzip")) |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
411 |
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
412 #get the directory name |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
413 suppressWarnings(filesInZip <- unzip(zipfile, list = T)) |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
414 directories <- unique(unlist(lapply(strsplit(filesInZip$Name, "/"), function(x) x[1]))) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
415 directories <- directories[!(directories %in% c("__MACOSX")) & file.info(directories)$isdir] |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
416 directory <- "." |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
417 if (length(directories) == 1) directory <- directories |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
418 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
419 cat("files_root_directory\t", directory, "\n") |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
420 |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
421 filepattern <- c("[Cc][Dd][Ff]", "[Nn][Cc]", "([Mm][Zz])?[Xx][Mm][Ll]", "[Mm][Zz][Dd][Aa][Tt][Aa]", "[Mm][Zz][Mm][Ll]") |
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
422 filepattern <- paste(paste("\\.", filepattern, "$", sep = ""), collapse = "|") |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
423 info <- file.info(directory) |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
424 listed <- list.files(directory[info$isdir], pattern = filepattern, recursive = TRUE, full.names = TRUE) |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
425 files <- c(directory[!info$isdir], listed) |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
426 exists <- file.exists(files) |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
427 files <- files[exists] |
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
428 |
8
4bfef820569b
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit a6f5f18b3d6130f7d7fbb9f2df856838c6217797
lecorguille
parents:
7
diff
changeset
|
429 } |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
430 return(list(zipfile = zipfile, singlefile = singlefile, files = files)) |
22
972e76582173
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit f1caf2a3bf23cf319a75dd12c86402555dd02617"
workflow4metabolomics
parents:
19
diff
changeset
|
431 |
2
54cc3edfe35c
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 83b80dcd96b379518c2e4ace992affc889d32ca6
lecorguille
parents:
diff
changeset
|
432 } |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
433 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
434 |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
435 # This function retrieve a xset like object |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
436 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
437 getxcmsSetObject <- function(xobject) { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
438 # XCMS 1.x |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
439 if (class(xobject) == "xcmsSet") |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
440 return(xobject) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
441 # XCMS 3.x |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
442 if (class(xobject) == "XCMSnExp") { |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
443 # Get the legacy xcmsSet object |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
444 suppressWarnings(xset <- as(xobject, "xcmsSet")) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
445 if (!is.null(xset@phenoData$sample_group)) |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
446 sampclass(xset) <- xset@phenoData$sample_group |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
447 else |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
448 sampclass(xset) <- "." |
25
aa252eec9229
"planemo upload for repository https://github.com/workflow4metabolomics/xcms commit dcc90f9cf76e6980c0a7d9698c89fab826e7adae"
workflow4metabolomics
parents:
24
diff
changeset
|
449 return(xset) |
10
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
450 } |
8828cba9aedd
planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 9f72e947d9c241d11221cad561f3525d27231857
lecorguille
parents:
8
diff
changeset
|
451 } |