Mercurial > repos > lecorguille > camera_annotate

annotate CAMERA_annotateDiffreport.r @ 47:1cd9574ee4f8 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 24d44ee26b7c23380c2b42fae2f7f6e58472100d
author workflow4metabolomics
date Sun, 24 Nov 2024 21:24:17 +0000
parents 67f0ef6236f1
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
1 #!/usr/bin/env Rscript
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
2
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
3 # ----- PACKAGE -----
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
4 cat("\tSESSION INFO\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
5
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
6 # Import the different functions
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
7 source_local <- function(fname) {
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
8 argv <- commandArgs(trailingOnly = FALSE)
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
9 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8))
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
10 source(paste(base_dir, fname, sep = "/"))
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
11 }
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
12 source_local("lib.r")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
13
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
14 pkgs <- c("CAMERA", "multtest", "batch")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
15 loadAndDisplayPackages(pkgs)
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
16 cat("\n\n")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
17 # ----- ARGUMENTS -----
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
18 cat("\tARGUMENTS INFO\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
19
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
20 args <- parseCommandArgs(evaluate = FALSE) # interpretation of arguments given in command line as an R list of objects
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
21 write.table(as.matrix(args), col.names = FALSE, quote = FALSE, sep = "\t")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
22
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
23 cat("\n\n")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
24 # ----- PROCESSING INFILE -----
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
25 cat("\tINFILE PROCESSING INFO\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
26
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
27 # image is an .RData file necessary to use xset variable given by previous tools
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
28 load(args$image)
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
29 args$image <- NULL
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
30
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
31 cat("\n\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
32
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
33
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
34 # ----- ARGUMENTS PROCESSING -----
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
35 cat("\tARGUMENTS PROCESSING INFO\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
36
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
37 # Save arguments to generate a report
46
67f0ef6236f1 planemo upload commit a923a432239a1f2ffee2dae1281fe2d8f41bac38
workflow4metabolomics
parents: 45
diff changeset
38 if (!exists("listOFlistArguments")) listOFlistArguments <- list()
67f0ef6236f1 planemo upload commit a923a432239a1f2ffee2dae1281fe2d8f41bac38
workflow4metabolomics
parents: 45
diff changeset
39 listOFlistArguments[[format(Sys.time(), "%y%m%d-%H:%M:%S_annotatediff")]] <- args
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
40
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
41 # We unzip automatically the chromatograms from the zip files.
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
42 if (!exists("zipfile")) zipfile <- NULL
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
43 if (!exists("singlefile")) singlefile <- NULL
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
44 rawFilePath <- getRawfilePathFromArguments(singlefile, zipfile, args)
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
45 zipfile <- rawFilePath$zipfile
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
46 singlefile <- rawFilePath$singlefile
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
47 args <- rawFilePath$args
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
48 directory <- retrieveRawfileInTheWorkingDir(singlefile, zipfile)
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
49
40
3d943f088143 "planemo upload commit 4fcbbcbc6d6b0a59e801870d31fe886a920ef429"
workflow4metabolomics
parents: 39
diff changeset
50
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
51 # Because so far CAMERA isn't compatible with the new XCMSnExp object
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
52 if (exists("xdata")) {
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
53 xset <- getxcmsSetObject(xdata)
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
54 }
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
55
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
56 cat("\n\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
57
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
58
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
59 # ----- PROCESSING INFO -----
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
60 cat("\tMAIN PROCESSING INFO\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
61
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
62 results_list <- annotatediff(xset = xset, args = args, variableMetadataOutput = "variableMetadata.tsv")
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
63 xa <- results_list$xa
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
64 diffrep <- results_list$diffrep
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
65 variableMetadata <- results_list$variableMetadata
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
66
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
67 cat("\n\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
68
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
69 # ----- EXPORT -----
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
70
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
71 cat("\tXSET OBJECT INFO\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
72 print(xa)
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
73 cat("\n\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
74
45
db77f2b67a7f planemo upload commit cfad09eb4dd6b1439b7de6a0852cd8fa22210f58
workflow4metabolomics
parents: 43
diff changeset
75 # saving R data in .Rdata file to save the variables used in the present tool
46
67f0ef6236f1 planemo upload commit a923a432239a1f2ffee2dae1281fe2d8f41bac38
workflow4metabolomics
parents: 45
diff changeset
76 objects2save <- c("xa", "variableMetadata", "diffrep", "cAnnot", "listOFlistArguments", "zipfile", "singlefile")
43
7f4ae4db39eb "planemo upload commit c1e276cf53b3c54d4702ab26d4f40a525a720998"
workflow4metabolomics
parents: 40
diff changeset
77 save(list = objects2save[objects2save %in% ls()], file = "annotatediff.RData")
39
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
78
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
79 cat("\n\n")
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
80
db9bc2c27a0a planemo upload commit d51a0d0a26833671b79711ee18d782e84f301e76
workflow4metabolomics
parents:
diff changeset
81 cat("\tDONE\n")