Mercurial > repos > tomnl > mspurity_createmsp
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
