外语网站建设目的,网络推广和网络运营,商标设计网图大全,网站费用预算1. 简介使用贝叶斯法构建进化树的软件有很多。在这里简要介绍MrBayes的安装和使用。以下介绍是对几种贝叶斯法构建进化树软件的简介#xff1a;MrBayes is a program for Bayesian inference and model choice across a wide range of phylogenetic and evolutionary models. …1. 简介使用贝叶斯法构建进化树的软件有很多。在这里简要介绍MrBayes的安装和使用。以下介绍是对几种贝叶斯法构建进化树软件的简介MrBayes is a program for Bayesian inference and model choice across a wide range of phylogenetic and evolutionary models. MrBayes uses Markov chain Monte Carlo (MCMC) methods to estimate the posterior distribution of model parameters.BAMBE A nice program by Bret Larget and Donald Simon for the Bayesian inference of phylogeny.Mac5 A program by Paul-Michael Agapow that deals with gaps as a fifth state.Beast BEAST, written by Alexei Drummond and Andrew Rambaut, is a cross-platform program for Bayesian MCMC analysis of molecular sequences. It is particularly good for molecular clock analyses.PHASE Paul Higgs is the author of Phase, designed specifically for use with RNA sequences that have a conserved secondary structure, e.g. rRNA and tRNA.2. MrBayes的安装通过MrByes官网http://mrbayes.sourceforge.net/来下载MrBayes软件并安装。软件包中有其PDF格式的Manual。在windows系统下的MrBayes不能支持多线程运行在Linux下则能很好地进行并行运算。MrBayes的安装过程需要注意其src文件夹的源码文件中有个名为CompileInstructions.txt的文件介绍了如何进行软件的安装。$ sudo yum install openmpi* mpi*$ wget http://sourceforge.net/projects/mrbayes/files/latest/download?sourcefiles$ tar zxf mrbayes-3.*.*.tar.gz$ cd mrbayes_3.*.*/src$ autoconf$ ./configure --with-beagleno --enable-mpiyes$ make -j 8$ sudo cp mb /usr/local/bin (optional)以下是使用MrBayes的指令单线程或多线程运行MrBayes.$ ./mb$ cat ~/.mpd.confMPD_SECRETWORDmr45-j9z$ chmod 600 ~/.mpd.conf$ mpd $ mpirun -np 8 ./mbMrBayes v3.2.1 x64(Bayesian Analysis of Phylogeny)(Parallel version)(24 processors available)Distributed under the GNU General Public LicenseType help or help for informationon the commands that are available.Type about for authorship and generalinformation about the program.MrBayes 附加使用心得1. 使用多线程版本得到的树状图和单线程版本的树状图完全不一样差别太大。多线程版本的树状图完全是所有的分支都集合到一个点上而单线程的就正常了。这可能是由于不会使用多线程运行MrBayes的原因 或 软件在多线程下的运算方法不好(可能性很小)2. 在使用MrBayes 3.2.1版本中发现默认下得出的tree文件中在treeview软件中显现不出后验概率而3.1.2版本有。3. 但是在64位的Linux系统中使用3.1.2版本总是会Crash (core dumped)。幸好在此网页中找到了解决方法Bioinformatics applications at University of Canterbury HPC。需要对Mrbayes安装包中多个文件进行修改方法就是打个补丁mb_64bit_safe.patch再以64位的参数来make。步骤如下$ wget http://sourceforge.net/projects/mrbayes/files/mrbayes/3.1.2/mrbayes-3.1.2.tar.gz$ tar zxf mrbayes-3.1.2.tar.gz$ cd mrbayes-3.1.2$ wget https://technical.bestgrid.org/images/7/73/Mb_64bit-safe.patch.txt$ patch -R -p 1 Mb_64bit-safe.patch.txt$ OBJECT_MODE64 make _64BITyes至此则运行MrBayes正常了。3. MrBayes的简单教程3.1 使用MrBayes来做一个典型的 Bayesian phylogenetic analysis包括4个步骤a. Read the Nexus data fileb. Set the evolutionary modelc. Run the analysisd. Summarize the samples3.2 MrBayes分步演示1. 导入nex文件.本案例使用多线程运行的演示,使用24个CPU运行程序。$ mpd $ mpirun -np 24 mbMrBayes execute example.nex2. 设置进化模型参数.本例中设定数据为DNA数据.MrBayes lset nst6 ratesinvgamma3.1 主程序运行。以下命令中nchains的值要 设置使用CPU数。在单线程运行的时候可以不需要设置,而在多线程运行的时候不设置则会报错ngen则是运行的长度默认1,000,000次samplefreq则是取样频率每隔多少次运行次数取一次样printfreq是打印频率即每运行多少次将打印一行结果到屏幕上默认为500diagnfreq则代表每运行多少次分析一次结果得出 Average standard deviation of split frequencies默认是5,000.运行时会在输出到屏幕的最后一列看到预测的程序剩余运行时间。MrBayes mcmc nchains24 ngen2000000 samplefreq1000 printfreq500 diagnfreq50003.2 如果在设定的代数运行完毕后给出的 Average standard deviation of split frequencies的值小于0.01则根据提示输入‘no来停止运行反之则输入yes继续运行直到满足其值小于0.01为止。If you are intersted mainly in the well-supported parts of the tree, a standard deviation below 0.05 may be adequate.4.1 使用sump来对参数值进行归纳。设置的burnin值为 (ngen / samplefreq) * 0.25 。程序给出一个概括的表要确保PSRF一列中的值接近 1.0否则需要运行该多的代数。MrBayes sump burnin5004.2 使用sumt来构树。burnin值和前一个相同MrBayes sumt burnin5004. 详细的MrBayes使用教程4.1 将数据导入到MrBayesMrBayes导入的数据为Nexus文件该文件可以有4中数据类型aligned nucleotide or amino acid sequences, morphological (standard) data, restriction site (binary) data。Nexus文件中可以混合有这4种数据。Nexus数据文件通常由其它程序产生比如 Mesquite。文件以 nex 为后缀。使用 execute fielename 或 exe filename将文件中的数据导入到MrBayes中。4.2 指定模型
