Question

我想运行此代码（如下）以为某些orthofinder数据生成损益树。我一直在遵循github指南 https://github.com/guyleonard/orthomcl_tools

，除了用于绘制数据的最后一部分以外，所有内容均已通过。

我敢肯定，您可以在ggplot2中使用多种方法来执行此操作，但是如果有人知道如何修复此代码以使subset参数不再存在，这将非常有帮助。

这是我收到的错误消息：

警告：忽略未知参数：子集警告：忽略未知参数：子集警告：忽略未知参数：子集警告：忽略未知参数：子集警告：忽略未知参数：子集警告信息：名称标准化后的重复美学：尺寸

我不确定这种美学错误是否会让我感到悲伤，但我知道R不再使用subset，所以这绝对是一个问题。

library(ggplot2)
library(ggtree)

user_report <- read.csv("outfile.phy_newstyle_report.txt", sep = "\t", header=TRUE, as.is=1, row.names=NULL)

user_tree <- read.tree("tree_4code_internal_labels.tree")

user_p <- ggplot(user_tree, aes(x, y), ladderize=TRUE) + geom_tree() + 
           theme_tree() + 
           geom_tiplab(size=3, align=TRUE, color="purple", x=13) + 
           xlab("") + ylab("") + 
           geom_text(aes(label=Gain, x=branch), size=3, color="springgreen4", vjust=-0.6, subset=.(!isTip)) + 
           geom_text(aes(label=Gain), size=3, color="springgreen4", hjust=0, subset=.(isTip), x=13.5) + 
           geom_text(aes(label=Loss, x=branch), size=3, color="firebrick3", vjust=1.3, subset=.(!isTip)) + 
           geom_text(aes(label=Loss), size=3, color="firebrick3", hjust=0, subset=.(isTip), x=14) + 
           geom_text(aes(label=node), size=2, hjust=-1.5, subset=.(!isTip), color="grey") + 
           scale_x_continuous(expand = c(.1, .2))

user_p <- user_p %<+% user_report

print(user_p)

如何使用子集参数修复旧的ggplot2代码

0 个答案: