需求是推動探索的重要動力！在學習的路上口渔，很多東西都是從需求出發(fā)的样屠！帶著目的去探索，能夠更快地找到與之對應的解決方法缺脉。上個周痪欲，同學問我能不能畫這樣一個圖：
靈魂畫手小藍哥.jpg
大概的意思就是先展示3個大類的百分比，再展示每個大類里面的小類的百分比枪向。
回到實驗室就Google一下勤揩，但是沒有現(xiàn)成的例子。也看了比較主流的幾款繪圖軟件秘蛔，貌似都沒有現(xiàn)成的解決辦法陨亡。就在要放棄的時候，突然想到Excel說不定可以呢深员，那就試試吧负蠕！
2.png
可是這也不對啊＞氤Ｕ谔恰！５汀Ｓ恕！０鸥拧Ｈ弧！
PPT啊罢洲，PPT可以拼圖疤吖省！Excel畫的圖在PPT里面肯定是可以編輯的惹苗，那就可以盡情的放大縮小換顏色暗罱稀！
1.png
那能不能用R語言來實現(xiàn)呢桩蓉？
真的是不知道怎么實現(xiàn)啊淋纲，嘗試了很多遍，還是不行院究！就在這時候帚戳，朋友圈一個老師分享了這個圖：
來自基因課張旭東老師朋友圈
簡直就是我想要的圖扮杌颉！
看到使用的是兩個數(shù)據(jù)集片任，猜想是把兩個數(shù)據(jù)集進行拼接得到的！簡單搜索一下蔬胯，發(fā)現(xiàn)使用ggplot2繪制餅圖的原理是先畫百分比柱狀圖再將直角坐標系轉換成極坐標对供！那也就意味著只需要將兩個柱狀圖挨著畫，然后轉換坐標系就??了氛濒！
最難的地方是標簽的位置产场！理清思路，水到渠成舞竿！
.png
萬事俱備京景，只欠旋轉！
6.png
不給代碼的教程都是在耍流氓骗奖！
library(ggplot2) # 繪圖
library(ggsci) # 配色
?
# 構建測試數(shù)據(jù)
dat = data.frame(x = rep('b',7),
y = rep('a',7),
z = rep('c',7),
cat1 = paste('c',1:7, sep = '_'),
cat2 = c('a1','a1','a2','a2','a3','a4','a4'),
value1 = 1:7,
value2 = 1:7)
?
# 分別求所占百分比
dat1 = aggregate(dat$value1, by = list(dat$cat1), FUN = sum)
dat1$per1 = dat1$x / sum(dat1$x)
?
# for循環(huán)構建標簽的相對位置
for (i in seq(nrow(dat1), 1)) {
if (i == nrow(dat1)) {
dat1$per.y1[i] = dat1$per1[i] / 2
}else{
dat1$per.y1[i] = sum(dat1$per1[(i + 1):nrow(dat1)]) + dat1$per1[i] / 2
}
}
?
# 構建標簽后合并數(shù)據(jù)
dat1$label1 = paste(dat1$Group.1,'(',round(dat1$per1*100, 2),'%',')', sep = '')
dat = merge(dat, dat1[,c(1,3,4,5)], by.x = 'cat1', by.y = 'Group.1')
?
# 重復操作
dat2 = aggregate(dat$value2, by = list(dat$cat2), FUN = sum)
dat2$per2 = dat2$x / sum(dat2$x)
?
for (i in seq(nrow(dat2), 1)) {
if (i == nrow(dat2)) {
dat2$per.y2[i] = dat2$per2[i] / 2
}else{
dat2$per.y2[i] = sum(dat2$per2[(i + 1):nrow(dat2)]) + dat2$per2[i] / 2
}
}
?
dat2$label2 = paste(dat2$Group.1,'(',round(dat2$per2*100, 2),'%',')', sep = '')
dat = merge(dat, dat2[,c(1,3,4,5)], by.x = 'cat2', by.y = 'Group.1')
?
# 繪圖
ggplot(dat) +
# 繪制柱狀圖
geom_bar(aes(y,
ifelse(cat2 == 'a3', per2, per2/2),
fill = cat2),
stat = 'identity', width = 1.3) +
# 添加標簽
geom_text(aes(1.25, as.numeric(per.y2),
label = label2),
size =2.5, color = 'black') +
# 繪制柱狀圖
geom_bar(aes(x, per1, fill = cat1),
stat = 'identity', width = .8, color = 'white') +
# 添加標簽
geom_text(aes(2, as.numeric(per.y1),label = label1),
size = 2.5, color = 'black') +
# 設置Y軸刻度
scale_y_continuous(labels = scales::percent) +
coord_polar(theta = "y") + # 轉換坐標軸
theme_void() +
scale_fill_igv() + # 設置填充色
theme(legend.position = 'none') # 隱藏圖例
