最近在看論文 Phased diploid genome assemblies and pan-genomes provide insights into the genetic history of apple domestication 今天的筆記記錄的是論文中Figure2圖b中的餅圖的畫法

 

ggplot2畫拼圖我個人認為相對是比較麻煩的，而且要實現上圖這種各個部分能夠分開的好像不太好實現。之前找餅狀圖的資料的時候發現了ggforce這個包，用他來做餅圖相對方便很多，今天的內容主要用這個包里的函數來實現

數據只需要兩列就可以了

install.packages("ggforce")

df1<-read.csv("pieplot/AA.csv",header = T)
df1
ggplot()+
  geom_arc_bar(data=df1,
               stat = "pie",
               aes(x0=0,y0=0,r0=0,r=2,
                   amount=value,fill=Var,
                   explode=c(0.05,0.05,0.05,0)),
               )

 

關鍵還可以把中間給弄成空心的

ggplot()+
  geom_arc_bar(data=df1,
               stat = "pie",
               aes(x0=0,y0=0,r0=1,r=2,
                   amount=value,fill=Var,
                   explode=c(0.05,0.05,0.05,0)),
               )

 

做餅狀圖用到的函數是geom_arc_bar()，這個函數里參數的作用之前錄制過視頻大家可以參考一下

library(ggforce)
df1<-read.csv("pieplot/AA.csv",header = T)
df1
df2<-edit(df1)
p1<-ggplot()+
  geom_arc_bar(data=df1,
               stat = "pie",
               aes(x0=0,y0=0,r0=0,r=2,
                   amount=value,fill=Var,
                   explode=c(0.05,0.05,0.05,0)),
               )+
  scale_fill_manual(values = c("#80c97f","#a68dc8",
                               "#ffc000","#c00000"))+
  annotate("text",x=0,y=-2.2,label="32.72%")+
  annotate("text",x=-1.6,y=1.5,label="36.76%",angle=50)+
  annotate("text",x=1.6,y=1.5,label="30.52%",angle=-50)+
  theme_void()+
  theme(legend.position = "none",
        plot.title = element_text(hjust = 0.5,face="italic"))+
  labs(title = "GDDH13")
p2<-ggplot()+
  geom_arc_bar(data=df1,
               stat = "pie",
               aes(x0=0,y0=0,r0=0,r=2,
                   amount=value,fill=Var,
                   explode=c(0.05,0.05,0.05,0)),
  )+
  scale_fill_manual(values = c("#80c97f","#a68dc8",
                               "#ffc000","#c00000"))+
  annotate("text",x=0,y=-2.2,label="32.72%")+
  annotate("text",x=-1.6,y=1.5,label="36.76%",angle=50)+
  annotate("text",x=1.6,y=1.5,label="30.52%",angle=-50)+
  theme_void()+
  theme(legend.position = "none",
        plot.title = element_text(hjust = 0.5,face="italic"))+
  labs(title = "GDDH13")
p3<-ggplot()+
  geom_arc_bar(data=df2,
               stat = "pie",
               aes(x0=0,y0=0,r0=0,r=2,
                   amount=value,fill=Var,
                   explode=c(0.05,0.05,0.05,0.05)),
  )+
  scale_fill_manual(values = c("#80c97f","#a68dc8",
                               "#ffc000","#c00000"))+
  annotate("text",x=0.5,y=-2.2,label="32.72%")+
  annotate("text",x=-1.6,y=1.5,label="16.76%",angle=50)+
  annotate("text",x=-2.2,y=0.5,label="30%",angle=80)+
  annotate("text",x=1.6,y=1.5,label="30.52%",angle=-50)+
  theme_void()+
  theme(legend.position = "none",
        plot.title = element_text(hjust = 0.5,face="italic"))+
  labs(title = "GDDH13")
p3
library(cowplot)
plot_grid(p1,p2,p3,ncol = 2,nrow=2)

 

這樣是沒有圖例的，我想到的辦法是重新做一個圖，拼接到右下角作為圖例

df4<-data.frame(x=1,y=c(1,2,3,4),label=c("A","B","C","D"))
p4<-ggplot(df4,aes(x=x,y=y))+
  geom_text(aes(x=x+0.2,y=y,label=label))+
  geom_point(aes(color=label),show.legend = F,size=4)+
  theme_void()+
  ylim(-8,8)+xlim(-1,3)+
  scale_color_manual(values = c("#80c97f","#a68dc8",
                               "#ffc000","#c00000"))
plot_grid(p1,p2,p3,p4,ncol = 2,nrow=2)

最終的效果如下