Mercurial > repos > lecorguille > camera_annotate
annotate lib.r @ 37:636e36a64d31 draft
planemo upload commit a5f400649acb64e1c66a92379027ac5abdb4f546
author | lecorguille |
---|---|
date | Tue, 09 Apr 2019 11:36:39 -0400 |
parents | b6a8143c397c |
children | db9bc2c27a0a |
rev | line source |
---|---|
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
1 # lib.r |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
2 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
3 # This function retrieve a xset like object |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
4 #@author Gildas Le Corguille lecorguille@sb-roscoff.fr |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
5 getxcmsSetObject <- function(xobject) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
6 # XCMS 1.x |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
7 if (class(xobject) == "xcmsSet") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
8 return (xobject) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
9 # XCMS 3.x |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
10 if (class(xobject) == "XCMSnExp") { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
11 # Get the legacy xcmsSet object |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
12 suppressWarnings(xset <- as(xobject, 'xcmsSet')) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
13 if (is.null(xset@phenoData$sample_group)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
14 sampclass(xset) = "." |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
15 else |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
16 sampclass(xset) <- xset@phenoData$sample_group |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
17 if (!is.null(xset@phenoData$sample_name)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
18 rownames(xset@phenoData) = xset@phenoData$sample_name |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
19 return (xset) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
20 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
21 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
22 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
23 #@author G. Le Corguille |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
24 #The function create a pdf from the different png generated by diffreport |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
25 diffreport_png2pdf <- function(filebase) { |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
26 dir.create("pdf") |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
27 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
28 pdfEicOutput = paste0("pdf/",filebase,"-eic_pdf.pdf") |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
29 pdfBoxOutput = paste0("pdf/",filebase,"-box_pdf.pdf") |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
30 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
31 system(paste0("gm convert ",filebase,"_eic/*.png ",pdfEicOutput)) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
32 system(paste0("gm convert ",filebase,"_box/*.png ",pdfBoxOutput)) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
33 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
34 } |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
35 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
36 #@author G. Le Corguille |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
37 #The function create a zip archive from the different png generated by diffreport |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
38 diffreport_png2zip <- function() { |
37
636e36a64d31
planemo upload commit a5f400649acb64e1c66a92379027ac5abdb4f546
lecorguille
parents:
36
diff
changeset
|
39 zip("eic.zip", dir(pattern="_eic"), zip=Sys.which("zip")) |
636e36a64d31
planemo upload commit a5f400649acb64e1c66a92379027ac5abdb4f546
lecorguille
parents:
36
diff
changeset
|
40 zip("box.zip", dir(pattern="_box"), zip=Sys.which("zip")) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
41 } |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
42 |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
43 #The function create a zip archive from the different tabular generated by diffreport |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
44 diffreport_tabular2zip <- function() { |
37
636e36a64d31
planemo upload commit a5f400649acb64e1c66a92379027ac5abdb4f546
lecorguille
parents:
36
diff
changeset
|
45 zip("tabular.zip", dir(pattern="tabular/*"), zip=Sys.which("zip")) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
46 } |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
47 |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
48 #@author G. Le Corguille |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
49 #This function convert if it is required the Retention Time in minutes |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
50 RTSecondToMinute <- function(variableMetadata, convertRTMinute) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
51 if (convertRTMinute){ |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
52 #converting the retention times (seconds) into minutes |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
53 print("converting the retention times into minutes in the variableMetadata") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
54 variableMetadata[,"rt"]=variableMetadata[,"rt"]/60 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
55 variableMetadata[,"rtmin"]=variableMetadata[,"rtmin"]/60 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
56 variableMetadata[,"rtmax"]=variableMetadata[,"rtmax"]/60 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
57 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
58 return (variableMetadata) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
59 } |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
60 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
61 #@author G. Le Corguille |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
62 #This function format ions identifiers |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
63 formatIonIdentifiers <- function(variableMetadata, numDigitsRT=0, numDigitsMZ=0) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
64 splitDeco = strsplit(as.character(variableMetadata$name),"_") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
65 idsDeco = sapply(splitDeco, function(x) { deco=unlist(x)[2]; if (is.na(deco)) return ("") else return(paste0("_",deco)) }) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
66 namecustom = make.unique(paste0("M",round(variableMetadata[,"mz"],numDigitsMZ),"T",round(variableMetadata[,"rt"],numDigitsRT),idsDeco)) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
67 variableMetadata=cbind(name=variableMetadata$name, namecustom=namecustom, variableMetadata[,!(colnames(variableMetadata) %in% c("name"))]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
68 return(variableMetadata) |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
69 } |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
70 |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
71 #The function annotateDiffreport without the corr function which bugs |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
72 annotatediff <- function(xset=xset, listArguments=listArguments, variableMetadataOutput="variableMetadata.tsv") { |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
73 # Resolve the bug with x11, with the function png |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
74 options(bitmapType='cairo') |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
75 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
76 #Check if the fillpeaks step has been done previously, if it hasn't, there is an error message and the execution is stopped. |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
77 res=try(is.null(xset@filled)) |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
78 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
79 # ------ annot ------- |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
80 listArguments[["calcCiS"]]=as.logical(listArguments[["calcCiS"]]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
81 listArguments[["calcIso"]]=as.logical(listArguments[["calcIso"]]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
82 listArguments[["calcCaS"]]=as.logical(listArguments[["calcCaS"]]) |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
83 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
84 # common parameters |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
85 listArguments4annotate = list(object=xset, |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
86 nSlaves=listArguments[["nSlaves"]],sigma=listArguments[["sigma"]],perfwhm=listArguments[["perfwhm"]], |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
87 maxcharge=listArguments[["maxcharge"]],maxiso=listArguments[["maxiso"]],minfrac=listArguments[["minfrac"]], |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
88 ppm=listArguments[["ppm"]],mzabs=listArguments[["mzabs"]],quick=listArguments[["quick"]], |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
89 polarity=listArguments[["polarity"]],max_peaks=listArguments[["max_peaks"]],intval=listArguments[["intval"]]) |
22
65e3e1a39e16
planemo upload commit 23ad0382ec59d66220fd94684a6d158fdef2c7f5-dirty
lecorguille
parents:
21
diff
changeset
|
90 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
91 # quick == FALSE |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
92 if(listArguments[["quick"]]==FALSE) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
93 listArguments4annotate = append(listArguments4annotate, |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
94 list(graphMethod=listArguments[["graphMethod"]],cor_eic_th=listArguments[["cor_eic_th"]],pval=listArguments[["pval"]], |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
95 calcCiS=listArguments[["calcCiS"]],calcIso=listArguments[["calcIso"]],calcCaS=listArguments[["calcCaS"]])) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
96 # no ruleset |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
97 if (!is.null(listArguments[["multiplier"]])) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
98 listArguments4annotate = append(listArguments4annotate, |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
99 list(multiplier=listArguments[["multiplier"]])) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
100 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
101 # ruleset |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
102 else { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
103 rulset=read.table(listArguments[["rules"]], h=T, sep=";") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
104 if (ncol(rulset) < 4) rulset=read.table(listArguments[["rules"]], h=T, sep="\t") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
105 if (ncol(rulset) < 4) rulset=read.table(listArguments[["rules"]], h=T, sep=",") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
106 if (ncol(rulset) < 4) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
107 error_message="Your ruleset file seems not well formatted. The column separators accepted are ; , and tabulation" |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
108 print(error_message) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
109 stop(error_message) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
110 } |
23
2a4a811c663d
planemo upload commit bf9996f1e380af6f783cc862f3e1d69f8c7b86e5
lecorguille
parents:
22
diff
changeset
|
111 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
112 listArguments4annotate = append(listArguments4annotate, |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
113 list(rules=rulset)) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
114 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
115 } |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
116 |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
117 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
118 # launch annotate |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
119 xa = do.call("annotate", listArguments4annotate) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
120 peakList=getPeaklist(xa,intval=listArguments[["intval"]]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
121 peakList=cbind(groupnames(xa@xcmsSet),peakList); colnames(peakList)[1] = c("name"); |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
122 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
123 # --- Multi condition : diffreport --- |
26
f769635d55f9
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa-dirty
lecorguille
parents:
25
diff
changeset
|
124 diffrepOri=NULL |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
125 if (!is.null(listArguments[["runDiffreport"]]) & nlevels(sampclass(xset))>=2) { |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
126 #Check if the fillpeaks step has been done previously, if it hasn't, there is an error message and the execution is stopped. |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
127 res=try(is.null(xset@filled)) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
128 classes=levels(sampclass(xset)) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
129 x=1:(length(classes)-1) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
130 for (i in seq(along=x) ) { |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
131 y=1:(length(classes)) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
132 for (n in seq(along=y)){ |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
133 if(i+n <= length(classes)){ |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
134 filebase=paste(classes[i],class2=classes[i+n],sep="-vs-") |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
135 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
136 diffrep=diffreport(object=xset,class1=classes[i],class2=classes[i+n],filebase=filebase,eicmax=listArguments[["eicmax"]],eicwidth=listArguments[["eicwidth"]],sortpval=TRUE,value=listArguments[["value"]],h=listArguments[["h"]],w=listArguments[["w"]],mzdec=listArguments[["mzdec"]],missing=0) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
137 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
138 diffrepOri = diffrep |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
139 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
140 # renamming of the column rtmed to rt to fit with camera peaklist function output |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
141 colnames(diffrep)[colnames(diffrep)=="rtmed"] <- "rt" |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
142 colnames(diffrep)[colnames(diffrep)=="mzmed"] <- "mz" |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
143 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
144 # combines results and reorder columns |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
145 diffrep = merge(peakList, diffrep[,c("name","fold","tstat","pvalue")], by.x="name", by.y="name", sort=F) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
146 diffrep = cbind(diffrep[,!(colnames(diffrep) %in% c(sampnames(xa@xcmsSet)))],diffrep[,(colnames(diffrep) %in% c(sampnames(xa@xcmsSet)))]) |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
147 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
148 diffrep = RTSecondToMinute(diffrep, listArguments[["convertRTMinute"]]) |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
149 diffrep = formatIonIdentifiers(diffrep, numDigitsRT=listArguments[["numDigitsRT"]], numDigitsMZ=listArguments[["numDigitsMZ"]]) |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
150 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
151 if(listArguments[["sortpval"]]){ |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
152 diffrep=diffrep[order(diffrep$pvalue), ] |
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
153 } |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
154 |
37
636e36a64d31
planemo upload commit a5f400649acb64e1c66a92379027ac5abdb4f546
lecorguille
parents:
36
diff
changeset
|
155 dir.create("tabular", showWarnings = FALSE) |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
156 write.table(diffrep, sep="\t", quote=FALSE, row.names=FALSE, file=paste("tabular/",filebase,"_tsv.tabular",sep="")) |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
157 |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
158 if (listArguments[["eicmax"]] != 0) { |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
159 if (listArguments[["png2"]] == "pdf") |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
160 diffreport_png2pdf(filebase) |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
161 } |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
162 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
163 } |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
164 } |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
165 if (listArguments[["png2"]] == "zip") |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
166 diffreport_png2zip() |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
167 if (listArguments[["tabular2"]] == "zip") |
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
168 diffreport_tabular2zip() |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
169 } |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
170 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
171 # --- variableMetadata --- |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
172 variableMetadata=peakList[,!(make.names(colnames(peakList)) %in% c(make.names(sampnames(xa@xcmsSet))))] |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
173 variableMetadata = RTSecondToMinute(variableMetadata, listArguments[["convertRTMinute"]]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
174 variableMetadata = formatIonIdentifiers(variableMetadata, numDigitsRT=listArguments[["numDigitsRT"]], numDigitsMZ=listArguments[["numDigitsMZ"]]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
175 # if we have 2 conditions, we keep stat of diffrep |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
176 if (!is.null(listArguments[["runDiffreport"]]) & nlevels(sampclass(xset))==2) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
177 variableMetadata = merge(variableMetadata, diffrep[,c("name","fold","tstat","pvalue")],by.x="name", by.y="name", sort=F) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
178 if(exists("listArguments[[\"sortpval\"]]")){ |
33
3b2b642542f6
planemo upload commit 896ef116e72008394779322f65137c632af5f4d2
lecorguille
parents:
30
diff
changeset
|
179 variableMetadata=variableMetadata[order(variableMetadata$pvalue), ] |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
180 } |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
181 } |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
182 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
183 variableMetadataOri=variableMetadata |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
184 write.table(variableMetadata, sep="\t", quote=FALSE, row.names=FALSE, file=variableMetadataOutput) |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
185 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
186 return(list("xa"=xa,"diffrep"=diffrepOri,"variableMetadata"=variableMetadataOri)); |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
187 |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
188 } |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
189 |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
190 |
30
fce12d692a41
planemo upload commit 2e3fcb07e9ff59974eaec8f868a9e78a63a89219
lecorguille
parents:
29
diff
changeset
|
191 combinexsAnnos_function <- function(xaP, xaN, listOFlistArgumentsP,listOFlistArgumentsN, diffrepP=NULL,diffrepN=NULL,pos=TRUE,tol=2,ruleset=NULL,keep_meta=TRUE, convertRTMinute=F, numDigitsMZ=0, numDigitsRT=0, variableMetadataOutput="variableMetadata.tsv"){ |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
192 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
193 #Load the two Rdata to extract the xset objects from positive and negative mode |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
194 cat("\tObject xset from positive mode\n") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
195 print(xaP) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
196 cat("\n") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
197 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
198 cat("\tObject xset from negative mode\n") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
199 print(xaN) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
200 cat("\n") |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
201 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
202 cat("\n") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
203 cat("\tCombining...\n") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
204 #Convert the string to numeric for creating matrix |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
205 row=as.numeric(strsplit(ruleset,",")[[1]][1]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
206 column=as.numeric(strsplit(ruleset,",")[[1]][2]) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
207 ruleset=cbind(row,column) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
208 #Test if the file comes from an older version tool |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
209 if ((!is.null(xaP)) & (!is.null(xaN))) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
210 #Launch the combinexsannos function from CAMERA |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
211 cAnnot=combinexsAnnos(xaP, xaN,pos=pos,tol=tol,ruleset=ruleset) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
212 } else { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
213 stop("You must relauch the CAMERA.annotate step with the lastest version.") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
214 } |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
215 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
216 if(pos){ |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
217 xa=xaP |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
218 listOFlistArgumentsP=listOFlistArguments |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
219 mode="neg. Mode" |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
220 } else { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
221 xa=xaN |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
222 listOFlistArgumentsN=listOFlistArguments |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
223 mode="pos. Mode" |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
224 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
225 |
30
fce12d692a41
planemo upload commit 2e3fcb07e9ff59974eaec8f868a9e78a63a89219
lecorguille
parents:
29
diff
changeset
|
226 peakList=getPeaklist(xa) |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
227 peakList=cbind(groupnames(xa@xcmsSet),peakList); colnames(peakList)[1] = c("name"); |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
228 variableMetadata=cbind(peakList, cAnnot[, c("isotopes", "adduct", "pcgroup",mode)]); |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
229 variableMetadata=variableMetadata[,!(colnames(variableMetadata) %in% c(sampnames(xa@xcmsSet)))] |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
230 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
231 #Test if there are more than two classes (conditions) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
232 if ( nlevels(sampclass(xaP@xcmsSet))==2 & (!is.null(diffrepN)) & (!is.null(diffrepP))) { |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
233 diffrepP = diffrepP[,c("name","fold","tstat","pvalue")]; colnames(diffrepP) = paste("P.",colnames(diffrepP),sep="") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
234 diffrepN = diffrepN[,c("name","fold","tstat","pvalue")]; colnames(diffrepN) = paste("N.",colnames(diffrepN),sep="") |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
235 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
236 variableMetadata = merge(variableMetadata, diffrepP, by.x="name", by.y="P.name") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
237 variableMetadata = merge(variableMetadata, diffrepN, by.x="name", by.y="N.name") |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
238 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
239 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
240 rownames(variableMetadata) = NULL |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
241 #TODO: checker |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
242 #colnames(variableMetadata)[1:2] = c("name","mz/rt"); |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
243 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
244 variableMetadata = RTSecondToMinute(variableMetadata, convertRTMinute) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
245 variableMetadata = formatIonIdentifiers(variableMetadata, numDigitsRT=numDigitsRT, numDigitsMZ=numDigitsMZ) |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
246 |
25
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
247 #If the user want to keep only the metabolites which match a difference |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
248 if(keep_meta){ |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
249 variableMetadata=variableMetadata[variableMetadata[,c(mode)]!="",] |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
250 } |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
251 |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
252 #Write the output into a tsv file |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
253 write.table(variableMetadata, sep="\t", quote=FALSE, row.names=FALSE, file=variableMetadataOutput) |
230f0bc9e792
planemo upload commit a634879c0e651eb0eb7b435a107ee40cf30524fa
lecorguille
parents:
23
diff
changeset
|
254 return(variableMetadata); |
9
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
255 |
c02d80efba80
planemo upload commit d8cc436fd91f5748dc396d0527a0a303d3221835
lecorguille
parents:
diff
changeset
|
256 } |
28
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
257 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
258 # This function get the raw file path from the arguments |
29
59bdfc862408
planemo upload commit 00bda51a5454373957a0b75024d9fc4282743c4d
lecorguille
parents:
28
diff
changeset
|
259 getRawfilePathFromArguments <- function(singlefile, zipfile, listArguments) { |
28
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
260 if (!is.null(listArguments[["zipfile"]])) zipfile = listArguments[["zipfile"]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
261 if (!is.null(listArguments[["zipfilePositive"]])) zipfile = listArguments[["zipfilePositive"]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
262 if (!is.null(listArguments[["zipfileNegative"]])) zipfile = listArguments[["zipfileNegative"]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
263 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
264 if (!is.null(listArguments[["singlefile_galaxyPath"]])) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
265 singlefile_galaxyPaths = listArguments[["singlefile_galaxyPath"]]; |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
266 singlefile_sampleNames = listArguments[["singlefile_sampleName"]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
267 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
268 if (!is.null(listArguments[["singlefile_galaxyPathPositive"]])) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
269 singlefile_galaxyPaths = listArguments[["singlefile_galaxyPathPositive"]]; |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
270 singlefile_sampleNames = listArguments[["singlefile_sampleNamePositive"]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
271 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
272 if (!is.null(listArguments[["singlefile_galaxyPathNegative"]])) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
273 singlefile_galaxyPaths = listArguments[["singlefile_galaxyPathNegative"]]; |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
274 singlefile_sampleNames = listArguments[["singlefile_sampleNameNegative"]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
275 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
276 if (exists("singlefile_galaxyPaths")){ |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
277 singlefile_galaxyPaths = unlist(strsplit(singlefile_galaxyPaths,",")) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
278 singlefile_sampleNames = unlist(strsplit(singlefile_sampleNames,",")) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
279 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
280 singlefile=NULL |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
281 for (singlefile_galaxyPath_i in seq(1:length(singlefile_galaxyPaths))) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
282 singlefile_galaxyPath=singlefile_galaxyPaths[singlefile_galaxyPath_i] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
283 singlefile_sampleName=singlefile_sampleNames[singlefile_galaxyPath_i] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
284 singlefile[[singlefile_sampleName]] = singlefile_galaxyPath |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
285 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
286 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
287 for (argument in c("zipfile","zipfilePositive","zipfileNegative","singlefile_galaxyPath","singlefile_sampleName","singlefile_galaxyPathPositive","singlefile_sampleNamePositive","singlefile_galaxyPathNegative","singlefile_sampleNameNegative")) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
288 listArguments[[argument]]=NULL |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
289 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
290 return(list(zipfile=zipfile, singlefile=singlefile, listArguments=listArguments)) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
291 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
292 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
293 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
294 # This function retrieve the raw file in the working directory |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
295 # - if zipfile: unzip the file with its directory tree |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
296 # - if singlefiles: set symlink with the good filename |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
297 retrieveRawfileInTheWorkingDirectory <- function(singlefile, zipfile) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
298 if(!is.null(singlefile) && (length("singlefile")>0)) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
299 for (singlefile_sampleName in names(singlefile)) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
300 singlefile_galaxyPath = singlefile[[singlefile_sampleName]] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
301 if(!file.exists(singlefile_galaxyPath)){ |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
302 error_message=paste("Cannot access the sample:",singlefile_sampleName,"located:",singlefile_galaxyPath,". Please, contact your administrator ... if you have one!") |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
303 print(error_message); stop(error_message) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
304 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
305 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
306 file.symlink(singlefile_galaxyPath,singlefile_sampleName) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
307 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
308 directory = "." |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
309 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
310 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
311 if(!is.null(zipfile) && (zipfile!="")) { |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
312 if(!file.exists(zipfile)){ |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
313 error_message=paste("Cannot access the Zip file:",zipfile,". Please, contact your administrator ... if you have one!") |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
314 print(error_message) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
315 stop(error_message) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
316 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
317 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
318 #list all file in the zip file |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
319 #zip_files=unzip(zipfile,list=T)[,"Name"] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
320 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
321 #unzip |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
322 suppressWarnings(unzip(zipfile, unzip="unzip")) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
323 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
324 #get the directory name |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
325 filesInZip=unzip(zipfile, list=T); |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
326 directories=unique(unlist(lapply(strsplit(filesInZip$Name,"/"), function(x) x[1]))); |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
327 directories=directories[!(directories %in% c("__MACOSX")) & file.info(directories)$isdir] |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
328 directory = "." |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
329 if (length(directories) == 1) directory = directories |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
330 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
331 cat("files_root_directory\t",directory,"\n") |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
332 |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
333 } |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
334 return (directory) |
591d26b9027e
planemo upload commit d708b85fc764e790fac071552d19cd7a89d02d31
lecorguille
parents:
26
diff
changeset
|
335 } |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
336 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
337 #@TODO: remove this function as soon as we can use xcms 3.x.x from Bioconductor 3.7 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
338 # https://github.com/sneumann/CAMERA/issues/33#issuecomment-405168524 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
339 # https://github.com/sneumann/xcms/commit/950a3fe794cdb6b0fda88696e31aab3d97a3b7dd |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
340 ############################################################ |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
341 ## getEIC |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
342 getEIC <- function(object, mzrange, rtrange = 200, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
343 groupidx, sampleidx = sampnames(object), |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
344 rt = c("corrected", "raw")) { |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
345 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
346 files <- filepaths(object) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
347 grp <- groups(object) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
348 samp <- sampnames(object) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
349 prof <- profinfo(object) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
350 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
351 rt <- match.arg(rt) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
352 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
353 if (is.numeric(sampleidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
354 sampleidx <- sampnames(object)[sampleidx] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
355 sampidx <- match(sampleidx, sampnames(object)) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
356 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
357 if (!missing(groupidx)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
358 if (is.numeric(groupidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
359 groupidx <- groupnames(object)[unique(as.integer(groupidx))] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
360 grpidx <- match(groupidx, groupnames(object, template = groupidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
361 } |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
362 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
363 if (missing(mzrange)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
364 if (missing(groupidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
365 stop("No m/z range or groups specified") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
366 if (any(is.na(groupval(object, value = "mz")))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
367 warning( |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
368 "`NA` values in xcmsSet. Use fillPeaks() on the object to fill", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
369 "-in missing peak values. Note however that this will also ", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
370 "insert intensities of 0 for peaks that can not be filled in.") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
371 mzmin <- apply(groupval(object, value = "mzmin"), 1, min, na.rm = TRUE) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
372 mzmax <- apply(groupval(object, value = "mzmax"), 1, max, na.rm = TRUE) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
373 mzrange <- matrix(c(mzmin[grpidx], mzmax[grpidx]), ncol = 2) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
374 ## if (any(is.na(groupval(object, value = "mz")))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
375 ## stop('Please use fillPeaks() to fill up NA values !') |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
376 ## mzmin <- -rowMax(-groupval(object, value = "mzmin")) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
377 ## mzmax <- rowMax(groupval(object, value = "mzmax")) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
378 ## mzrange <- matrix(c(mzmin[grpidx], mzmax[grpidx]), ncol = 2) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
379 } else if (all(c("mzmin","mzmax") %in% colnames(mzrange))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
380 mzrange <- mzrange[,c("mzmin", "mzmax"),drop=FALSE] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
381 else if (is.null(dim(mzrange))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
382 stop("mzrange must be a matrix") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
383 colnames(mzrange) <- c("mzmin", "mzmax") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
384 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
385 if (length(rtrange) == 1) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
386 if (missing(groupidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
387 rtrange <- matrix(rep(range(object@rt[[rt]][sampidx]), nrow(mzrange)), |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
388 ncol = 2, byrow = TRUE) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
389 else { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
390 rtrange <- retexp(grp[grpidx,c("rtmin","rtmax"),drop=FALSE], rtrange) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
391 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
392 } else if (is.null(dim(rtrange))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
393 stop("rtrange must be a matrix or single number") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
394 colnames(rtrange) <- c("rtmin", "rtmax") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
395 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
396 ## Ensure that we've got corrected retention time if requested. |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
397 if (is.null(object@rt[[rt]])) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
398 stop(rt, " retention times not present in 'object'!") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
399 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
400 ## Ensure that the defined retention time range is within the rtrange of the |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
401 ## object: we're using the max minimal rt of all files and the min maximal rt |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
402 rtrs <- lapply(object@rt[[rt]], range) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
403 rtr <- c(max(unlist(lapply(rtrs, "[", 1))), |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
404 min(unlist(lapply(rtrs, "[", 2)))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
405 ## Check if we've got a range which is completely off: |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
406 if (any(rtrange[, "rtmin"] >= rtr[2] | rtrange[, "rtmax"] <= rtr[1])) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
407 outs <- which(rtrange[, "rtmin"] >= rtr[2] | |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
408 rtrange[, "rtmax"] <= rtr[1]) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
409 stop(length(outs), " of the specified 'rtrange' are completely outside ", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
410 "of the retention time range of 'object' which is (", rtr[1], ", ", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
411 rtr[2], "). The first was: (", rtrange[outs[1], "rtmin"], ", ", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
412 rtrange[outs[1], "rtmax"], "!") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
413 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
414 lower_rt_outside <- rtrange[, "rtmin"] < rtr[1] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
415 upper_rt_outside <- rtrange[, "rtmax"] > rtr[2] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
416 if (any(lower_rt_outside) | any(upper_rt_outside)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
417 ## Silently fix these ranges. |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
418 rtrange[lower_rt_outside, "rtmin"] <- rtr[1] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
419 rtrange[upper_rt_outside, "rtmax"] <- rtr[2] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
420 } |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
421 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
422 if (missing(groupidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
423 gnames <- character(0) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
424 else |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
425 gnames <- groupidx |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
426 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
427 eic <- vector("list", length(sampleidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
428 names(eic) <- sampleidx |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
429 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
430 for (i in seq(along = sampidx)) { |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
431 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
432 ## cat(sampleidx[i], "") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
433 flush.console() |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
434 ## getXcmsRaw takes care of rt correction, susetting to scanrage and other |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
435 ## stuff. |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
436 lcraw <- getXcmsRaw(object, sampleidx = sampidx[i], rt=rt) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
437 currenteic <- xcms::getEIC(lcraw, mzrange, rtrange, step = prof$step) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
438 eic[[i]] <- currenteic@eic[[1]] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
439 rm(lcraw) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
440 gc() |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
441 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
442 ## cat("\n") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
443 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
444 invisible(new("xcmsEIC", eic = eic, mzrange = mzrange, rtrange = rtrange, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
445 rt = rt, groupnames = gnames)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
446 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
447 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
448 #@TODO: remove this function as soon as we can use xcms 3.x.x from Bioconductor 3.7 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
449 # https://github.com/sneumann/CAMERA/issues/33#issuecomment-405168524 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
450 # https://github.com/sneumann/xcms/commit/950a3fe794cdb6b0fda88696e31aab3d97a3b7dd |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
451 ############################################################ |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
452 ## diffreport |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
453 diffreport = function(object, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
454 class1 = levels(sampclass(object))[1], |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
455 class2 = levels(sampclass(object))[2], |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
456 filebase = character(), |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
457 eicmax = 0, eicwidth = 200, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
458 sortpval = TRUE, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
459 classeic = c(class1,class2), |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
460 value = c("into","maxo","intb"), |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
461 metlin = FALSE, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
462 h = 480, w = 640, mzdec=2, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
463 missing = numeric(), ...) { |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
464 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
465 if ( nrow(object@groups)<1 || length(object@groupidx) <1) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
466 stop("No group information. Use group().") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
467 } |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
468 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
469 if (!is.numeric(w) || !is.numeric(h)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
470 stop("'h' and 'w' have to be numeric") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
471 ## require(multtest) || stop("Couldn't load multtest") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
472 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
473 value <- match.arg(value) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
474 groupmat <- groups(object) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
475 if (length(groupmat) == 0) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
476 stop("No group information found") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
477 samples <- sampnames(object) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
478 n <- length(samples) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
479 classlabel <- sampclass(object) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
480 classlabel <- levels(classlabel)[as.vector(unclass(classlabel))] |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
481 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
482 values <- groupval(object, "medret", value=value) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
483 indecies <- groupval(object, "medret", value = "index") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
484 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
485 if (!all(c(class1,class2) %in% classlabel)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
486 stop("Incorrect Class Labels") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
487 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
488 ## c1 and c2 are column indices of class1 and class2 resp. |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
489 c1 <- which(classlabel %in% class1) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
490 c2 <- which(classlabel %in% class2) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
491 ceic <- which(classlabel %in% classeic) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
492 if (length(intersect(c1, c2)) > 0) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
493 stop("Intersecting Classes") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
494 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
495 ## Optionally replace NA values with the value provided with missing |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
496 if (length(missing)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
497 if (is.numeric(missing)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
498 ## handles NA, Inf and -Inf |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
499 values[, c(c1, c2)][!is.finite(values[, c(c1, c2)])] <- missing[1] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
500 } else |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
501 stop("'missing' should be numeric") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
502 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
503 ## Check against missing Values |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
504 if (any(is.na(values[, c(c1, c2)]))) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
505 warning("`NA` values in xcmsSet. Use fillPeaks() on the object to fill", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
506 "-in missing peak values. Note however that this will also ", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
507 "insert intensities of 0 for peaks that can not be filled in.") |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
508 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
509 mean1 <- rowMeans(values[,c1,drop=FALSE], na.rm = TRUE) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
510 mean2 <- rowMeans(values[,c2,drop=FALSE], na.rm = TRUE) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
511 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
512 ## Calculate fold change. |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
513 ## For foldchange <1 set fold to 1/fold |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
514 ## See tstat to check which was higher |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
515 fold <- mean2 / mean1 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
516 fold[!is.na(fold) & fold < 1] <- 1/fold[!is.na(fold) & fold < 1] |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
517 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
518 testval <- values[,c(c1,c2)] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
519 ## Replace eventual infinite values with NA (CAMERA issue #33) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
520 testval[is.infinite(testval)] <- NA |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
521 testclab <- c(rep(0,length(c1)),rep(1,length(c2))) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
522 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
523 if (min(length(c1), length(c2)) >= 2) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
524 tstat <- mt.teststat(testval, testclab, ...) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
525 pvalue <- xcms:::pval(testval, testclab, tstat) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
526 } else { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
527 message("Too few samples per class, skipping t-test.") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
528 tstat <- pvalue <- rep(NA,nrow(testval)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
529 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
530 stat <- data.frame(fold = fold, tstat = tstat, pvalue = pvalue) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
531 if (length(levels(sampclass(object))) >2) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
532 pvalAnova<-c() |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
533 for(i in 1:nrow(values)){ |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
534 var<-as.numeric(values[i,]) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
535 ano<-summary(aov(var ~ sampclass(object)) ) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
536 pvalAnova<-append(pvalAnova, unlist(ano)["Pr(>F)1"]) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
537 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
538 stat<-cbind(stat, anova= pvalAnova) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
539 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
540 if (metlin) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
541 neutralmass <- groupmat[,"mzmed"] + ifelse(metlin < 0, 1, -1) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
542 metlin <- abs(metlin) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
543 digits <- ceiling(-log10(metlin))+1 |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
544 metlinurl <- |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
545 paste("http://metlin.scripps.edu/simple_search_result.php?mass_min=", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
546 round(neutralmass - metlin, digits), "&mass_max=", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
547 round(neutralmass + metlin, digits), sep="") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
548 values <- cbind(metlin = metlinurl, values) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
549 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
550 twosamp <- cbind(name = groupnames(object), stat, groupmat, values) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
551 if (sortpval) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
552 tsidx <- order(twosamp[,"pvalue"]) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
553 twosamp <- twosamp[tsidx,] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
554 rownames(twosamp) <- 1:nrow(twosamp) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
555 values<-values[tsidx,] |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
556 } else |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
557 tsidx <- 1:nrow(values) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
558 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
559 if (length(filebase)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
560 write.table(twosamp, paste(filebase, ".tsv", sep = ""), quote = FALSE, sep = "\t", col.names = NA) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
561 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
562 if (eicmax > 0) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
563 if (length(unique(peaks(object)[,"rt"])) > 1) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
564 ## This looks like "normal" LC data |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
565 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
566 eicmax <- min(eicmax, length(tsidx)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
567 eics <- getEIC(object, rtrange = eicwidth*1.1, sampleidx = ceic, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
568 groupidx = tsidx[seq(length = eicmax)]) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
569 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
570 if (length(filebase)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
571 eicdir <- paste(filebase, "_eic", sep="") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
572 boxdir <- paste(filebase, "_box", sep="") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
573 dir.create(eicdir) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
574 dir.create(boxdir) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
575 if (capabilities("png")){ |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
576 xcms:::xcmsBoxPlot(values[seq(length = eicmax),], |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
577 sampclass(object), dirpath=boxdir, pic="png", width=w, height=h) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
578 png(file.path(eicdir, "%003d.png"), width = w, height = h) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
579 } else { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
580 xcms:::xcmsBoxPlot(values[seq(length = eicmax),], |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
581 sampclass(object), dirpath=boxdir, pic="pdf", width=w, height=h) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
582 pdf(file.path(eicdir, "%003d.pdf"), width = w/72, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
583 height = h/72, onefile = FALSE) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
584 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
585 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
586 plot(eics, object, rtrange = eicwidth, mzdec=mzdec) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
587 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
588 if (length(filebase)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
589 dev.off() |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
590 } else { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
591 ## This looks like a direct-infusion single spectrum |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
592 if (length(filebase)) { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
593 eicdir <- paste(filebase, "_eic", sep="") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
594 boxdir <- paste(filebase, "_box", sep="") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
595 dir.create(eicdir) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
596 dir.create(boxdir) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
597 if (capabilities("png")){ |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
598 xcmsBoxPlot(values[seq(length = eicmax),], |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
599 sampclass(object), dirpath=boxdir, pic="png", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
600 width=w, height=h) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
601 png(file.path(eicdir, "%003d.png"), width = w, height = h, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
602 units = "px") |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
603 } else { |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
604 xcmsBoxPlot(values[seq(length = eicmax),], |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
605 sampclass(object), dirpath=boxdir, pic="pdf", |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
606 width=w, height=h) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
607 pdf(file.path(eicdir, "%003d.pdf"), width = w/72, |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
608 height = h/72, onefile = FALSE) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
609 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
610 } |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
611 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
612 plotSpecWindow(object, gidxs = tsidx[seq(length = eicmax)], borderwidth=1) |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
613 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
614 if (length(filebase)) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
615 dev.off() |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
616 } |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
617 } |
36
b6a8143c397c
planemo upload commit aa538ae477bc2f901c95d916e11d70eab75b5e00
lecorguille
parents:
35
diff
changeset
|
618 |
35
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
619 invisible(twosamp) |
b02797fbead1
planemo upload commit 688b63e161c4df4364d248bf78536a7649564db4
lecorguille
parents:
33
diff
changeset
|
620 } |