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