Mercurial > repos > ethevenot > heatmap

annotate heatmap_script.R @ 0:81ffd91ba495 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
author ethevenot
date Tue, 24 Oct 2017 09:32:23 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
1 ## Etienne Thevenot
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
2 ## CEA, MetaboHUB
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
3 ## W4M Core Development Team
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
4 ## etienne.thevenot@cea.fr
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
5 ## 2015-05-30
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
6
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
7
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
8 heatmapF <- function(proMN,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
9 obsDF,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
10 feaDF,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
11 disC, ## dissimilarity
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
12 cutSamN, ## number of sample clusters
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
13 cutVarN, ## number of variable clusters
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
14 fig.pdfC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
15 corMetC, ## correlation method
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
16 aggMetC, ## agglomeration method
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
17 colC, ## color scale
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
18 scaL,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
19 cexN) {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
20
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
21 ncaN <- 14 ## Sample and variable name truncature for display
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
22
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
23 if(aggMetC == "ward") {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
24
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
25 rvsLs <- R.Version()
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
26 aggMetC <- paste0(aggMetC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
27 ifelse(as.numeric(rvsLs[["major"]]) > 3 ||
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
28 as.numeric(rvsLs[["major"]]) == 3 && as.numeric(rvsLs[["minor"]]) > 0.3,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
29 ".D",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
30 ""))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
31
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
32 }
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
33
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
34 if(disC %in% c("euclidean", "maximum", "manhattan", "canberra", "binary", "minkowski")) {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
35
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
36 obsHcl <- hclust(dist(proMN, method = disC),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
37 method = aggMetC)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
38
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
39 feaHcl <- hclust(dist(t(proMN), method = disC),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
40 method = aggMetC)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
41
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
42 } else if(disC == "1-cor") {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
43
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
44 obsHcl <- hclust(as.dist(1-cor(t(proMN),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
45 method = corMetC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
46 use = "pairwise.complete.obs")),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
47 method = aggMetC)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
48
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
49 feaHcl <- hclust(as.dist(1-cor(proMN,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
50 method = corMetC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
51 use = "pairwise.complete.obs")),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
52 method = aggMetC)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
53
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
54 } else if(disC == "1-abs(cor)") {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
55
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
56 obsHcl <- hclust(as.dist(1-abs(cor(t(proMN),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
57 method = corMetC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
58 use = "pairwise.complete.obs"))),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
59 method = aggMetC)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
60
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
61 feaHcl <- hclust(as.dist(1-abs(cor(proMN,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
62 method = corMetC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
63 use = "pairwise.complete.obs"))),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
64 method = aggMetC)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
65
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
66 }
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
67
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
68 heaMN <- proMN <- proMN[obsHcl[["order"]], feaHcl[["order"]]]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
69
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
70 if(scaL)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
71 heaMN <- scale(heaMN)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
72
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
73 heaMN <- heaMN[, rev(1:ncol(heaMN)), drop = FALSE]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
74
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
75 switch(colC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
76 blueOrangeRed = {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
77 imaPalVn <- colorRampPalette(c("blue", "orange", "red"),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
78 space = "rgb")(5)[1:5]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
79 },
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
80 redBlackGreen = {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
81 imaPalVn <- colorRampPalette(c("red", "black", "green"),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
82 space = "rgb")(5)[1:5]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
83 })
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
84
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
85
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
86 ## figure
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
87 ##-------
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
88
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
89 pdf(fig.pdfC,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
90 width = 14,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
91 height = 14)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
92
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
93 layout(matrix(1:4, nrow = 2),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
94 widths = c(1, 9), heights = c(1, 9))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
95
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
96 ## Color scale
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
97
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
98 scaN <- length(imaPalVn)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
99
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
100 par(mar = c(0.6, 0.6, 0.6, 4.1))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
101
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
102 ylimVn <- c(0, scaN)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
103 ybottomVn <- 0:(scaN - 1)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
104 ytopVn <- 1:scaN
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
105
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
106 plot(x = 0,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
107 y = 0,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
108 bty = "n",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
109 font.axis = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
110 font.lab = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
111 type = "n",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
112 xlim = c(0, 1),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
113 ylim = ylimVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
114 xlab = "",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
115 ylab = "",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
116 xaxs = "i",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
117 yaxs = "i",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
118 xaxt = "n",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
119 yaxt = "n")
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
120
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
121 rect(xleft = 0.8,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
122 ybottom = ybottomVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
123 xright = 1,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
124 ytop = ytopVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
125 col = imaPalVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
126 border = NA)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
127
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
128 prtVn <- pretty(range(heaMN, na.rm = TRUE))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
129 axis(at = scaN / diff(range(prtVn)) * (prtVn - min(prtVn)),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
130 font = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
131 font.axis = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
132 labels = prtVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
133 las = 1,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
134 lwd = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
135 lwd.ticks = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
136 side = 4,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
137 xpd = TRUE)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
138
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
139 arrows(par("usr")[2],
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
140 par("usr")[4],
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
141 par("usr")[2],
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
142 par("usr")[3],
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
143 code = 0,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
144 lwd = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
145 xpd = TRUE)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
146
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
147 ## Feature dendrogram
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
148
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
149 par(mar = c(7.1, 0.6, 0, 0.1),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
150 lwd = 2)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
151
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
152 plot(rev(as.dendrogram(feaHcl)), horiz = TRUE,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
153 leaflab = "none",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
154 main = "", xaxs = "i", yaxs = "i",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
155 xaxt = "n", yaxt = "n", xlab = "", ylab = "")
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
156
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
157 revFeaHcl <- list(merge = cbind(feaHcl[["merge"]][, 2], feaHcl[["merge"]][, 1]),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
158 height = feaHcl[["height"]],
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
159 order = rev(feaHcl[["order"]]),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
160 labels = feaHcl[["labels"]])
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
161
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
162 if(cutVarN > 1) {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
163 cluFeaVn <- cutree(revFeaHcl, k = cutVarN)[revFeaHcl[["order"]]]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
164 cutFeaVn <- which(abs(diff(cluFeaVn)) > 0)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
165 cutFeaTxtVn <- c(cutFeaVn[1] / 2, cutFeaVn + diff(c(cutFeaVn, length(cluFeaVn))) / 2) + 0.5
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
166 cutFeaLinVn <- cutFeaVn + 0.5
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
167 text(par("usr")[1] + 0.2 * diff(par("usr")[1:2]),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
168 cutFeaTxtVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
169 labels = unique(cluFeaVn),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
170 cex = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
171 font = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
172 las = 2)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
173 }
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
174
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
175 ## Observation dendrogram
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
176
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
177 par(mar = c(0.1, 0, 0.6, 7.1),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
178 lwd = 2)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
179
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
180 plot(as.dendrogram(obsHcl), leaflab = "none",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
181 main = "", xaxs = "i", yaxs = "i",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
182 yaxt = "n", xlab = "", ylab = "")
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
183
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
184 if(cutSamN > 1) {
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
185 cluObsVn <- cutree(obsHcl, k = cutSamN)[obsHcl[["order"]]]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
186 cutObsVn <- which(abs(diff(cluObsVn)) > 0)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
187 cutObsTxtVn <- c(cutObsVn[1] / 2, cutObsVn + diff(c(cutObsVn, length(cluObsVn))) / 2) + 0.5
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
188 cutObsLinVn <- cutObsVn + 0.5
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
189 text(cutObsTxtVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
190 0.8 * par("usr")[4],
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
191 labels = unique(cluObsVn),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
192 cex = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
193 font = 2)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
194 }
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
195
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
196 ## Heatmap
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
197
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
198 par(mar = c(7.1, 0, 0, 7.1))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
199
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
200 image(x = 1:nrow(heaMN),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
201 y = 1:ncol(heaMN),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
202 z = round(heaMN),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
203 col = imaPalVn,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
204 font.axis = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
205 font.lab = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
206 xaxt = "n",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
207 yaxt = "n",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
208 xlab = "",
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
209 ylab = "")
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
210
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
211 obsOrdVc <- obsHcl[["labels"]][obsHcl[["order"]]]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
212 obsOrdLenVn <- sapply(obsOrdVc, nchar)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
213 obsOrdVc <- substr(obsOrdVc, 1, ncaN)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
214 obsOrdVc <- paste0(obsOrdVc, ifelse(obsOrdLenVn > ncaN, ".", ""), " ")
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
215
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
216 mtext(obsOrdVc,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
217 at = 1:nrow(heaMN),
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
218 cex = cexN,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
219 las = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
220 side = 1)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
221
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
222 feaOrdVc <- feaHcl[["labels"]][feaHcl[["order"]]]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
223 feaOrdLenVn <- sapply(feaOrdVc, nchar)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
224 feaOrdVc <- substr(feaOrdVc, 1, ncaN)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
225 feaOrdVc <- paste0(" ", feaOrdVc, ifelse(feaOrdLenVn > ncaN, ".", ""))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
226
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
227 mtext(feaOrdVc,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
228 at = ncol(heaMN):1,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
229 cex = cexN,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
230 las = 2,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
231 side = 4)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
232
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
233 if(cutVarN > 1)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
234 abline(h = cutFeaLinVn)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
235 if(cutSamN > 1)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
236 abline(v = cutObsLinVn)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
237
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
238 box()
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
239
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
240 dev.off()
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
241
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
242 ## Returning
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
243 ##----------
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
244
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
245 if(cutSamN > 1)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
246 obsDF[, "heat_clust"] <- cutree(obsHcl, k = cutSamN)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
247 obsDF <- obsDF[obsHcl[["order"]], , drop = FALSE]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
248
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
249 if(cutVarN > 1)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
250 feaDF[, "heat_clust"] <- cutree(feaHcl, k = cutVarN)
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
251 feaDF <- feaDF[feaHcl[["order"]], , drop = FALSE]
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
252
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
253 return(invisible(list(proMN = proMN,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
254 obsDF = obsDF,
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
255 feaDF = feaDF)))
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
256
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
257
81ffd91ba495 planemo upload for repository https://github.com/workflow4metabolomics/heatmap.git commit bbfc13f2e4fa9e7e5b562c96d0e570318d3482d9
ethevenot
parents:
diff changeset
258 } ## end of heatmapF