comparison test-data/out.test.R @ 0:eab37e8fea75 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mageck commit 71cef018eec5ee7ff7f3853599c027e80e2637fe
author iuc
date Wed, 14 Feb 2018 06:42:36 -0500
parents
children 9cd937788131
comparison
equal deleted inserted replaced
-1:000000000000 0:eab37e8fea75
1 pdf(file='output.pdf',width=4.5,height=4.5);
2 gstable=read.table('output.gene_summary.txt',header=T)
3 #
4 #
5 # parameters
6 # Do not modify the variables beginning with "__"
7
8 # gstablename='__GENE_SUMMARY_FILE__'
9 startindex=3
10 # outputfile='__OUTPUT_FILE__'
11 targetgenelist=c("ACIN1","ACTR8","AHCY","ACLY","AATF","AGBL5","AHCTF1","ABT1","ADIRF","ABCF1")
12 # samplelabel=sub('.\\w+.\\w+$','',colnames(gstable)[startindex]);
13 samplelabel='HL60.final,KBM7.final_vs_HL60.initial,KBM7.initial neg.'
14
15
16 # You need to write some codes in front of this code:
17 # gstable=read.table(gstablename,header=T)
18 # pdf(file=outputfile,width=6,height=6)
19
20
21 # set up color using RColorBrewer
22 #library(RColorBrewer)
23 #colors <- brewer.pal(length(targetgenelist), "Set1")
24
25 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
26 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
27 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
28 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
29
30 ######
31 # function definition
32
33 plotrankedvalues<-function(val, tglist, ...){
34
35 plot(val,log='y',ylim=c(max(val),min(val)),type='l',lwd=2, ...)
36 if(length(tglist)>0){
37 for(i in 1:length(tglist)){
38 targetgene=tglist[i];
39 tx=which(names(val)==targetgene);ty=val[targetgene];
40 points(tx,ty,col=colors[(i %% length(colors)) ],cex=2,pch=20)
41 # text(tx+50,ty,targetgene,col=colors[i])
42 }
43 legend('topright',tglist,pch=20,pt.cex = 2,cex=1,col=colors)
44 }
45 }
46
47
48
49 plotrandvalues<-function(val,targetgenelist, ...){
50 # choose the one with the best distance distribution
51
52 mindiffvalue=0;
53 randval=val;
54 for(i in 1:20){
55 randval0=sample(val)
56 vindex=sort(which(names(randval0) %in% targetgenelist))
57 if(max(vindex)>0.9*length(val)){
58 # print('pass...')
59 next;
60 }
61 mindiffind=min(diff(vindex));
62 if (mindiffind > mindiffvalue){
63 mindiffvalue=mindiffind;
64 randval=randval0;
65 # print(paste('Diff: ',mindiffvalue))
66 }
67 }
68 plot(randval,log='y',ylim=c(max(randval),min(randval)),pch=20,col='grey', ...)
69
70 if(length(targetgenelist)>0){
71 for(i in 1:length(targetgenelist)){
72 targetgene=targetgenelist[i];
73 tx=which(names(randval)==targetgene);ty=randval[targetgene];
74 points(tx,ty,col=colors[(i %% length(colors)) ],cex=2,pch=20)
75 text(tx+50,ty,targetgene,col=colors[i])
76 }
77 }
78
79 }
80
81
82
83
84 # set.seed(1235)
85
86
87
88 pvec=gstable[,startindex]
89 names(pvec)=gstable[,'id']
90 pvec=sort(pvec);
91
92 plotrankedvalues(pvec,targetgenelist,xlab='Genes',ylab='RRA score',main=paste('Distribution of RRA scores in \\n',samplelabel))
93
94 # plotrandvalues(pvec,targetgenelist,xlab='Genes',ylab='RRA score',main=paste('Distribution of RRA scores in \\n',samplelabel))
95
96
97 pvec=gstable[,startindex+1]
98 names(pvec)=gstable[,'id']
99 pvec=sort(pvec);
100
101 plotrankedvalues(pvec,targetgenelist,xlab='Genes',ylab='p value',main=paste('Distribution of p values in \\n',samplelabel))
102
103 # plotrandvalues(pvec,targetgenelist,xlab='Genes',ylab='p value',main=paste('Distribution of p values in \\n',samplelabel))
104
105
106
107 # you need to write after this code:
108 # dev.off()
109
110
111
112
113
114
115 # parameters
116 # Do not modify the variables beginning with "__"
117 targetmat=list(c(561.4907165816957,824.0396348113272,428.37415340969943,579.047491896501),c(3424.79939695118,3818.2871009576584,1992.3498917052,690.0506672205338),c(846.6456878299913,985.6508562937211,335.0024675413113,415.97581680707134),c(2432.636481525409,2122.257249136931,1067.465489792653,155.6333179800872),c(1308.1851773762019,2186.1913587343615,1482.5909580453515,997.3120339679854),c(405.68439208520414,268.16807081144486,170.34023773287015,109.85881269182627),c(640.8637498157573,559.4234589775174,711.6436598617687,632.2603542941043),c(946.5969148654764,470.6260845366416,663.0651476194316,457.74505288260946),c(246.9383256170808,177.59474888175154,28.39003962214503,0.0),c(568.8400715107754,612.7018836420428,564.0154538266146,270.64176251684285))
118 targetgene="ACIN1"
119 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
120
121 # set up color using RColorBrewer
122 #library(RColorBrewer)
123 #colors <- brewer.pal(length(targetgenelist), "Set1")
124
125 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
126 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
127 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
128 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
129
130
131 ## code
132
133 targetmatvec=unlist(targetmat)+1
134 yrange=range(targetmatvec[targetmatvec>0]);
135 # yrange[1]=1; # set the minimum value to 1
136 for(i in 1:length(targetmat)){
137 vali=targetmat[[i]]+1;
138 if(i==1){
139 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
140 axis(1,at=1:length(vali),labels=(collabel),las=2)
141 # lines(0:100,rep(1,101),col='black');
142 }else{
143 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
144 }
145 }
146
147
148
149
150 # parameters
151 # Do not modify the variables beginning with "__"
152 targetmat=list(c(2484.0819660289676,2349.578527705573,2172.7843657481662,910.9126552363929),c(992.1629154257711,1005.1862786707138,743.8190381001997,200.26346063614164),c(1267.0287897733551,1156.1418152202027,251.09412821363824,42.34141739164138),c(1500.738276518092,1315.977089213779,800.5991173444897,1476.2277955464156),c(1925.5309914189038,2054.7712445618654,194.94493873872918,235.16652091844063),c(351.29916561001374,781.4168950797068,227.75120674654121,624.2498158686586),c(1719.74905340467,1006.9622261595313,356.45271970026533,222.0063506480656),c(903.9706562768137,1445.6212558974576,1482.5909580453515,1055.1023468944147),c(651.152846716469,1081.552020689867,576.0023594448536,1072.2677863775127),c(285.1549712482957,408.46792242802854,99.0496937928171,44.630142656054424))
153 targetgene="ACTR8"
154 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
155
156 # set up color using RColorBrewer
157 #library(RColorBrewer)
158 #colors <- brewer.pal(length(targetgenelist), "Set1")
159
160 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
161 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
162 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
163 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
164
165
166 ## code
167
168 targetmatvec=unlist(targetmat)+1
169 yrange=range(targetmatvec[targetmatvec>0]);
170 # yrange[1]=1; # set the minimum value to 1
171 for(i in 1:length(targetmat)){
172 vali=targetmat[[i]]+1;
173 if(i==1){
174 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
175 axis(1,at=1:length(vali),labels=(collabel),las=2)
176 # lines(0:100,rep(1,101),col='black');
177 }else{
178 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
179 }
180 }
181
182
183
184
185 # parameters
186 # Do not modify the variables beginning with "__"
187 targetmat=list(c(301.3235520922712,657.1005708624807,228.38209651592223,137.32351586478285),c(1142.0897559789987,1099.311495578042,112.92926871919911,100.70391163417409),c(789.3207193831689,671.3081507730209,723.6305654800077,588.7745742702564),c(392.45555321286054,412.0198174056636,334.37157777193033,213.99581222261992),c(2009.3136376104133,2235.917888421252,2437.1271791188055,1937.9781176417478),c(1071.5359486598327,406.69197493921104,645.4002340767636,349.602784139093),c(61.7345814042702,218.44154112455442,614.4866353770946,452.5954210376801),c(651.152846716469,879.0940069646701,237.21455328725622,18.88198343140764),c(1625.6773103124485,1410.1023061211074,2146.286995434164,1986.613529510525),c(1053.8974968300413,882.6459019423052,106.6203710253891,105.85354347910344))
188 targetgene="AHCY"
189 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
190
191 # set up color using RColorBrewer
192 #library(RColorBrewer)
193 #colors <- brewer.pal(length(targetgenelist), "Set1")
194
195 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
196 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
197 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
198 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
199
200
201 ## code
202
203 targetmatvec=unlist(targetmat)+1
204 yrange=range(targetmatvec[targetmatvec>0]);
205 # yrange[1]=1; # set the minimum value to 1
206 for(i in 1:length(targetmat)){
207 vali=targetmat[[i]]+1;
208 if(i==1){
209 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
210 axis(1,at=1:length(vali),labels=(collabel),las=2)
211 # lines(0:100,rep(1,101),col='black');
212 }else{
213 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
214 }
215 }
216
217
218
219
220 # parameters
221 # Do not modify the variables beginning with "__"
222 targetmat=list(c(1268.498660759171,1411.8782536099247,1136.2324746551822,603.6512884889412),c(327.78122983695846,454.642557137284,51.73296108924205,24.031615276336996),c(132.28838872343613,241.5288584791821,123.02350502929512,65.80085135187511),c(495.34652221997754,586.0626713097802,279.4841678357833,243.74924065998954),c(1009.8013672555626,1102.8633905556771,1237.174837756142,1004.7503910773278),c(877.5129785321263,715.7068379934587,538.1489732819936,594.496387431289),c(1594.8100196103135,1108.1912330221296,605.6541786057605,127.59643349102738),c(314.5523909646148,252.1845434120872,88.95545748272109,359.9020478289517),c(512.984974049769,269.94401830026237,205.67006481820619,126.45207085882086),c(761.3931706526657,475.9539270030942,559.5992254409475,596.7851126957021))
223 targetgene="ACLY"
224 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
225
226 # set up color using RColorBrewer
227 #library(RColorBrewer)
228 #colors <- brewer.pal(length(targetgenelist), "Set1")
229
230 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
231 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
232 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
233 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
234
235
236 ## code
237
238 targetmatvec=unlist(targetmat)+1
239 yrange=range(targetmatvec[targetmatvec>0]);
240 # yrange[1]=1; # set the minimum value to 1
241 for(i in 1:length(targetmat)){
242 vali=targetmat[[i]]+1;
243 if(i==1){
244 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
245 axis(1,at=1:length(vali),labels=(collabel),las=2)
246 # lines(0:100,rep(1,101),col='black');
247 }else{
248 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
249 }
250 }
251
252
253
254
255 # parameters
256 # Do not modify the variables beginning with "__"
257 targetmat=list(c(659.9720726313648,809.832054900787,880.7221180558769,802.1982051767731),c(724.6463960072668,1086.8798631563195,695.2405258578626,307.26136674745163),c(836.3565909292796,1289.3378768815162,468.75109865008346,177.94838930811443),c(367.46774645398926,571.85509139924,300.30353022535627,116.72498848506541),c(518.8644579930328,632.2373060190355,627.7353205340956,308.9779106957614),c(405.68439208520414,259.28833336735727,324.27734146183434,166.5047629860492),c(2096.0360257735547,1960.6460276545372,1573.4390848362154,629.9716290296913),c(277.8056163192159,435.1071347602913,182.32714335110919,0.0),c(995.1026573974029,477.7298744919117,728.0467938656747,275.21921304566894),c(2185.6981559083283,1482.9161531626255,1741.8866532609427,1862.4501839161173))
258 targetgene="AATF"
259 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
260
261 # set up color using RColorBrewer
262 #library(RColorBrewer)
263 #colors <- brewer.pal(length(targetgenelist), "Set1")
264
265 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
266 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
267 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
268 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
269
270
271 ## code
272
273 targetmatvec=unlist(targetmat)+1
274 yrange=range(targetmatvec[targetmatvec>0]);
275 # yrange[1]=1; # set the minimum value to 1
276 for(i in 1:length(targetmat)){
277 vali=targetmat[[i]]+1;
278 if(i==1){
279 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
280 axis(1,at=1:length(vali),labels=(collabel),las=2)
281 # lines(0:100,rep(1,101),col='black');
282 }else{
283 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
284 }
285 }
286
287
288
289
290 # parameters
291 # Do not modify the variables beginning with "__"
292 targetmat=list(c(640.8637498157573,602.0461987091378,307.2433176885473,192.82510352679924),c(354.23890758164566,280.5997032331675,204.4082852794442,275.79139436177223),c(779.0316224824572,932.3724316291956,778.5179754161547,905.1908420753603),c(624.6951689717818,554.0956165110648,370.96318439602834,558.4489645167836),c(1133.270530064103,1394.1187787217498,639.0913363829536,1131.2024619361487),c(423.32284391499564,412.0198174056636,224.59675789963623,426.84726181303336),c(296.91393913482335,829.3674772777797,489.5704610396565,1233.0507362025292),c(684.959879390236,546.9918265557948,394.30610586312537,566.4595029422292),c(440.96129574478715,630.461358530218,434.6830511035094,457.1728715665062),c(1108.2827233052317,1969.5257650986248,1066.2037102538911,1333.7546478367033))
293 targetgene="AGBL5"
294 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
295
296 # set up color using RColorBrewer
297 #library(RColorBrewer)
298 #colors <- brewer.pal(length(targetgenelist), "Set1")
299
300 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
301 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
302 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
303 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
304
305
306 ## code
307
308 targetmatvec=unlist(targetmat)+1
309 yrange=range(targetmatvec[targetmatvec>0]);
310 # yrange[1]=1; # set the minimum value to 1
311 for(i in 1:length(targetmat)){
312 vali=targetmat[[i]]+1;
313 if(i==1){
314 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
315 axis(1,at=1:length(vali),labels=(collabel),las=2)
316 # lines(0:100,rep(1,101),col='black');
317 }else{
318 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
319 }
320 }
321
322
323
324
325 # parameters
326 # Do not modify the variables beginning with "__"
327 targetmat=list(c(196.96271209933826,301.9110730989776,423.9579250240324,34.33087896619571),c(1106.8128523194157,1056.6887558464218,1743.1484327997048,807.3478370217025),c(748.1643317803222,488.3855594248168,239.73811236478022,477.77139894622366),c(1095.053884432888,882.6459019423052,837.8216137379688,365.05167967388104),c(677.6105244611563,316.11865300951774,613.8557456077136,819.3636446598709),c(1078.8853035889126,1609.008424868669,348.88204246769334,193.96946615900578),c(1437.533824128006,1095.759600600407,320.4920028455483,161.35513114111984),c(845.1758168441753,660.6524658401157,541.3034221288985,640.8430740356532),c(551.2016196809839,740.570102836904,1103.42620664737,622.5332719203489),c(601.1772331987264,900.4053768304803,735.6174710982467,754.1349746240991))
328 targetgene="AHCTF1"
329 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
330
331 # set up color using RColorBrewer
332 #library(RColorBrewer)
333 #colors <- brewer.pal(length(targetgenelist), "Set1")
334
335 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
336 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
337 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
338 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
339
340
341 ## code
342
343 targetmatvec=unlist(targetmat)+1
344 yrange=range(targetmatvec[targetmatvec>0]);
345 # yrange[1]=1; # set the minimum value to 1
346 for(i in 1:length(targetmat)){
347 vali=targetmat[[i]]+1;
348 if(i==1){
349 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
350 axis(1,at=1:length(vali),labels=(collabel),las=2)
351 # lines(0:100,rep(1,101),col='black');
352 }else{
353 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
354 }
355 }
356
357
358
359
360 # parameters
361 # Do not modify the variables beginning with "__"
362 targetmat=list(c(487.9971672908978,367.6211301852257,312.2904358435953,441.15179471561487),c(358.6485205390935,394.2603425174884,593.0363832181406,268.35303725242983),c(1743.266989177725,1980.1814500315297,837.1907239685878,281.5132075228048),c(1597.7497615819454,1465.1566782744503,1065.57282048451,992.7345834391593),c(119.05954985109253,378.2768151181308,185.48159219801417,128.7407961232339),c(986.2834314825072,745.8979453033566,328.0626800781203,302.11173490252224),c(523.2740709504807,694.3954681276485,336.89513684945433,597.9294753279087),c(1562.4728579223624,763.6574201915316,422.0652557158894,220.28980669975581),c(30.8672907021351,179.37069637056908,238.47633282601822,184.81456510135357),c(339.5401977234861,447.5387671820139,310.3977665354523,205.98527379717427))
363 targetgene="ABT1"
364 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
365
366 # set up color using RColorBrewer
367 #library(RColorBrewer)
368 #colors <- brewer.pal(length(targetgenelist), "Set1")
369
370 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
371 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
372 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
373 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
374
375
376 ## code
377
378 targetmatvec=unlist(targetmat)+1
379 yrange=range(targetmatvec[targetmatvec>0]);
380 # yrange[1]=1; # set the minimum value to 1
381 for(i in 1:length(targetmat)){
382 vali=targetmat[[i]]+1;
383 if(i==1){
384 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
385 axis(1,at=1:length(vali),labels=(collabel),las=2)
386 # lines(0:100,rep(1,101),col='black');
387 }else{
388 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
389 }
390 }
391
392
393
394
395 # parameters
396 # Do not modify the variables beginning with "__"
397 targetmat=list(c(492.4067802483456,221.99343610218943,309.7668767660713,102.99263689858714),c(243.9985836454489,239.7529109903646,130.59418226186713,174.51530141149487),c(734.9354929079785,673.0840982618383,620.7955330709046,470.9052231529845),c(1074.4756906314647,950.1319065173708,1100.902647569846,743.8357109342404),c(702.5983312200275,1010.5141211371663,1291.4313579229083,1017.3383800315995),c(1647.7253750996879,760.1055252138966,685.7771793171477,608.2287390177673),c(951.0065278229242,864.8864270541301,606.9159581445226,769.0116888427839),c(435.0818118015233,435.1071347602913,275.69882921949727,339.8757017653375),c(89.66213013477338,209.56180368046682,208.8245136651112,304.4004601669353),c(1328.7633711776252,1571.7135276035012,1122.983789498181,1356.6419004808338))
398 targetgene="ADIRF"
399 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
400
401 # set up color using RColorBrewer
402 #library(RColorBrewer)
403 #colors <- brewer.pal(length(targetgenelist), "Set1")
404
405 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
406 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
407 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
408 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
409
410
411 ## code
412
413 targetmatvec=unlist(targetmat)+1
414 yrange=range(targetmatvec[targetmatvec>0]);
415 # yrange[1]=1; # set the minimum value to 1
416 for(i in 1:length(targetmat)){
417 vali=targetmat[[i]]+1;
418 if(i==1){
419 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
420 axis(1,at=1:length(vali),labels=(collabel),las=2)
421 # lines(0:100,rep(1,101),col='black');
422 }else{
423 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
424 }
425 }
426
427
428
429
430 # parameters
431 # Do not modify the variables beginning with "__"
432 targetmat=list(c(216.0710349149457,289.479440677255,192.42137966120518,498.36992632594104),c(1127.391046120839,1198.764554951823,371.5940741654094,370.2013115188104),c(1111.2224652768637,1038.9292809582466,948.227323379644,922.3562815584581),c(1164.137820766238,1204.0923974182756,1686.9992433247955,2089.033985093009),c(48.505742531926586,248.63264843445216,665.5887066969557,248.8988725049189),c(501.2260061632414,387.1565525622184,436.5757204116524,314.69972385679404),c(1975.5066049366465,1797.2588586833258,1628.3264947723626,1289.6966864967521),c(213.13129294331378,376.5008676293133,404.4003421732214,482.921030791153),c(2012.2533795820452,1989.0611874756173,1064.3110409457481,431.9968936579627),c(264.57677744687226,353.4135502746856,442.25372833608145,191.6807408945927))
433 targetgene="ABCF1"
434 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
435
436 # set up color using RColorBrewer
437 #library(RColorBrewer)
438 #colors <- brewer.pal(length(targetgenelist), "Set1")
439
440 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
441 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
442 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
443 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
444
445
446 ## code
447
448 targetmatvec=unlist(targetmat)+1
449 yrange=range(targetmatvec[targetmatvec>0]);
450 # yrange[1]=1; # set the minimum value to 1
451 for(i in 1:length(targetmat)){
452 vali=targetmat[[i]]+1;
453 if(i==1){
454 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
455 axis(1,at=1:length(vali),labels=(collabel),las=2)
456 # lines(0:100,rep(1,101),col='black');
457 }else{
458 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
459 }
460 }
461
462
463
464 #
465 #
466 # parameters
467 # Do not modify the variables beginning with "__"
468
469 # gstablename='__GENE_SUMMARY_FILE__'
470 startindex=9
471 # outputfile='__OUTPUT_FILE__'
472 targetgenelist=c("ACRC","AGAP3","ADCK4","AHRR","ADRBK1","ADK","ADCK1","ADARB2","ACSS2","ADNP")
473 # samplelabel=sub('.\\w+.\\w+$','',colnames(gstable)[startindex]);
474 samplelabel='HL60.final,KBM7.final_vs_HL60.initial,KBM7.initial pos.'
475
476
477 # You need to write some codes in front of this code:
478 # gstable=read.table(gstablename,header=T)
479 # pdf(file=outputfile,width=6,height=6)
480
481
482 # set up color using RColorBrewer
483 #library(RColorBrewer)
484 #colors <- brewer.pal(length(targetgenelist), "Set1")
485
486 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
487 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
488 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
489 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
490
491 ######
492 # function definition
493
494 plotrankedvalues<-function(val, tglist, ...){
495
496 plot(val,log='y',ylim=c(max(val),min(val)),type='l',lwd=2, ...)
497 if(length(tglist)>0){
498 for(i in 1:length(tglist)){
499 targetgene=tglist[i];
500 tx=which(names(val)==targetgene);ty=val[targetgene];
501 points(tx,ty,col=colors[(i %% length(colors)) ],cex=2,pch=20)
502 # text(tx+50,ty,targetgene,col=colors[i])
503 }
504 legend('topright',tglist,pch=20,pt.cex = 2,cex=1,col=colors)
505 }
506 }
507
508
509
510 plotrandvalues<-function(val,targetgenelist, ...){
511 # choose the one with the best distance distribution
512
513 mindiffvalue=0;
514 randval=val;
515 for(i in 1:20){
516 randval0=sample(val)
517 vindex=sort(which(names(randval0) %in% targetgenelist))
518 if(max(vindex)>0.9*length(val)){
519 # print('pass...')
520 next;
521 }
522 mindiffind=min(diff(vindex));
523 if (mindiffind > mindiffvalue){
524 mindiffvalue=mindiffind;
525 randval=randval0;
526 # print(paste('Diff: ',mindiffvalue))
527 }
528 }
529 plot(randval,log='y',ylim=c(max(randval),min(randval)),pch=20,col='grey', ...)
530
531 if(length(targetgenelist)>0){
532 for(i in 1:length(targetgenelist)){
533 targetgene=targetgenelist[i];
534 tx=which(names(randval)==targetgene);ty=randval[targetgene];
535 points(tx,ty,col=colors[(i %% length(colors)) ],cex=2,pch=20)
536 text(tx+50,ty,targetgene,col=colors[i])
537 }
538 }
539
540 }
541
542
543
544
545 # set.seed(1235)
546
547
548
549 pvec=gstable[,startindex]
550 names(pvec)=gstable[,'id']
551 pvec=sort(pvec);
552
553 plotrankedvalues(pvec,targetgenelist,xlab='Genes',ylab='RRA score',main=paste('Distribution of RRA scores in \\n',samplelabel))
554
555 # plotrandvalues(pvec,targetgenelist,xlab='Genes',ylab='RRA score',main=paste('Distribution of RRA scores in \\n',samplelabel))
556
557
558 pvec=gstable[,startindex+1]
559 names(pvec)=gstable[,'id']
560 pvec=sort(pvec);
561
562 plotrankedvalues(pvec,targetgenelist,xlab='Genes',ylab='p value',main=paste('Distribution of p values in \\n',samplelabel))
563
564 # plotrandvalues(pvec,targetgenelist,xlab='Genes',ylab='p value',main=paste('Distribution of p values in \\n',samplelabel))
565
566
567
568 # you need to write after this code:
569 # dev.off()
570
571
572
573
574
575
576 # parameters
577 # Do not modify the variables beginning with "__"
578 targetmat=list(c(461.5394895462105,502.5931393353569,445.40817718298644,889.1697652244688),c(76.43329126242978,90.5733219296933,447.30084649112945,357.0411412484354),c(258.6972935036084,685.515730683561,533.7327448963265,560.7376897811967),c(232.23961575892122,681.9638357059259,275.69882921949727,467.47213525636494),c(1393.4376945535273,1472.2604682297203,1039.706339939889,532.7008052921368),c(2395.88970688001,2441.927797124084,2462.9936596634266,2461.5240218762324),c(495.34652221997754,605.5980936867728,1159.575396122279,1617.5565806239213),c(682.0201374186041,822.2636873225097,1572.1773052974536,1333.7546478367033),c(961.2956247236359,1097.5355480892247,959.5833392285019,905.1908420753603),c(1940.2297012770634,1289.3378768815162,942.5493154552149,1103.737758763192))
579 targetgene="ACRC"
580 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
581
582 # set up color using RColorBrewer
583 #library(RColorBrewer)
584 #colors <- brewer.pal(length(targetgenelist), "Set1")
585
586 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
587 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
588 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
589 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
590
591
592 ## code
593
594 targetmatvec=unlist(targetmat)+1
595 yrange=range(targetmatvec[targetmatvec>0]);
596 # yrange[1]=1; # set the minimum value to 1
597 for(i in 1:length(targetmat)){
598 vali=targetmat[[i]]+1;
599 if(i==1){
600 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
601 axis(1,at=1:length(vali),labels=(collabel),las=2)
602 # lines(0:100,rep(1,101),col='black');
603 }else{
604 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
605 }
606 }
607
608
609
610
611 # parameters
612 # Do not modify the variables beginning with "__"
613 targetmat=list(c(1387.5582106102636,1120.6228654438523,1214.4628060584262,1111.1761158725344),c(388.0459402554127,509.69692929062694,933.0859689144999,750.1297054113762),c(326.3113588511425,635.7892009966705,960.8451187672639,615.6670961271097),c(1328.7633711776252,1038.9292809582466,1346.3187678590552,1596.3858719281006),c(352.7690365958297,234.42506852391205,310.3977665354523,429.1359870774464),c(693.7791053051318,678.4119407282909,784.1959833405838,895.4637597016048),c(837.8264619150956,719.2587329710938,374.74852301231437,993.8789460713658),c(365.99787546817333,369.3970776740432,333.74068800254935,746.6966175147567),c(707.0079441774753,635.7892009966705,837.1907239685878,1465.3563505404536),c(486.5272963050818,673.0840982618383,784.8268731099647,734.6808098765882))
614 targetgene="AGAP3"
615 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
616
617 # set up color using RColorBrewer
618 #library(RColorBrewer)
619 #colors <- brewer.pal(length(targetgenelist), "Set1")
620
621 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
622 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
623 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
624 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
625
626
627 ## code
628
629 targetmatvec=unlist(targetmat)+1
630 yrange=range(targetmatvec[targetmatvec>0]);
631 # yrange[1]=1; # set the minimum value to 1
632 for(i in 1:length(targetmat)){
633 vali=targetmat[[i]]+1;
634 if(i==1){
635 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
636 axis(1,at=1:length(vali),labels=(collabel),las=2)
637 # lines(0:100,rep(1,101),col='black');
638 }else{
639 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
640 }
641 }
642
643
644
645
646 # parameters
647 # Do not modify the variables beginning with "__"
648 targetmat=list(c(830.4771069860158,864.8864270541301,1349.4732167059603,740.974804353724),c(1481.6299537024847,1994.38902994207,2044.082852794442,1810.9538654668238),c(1234.6916280854039,1299.9935618144214,1357.6747837079133,2232.6514954349277),c(224.89026082984142,188.25043381465665,700.2876440129107,81.24974688666317),c(812.8386551562243,845.3510046771374,946.334654071501,999.6007592323984),c(1978.4463469082782,1751.0842239740703,2659.2003779409174,2851.1794981425537),c(565.9003295391435,776.0890526132542,878.1985589783528,445.72924524444096),c(680.5502664327881,534.5601941340722,550.7667686696135,1025.9210997731484),c(161.68580843975528,333.87812789769293,275.0679394501163,465.18340999195186),c(2523.768482645998,2445.4796921017187,2153.226782897355,1516.8526689897471))
649 targetgene="ADCK4"
650 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
651
652 # set up color using RColorBrewer
653 #library(RColorBrewer)
654 #colors <- brewer.pal(length(targetgenelist), "Set1")
655
656 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
657 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
658 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
659 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
660
661
662 ## code
663
664 targetmatvec=unlist(targetmat)+1
665 yrange=range(targetmatvec[targetmatvec>0]);
666 # yrange[1]=1; # set the minimum value to 1
667 for(i in 1:length(targetmat)){
668 vali=targetmat[[i]]+1;
669 if(i==1){
670 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
671 axis(1,at=1:length(vali),labels=(collabel),las=2)
672 # lines(0:100,rep(1,101),col='black');
673 }else{
674 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
675 }
676 }
677
678
679
680
681 # parameters
682 # Do not modify the variables beginning with "__"
683 targetmat=list(c(345.4196816667499,163.38716897121142,474.42910657451245,481.2044868428432),c(415.9734889859159,372.9489726516783,212.6098522813972,349.03060282298975),c(1.469870985815957,83.46953197442323,0.0,62.9399447713588),c(351.29916561001374,150.9555365494888,288.9475143764983,416.54799812317464),c(561.4907165816957,170.49095892648148,199.3611671243962,411.97054759434855),c(251.34793857452865,221.99343610218943,1564.6066280648815,1502.5481360871656),c(736.4053638937945,893.3015868752103,1114.782222496228,459.46159683091923),c(338.07032673767014,607.3740411755903,378.5338616286004,65.22867003577186),c(1230.2820151279561,525.6804566899846,837.1907239685878,945.2435342025885))
684 targetgene="AHRR"
685 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
686
687 # set up color using RColorBrewer
688 #library(RColorBrewer)
689 #colors <- brewer.pal(length(targetgenelist), "Set1")
690
691 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
692 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
693 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
694 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
695
696
697 ## code
698
699 targetmatvec=unlist(targetmat)+1
700 yrange=range(targetmatvec[targetmatvec>0]);
701 # yrange[1]=1; # set the minimum value to 1
702 for(i in 1:length(targetmat)){
703 vali=targetmat[[i]]+1;
704 if(i==1){
705 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
706 axis(1,at=1:length(vali),labels=(collabel),las=2)
707 # lines(0:100,rep(1,101),col='black');
708 }else{
709 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
710 }
711 }
712
713
714
715
716 # parameters
717 # Do not modify the variables beginning with "__"
718 targetmat=list(c(371.87735941143717,877.3180594758527,2395.4884543396593,1564.9158995424211),c(1109.7525942910477,1138.3823403320275,970.308465307979,999.0285779162951),c(1462.5216308868773,1209.420239884728,1537.4783679814984,1519.14139425416),c(586.4785233405669,987.4268037825386,743.8190381001997,1312.0117578247794),c(1018.6205931704583,717.4827854822763,1070.619938639558,1144.3626322065236),c(1269.9685317449869,1212.9721348623632,1591.1039983788835,1624.9949377332637),c(1321.4140162485455,1795.4829111945082,1478.8056194290655,1237.056005415252),c(908.3802692342615,832.9193722554148,1639.6825106212207,1268.5259778009315),c(923.078979092421,758.3295777250792,1479.4365091984464,1275.964334910274),c(680.5502664327881,634.013253507853,318.5993335374053,631.1159916618979))
719 targetgene="ADRBK1"
720 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
721
722 # set up color using RColorBrewer
723 #library(RColorBrewer)
724 #colors <- brewer.pal(length(targetgenelist), "Set1")
725
726 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
727 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
728 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
729 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
730
731
732 ## code
733
734 targetmatvec=unlist(targetmat)+1
735 yrange=range(targetmatvec[targetmatvec>0]);
736 # yrange[1]=1; # set the minimum value to 1
737 for(i in 1:length(targetmat)){
738 vali=targetmat[[i]]+1;
739 if(i==1){
740 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
741 axis(1,at=1:length(vali),labels=(collabel),las=2)
742 # lines(0:100,rep(1,101),col='black');
743 }else{
744 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
745 }
746 }
747
748
749
750
751 # parameters
752 # Do not modify the variables beginning with "__"
753 targetmat=list(c(1472.810727787589,1829.225913482041,1263.0413183007631,1315.444845721399),c(208.7216799858659,65.71005708624807,292.1019632234033,350.17496545519623),c(1011.2712382413785,1166.7975001531076,652.9709113093356,860.5606994193058),c(557.0811036242477,685.515730683561,875.0441101314478,1019.6271052960126),c(363.0581334965414,825.8155823001447,736.8792506370087,349.602784139093),c(1505.14788947554,451.09066215964896,653.6018010787167,991.0180394908496),c(198.43258308515422,28.41515982108025,249.83234867487624,114.43626322065236),c(438.02155377315523,74.58979453033565,254.87946682992424,231.16125170571777),c(804.0194292413286,472.4020320254591,1336.2245315489592,1203.2973077651598),c(454.19013461713075,490.1615069136343,896.4943622904019,685.4732166917076))
754 targetgene="ADK"
755 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
756
757 # set up color using RColorBrewer
758 #library(RColorBrewer)
759 #colors <- brewer.pal(length(targetgenelist), "Set1")
760
761 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
762 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
763 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
764 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
765
766
767 ## code
768
769 targetmatvec=unlist(targetmat)+1
770 yrange=range(targetmatvec[targetmatvec>0]);
771 # yrange[1]=1; # set the minimum value to 1
772 for(i in 1:length(targetmat)){
773 vali=targetmat[[i]]+1;
774 if(i==1){
775 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
776 axis(1,at=1:length(vali),labels=(collabel),las=2)
777 # lines(0:100,rep(1,101),col='black');
778 }else{
779 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
780 }
781 }
782
783
784
785
786 # parameters
787 # Do not modify the variables beginning with "__"
788 targetmat=list(c(662.9118146029966,1008.7381736483488,1101.533537339227,1694.8010582978616),c(1547.7741480642028,1965.9738701209897,1869.9572764452857,2353.9539344488194),c(1459.5818889152454,1179.2291325748304,1296.4784760779562,1222.1792911965672),c(1193.5352404825571,1355.0479339677643,1622.0175970785526,1905.9359639399652),c(868.6937526172306,701.4992580829187,720.4761166331027,603.6512884889412),c(798.1399452980647,768.9852626579842,1478.8056194290655,1756.0244591209105),c(1168.5474337236858,907.5091667857504,879.4603385171149,977.8578692204745),c(809.8989131845924,687.2916781723785,678.8373918539567,865.7103312642352),c(1246.4505959719315,753.0017352586266,1301.5255942330043,1264.5207085882087),c(826.0674940285679,797.4004224790644,977.8791425405509,2066.7189137649816))
789 targetgene="ADCK1"
790 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
791
792 # set up color using RColorBrewer
793 #library(RColorBrewer)
794 #colors <- brewer.pal(length(targetgenelist), "Set1")
795
796 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
797 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
798 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
799 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
800
801
802 ## code
803
804 targetmatvec=unlist(targetmat)+1
805 yrange=range(targetmatvec[targetmatvec>0]);
806 # yrange[1]=1; # set the minimum value to 1
807 for(i in 1:length(targetmat)){
808 vali=targetmat[[i]]+1;
809 if(i==1){
810 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
811 axis(1,at=1:length(vali),labels=(collabel),las=2)
812 # lines(0:100,rep(1,101),col='black');
813 }else{
814 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
815 }
816 }
817
818
819
820
821 # parameters
822 # Do not modify the variables beginning with "__"
823 targetmat=list(c(1863.7964100146337,1585.9211075140413,1761.4442361117538,1464.211987908247),c(742.2848478370584,598.4943037315028,943.8110949939769,820.5080072920774),c(1568.3523418656262,2083.1864043829455,1810.6536381234716,1887.6261618246608),c(1018.6205931704583,513.248824268262,679.4682816233377,824.5132765048003),c(1140.6198849931827,1191.6607649965529,880.0912282864958,977.8578692204745),c(135.22813069506805,118.98848175077354,351.40560154521734,399.95473995618005),c(665.8515565746286,701.4992580829187,986.7115993118849,746.6966175147567),c(418.9132309575478,300.1351256101601,376.6411923204574,645.4205245644794),c(561.4907165816957,543.4399315781598,881.9838975946388,580.7640358448108),c(442.4311667306031,229.0972260574595,395.5678854018874,651.142337725512))
824 targetgene="ADARB2"
825 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
826
827 # set up color using RColorBrewer
828 #library(RColorBrewer)
829 #colors <- brewer.pal(length(targetgenelist), "Set1")
830
831 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
832 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
833 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
834 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
835
836
837 ## code
838
839 targetmatvec=unlist(targetmat)+1
840 yrange=range(targetmatvec[targetmatvec>0]);
841 # yrange[1]=1; # set the minimum value to 1
842 for(i in 1:length(targetmat)){
843 vali=targetmat[[i]]+1;
844 if(i==1){
845 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
846 axis(1,at=1:length(vali),labels=(collabel),las=2)
847 # lines(0:100,rep(1,101),col='black');
848 }else{
849 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
850 }
851 }
852
853
854
855
856 # parameters
857 # Do not modify the variables beginning with "__"
858 targetmat=list(c(734.9354929079785,358.74139274113816,541.9343118982795,378.7840312603593),c(595.2977492554626,591.3905137762326,1061.787481868224,887.4532212761591),c(1655.0747300287676,943.0281165621008,1069.358159100796,2038.1098479598186),c(626.1650399575977,884.4218494311227,517.3296108924205,858.2719741548927),c(680.5502664327881,747.673892792174,533.1018551269456,1016.194017399393),c(662.9118146029966,777.8650001020718,864.9498738213518,787.3214909580882),c(880.4527205037583,621.5816210861304,671.8976043907657,1040.7978139918332),c(94.07174309222125,447.5387671820139,711.6436598617687,927.5059134033875),c(399.80490814194036,806.280159923152,1147.58849050404,1059.1076161071376),c(698.1887182625796,531.0082991564371,504.0809257354195,347.8862401907832))
859 targetgene="ACSS2"
860 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
861
862 # set up color using RColorBrewer
863 #library(RColorBrewer)
864 #colors <- brewer.pal(length(targetgenelist), "Set1")
865
866 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
867 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
868 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
869 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
870
871
872 ## code
873
874 targetmatvec=unlist(targetmat)+1
875 yrange=range(targetmatvec[targetmatvec>0]);
876 # yrange[1]=1; # set the minimum value to 1
877 for(i in 1:length(targetmat)){
878 vali=targetmat[[i]]+1;
879 if(i==1){
880 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
881 axis(1,at=1:length(vali),labels=(collabel),las=2)
882 # lines(0:100,rep(1,101),col='black');
883 }else{
884 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
885 }
886 }
887
888
889
890
891 # parameters
892 # Do not modify the variables beginning with "__"
893 targetmat=list(c(408.62413405683606,523.9045092011671,483.89245311522745,701.494293542599),c(1805.0015705819953,1434.9655709645526,1712.2348341000356,2152.546111180471),c(3017.64513388016,2642.609863360463,1834.6274493599499,3573.2723190648703),c(1649.1952460855039,783.1928425685244,773.4708572611067,1332.0381038883936),c(959.82575373782,1397.6706736993847,1429.5962174173474,2811.126806015325),c(495.34652221997754,301.9110730989776,336.89513684945433,555.015876620164),c(1491.9190506031964,1331.9606166131366,2087.614246881731,1983.1804416139055),c(429.2023278582595,889.7496918975753,567.8007924429005,1132.9190058844583),c(427.7324568724435,573.6310388880576,655.4944703868597,899.4690289143276),c(690.8393633334998,767.2093151691668,1040.33722970927,993.3067647552625))
894 targetgene="ADNP"
895 collabel=c("HL60.initial","KBM7.initial","HL60.final","KBM7.final")
896
897 # set up color using RColorBrewer
898 #library(RColorBrewer)
899 #colors <- brewer.pal(length(targetgenelist), "Set1")
900
901 colors=c( "#E41A1C", "#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#A65628", "#F781BF",
902 "#999999", "#66C2A5", "#FC8D62", "#8DA0CB", "#E78AC3", "#A6D854", "#FFD92F", "#E5C494", "#B3B3B3",
903 "#8DD3C7", "#FFFFB3", "#BEBADA", "#FB8072", "#80B1D3", "#FDB462", "#B3DE69", "#FCCDE5",
904 "#D9D9D9", "#BC80BD", "#CCEBC5", "#FFED6F")
905
906
907 ## code
908
909 targetmatvec=unlist(targetmat)+1
910 yrange=range(targetmatvec[targetmatvec>0]);
911 # yrange[1]=1; # set the minimum value to 1
912 for(i in 1:length(targetmat)){
913 vali=targetmat[[i]]+1;
914 if(i==1){
915 plot(1:length(vali),vali,type='b',las=1,pch=20,main=paste('sgRNAs in',targetgene),ylab='Read counts',xlab='Samples',xlim=c(0.7,length(vali)+0.3),ylim = yrange,col=colors[(i %% length(colors))],xaxt='n',log='y')
916 axis(1,at=1:length(vali),labels=(collabel),las=2)
917 # lines(0:100,rep(1,101),col='black');
918 }else{
919 lines(1:length(vali),vali,type='b',pch=20,col=colors[(i %% length(colors))])
920 }
921 }
922
923
924
925 dev.off()
926 Sweave("output_summary.Rnw");
927 library(tools);
928
929 texi2dvi("output_summary.tex",pdf=TRUE);
930