跟着Nature学作图：R语言ggtree给进化树的节点添加饼状图

library(ggtree)
library(scatterpie)


tree<-read.tree("data/20220930/fig4a.nwk")

ggtree(tree)+
  geom_tiplab()+
  xlim(NA,15)

ggtree(tree)+
  geom_tiplab(aes(image=paste0("data/20220930/",label,".png")),
              geom = "image",size=0.2,offset = 0.2)+
  xlim(NA,16)+
  geom_tiplab(offset = 4)

ggtree(tree)+
  geom_nodelab(aes(label=node))+
  geom_tiplab(aes(label=node)) -> p

library(tidyverse)



list(ggplot2::ggplot_build(p)$data[[3]],
     ggplot2::ggplot_build(p)$data[[4]])%>% 
  bind_rows() %>% 
  select(y,x,node) %>% 
  arrange(node) %>% 
  write_csv("data/20220930/pie_df1.csv")

pie.df<-read_csv("data/20220930/pie_df1.csv")

pie.df

ggtree(tree,size=1)+
  #geom_tiplab(offset=2)+
  geom_scatterpie(data=pie.df,
                  aes(x=x,y=y,group=node,
                      r=0.4),
                  cols = c("Absent","Present"))+
  scale_fill_manual(values = c("Absent"="white",
                               "Present"="#55cc81"))+
  geom_tiplab(aes(image=paste0("data/20220930/",label,".png")),
              geom = "image",size=0.2,offset = 0.4)+
  xlim(NA,16)+
  geom_tiplab(offset = 4)+
  theme(legend.position = c(0.2,0.8))+
  coord_fixed()