Mercurial > repos > workflow4metabolomics > xcms_refine
diff xcms_refine.r @ 3:6c2798bce3eb draft default tip
planemo upload for repository https://github.com/workflow4metabolomics/tools-metabolomics/ commit 95721ced8347c09e79340e6d67ecb41c5cc64163
author | workflow4metabolomics |
---|---|
date | Mon, 03 Feb 2025 14:46:55 +0000 |
parents | 20f8ebc3a391 |
children |
line wrap: on
line diff
--- a/xcms_refine.r Mon Jul 15 16:00:52 2024 +0000 +++ b/xcms_refine.r Mon Feb 03 14:46:55 2025 +0000 @@ -1,7 +1,7 @@ #!/usr/bin/env Rscript # ----- LOG FILE ----- -log_file <- file("log.txt", open = "wt") +log_file <- file("log.txt", open = "wt") sink(log_file) sink(log_file, type = "output") @@ -9,11 +9,11 @@ # ----- PACKAGE ----- cat("\tSESSION INFO\n") -#Import the different functions +# Import the different functions source_local <- function(fname) { - argv <- commandArgs(trailingOnly = FALSE) - base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)) - source(paste(base_dir, fname, sep = "/")) + argv <- commandArgs(trailingOnly = FALSE) + base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)) + source(paste(base_dir, fname, sep = "/")) } source_local("lib.r") @@ -32,23 +32,23 @@ # ----- PROCESSING INFILE ----- cat("\tARGUMENTS PROCESSING INFO\n") -#saving the specific parameters -args_method <- args$method -args_image <- args$image +# saving the specific parameters +args_method <- args$method +args_image <- args$image args_msLevel <- args$msLevel param_args <- list() if (args_method == "CleanPeaks") { - param_args$maxPeakwidth <- args$maxPeakwidth + param_args$maxPeakwidth <- args$maxPeakwidth } else if (args_method == "FilterIntensity") { - param_args$threshold <- args$threshold - param_args$value <- args$value - param_args$nValues <- args$nValues + param_args$threshold <- args$threshold + param_args$value <- args$value + param_args$nValues <- args$nValues } else if (args_method == "MergeNeighboringPeaks") { - param_args$expandRt <- args$expandRt - param_args$expandMz <- args$expandMz - param_args$ppm <- args$ppm - param_args$minProp <- args$minProp + param_args$expandRt <- args$expandRt + param_args$expandMz <- args$expandMz + param_args$ppm <- args$ppm + param_args$minProp <- args$minProp } cat("\n\n") @@ -57,7 +57,7 @@ # ----- ARGUMENTS PROCESSING ----- cat("\tINFILE PROCESSING INFO\n") -#image is an .RData file necessary to use xset variable given by previous tools +# image is an .RData file necessary to use xset variable given by previous tools load(args_image) if (!exists("xdata")) stop("\n\nERROR: The RData doesn't contain any object called 'xdata'. Such RData as this might have been created by an old version of XMCS 2.*") @@ -78,20 +78,20 @@ cat("\t\tPREPARE PARAMETERS\n\n") if (args_method == "CleanPeaks") { - refineChromPeaksParam <- CleanPeaksParam(maxPeakwidth = param_args$maxPeakwidth) + refineChromPeaksParam <- CleanPeaksParam(maxPeakwidth = param_args$maxPeakwidth) } else if (args_method == "FilterIntensity") { - refineChromPeaksParam <- FilterIntensityParam( - threshold = param_args$threshold, - nValues = param_args$nValues, - value = param_args$value - ) + refineChromPeaksParam <- FilterIntensityParam( + threshold = param_args$threshold, + nValues = param_args$nValues, + value = param_args$value + ) } else if (args_method == "MergeNeighboringPeaks") { - refineChromPeaksParam <- MergeNeighboringPeaksParam( - expandRt = param_args$expandRt, - expandMz = param_args$expandMz, - ppm = param_args$ppm, - minProp = param_args$minProp - ) + refineChromPeaksParam <- MergeNeighboringPeaksParam( + expandRt = param_args$expandRt, + expandMz = param_args$expandMz, + ppm = param_args$ppm, + minProp = param_args$minProp + ) } cat(str(refineChromPeaksParam)) @@ -116,7 +116,7 @@ print(xset) cat("\n\n") -#saving R data in .Rdata file to save the variables used in the present tool +# saving R data in .Rdata file to save the variables used in the present tool objects2save <- c("xdata", "zipfile", "singlefile", "md5sumList", "sampleNamesList") save(list = objects2save[objects2save %in% ls()], file = "xcmsSet.RData")