贝叶斯系统发育推断
以贝叶斯统计框架估计进化树,输出的分支支持度是后验概率,与最大似然法给出的自举支持度口径不同,解释时要留意。
拖动滑块完成拼图,通过后发送验证码
独立第三方 AI 评测平台 · 只用实测数据
MrBayes 是一款做系统发育分析的命令行程序,核心方法是用贝叶斯推断从分子序列数据里估计进化树,并通过马尔可夫链蒙特卡洛(MCMC)采样给出后验概率作为分支支持度。它支持多基因、多分区数据的联合分析,可同时估计树的拓扑结构与分支长度。使用者需要自己在命令行里写 nexus 格式的输入文件、设置模型和运行参数,没有图形界面,入门门槛偏高,输出结果要靠 Tracer、FigTree 这类外部工具再做诊断和绘图。
MrBayes 是一个用于系统发育推断的命令行程序,采用贝叶斯统计方法,从 DNA、RNA 或蛋白质序列比对数据中估计物种或基因之间的进化关系。它由学术团队开发并长期维护,主要面向做分子系统学、进化生物学和群体遗传学研究的人员。
它的工作方式是构建一个包含树拓扑、分支长度和核苷酸替换模型参数的联合后验分布,再用 MCMC 方法对这个分布采样。软件内置了多种核苷酸与氨基酸替换模型,支持把不同基因、不同密码子位置拆成分区分别指定模型,还支持混合模型来处理异质性。跑完后会输出参数采样日志和树文件,需要用户自己判断链是否收敛、采样是否充分。
普通研究者通常用它做多基因联合建树或物种界定分析;有一定计算基础的用户会把它挂到集群或服务器上长时间跑;不熟悉命令行和统计推断的人,往往需要先借用带图形界面的替代工具,或请人代跑。
以贝叶斯统计框架估计进化树,输出的分支支持度是后验概率,与最大似然法给出的自举支持度口径不同,解释时要留意。
通过马尔可夫链蒙特卡洛对后验分布采样,可并行跑多条链并比较链间差异,收敛情况需要用户用外部工具自行判断。
可以把多个基因或不同密码子位置拆成分区,各自指定替换模型和参数,也支持混合模型来刻画位点间的速率异质性。
内置常见的核苷酸和氨基酸替换模型,覆盖从简单到含速率异质性的多种设定,用户按数据特点在配置里选用。
通过配置文件加命令的方式启动分析,适合写脚本批量跑多个数据集,也能投到集群上长时间计算,但需要自己管理任务。
输入采用 NEXUS 格式的序列比对文件,数据块、分区和模型设定都写在同一个文件里,格式写错往往会导致运行直接报错。
当你要做多基因联合的贝叶斯建树、并且愿意花时间调参和看收敛诊断时,MrBayes 是这一路线里被广泛引用的选择。当你的数据量不大、只想快速得到一棵树,或者你完全不熟悉命令行和 MCMC 诊断时,换用带图形界面的建树工具会更省事。它本身不负责比对、不负责剪接、也不负责画图,这些环节都要靠别的软件补齐,把它当成建树流水线里的一环而不是整套方案,用起来才顺手。
手里有多基因比对数据、需要联合建树并报告分支支持度的科研人员,这是它最对口的用法。
课题组里有服务器资源、有人会写脚本,把 MrBayes 当成建树流水线的一环来调用,效率比较合适。
课程作业只要一棵树、又不打算深入统计推断细节时,用图形界面的替代工具更容易在期限内完成。
不熟悉配置文件和终端操作的话,光是让程序跑起来就会卡住,先补命令行基础再考虑使用。
它只负责建树,比对、剪接、结果绘图都要另找工具,期望一体化流程的话不适合。
MrBayes 走的是贝叶斯路线,输出的是后验概率,你看到的分支数值含义是「在该模型和数据下这棵分支成立的概率」。最大似然工具给出的通常是自举支持度,含义是重采样下的稳定程度。两者数值不能直接横向比较,报告结果时要写清用的是哪一种。
取决于序列长度、分类群数量和模型复杂度,没有统一答案。实际操作中通常先跑一个短链看日志是否稳定,再决定要不要延长。建议一开始就把任务提交到服务器或集群上,避免在个人电脑上长时间占用。
不要只看树文件里的数值。先把参数采样日志导入诊断工具,确认链已收敛、有效样本量足够,再决定舍弃多少早期样本,最后看树的拓扑是否稳定。如果多条独立运行的链给出的拓扑差异很大,说明采样还不够,需要继续跑。
可以。除了核苷酸序列,它也内置了氨基酸替换模型,用蛋白质比对数据建树时按相应模型配置即可。具体支持哪些模型、怎么写进配置文件,需要对照官方文档里的说明逐项确认。
MrBayes 自己不做绘图。它输出的是树文件,需要另外用树可视化工具打开、查看拓扑并出图。如果你要的是一张能直接放进论文的图,要提前准备好绘图环节的工具。
| 工具 | 编辑评分 | 价格 | 核心功能 | 平台 |
|---|---|---|---|---|
| MrBayes本工具 | 4.1 | 未标价 | 贝叶斯系统发育推断;MCMC 采样与收敛诊断;分区与混合模型;多种替换模型 | — |
| PAML | — | 未标价 | 一款利用最大似然法对 DNA 或蛋白质序列进行系统发育分析的软件包 | — |
| FULL STACK AI | — | 未标价 | 使用 AI 提示构建全栈应用 | CLI |
| R.O.B. | — | 未标价 | AI 网站构建工具 | Web |
| ReactAI Components | — | 未标价 | 使用 AI 快速构建 React 组件 | — |
| Go to Fluent 1 | — | 未标价 | 流畅设计系统 | Web / iOS / Android |