【R語言】多層餅圖

需求是推動探索的重要動力!在學習的路上,很多東西都是從需求出發的!帶着目的去探索,能夠更快地找到與之對應的解決方法。上個周,同學問我能不能畫這樣一個圖:

大概的意思就是先展示3個大類的百分比,再展示每個大類裏面的小類的百分比。


回到實驗室就Google一下,但是沒有現成的例子。也看了比較主流的幾款繪圖軟件,貌似都沒有現成的解決辦法。就在要放棄的時候,突然想到Excel說不定可以呢,那就試試吧!

可是這也不對啊!!!!!!!!!

PPT啊,PPT可以拼圖啊!Excel畫的圖在PPT裏面肯定是可以編輯的,那就可以盡情的放大縮小換顏色啊!


那能不能用R語言來實現呢?

真的是不知道怎麼實現啊,嘗試了很多遍,還是不行!就在這時候,朋友圈一個老師分享了這個圖:

簡直就是我想要的圖啊!

看到使用的是兩個數據集,猜想是把兩個數據集進行拼接得到的!簡單搜索一下,發現使用ggplot2繪製餅圖的原理是先畫百分比柱狀圖再將直角座標系轉換成極座標!那也就意味着只需要將兩個柱狀圖挨着畫,然後轉換座標系就🆗了!

最難的地方是標籤的位置!理清思路,水到渠成!

萬事俱備,只欠旋轉!


不給代碼的教程都是在耍流氓!

 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') # 隱藏圖例
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章