comparison lib.r @ 36:e309e6af6744 draft

planemo upload for repository https://github.com/workflow4metabolomics/xcms commit 3922356e5edae763800a4311146e2fdeb52c296e
author lecorguille
date Fri, 31 Aug 2018 09:02:01 -0400
parents 2b0a4c7a4a48
children 35a20d7c9f33
comparison
equal deleted inserted replaced
35:2b0a4c7a4a48 36:e309e6af6744
389 } 389 }
390 390
391 391
392 # This function get the raw file path from the arguments 392 # This function get the raw file path from the arguments
393 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr 393 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr
394 getRawfilePathFromArguments <- function(singlefile, zipfile, args) { 394 getRawfilePathFromArguments <- function(singlefile, zipfile, args, prefix="") {
395 if (!is.null(args$zipfile)) zipfile <- args$zipfile 395 if (!(prefix %in% c("","Positive","Negative","MS1","MS2"))) stop("prefix must be either '', 'Positive', 'Negative', 'MS1' or 'MS2'")
396 if (!is.null(args$zipfilePositive)) zipfile <- args$zipfilePositive 396
397 if (!is.null(args$zipfileNegative)) zipfile <- args$zipfileNegative 397 if (!is.null(args[[paste0("zipfile",prefix)]])) zipfile <- args[[paste0("zipfile",prefix)]]
398 398
399 if (!is.null(args$singlefile_galaxyPath)) { 399 if (!is.null(args[[paste0("singlefile_galaxyPath",prefix)]])) {
400 singlefile_galaxyPaths <- args$singlefile_galaxyPath; 400 singlefile_galaxyPaths <- args[[paste0("singlefile_galaxyPath",prefix)]]
401 singlefile_sampleNames <- args$singlefile_sampleName 401 singlefile_sampleNames <- args[[paste0("singlefile_sampleName",prefix)]]
402 } 402 }
403 if (!is.null(args$singlefile_galaxyPathPositive)) { 403 if (exists("singlefile_galaxyPaths")){
404 singlefile_galaxyPaths <- args$singlefile_galaxyPathPositive; 404 singlefile_galaxyPaths <- unlist(strsplit(singlefile_galaxyPaths,"\\|"))
405 singlefile_sampleNames <- args$singlefile_sampleNamePositive 405 singlefile_sampleNames <- unlist(strsplit(singlefile_sampleNames,"\\|"))
406 } 406
407 if (!is.null(args$singlefile_galaxyPathNegative)) { 407 singlefile <- NULL
408 singlefile_galaxyPaths <- args$singlefile_galaxyPathNegative; 408 for (singlefile_galaxyPath_i in seq(1:length(singlefile_galaxyPaths))) {
409 singlefile_sampleNames <- args$singlefile_sampleNameNegative 409 singlefile_galaxyPath <- singlefile_galaxyPaths[singlefile_galaxyPath_i]
410 } 410 singlefile_sampleName <- singlefile_sampleNames[singlefile_galaxyPath_i]
411 if (exists("singlefile_galaxyPaths")){ 411 # In case, an url is used to import data within Galaxy
412 singlefile_galaxyPaths <- unlist(strsplit(singlefile_galaxyPaths,"\\|")) 412 singlefile_sampleName <- tail(unlist(strsplit(singlefile_sampleName,"/")), n=1)
413 singlefile_sampleNames <- unlist(strsplit(singlefile_sampleNames,"\\|")) 413 singlefile[[singlefile_sampleName]] <- singlefile_galaxyPath
414 414 }
415 singlefile <- NULL 415 }
416 for (singlefile_galaxyPath_i in seq(1:length(singlefile_galaxyPaths))) { 416 return(list(zipfile=zipfile, singlefile=singlefile))
417 singlefile_galaxyPath <- singlefile_galaxyPaths[singlefile_galaxyPath_i] 417 }
418 singlefile_sampleName <- singlefile_sampleNames[singlefile_galaxyPath_i]
419 singlefile[[singlefile_sampleName]] <- singlefile_galaxyPath
420 }
421 }
422 return(list(zipfile=zipfile, singlefile=singlefile))
423 }
424
425 418
426 # This function retrieve the raw file in the working directory 419 # This function retrieve the raw file in the working directory
427 # - if zipfile: unzip the file with its directory tree 420 # - if zipfile: unzip the file with its directory tree
428 # - if singlefiles: set symlink with the good filename 421 # - if singlefiles: set symlink with the good filename
429 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr 422 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr
478 return (xobject) 471 return (xobject)
479 # XCMS 3.x 472 # XCMS 3.x
480 if (class(xobject) == "XCMSnExp") { 473 if (class(xobject) == "XCMSnExp") {
481 # Get the legacy xcmsSet object 474 # Get the legacy xcmsSet object
482 suppressWarnings(xset <- as(xobject, 'xcmsSet')) 475 suppressWarnings(xset <- as(xobject, 'xcmsSet'))
483 sampclass(xset) <- xset@phenoData$sample_group 476 if (!is.null(xset@phenoData$sample_group))
477 sampclass(xset) <- xset@phenoData$sample_group
478 else
479 sampclass(xset) <- "."
484 return (xset) 480 return (xset)
485 } 481 }
486 } 482 }
487 483
488 484