建树是分子系统发育与进化研究中绕不过去的一道坎，本文就本人实际经验介绍贝叶斯树的实际操作：
基于不同原理的系统发育树构建方法的比较
Bayesian Inference (BI）建树
需要使用的软件包括: PAUP、SequenceMatrix、MrModeltest、Mrbayes软件
数据准备
1.1 数据格式转换
将fasta数据格式转为nex文件格式，可使用MEGA等软件。
Paup中fasta转换nex格式代码：
export file=newfile.nex format=nexus interleaved=no;
1.2 为nex格式文件添加MrModelblock命令文本
2. 使用Paup计算Score值
直接使用windows-PAUP打开上一步粘贴有MrModelblock的nex文件。
3. 使用mrmodeltest计算推荐模型 （MrMTgui软件已失效）
本文使用的是linux版本的mrmodeltest，命令行如下；
4. 准备mrbayes输入文件
4.1 打开界面版windows-PAUP，运行，open， execute ----选择.nex格式文件，
转换至标准NEXUS。
Paup中fasta转换nex格式代码：
export file=newfile.nex format=nexus interleaved=no;
4.2 将3步骤得到的best-fit model粘贴到到Mrbayes block参考脚本中，复制粘贴到转换好的nexus标准文件末尾中。 请根据自己的实际情况更改Mrbayes block。
5. 运行mrBayes软件
双击软件中mrbayes.3.2.7-win64.exe，输入带有Mrbayes block的nex文件。命令如下：
From Rindy
6. 运算结束后，可进行后续可视化操作（figtree）
mrbayes运行结束后输出结果
网友评论