comparison xcms_xcmsSet.r @ 34:efd23113d5f4 draft

planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 49203f8a5271fa5e6bb889e907df71ebf7757309
author lecorguille
date Thu, 08 Mar 2018 05:51:32 -0500
parents c363b9f1caef
children 7f636ad1162b
comparison
equal deleted inserted replaced
33:c363b9f1caef 34:efd23113d5f4
113 xdata@processingData@files <- sub(paste(getwd(), "/", sep="") , "", xdata@processingData@files) 113 xdata@processingData@files <- sub(paste(getwd(), "/", sep="") , "", xdata@processingData@files)
114 114
115 # Create a sampleMetada file 115 # Create a sampleMetada file
116 sampleNamesList <- getSampleMetadata(xdata=xdata, sampleMetadataOutput="sampleMetadata.tsv") 116 sampleNamesList <- getSampleMetadata(xdata=xdata, sampleMetadataOutput="sampleMetadata.tsv")
117 117
118 # Get the legacy xcmsSet object
119 xset <- getxcmsSetObject(xdata)
120
121 cat("\n\n")
122
123
124 # -- TIC --
125 cat("\t\tGET TIC GRAPH\n")
126 #@TODO: one day, use xdata instead of xset to draw the TICs and BPC or a complete other method
127 getTICs(xcmsSet=xset, rt="raw", pdfname="TICs.pdf")
128 getBPCs(xcmsSet=xset, rt="raw", pdfname="BICs.pdf")
129
130 cat("\n\n") 118 cat("\n\n")
131 119
132 # ----- EXPORT ----- 120 # ----- EXPORT -----
133 121
134 cat("\tXCMSnExp OBJECT INFO\n") 122 cat("\tXCMSnExp OBJECT INFO\n")
135 print(xdata) 123 print(xdata)
136 cat("\n\n") 124 cat("\n\n")
137 125
138 cat("\txcmsSet OBJECT INFO\n") 126 cat("\txcmsSet OBJECT INFO\n")
127 # Get the legacy xcmsSet object
128 xset <- getxcmsSetObject(xdata)
139 print(xset) 129 print(xset)
140 cat("\n\n") 130 cat("\n\n")
141 131
142 #saving R data in .Rdata file to save the variables used in the present tool 132 #saving R data in .Rdata file to save the variables used in the present tool
143 objects2save <- c("xdata", "zipfile", "singlefile", "md5sumList", "sampleNamesList") 133 objects2save <- c("xdata", "zipfile", "singlefile", "md5sumList", "sampleNamesList")