R 如何将饼图添加到ggplot系统发育树的顶端？

我正试图在树的顶端用饼图绘制一个系统发育树，但我真的被如何做到这一点所困扰

这里有一个简单的虚构的例子

library("treeio")
library("ggtree")

nwk <- system.file("extdata", "sample.nwk", package="treeio")
tree <- read.tree(nwk)

ggplot(tree, aes(x, y)) + geom_tree() + theme_tree() + geom_tiplab()

library（“treeio”）
图书馆（“ggtree”）
nwk此代码应该可以使用ggplot和ggtree

library(ggtree)
library(reshape2)
library(ggplot2)

set.seed(2020-12-18)

## generate random data with 15 tips
num_tips = 15
tr <- rtree(num_tips)
p <- ggtree(tr)

a <- runif(num_tips, 0, 0.33)
b <- runif(num_tips, 0, 0.33)
c <- runif(num_tips, 0, 0.33)
d <- 1 - a - b - c
dat <- data.frame(a=a, b=b, c=c, d=d)

## melt the dataframe to create pie-charts
dat.m = melt(dat, id.vars = "tip")

## create list to hold all pie charts
pies = list()
for (i in 1:num_tips) {
  curr_dat = melt(dat[i,])
  ## create a ggplot object for each pie chart
  pies[[i]] =  ggplot(curr_dat, aes(y = value, fill = variable, x="")) + 
    geom_bar(stat = "identity") +
    coord_polar("y", start=0) +
    theme_void() + scale_fill_brewer(palette = "Set1", guide = F)
}
# give them the appropriate names and plot on tree
names(pies) = 1:15
inset(p, pies, width=0.1, height=0.1, hjust=0)

库（ggtree）
图书馆（E2）
图书馆（GG2）
种子集（2020-12-18）
##生成带有15个提示的随机数据
num_tips=15
tr这段代码应该可以使用ggplot和ggtree

library(ggtree)
library(reshape2)
library(ggplot2)

set.seed(2020-12-18)

## generate random data with 15 tips
num_tips = 15
tr <- rtree(num_tips)
p <- ggtree(tr)

a <- runif(num_tips, 0, 0.33)
b <- runif(num_tips, 0, 0.33)
c <- runif(num_tips, 0, 0.33)
d <- 1 - a - b - c
dat <- data.frame(a=a, b=b, c=c, d=d)

## melt the dataframe to create pie-charts
dat.m = melt(dat, id.vars = "tip")

## create list to hold all pie charts
pies = list()
for (i in 1:num_tips) {
  curr_dat = melt(dat[i,])
  ## create a ggplot object for each pie chart
  pies[[i]] =  ggplot(curr_dat, aes(y = value, fill = variable, x="")) + 
    geom_bar(stat = "identity") +
    coord_polar("y", start=0) +
    theme_void() + scale_fill_brewer(palette = "Set1", guide = F)
}
# give them the appropriate names and plot on tree
names(pies) = 1:15
inset(p, pies, width=0.1, height=0.1, hjust=0)

库（ggtree）
图书馆（E2）
图书馆（GG2）
种子集（2020-12-18）
##生成带有15个提示的随机数据
num_tips=15
tr若要包括尖端标签，请将最后一行更改为：“插入（p+geom_tiplab（），pies，宽度=0.1，高度=0.1，hjust=-0.2）”若要包括尖端标签，请将最后一行更改为：“插入（p+geom_tiplab（），pies，宽度=0.1，高度=0.1，hjust=-0.2）”