软件的贝叶斯检验使用课件

2023/6/51MrBayes软件2023/6/521.贝叶斯算法贝叶斯树：指定树的结构和进化模型，计算序列组成的概率，从而推断出对应的进化树。理论基础：最大后验概率原理，根据先验知识求出后验分布，用于刻画所推断发育树的为真的可能性，选择具有最大后验概率的那棵树作为系统发育树的最佳估计。MrBayes(Bayesianinferenceofphylogeny)是由JohnHuelsenbeck等编写，使用马儿可夫链方法来估计参数模型的后验概率分布。优点：能以很高的计算速度处理大型数据集，同时还使用后验概率衡量树的置信度。缺点：结果有效性受模型选择的影响很多、耗时长。2023/6/532.MrBayes软件简介MrBayes是一个采用贝叶斯法进行系统发育分析的免费软件。MrBayes对输入文件格式的要求是“*.nex”格式数据可以是核酸或氨基酸序列,也可以是限制性位点或以0、1表示的形态数据。可以执行批处理或逐条执行命令。该格式可从MEGA，ClustalW中获得。2023/6/543.MrBayes主要命令说明使用“Help”命令可以显示详细的命令列表：在MrBayes命令提示符“>”后面键入“help”命令后按“Enter”键即可。使用“help”命令后跟命令名可以显示该命令的详细说明：在MrBayes命令提示符“>”后面键入“help[命令名]”后按“Enter”键即可。2023/6/554.MrBayes的主要参数：“outgroup”命令该命令用于指定外群（outgroup）。“lset”命令该命令用于设置似然模型的参数。“Prset”命令该命令用于设置系统发育模型的先验概率。“mcmc”命令该命令用于设置执行马尔可夫链蒙特卡罗分析的参数和启动分析。2023/6/56（1）执行文件（2）选择外群（3）选择模型（4）设置先验参数（5）分析及设置（6）运行分析（7）停止分析（8）总结样品替代模型参数（9）总结样品树和枝长5.MrBayes的运行步骤：2023/6/57（1）执行文件--execute<filename>（2）选择外群--outgroup2023/6/58（3）选择模型--“lset”命令用于指定DNA模型的一般类型用于设置替换类型数用于设置位点间速率变异模型2023/6/59（4）设置先验参数--“Prset”命令Revmatpr：该参数用于设置核酸数据GTR模型的替换率的先验概率。Statefreqpr：该参数用于指定状态频率的先验分布概率。Shapepr：该参数用于指定位点间速率变异的gamma形状参数的先验概率。Privarpr：该参数用于设置不变位点速率。2023/6/510（5）分析及设置--“mcmc”命令Ngen：设置分析要跑的代数。Samplefreq：对链取样的频率。默认状态下，每第100代，对链取样一次。如果分析量较小，我们也许想尽快使其收敛，可设置为每10代取样一次。Printfreq：定义链的状态输出到屏幕上的频率。默认为每100代输出一次。2023/6/511（6）运行分析然后分析就开始运行，窗口会输出每100代链的状态信息。其中第1栏为代数，2－5为其中一个分析的4个链的loglikelihoodvalues，中括号为冷链。最后一栏为运行剩余时间，在运行初始，该值可能偏大，逐渐平稳而代表真实所剩余时间。2023/6/5122023/6/513当要求的代数已经运行完毕，窗口会提示询问是否继续运行，如果回答yes，会要求输入继续运行的代数。在回答之前，我们一般要先检查theaveragestandarddeviationofsplitfrequencies的值，该值代表两个独立分析当前的相似性程度，越接近0越好。（7）停止分析2023/6/514平衡后将会出现：四条链之间的交换频率在0.1-0.8之间，认为是合理的，则进行下一项。默认状态下，bayes会同时运行两个（Nruns=2）完全独立的但由不同的随机树开始的分析。2023/6/515（8）总结样品替代模型参数图看起来很平稳，没有上升或者下降的趋势。如果有任何上升或者下降的趋势，可能需要延长分析时间以获得充分的后掩盖率分布取样。2023/6/5162023/6/517（9）总结样品树和枝长sumtburnin=2500cladecredibilitytree给出每一分枝的可信度Phylogram系统演化树sumt命令还产生3个附加文件：filename.parts文件:包含了两分列表thelistoftaxonbipartitions与其后验概率posteriorprobability(theproportionofsampledtreescontainingthem),以及与之有关的枝长thebranchlengthsassociatedwiththem。枝长值是基于包含相关两分枝的树的。Filename.con文件:包含两棵一致树，第1棵同时包含了枝的后验概率（以内部结点标签的形式）和枝长，可由treeview读取。第2棵仅包含枝长，可由多种软件读取，如MacCladeMesquite等。filename.trprobs文件:包含了mcmc搜索过程中找到的树，由后验概率分类。2023/6/518Treeview输出贝叶斯树