科学网

 找回密码
  注册

tag 标签: 节点支持率

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

R语言设置系统发育树节点支持率的位置
Bearjazz 2014-5-2 16:32
R 语言设置系统发育树节点支持率的位置 # 作者信息 熊荣川 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz 生物信息学上构建系统发育树的方法有很多,如邻接法、最大简约法、最大似然法、贝叶斯法等等。在相当多的案例里面,针对同一个数据集而使用不同方法构建的系统发育树的拓扑结构通常大致一致。于是在发表结果时,为了节约篇幅或是更直观的对比不同方法得到的支持率差异,常常只展现某一个方法构建的系统树,而在节点处同时列出各种方法得到的支持率。其中一种方法就是,将不同的支持率展现在节点(或支 branch )的不同位置,下面是来自 R 语言 ape 包的例子。 intree = exampleTree.nwk Otr - read.tree(intree) # 读入系统发育树 plot(Otr,show.node.label = T,use.edge.length = F) bs.pars = c(NA,72,86,91) nodelabels(bs.pars, adj = c(-0.2, -0.1), frame = n, + cex = 0.8, font = 2) # 标在后 bs.nj = c(NA,81,87,91) nodelabels(bs.nj, adj = c(1.2, -0.5), frame = n, + cex = 0.8, font = 3) # 标在上 bs.ml = c(NA,77,86,90) nodelabels(bs.ml, adj = c(1.2, 1.5), frame = n, cex = 0.8) # 标在下
个人分类: 我的研究|6398 次阅读|0 个评论
分子系统发育树的节点支持率多少才可信?
Bearjazz 2014-3-15 13:27
分子系统发育树的节点支持率多少才可信? 熊荣川 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz 同时使用最大似然法( Maximum likelihood , ML )和贝叶斯方法( Bayesianinference )构建系统发育树,两种方法分析分别产生的非参数自举检验值和贝叶斯后验概率代表节点的支持率。因为非参数自举检验方法偏向低估,而贝叶斯方法则偏向高估节点支持率 ( Suzuki etal., 2002 ) ,最大似然率大于 70% 表明该支系关系得到充分解决,在 50% - 70% 之间为弱支持,否则视为未解决 ( Huelsenbecket al., 1993 ) 。 贝叶斯后验概率大于等于 95% 说明其支系支持率得到充分解决 ( Leaché et al., 2002 ) 。 Huelsenbeck,J. P., Hillis, D. M. 1993. Success of phylogenetic methods in the four-taxoncase. Systematic Biology , 42 (3): 247-264. Leaché, A. D., Reeder, T. W. 2002. Molecularsystematics of the eastern fence lizard (Sceloporus undulatus): a comparison ofparsimony, likelihood, and Bayesian approaches. Systematic biology , 51 (1): 44-68. Suzuki, Y., Glazko, G. V., Nei, M. 2002.Overcredibility of molecular phylogenies obtained by Bayesian phylogenetics. Proceedings of the National Academy ofSciences , 99 (25): 16138-16143.
个人分类: 我的研究|15556 次阅读|0 个评论
在R语言中有选择性地显示系统发育树的节点支持率
Bearjazz 2013-3-23 15:01
在 R 语言中有选择性地显示系统发育树的节点支持率 熊荣川 xiongrongchuan 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz 使用生物信息学构建的系统发育树通常应用节点支持率来讨论一定的生物学问题,然而支持率太低,其生物学意义就不到了。所以很多科学文献中通常会有选择性的在系统发育树上显示节点支持率——显示较高的,省去较低的。许多软件里面通常有这样的功能,如 mega5 之类,下面讲讲如何在 R 语言中如何实现 library(ape) Otr-read.tree(chou.nwk)# 读入系统发育树 plot(Otr,show.node.label=T)# 输出系统发育树,并显示节点信息 # 如上图,我们希望小于 0.5 的节点支持率不被显示 Otr$node.label# 查看中间节点信息 1.00000.41101.00000.7520 index=Otr$node.label0.5# 找到支持率小于 0.5 的节点位置 Otr$node.label =# 节点支持率小于 0.5 ,则直接复制为空 plot(Otr,show.node.label=T)# 输出系统发育树,并显示重新配置后的节点信息
个人分类: 我的研究|6004 次阅读|0 个评论
R语言中提取系统发育树的各种信息
Bearjazz 2013-3-23 14:42
R 语言中提取系统发育树的各种信息 熊荣川 xiongrongchuan 六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz library(ape) Otr-read.tree(chou.nwk)# 读入系统发育树 plot(Otr,show.node.label=T)# 输出系统发育树,并显示节点信息 Otr$tip.label# 末梢信息 DQ650503DQ650501DQ650500DQ650545DQ650547DQ650546 Otr$edge# 支系节点信息 78 89 91 92 83 710 104 1011 115 116 Otr$edge.length# 支长 0.076593150.000000000.001899570.000000000.000000000.07660045 0.000000000.000946960.000945320.00094532 Otr$node.label# 中间节点信息 1.00000.41101.00000.7520 Otr$root.edge# 根部支系长度,这里是无根树,所以结果 NULL NULL
个人分类: 我的研究|9398 次阅读|0 个评论
用RAxML构建系统发育树并计算节点支持率
热度 2 Bearjazz 2011-12-12 16:57
用 RAxML 构建系统发育树并计算节点支持率 熊荣川 中国科学院成都生物研究所 xiongrongchuan@126.com 刘国华 中国农业科学院兰州兽医研究所 RAxML 是用极大似然法建立进化树的软件之一,可以处理超大规模的序列数据,包括上千至上万个物种,几百至上万个已经比对好的碱基序列。作者是德国慕尼黑大学的 A.Stamatak 博士。 关于如何开始用 RAxML 构建系统发育树之前已有博文提及“ 最大似然法建立进化树的软件之一 RAxML ” http://bbs.sciencenet.cn/home.php?mod=spaceuid=508298do=blogid=467160 该博文为转载,内容翔实但没有明确如何构建带有节点支持率的系统发育树 笔者总结 bat 文件指令如下 raxmlHPC -f a -s sequence.phy -n boot2 -m GTRGAMMA -x 1234 -# 50 -n outname.trees 将以上指令复制到 txt 文件中,另存为 bat 格式,和 phy 格式的序列文件 sequence.phy 一起保存于 RAxML 程序所在的文件夹中,点击运行 bat 文件及开始运行。outname.trees表示所有输出的后缀名。以“.trees”结尾便于结果的查看。 -m GTRGAMMA 表示使用 GTRGAMMA 模型 -# 50 表示 bootstrap 50 次 生成文件中 RAxML_bipartitions.boot2 即位带节点支持率的系统发育树,可以把它改成 trees 格式即可用 treeview 等软件查看该系统发育树。 以上为适用 DNA 序列的指令 如果是氨基酸序列, bat 文件指令如下 raxmlHPC -f a -s sequence.phy -n boot2 -m PROTCATJTT -x 1234 -# 50 -m PROTCATJTT 为氨基酸替换模型 -# 50 表示 bootstrap 50 次 生成文件中 RAxML_bipartitions.boot2 即位带节点支持率的系统发育树,可以把它改成 trees 格式即可用 treeview 等软件查看该系统发育树。 另外,鉴于手工输入各种命令较为繁杂,且容易出错,笔者推荐一个服务器地址 http://phylobench.vital-it.ch/raxml-bb/index.php 可以方便的使用 RAxML 构建系统发育树。 附连续无间断构建两棵系统发育树 raxmlHPC -f a -s sequence1.phy -n boot2 -m GTRGAMMA -x 1234 -# 50 -n outname2.trees raxmlHPC -f a -s sequence2.phy -n boot2 -m GTRGAMMA -x 1234 -# 50 -n outname2.trees
个人分类: 我的研究|16716 次阅读|4 个评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-23 12:32

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部