Mercurial > repos > lecorguille > xcms_xcmsset
comparison xcms_xcmsSet.r @ 33:f5d51091cf84 draft
planemo upload for repository https://github.com/workflow4metabolomics/tools-metabolomics/ commit 2cb157bd9a8701a3d6874e084032cbd050b8953e
| author | workflow4metabolomics |
|---|---|
| date | Mon, 11 Sep 2023 09:22:36 +0000 |
| parents | b02d1992a43a |
| children | 3990a65b45a4 |
comparison
equal
deleted
inserted
replaced
| 32:b02d1992a43a | 33:f5d51091cf84 |
|---|---|
| 9 # ----- PACKAGE ----- | 9 # ----- PACKAGE ----- |
| 10 cat("\tSESSION INFO\n") | 10 cat("\tSESSION INFO\n") |
| 11 | 11 |
| 12 #Import the different functions | 12 #Import the different functions |
| 13 source_local <- function(fname) { | 13 source_local <- function(fname) { |
| 14 argv <- commandArgs(trailingOnly = FALSE); base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)); source(paste(base_dir, fname, sep = "/")) | 14 argv <- commandArgs(trailingOnly = FALSE) |
| 15 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)) | |
| 16 source(paste(base_dir, fname, sep = "/")) | |
| 15 } | 17 } |
| 16 source_local("lib.r") | 18 source_local("lib.r") |
| 17 | 19 |
| 18 pkgs <- c("xcms", "batch") | 20 pkgs <- c("xcms", "batch") |
| 19 loadAndDisplayPackages(pkgs) | 21 loadAndDisplayPackages(pkgs) |
| 20 cat("\n\n"); | 22 cat("\n\n") |
| 21 | 23 |
| 22 | 24 |
| 23 # ----- ARGUMENTS ----- | 25 # ----- ARGUMENTS ----- |
| 24 cat("\tARGUMENTS INFO\n") | 26 cat("\tARGUMENTS INFO\n") |
| 25 args <- parseCommandArgs(evaluate = FALSE) #interpretation of arguments given in command line as an R list of objects | 27 args <- parseCommandArgs(evaluate = FALSE) #interpretation of arguments given in command line as an R list of objects |
| 26 write.table(as.matrix(args), col.names = F, quote = F, sep = "\t") | 28 write.table(as.matrix(args), col.names = FALSE, quote = FALSE, sep = "\t") |
| 27 | 29 |
| 28 cat("\n\n") | 30 cat("\n\n") |
| 29 | 31 |
| 30 | 32 |
| 31 # ----- PROCESSING INFILE ----- | 33 # ----- PROCESSING INFILE ----- |
| 32 cat("\tARGUMENTS PROCESSING INFO\n") | 34 cat("\tARGUMENTS PROCESSING INFO\n") |
| 33 | 35 |
| 34 #saving the commun parameters | 36 #saving the commun parameters |
| 35 BPPARAM <- MulticoreParam(1) | 37 BPPARAM <- MulticoreParam(1) |
| 36 if (!is.null(args$BPPARAM)) { | 38 if (!is.null(args$BPPARAM)) { |
| 37 BPPARAM <- MulticoreParam(args$BPPARAM) | 39 BPPARAM <- MulticoreParam(args$BPPARAM) |
| 38 } | 40 } |
| 39 register(BPPARAM) | 41 register(BPPARAM) |
| 40 | 42 |
| 41 #saving the specific parameters | 43 #saving the specific parameters |
| 42 if (!is.null(args$filterAcquisitionNum)) filterAcquisitionNumParam <- args$filterAcquisitionNum | 44 if (!is.null(args$filterAcquisitionNum)) filterAcquisitionNumParam <- args$filterAcquisitionNum |
| 45 if (!is.null(args$peaklist)) peaklistParam <- args$peaklist | 47 if (!is.null(args$peaklist)) peaklistParam <- args$peaklist |
| 46 | 48 |
| 47 method <- args$method | 49 method <- args$method |
| 48 | 50 |
| 49 if (!is.null(args$roiList)) { | 51 if (!is.null(args$roiList)) { |
| 50 cat("\t\troiList provided\n") | 52 cat("\t\troiList provided\n") |
| 51 args$roiList <- list(getDataFrameFromFile(args$roiList)) | 53 args$roiList <- list(getDataFrameFromFile(args$roiList)) |
| 52 print(args$roiList) | 54 print(args$roiList) |
| 53 } | 55 } |
| 54 | 56 |
| 55 cat("\n\n") | 57 cat("\n\n") |
| 56 | 58 |
| 57 # ----- INFILE PROCESSING ----- | 59 # ----- INFILE PROCESSING ----- |
| 101 | 103 |
| 102 # Create a chromPeaks table if required | 104 # Create a chromPeaks table if required |
| 103 if (exists("peaklistParam")) { | 105 if (exists("peaklistParam")) { |
| 104 if (peaklistParam) { | 106 if (peaklistParam) { |
| 105 cat("\nCreating the chromatographic peaks' table...\n") | 107 cat("\nCreating the chromatographic peaks' table...\n") |
| 106 write.table(chromPeaks(xdata), file = "chromPeak_table.tsv", sep = "\t", quote = F, row.names = F) | 108 write.table(chromPeaks(xdata), file = "chromPeak_table.tsv", sep = "\t", quote = FALSE, row.names = FALSE) |
| 107 } | 109 } |
| 108 } | 110 } |
| 109 | 111 |
| 110 cat("\n\n") | 112 cat("\n\n") |
| 111 | 113 |
