【R语言】多层饼图

需求是推动探索的重要动力！在学习的路上，很多东西都是从需求出发的！带着目的去探索，能够更快地找到与之对应的解决方法。上个周，同学问我能不能画这样一个图：

灵魂画手小蓝哥.jpg

大概的意思就是先展示3个大类的百分比，再展示每个大类里面的小类的百分比。

---

回到实验室就Google一下，但是没有现成的例子。也看了比较主流的几款绘图软件，貌似都没有现成的解决办法。就在要放弃的时候，突然想到Excel说不定可以呢，那就试试吧！

2.png

可是这也不对啊！！！！！！！！！

PPT啊，PPT可以拼图啊！Excel画的图在PPT里面肯定是可以编辑的，那就可以尽情的放大缩小换颜色啊！

1.png

---

那能不能用R语言来实现呢？

真的是不知道怎么实现啊，尝试了很多遍，还是不行！就在这时候，朋友圈一个老师分享了这个图：

来自基因课张旭东老师朋友圈

简直就是我想要的图啊！

看到使用的是两个数据集，猜想是把两个数据集进行拼接得到的！简单搜索一下，发现使用ggplot2绘制饼图的原理是先画百分比柱状图再将直角坐标系转换成极坐标！那也就意味着只需要将两个柱状图挨着画，然后转换坐标系就🆗了！

最难的地方是标签的位置！理清思路，水到渠成！

.png

万事俱备，只欠旋转！

6.png

---

不给代码的教程都是在耍流氓！

 library(ggplot2) # 绘图
 library(ggsci) # 配色
 ​
 # 构建测试数据
 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循环构建标签的相对位置
 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
  }
 }
 ​
 # 构建标签后合并数据
 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') # 隐藏图例