changeset 3:5dd59c833c34 draft

planemo upload for repository https://github.com/computational-metabolomics/mspurity-galaxy commit 7c8cb8294e711ca45f19de7bd9564157527f98db
author tomnl
date Fri, 14 Jun 2019 09:00:19 -0400
parents 7e7223015600
children 13985cdcf0ba
files flagRemove.R
diffstat 1 files changed, 20 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/flagRemove.R	Fri Jun 14 07:05:18 2019 -0400
+++ b/flagRemove.R	Fri Jun 14 09:00:19 2019 -0400
@@ -71,11 +71,9 @@
               help="Assign True if files for each step saved (for testing purposes) [default = %default]"
   ),
 
-  make_option("--samplelist",  type="character", help="Sample list to determine the blank class"),
+  make_option("--samplelist",  type="character", help="Sample list to determine the blank class")
 
-  make_option("--xset_name",  default="xset",
-              help="Name of the xcmsSet object within the RData file [default = %default]"
-  )
+
 
 
 
@@ -108,13 +106,28 @@
 
 print(opt)
 
-loadRData <- function(rdata_path, xset_name){
+getxcmsSetObject <- function(xobject) {
+    # XCMS 1.x
+    if (class(xobject) == "xcmsSet")
+        return (xobject)
+    # XCMS 3.x
+    if (class(xobject) == "XCMSnExp") {
+        # Get the legacy xcmsSet object
+        suppressWarnings(xset <- as(xobject, 'xcmsSet'))
+        sampclass(xset) <- xset@phenoData$sample_group
+        return (xset)
+    }
+}
+
+
+loadRData <- function(rdata_path, name){
 #loads an RData file, and returns the named xset object if it is there
     load(rdata_path)
-    return(get(ls()[ls() == xset_name]))
+    return(get(ls()[ls() %in% name]))
 }
 
-xset<-loadRData(opt$xset_path, opt$xset_name)
+xset <- loadRData(opt$xset_path, c('xset','xdata'))
+
 print(xset)
 if (is.null(opt$samplelist)){
     blank_class <- opt$blank_class