-
Notifications
You must be signed in to change notification settings - Fork 1
/
reads_multilayer_pie.R
38 lines (38 loc) · 3.16 KB
/
reads_multilayer_pie.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
library('plotrix')
total = 3220928 # ICC #3558692 # DS2_2 #3718896 # DS2_1
assigned = 2211826 #2484110 #2560484
unassigned = 1009102 #1074582 #1158412
genus= 1741930 #1970997 #2028347
above_genus= 469896 #513113 #532137
cellular = 2199641 #2471315 #2547598
virus = 12185 #12795 #12886
archaea = 2021856 #2278316 #2344307
bacteria = 138349 #150062 #159313
eukaryota = 395 #482 #460
cellular_root = 39041 #42455 #43518
ahs = 1984877 #2238027 #2302471
anhs = 36979 #40289 #41836
bhs = 118858 #129382 #137915
bnhs = 19491 #20680 #21398
ahsgenus = 1612629 #1830660 #1879348
ahsnongenus = 372248 #407367 #423123
bhsgenus = 118858 #129382 #137915 # Salinbacter
bhsnongenus = 0 #0 #0
colors=list(NO='white',total='black',unassigned='gray84',assigned='indianred4',genus='lightsalmon4',above_genus='navajowhite4',
cellular='lightcoral',virus='dodgerblue4', archaea='lightpink3',bacteria='lightgoldenrod4',eukaryota='hotpink4',
bhypersaline='lightgoldenrod3',bnonhypersaline='steelblue2',ahypersaline ='rosybrown2', anonhypersaline='lightcyan3',
bhsgenus='lightgoldenrod2', bnongenus='dodgerblue2', ahsgenus='lightpink', ahsnongenus='thistle2')
iniR=0.15
par(mar=c(2,2,2,2))
png('/Users/Xin/Desktop/IC_project/output/01142016/ICC_assignment_summary.png', width = 2000, height = 2000)
pie(1, radius=iniR,col=c('white'), border = NA, labels='')
title(main = "ICC reads assignment summary", sub = 'Total reads: 3220928', line = -2, cex.main = 4, cex.sub = 3)
floating.pie(0,0,c(unassigned, virus, anhs,ahsnongenus,ahsgenus,bhsgenus,bhsnongenus,bnhs, eukaryota, cellular_root), radius=6*iniR, startpos=pi/2, col=as.character(colors[c('NO', 'NO', 'NO','ahsnongenus','ahsgenus','bhsgenus','NO','NO','NO', 'NO')]),border='white')
floating.pie(0,0,c(unassigned, virus, anhs,ahs,bhs,bnhs, eukaryota, cellular_root), radius=5*iniR, startpos=pi/2, col=as.character(colors[c('NO', 'NO', 'anonhypersaline','ahypersaline','bhypersaline','bnonhypersaline','NO', 'NO')]),border='white')
floating.pie(0,0,c(unassigned, virus, archaea,bacteria, eukaryota, cellular_root), radius=4*iniR, startpos=pi/2, col=as.character(colors[c('NO', 'NO', 'archaea','bacteria','eukaryota', 'NO')]),border='white')
floating.pie(0,0,c(unassigned,virus,cellular),radius=3*iniR, startpos=pi/2, col=as.character(colors[c('NO', 'virus','cellular')]),border='white')
floating.pie(0,0,c(unassigned, genus, above_genus), radius=2*iniR,startpos=pi/2, col=as.character(colors[c('unassigned','genus', 'above_genus')]), border = 'white')
floating.pie(0,0,c(unassigned, assigned), radius=1.2*iniR,startpos=pi/2, col=as.character(colors[c('unassigned','assigned')]), border = NA)
labelnames=c('Unassigned', 'Assigned', 'Assigned at genus level','Assigned above genus level','Cellular organisms', 'Virus_related', 'Archaea', 'Bacteria', 'Eukaryota', 'Hypersaline bacteria', 'Non-hypersaline bacteria', 'Hypersaline archaea', 'Non-hypersaline archaea', 'Hypersaline bacteria assigned at genus level', 'Hypersaline bacteria assigned above genus level', 'Hypersaline archaea assigned at genus level', 'Hypersaline archaea assigned above genus level')
legend(-1,0.9, labelnames, col=as.character(colors[-1:-2]), pch=19,bty='n', cex = 2)
dev.off()