从品牌网站建设到网络营销策划,从策略到执行的一站式服务
R语言怎样读入比对好的fasta文件做NJ树并做boostrap检验,很多新手对此不是很清楚,为了帮助大家解决这个难题,下面小编将为大家详细讲解,有这方面需求的人可以来学习下,希望你能有所收获。
在巩义等地区,都构建了全面的区域性战略布局,加强发展的系统性、市场前瞻性、产品创新能力,以专注、极致的服务理念,为客户提供成都网站制作、网站建设、外贸网站建设 网站设计制作按需求定制网站,公司网站建设,企业网站建设,品牌网站设计,网络营销推广,成都外贸网站制作,巩义网站建设费用合理。
今天推文用到的示例数据是参考链接2中提供的
usflu.fasta
,fasta文件已经比对好,R语言里读入fasta格式的数据可以使用adegenet
包中的fasta2DNAbin
函数
#install.packages("adegenet")
library(adegenet)
dna<-fasta2DNAbin(file = "usflu.fasta")
dna
library(ape)
dd<-dist.dna(dna)
用到的是
ape
包中的dist.dna()
函数
tree<-nj(dd)
用到的是
ape
包中的nj()
函数
library(ggtree)
ggtree(tree)+
geom_tiplab(size=2)
bs.tree<-boot.phylo(tree,dna,
function(x)nj(dist.dna(x)),1000,
trees = TRUE)
tree$node.label<-bs.tree$BP
这一步不知道对不对,好像是有问题,暂时还不知道如何验证
ggtree(tree)+
geom_tiplab(size=2)+
geom_nodelab(hjust=-0.2,size=2)
看完上述内容是否对您有帮助呢?如果还想对相关知识有进一步的了解或阅读更多相关文章,请关注创新互联行业资讯频道,感谢您对创新互联的支持。
成都网站建设公司地址:成都市青羊区太升南路288号锦天国际A座10层 建设咨询028-86922220
成都快上网科技有限公司-四川网站建设设计公司 | 蜀ICP备19037934号 Copyright 2020,ALL Rights Reserved cdkjz.cn | 成都网站建设 | © Copyright 2020版权所有.
专家团队为您提供成都网站建设,成都网站设计,成都品牌网站设计,成都营销型网站制作等服务,成都建网站就找快上网! | 成都网站建设哪家好? | 网站建设地图