Mercurial > repos > alvarofaure > bitlab

comparison chromeister/bin/plot_diags.R @ 0:7fdf47a0bae8 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author alvarofaure
date Wed, 12 Dec 2018 07:18:40 -0500
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:7fdf47a0bae8
1 #!/usr/bin/env Rscript
2 args = commandArgs(trailingOnly=TRUE)
3
4 if(length(args) < 1){
5 stop("USE: Rscript --vanilla plot.R <matrix>")
6 }
7
8
9 path_mat = args[1]
10
11
12
13
14 data <- as.matrix(read.csv(path_mat, sep = " "))
15 percentage <- 0.1
16
17
18
19
20 d <- col(data) - row(data)
21 groups <- split(data, d)
22 acum <- c()
23 for(i in 1:length(groups)){
24 acum <- c(acum, sum(unlist(groups[i])))
25 }
26
27
28 # reverse
29 data_r <- data
30 for(i in 1:length(data_r[,1])){
31 data_r[i,] <- rev(data_r[i,])
32 }
33 d <- col(data_r) - row(data_r)
34 groups <- split(data_r, d)
35 acum_r <- c()
36 for(i in 1:length(groups)){
37 acum_r <- c(acum, sum(unlist(groups[i]))/length(unlist(groups[i])) )
38 }
39
40 indexes <- matrix(0, nrow=length(acum), ncol = 3)
41 indexes[,1] <- c(1:length(acum))
42 indexes[,2] <- acum
43
44 #sort
45 indexes[order(indexes[,2], decreasing=TRUE),]
46
47 #crop by percentage
48 n_percent <- percentage * length(acum)
49
50 for(i in 1:n_percent){
51 indexes[i,3] <- 1
52 }
53
54 # Now resort based on first column to have them again in order
55 indexes[order(indexes[,1], decreasing=FALSE),]
56
57
58 # same in reverse
59 #put into indexes to sort
60 indexes_r <- matrix(0, nrow=length(acum_r), ncol = 3)
61 indexes_r[,1] <- c(1:length(acum_r))
62 indexes_r[,2] <- acum_r
63
64 #sort
65 indexes_r[order(indexes_r[,2], decreasing=TRUE),]
66
67 #crop by percentage
68 n_percent <- percentage * length(acum_r)
69
70 for(i in 1:n_percent){
71 indexes_r[i,3] <- 1
72 }
73
74 # Now resort based on first column to have them again in order
75 indexes_r[order(indexes_r[,1], decreasing=FALSE),]
76
77
78
79
80 finalmat <- matrix(0, nrow=length(data[,1]), ncol=length(data[1,]))
81
82
83 for(i in length(data[1,]):1){
84 x <- i
85 y <- 1
86
87 if(indexes[i,3] == 1){
88 while(x < length(data[1,]) && y <= length(data[,1])){
89
90
91 if(x < length(data[1,]) && y <= length(data[,1])){
92 finalmat[x,y] <- 1
93 }
94
95 y <- y + 1
96 x <- x + 1
97
98
99 }
100 }
101
102 }
103
104 for(i in 2:length(data[,1])){
105 x <- 1
106 y <- i
107
108 if(indexes[i,3] == 1){
109
110 while(x < length(data[1,]) && y <= length(data[,1])){
111
112 if(x < length(data[1,]) && y <= length(data[,1])){
113 finalmat[x,y] <- 1
114 }
115 y <- y + 1
116 x <- x + 1
117
118 }
119 }
120 }
121
122
123 #same for reverse
124 for(i in length(data_r[1,]):1){
125 x <- i
126 y <- 1
127
128 if(indexes[i,3] == 1){
129 while(x < length(data_r[1,]) && y <= length(data_r[,1])){
130
131
132 if(x < length(data_r[1,]) && y <= length(data_r[,1])){
133 finalmat[x,length(data_r[,1])-y] <- 1
134 }
135
136 y <- y + 1
137 x <- x + 1
138
139
140 }
141 }
142
143 }
144
145 for(i in 2:length(data_r[,1])){
146 x <- 1
147 y <- i
148
149 if(indexes[i,3] == 1){
150
151 while(x < length(data_r[1,]) && y <= length(data_r[,1])){
152
153 if(x < length(data_r[1,]) && y <= length(data_r[,1])){
154 finalmat[x,length(data_r[,1])-y] <- 1
155 }
156 y <- y + 1
157 x <- x + 1
158
159 }
160 }
161 }
162
163 # fully write last
164
165 for(i in 1:length(data[,1])){
166 for(j in 1:length(data[1,])){
167 if(finalmat[i,j] == 0){
168
169 data[i,j] <- 0
170 }
171 }
172 }
173
174 png(paste(path_mat, ".png", sep=""), width = length(data[,1]), height = length(data[,1]))
175 image((data), col = grey(seq(1, 0, length = 256)), xaxt='n', yaxt='n')
176 dev.off()
177
178
179