千家信息网

请输入关键字词

热门搜索排行

最新搜索排行

导航： 首页 > 互联网科技 >

R语言ggtree如何将进化树中的序列id改成物种名称

发表于：2024-11-11 作者：千家信息网编辑
千家信息网最后更新 2024年11月11日，R语言ggtree如何将进化树中的序列id改成物种名称，很多新手对此不是很清楚，为了帮助大家解决这个难题，下面小编将为大家详细讲解，有这方面需求的人可以来学习下，希望你能有所收获。通常我们会使用比对好
千家信息网最后更新 2024年11月11日R语言ggtree如何将进化树中的序列id改成物种名称

R语言ggtree如何将进化树中的序列id改成物种名称，很多新手对此不是很清楚，为了帮助大家解决这个难题，下面小编将为大家详细讲解，有这方面需求的人可以来学习下，希望你能有所收获。

通常我们会使用比对好的fasta文件构建进化树，fasta文件中大于号后的内容就是最终进化树上的文字标签。如果拿到进化树文件后你想替换掉其中的一些内容，那该怎么办呢？本篇推文介绍一下使用R语言的ggtree包实现这个目的

这个问题是来源于公众号的一位读者的提问

大家可以关注我的公众号 小明的数据分析笔记本 留言相关问题，如果我恰巧会的话，我会抽出时间介绍对应的解决办法

首先你已经有了构建好的进化树文件
(Synergus:0.1976902387,(((((Periclistus:0.1403183720,Synophromorpha:0.0325185390)93:0.0313182375,(Xestophanes:0.0275715134,(Diastrophus:0.0456139475,Gonaspis:0.1146402107)97:0.0603746476)86:0.0275523221)91:0.0396704245,Ibalia:0.1295291852)93:0.0678466304,(((Liposthenes_ker:0.0568838340,Rhodus:0.4243267334)73:0.0825510697,Plagiotrochus:0.0778290252)71:0.0457931797,Phanacis_2:0.1416544135)42:0.0142517743)48:0.0209026386,(((Liposthenes_gle:0.1641119081,((((Antistrophus:0.1098867540,Hedickiana:0.2313789580)73:0.0566918206,Neaylax:0.1747090949)53:0.0027850349,(Isocolus:0.0980216531,Aulacidea:0.1315344980)40:0.0147148853)54:0.0123010924,((Andricus:0.0479556214,Neuroterus:0.0392025403)95:0.0395094917,Biorhiza:0.0640188941)87:0.0159496082)20:0.0000025961)50:0.0194234721,((((Panteliella:0.0792235900,Diplolepis:0.3184402599)84:0.0461941800,Phanacis_1:0.1153410113)66:0.0099961323,(Eschatocerus:0.2548694740,Parnips:0.0000022831)64:0.0802390069)34:0.0241704495,((Barbotinia:0.0731026287,Aylax:0.0957869567)87:0.0269932737,Iraella:0.0390833327)95:0.0797807340)18:0.0000021284)23:0.0095262346,Timaspis:0.0585073936)19:0.0170106400)57:0.0526944283,(Ceroptres:0.1057541047,(Pediaspis:0.1932340906,Paramblynotus:0.1711455809)28:0.0000021043)48:0.0416999011);
也准备好了需要替换的数据
image.png
  • 第一列x就是进化树中原本的序列名称
  • 第二列y是想要替换成的id名称
读入进化树文件
library(treeio)

tree<-read.newick("ggtree_practice_aligned.fasta.treefile",
                  node.label = "support")
使用ggtree进行可视化展示
ggtree(tree)+
  geom_tiplab()+
  xlim(NA,0.8)
读入已经准备好打算替换内容
df<-read.csv("pra.csv",header=T)
替换内容
df<-read.csv("pra.csv",header=T)
tree1<-tree
tree1@phylo$tip.label<-
  df[match(tree1@phylo$tip.label,df$x),]$y

这样就替换过来了

接下来可视化展示一下新的进化树
ggtree(tree1)+
  geom_tiplab()
image.png
把这个新的进化树写出到文件里
write.tree(tree1@phylo,file = "pra.nwk")

这样就达成目的了

这里导出的进化树文件没有了最初的支持率的信息，我们再通过一行代码给他加上就好了

tree1@phylo$node.label<-tree1@data$support
write.tree(tree1@phylo,file = "pra.nwk")

这样就没有问题了。

看完上述内容是否对您有帮助呢？如果还想对相关知识有进一步的了解或阅读更多相关文章，请关注行业资讯频道，感谢您对的支持。

很赞哦！
进化 文件 内容 名称 问题 序列 语言 公众 就是 数据 目的 准备 可视化 帮助 支持 物种 清楚 接下来 一行 代码 数据库的安全要保护哪些东西 数据库安全各自的含义是什么 生产安全数据库录入 数据库的安全性及管理 数据库安全策略包含哪些 海淀数据库安全审计系统 建立农村房屋安全信息数据库 易用的数据库客户端支持安全管理 连接数据库失败ssl安全错误 数据库的锁怎样保障安全 互联网金融科技元年 如何预防网络安全威胁 安卓软件开发与小程序 金蝶专用版没有加密服务器 山西网络技术转让包括什么 成广平全区网络安全和信息化 如何抵御黑客攻击服务器 2台电脑数据库数据互通 企业诚信积分管理软件开发 签订网络安全承诺书新闻 数据库基本操作视频 服务器内存 家用 路由器做打印服务器 滨州招聘hr软件开发公司 宁波海曙区hp服务器 网络技术工程师证书是什么样的 软件开发在鞍山的发展地位 什么因素对服务器安全 数据库中用户头像表 文件备份服务器管理制度 VS金仓数据库 网络安全法的规定国家采取 广州天翔软件开发有限公司 美国应用软件开发 数据库实验二思考题 黄浦区电商软件开发市场报价 安全狗拦截 数据库修改 警钟长鸣网络安全漫画 网络安全高新技术岗位 数字网络安全签名

相关文章

0