changeset 3:60b156e58529 draft

planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 5e066c50d2e4ac6d5796b53331cbd3198ab8c4f9-dirty

--- a/lib.r	Tue Apr 03 11:37:40 2018 -0400
+++ b/lib.r	Thu Apr 05 18:06:32 2018 -0400
@@ -28,8 +28,18 @@
 }
 
 #@author G. Le Corguille
+# This function merge several chromBPI or chromTIC into one.
+mergeChrom <- function(chromTIC_merged, chromTIC, xdata_merged) {
+    if (is.null(chromTIC_merged)) return(NULL)
+    chromTIC_merged@.Data <- cbind(chromTIC_merged@.Data, chromTIC@.Data)
+    chromTIC_merged@phenoData <- xdata_merged@phenoData
+    return(chromTIC_merged)
+}
+
+#@author G. Le Corguille
 # This function merge several xdata into one.
 mergeXData <- function(args) {
+    chromTIC <- NULL; chromBPI <- NULL
     for(image in args$images) {
         load(image)
         # Handle infiles
@@ -47,6 +57,8 @@
             singlefile_merged <- singlefile
             md5sumList_merged <- md5sumList
             sampleNamesList_merged <- sampleNamesList
+            chromTIC_merged <- chromTIC
+            chromBPI_merged <- chromBPI
         } else {
             if (is(xdata, "XCMSnExp")) xdata_merged <- c(xdata_merged,xdata)
             else if (is(xdata, "OnDiskMSnExp")) xdata_merged <- .concatenate_OnDiskMSnExp(xdata_merged,xdata)
@@ -55,6 +67,8 @@
             md5sumList_merged$origin <- rbind(md5sumList_merged$origin,md5sumList$origin)
             sampleNamesList_merged$sampleNamesOrigin <- c(sampleNamesList_merged$sampleNamesOrigin,sampleNamesList$sampleNamesOrigin)
             sampleNamesList_merged$sampleNamesMakeNames <- c(sampleNamesList_merged$sampleNamesMakeNames,sampleNamesList$sampleNamesMakeNames)
+            chromTIC_merged <- mergeChrom(chromTIC_merged, chromTIC, xdata_merged)
+            chromBPI_merged <- mergeChrom(chromBPI_merged, chromBPI, xdata_merged)
         }
     }
     rm(image)
@@ -62,6 +76,8 @@
     singlefile <- singlefile_merged; rm(singlefile_merged)
     md5sumList <- md5sumList_merged; rm(md5sumList_merged)
     sampleNamesList <- sampleNamesList_merged; rm(sampleNamesList_merged)
+    chromTIC <- chromTIC_merged; rm(chromTIC_merged)
+    chromBPI <- chromBPI_merged; rm(chromBPI_merged)
 
     if (!is.null(args$sampleMetadata)) {
         cat("\tXSET PHENODATA SETTING...\n")
@@ -76,7 +92,7 @@
             stop(error_message)
         }
     }
-    return(list("xdata"=xdata, "singlefile"=singlefile, "md5sumList"=md5sumList,"sampleNamesList"=sampleNamesList))
+    return(list("xdata"=xdata, "singlefile"=singlefile, "md5sumList"=md5sumList,"sampleNamesList"=sampleNamesList, "chromTIC"=chromTIC, "chromBPI"=chromBPI))
 }
 
 #@author G. Le Corguille
@@ -176,9 +192,8 @@
     return(myDataFrame)
 }
 
-getPlotChromatogram <- function(xdata, pdfname="Chromatogram.pdf", aggregationFun = "max") {
+getPlotChromatogram <- function(chrom, xdata, pdfname="Chromatogram.pdf", aggregationFun = "max") {
 
-    chrom <- chromatogram(xdata, aggregationFun = aggregationFun)
     if (aggregationFun == "sum")
         type="Total Ion Chromatograms"
     else
@@ -207,16 +222,6 @@
     dev.off()
 }
 
-#@author G. Le Corguille
-getPlotTICs <- function(xdata, pdfname="TICs.pdf") {
-    getPlotChromatogram(xdata, pdfname, aggregationFun = "sum")
-}
-
-#@author G. Le Corguille
-getPlotBPIs <- function(xdata, pdfname="BPIs.pdf") {
-    getPlotChromatogram(xdata, pdfname, aggregationFun = "max")
-}
-
 
 # Get the polarities from all the samples of a condition
 #@author Misharl Monsoor [e-mail] ABiMS TEAM

Binary file test-data/ko15-raw.Rdata has changed

Binary file test-data/ko16-raw.Rdata has changed

Binary file test-data/wt15-raw.Rdata has changed

Binary file test-data/wt16-raw.Rdata has changed

--- a/xcms_plot_chromatogram.r	Tue Apr 03 11:37:40 2018 -0400
+++ b/xcms_plot_chromatogram.r	Thu Apr 05 18:06:32 2018 -0400
@@ -39,6 +39,8 @@
 singlefile <- mergeXDataReturn$singlefile
 md5sumList <- mergeXDataReturn$md5sumList
 sampleNamesList <- mergeXDataReturn$sampleNamesList
+chromTIC <- mergeXDataReturn$chromTIC
+chromBPI <- mergeXDataReturn$chromBPI
 
 cat("\n\n")
 
@@ -49,8 +51,11 @@
 
 cat("\t\tDRAW GRAPHICS\n")
 
-getPlotTICs(xdata, pdfname="TICs.pdf")
-getPlotBPIs(xdata, pdfname="BPIs.pdf")
+if (!exists("chromTIC") || is.null(chromTIC)) { cat("\t\t\tCompute TIC\n"); chromTIC <- chromatogram(xdata, aggregationFun = "sum") }
+if (!exists("chromBPI") || is.null(chromBPI)) { cat("\t\t\tCompute BPI\n"); chromBPI <- chromatogram(xdata, aggregationFun = "max") }
+
+getPlotChromatogram(chromTIC, xdata, pdfname="TICs.pdf", aggregationFun = "sum")
+getPlotChromatogram(chromBPI, xdata, pdfname="BPIs.pdf", aggregationFun = "max")
 
 cat("\n\n")