Mercurial > repos > mmonsoor > camera_combinexsannos

annotate CAMERA.r @ 5:d245f818f20a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 9d47e3b467dbbe0af0d90a937c5e9f2c4b958c4b
author lecorguille
date Mon, 25 Apr 2016 11:08:14 -0400
parents 9a1a72ee9418
children 972152941964
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
1 #!/usr/bin/env Rscript
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
2 # CAMERA.r version="2.2.1"
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
3
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
4
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
5
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
6 # ----- PACKAGE -----
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
7 cat("\tPACKAGE INFO\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
8
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
9 setRepositories(graphics=F, ind=31)
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
10
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
11 #pkgs=c("xcms","batch")
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
12 pkgs=c("parallel","BiocGenerics", "Biobase", "Rcpp", "mzR", "xcms","snow","igraph","CAMERA","multtest","batch")
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
13 for(p in pkgs) {
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
14 suppressPackageStartupMessages(suppressWarnings(library(p, quietly=TRUE, logical.return=TRUE, character.only=TRUE)))
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
15 cat(p,"\t",as.character(packageVersion(p)),"\n",sep="")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
16 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
17 source_local <- function(fname){ argv <- commandArgs(trailingOnly = FALSE); base_dir <- dirname(substring(argv[grep("--file=", argv)], 8)); source(paste(base_dir, fname, sep="/")) }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
18
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
19 cat("\n\n");
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
20
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
21
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
22
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
23 # ----- ARGUMENTS -----
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
24 cat("\tARGUMENTS INFO\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
25
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
26 listArguments = parseCommandArgs(evaluate=FALSE) #interpretation of arguments given in command line as an R list of objects
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
27 write.table(as.matrix(listArguments), col.names=F, quote=F, sep='\t')
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
28
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
29 cat("\n\n");
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
30
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
31
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
32 # ----- PROCESSING INFILE -----
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
33 cat("\tINFILE PROCESSING INFO\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
34
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
35 #image is an .RData file necessary to use xset variable given by previous tools
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
36 if (!is.null(listArguments[["image"]])){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
37 load(listArguments[["image"]]); listArguments[["image"]]=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
38 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
39
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
40 if (listArguments[["xfunction"]] %in% c("combinexsAnnos")) {
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
41 load(listArguments[["image_pos"]])
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
42 xaP=xa
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
43 listOFlistArgumentsP=listOFlistArguments
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
44 if (exists("xsAnnotate_object")) xaP=xsAnnotate_object
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
45
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
46 diffrepP=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
47 if (exists("diffrep")) diffrepP=diffrep
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
48
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
49 load(listArguments[["image_neg"]])
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
50 xaN=xa
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
51 listOFlistArgumentsN=listOFlistArguments
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
52 if (exists("xsAnnotate_object")) xaN=xsAnnotate_object
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
53
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
54 diffrepN=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
55 if (exists("diffrep")) diffrepN=diffrep
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
56 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
57
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
58
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
59 cat("\n\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
60
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
61
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
62 # ----- ARGUMENTS PROCESSING -----
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
63 cat("\tARGUMENTS PROCESSING INFO\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
64
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
65 # Save arguments to generate a report
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
66 if (!exists("listOFlistArguments")) listOFlistArguments=list()
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
67 listOFlistArguments[[paste(format(Sys.time(), "%y%m%d-%H:%M:%S_"),listArguments[["xfunction"]],sep="")]] = listArguments
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
68
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
69
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
70 #saving the commun parameters
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
71 thefunction = listArguments[["xfunction"]]
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
72 listArguments[["xfunction"]]=NULL #delete from the list of arguments
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
73
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
74 xsetRdataOutput = paste(thefunction,"RData",sep=".")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
75 if (!is.null(listArguments[["xsetRdataOutput"]])){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
76 xsetRdataOutput = listArguments[["xsetRdataOutput"]]; listArguments[["xsetRdataOutput"]]=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
77 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
78
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
79 rplotspdf = "Rplots.pdf"
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
80 if (!is.null(listArguments[["rplotspdf"]])){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
81 rplotspdf = listArguments[["rplotspdf"]]; listArguments[["rplotspdf"]]=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
82 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
83
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
84 dataMatrixOutput = "dataMatrix.tsv"
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
85 if (!is.null(listArguments[["dataMatrixOutput"]])){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
86 dataMatrixOutput = listArguments[["dataMatrixOutput"]]; listArguments[["dataMatrixOutput"]]=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
87 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
88
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
89 variableMetadataOutput = "variableMetadata.tsv"
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
90 if (!is.null(listArguments[["variableMetadataOutput"]])){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
91 variableMetadataOutput = listArguments[["variableMetadataOutput"]]; listArguments[["variableMetadataOutput"]]=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
92 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
93
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
94 if (!is.null(listArguments[["new_file_path"]])){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
95 new_file_path = listArguments[["new_file_path"]]; listArguments[["new_file_path"]]=NULL
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
96 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
97
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
98 #Import the different functions
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
99 source_local("lib.r")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
100
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
101 #necessary to unzip .zip file uploaded to Galaxy
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
102 #thanks to .zip file it's possible to upload many file as the same time conserving the tree hierarchy of directories
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
103
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
104
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
105 if (!is.null(listArguments[["zipfile"]])){
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
106 zipfile= listArguments[["zipfile"]]; listArguments[["zipfile"]]=NULL
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
107 }
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
108
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
109 # We unzip automatically the chromatograms from the zip files.
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
110 if (thefunction %in% c("annotatediff")) {
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
111 if(exists("zipfile") && (zipfile!="")) {
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
112 if(!file.exists(zipfile)){
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
113 error_message=paste("Cannot access the Zip file:",zipfile,". Please, contact your administrator ... if you have one!")
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
114 print(error_message)
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
115 stop(error_message)
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
116 }
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
117
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
118 #unzip
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
119 suppressWarnings(unzip(zipfile, unzip="unzip"))
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
120
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
121 #get the directory name
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
122 filesInZip=unzip(zipfile, list=T);
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
123 directories=unique(unlist(lapply(strsplit(filesInZip$Name,"/"), function(x) x[1])));
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
124 directories=directories[!(directories %in% c("__MACOSX")) & file.info(directories)$isdir]
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
125 directory = "."
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
126 if (length(directories) == 1) directory = directories
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
127
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
128 cat("files_root_directory\t",directory,"\n")
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
129 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
130 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
131
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
132
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
133
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
134 #addition of xset object to the list of arguments in the first position
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
135 if (exists("xset") != 0){
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
136 listArguments=append(list(xset), listArguments)
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
137 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
138
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
139 cat("\n\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
140
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
141
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
142
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
143
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
144 # ----- PROCESSING INFO -----
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
145 cat("\tMAIN PROCESSING INFO\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
146
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
147 #change the default display settings
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
148 pdf(file=rplotspdf, width=16, height=12)
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
149
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
150
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
151 if (thefunction %in% c("annotatediff")) {
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
152 results_list=annotatediff(xset=xset,listArguments=listArguments,variableMetadataOutput=variableMetadataOutput,dataMatrixOutput=dataMatrixOutput,new_file_path=new_file_path)
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
153 xa=results_list[["xa"]]
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
154 diffrep=results_list[["diffrep"]]
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
155 variableMetadata=results_list[["variableMetadata"]]
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
156
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
157 cat("\n\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
158 cat("\tXSET OBJECT INFO\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
159 print(xa)
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
160 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
161
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
162 if (thefunction %in% c("combinexsAnnos")) {
4
9a1a72ee9418 planemo upload commit 2761472c8afc41d54cecd54cabd1c5a5713514a3-dirty
lecorguille
parents: 2
diff changeset
163 cAnnot=combinexsAnnos_function(xaP=xaP,xaN=xaN,listOFlistArgumentsP=listOFlistArgumentsP,listOFlistArgumentsN=listOFlistArgumentsN,diffrepP=diffrepP,diffrepN=diffrepN,convert_param=listArguments[["convert_param"]],pos=listArguments[["pos"]],tol=listArguments[["tol"]],ruleset=listArguments[["ruleset"]],keep_meta=listArguments[["keep_meta"]],variableMetadataOutput=variableMetadataOutput)
2
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
164 }
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
165
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
166 dev.off()
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
167
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
168
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
169 #saving R data in .Rdata file to save the variables used in the present tool
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
170 objects2save = c("xa","variableMetadata","diffrep","cAnnot","listOFlistArguments","zipfile")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
171 save(list=objects2save[objects2save %in% ls()], file=xsetRdataOutput)
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
172
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
173 cat("\n\n")
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
174
2591056df814 planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff changeset
175 cat("\tDONE\n")