使用 PhyML 构建进化树
使用 PhyML 构建进化树
使用 PhyML 构建进化树
PhyML 简介
使用 PhyML 构建最大似然树。
参考文献:New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0PhyML 的下载和安装
$ wget http://www.atgc-montpellier.fr/download/binaries/phyml/PhyML-3.1.zip
$ unzip PhyML-3.1.zip
$ mv PhyML-3.1 /opt/biosoft/
$ ln -s /opt/biosoft/PhyML-3.1/PhyML-3.1_linux64 /opt/biosoft/PhyML-3.1/PhyML
$ echo ‘PATH=$PATH:/opt/biosoft/PhyML-3.1/’ >> ~/.
$ source ~/.bashrcPhyML 的使用
PhyML 的输入文件为 phylip 格式。
常用例子:
$ PhyML -i proteins.phy -d aa -b 1000 -m LG -f m -v e -a e -o tlr
常用参数:
-i seq_file_name
输入文件,phylip 格式的多序列比对结果。
-d data_type default:nt
该参数的值为 nt, aa 或 generic。
-b int
设置 bootstrap 次数。
-m model
设置替代模型。 核酸的模型有: HKY85(默认的), JC69, K80, F81, TN93, GTR ; 氨基酸的模型有:LG (默认的), WAG, JTT, MtREV, Dayhoff, DCMut, RtREV, CpREV, VT, Blosum62, MtMam, HIVw, HIVb 。
-f e,m or fA,fC,fG,fT
设置频率计算的方法。 e 表示使用比对结果中不同氨基酸或碱基出现的频率来计算; m 表示使用最大似然法计算碱基频率,或使用替换模型计算氨基酸频率; fA,fC,fG,fT 则是 4 个浮点数,表示 4 中碱基的频率,仅适合核酸序列。
-v prop_invar
设置不变位点的比例,是一个[0,1]区间的值。或者使用 e 表示程序获得其最大似然估计值。
-a gamma
gamma 分布的参数。此参数值是个正数,或者使用 e 表示程序获得其最大似然估计值。在 ProtTest 软件给出的最优模型中含有 G 时,使用该参数。
-o params
参数优化的选项。t 表示对 tree topology 进行优化; l 表示对 branch length 进行优化; r 表示对 rate parameters 优化。
params=tlr 这表示对 3 者都进行优化。 params=n 表示不进行优化。
4. PhyML 结果
PhyML 的输出结果为:
proteins.phy_phyml_tree.txt : 最大似然法构建的进化树
proteins.phy_phyml_boot_stats.txt : bootstrap 的统计信息
proteins.phy_phyml_boot_trees.txt : bootstrap 树
proteins.phy_phyml_stats.txt : 程序运行的中的参数和结果统计# 欢迎使用Markdown编辑器
原文来自:
http://www.chenlianfu.com/?p=2221
https://www.plob.org/article/9891.html
使用 PhyML 构建进化树相关推荐
- phyml:基于最大似然法构建进化树
欢迎关注"生信修炼手册"! phyml 是基于最大似然法原理构建系统发生树的软件,官网如下 http://www.atgc-montpellier.fr/phyml/ 官网提供了在 ...
- linux中mpich的运行线程,贝叶斯法构建进化树:MrBayes
1. 简介 使用贝叶斯法构建进化树的软件有很多.在这里简要介绍MrBayes的安装和使用.以下介绍是对几种贝叶斯法构建进化树软件的简介: MrBayes is a program for Bayesi ...
- 群体结构分析:用 phylip 构建进化树
用 phylip 构N-J树 在 linux 系统,可以用 conda 安装 phylip 软件 1. 将 SNP 文件转换为 phylip 格式 用 tassel 的格式转换功能将 plink 格式 ...
- linux系统发育树的构建步骤,megacc构建进化树.mao文件生成方法
基于linux版megacc构建系统进化树命令如下: /biosoft/MEGA/megacc -a .mao -d fasta -o ./ -a输入的是构建进化树的参数 .mao文件 -d输入的是多 ...
- 论文模型构建的步骤_构建进化树---核苷酸替换模型的选择
在进行系统发育分析时需要通过构建系统发育树来描述不同物种之间的进化关系,通过同源DNA的核苷酸序列或者同源蛋白质分子的氨基酸序列可以实现构建进化树.在系统发育分析中,最大似然法(ML)和贝叶斯法(BI ...
- 构建进化树常用方法:生物进化距离(NJ)、统计特征(ML)和离散特征(MP)
常用方法列表 方法名 方法名 ML,Maximum likelihood 最大似然法 NJ,Neighbor-Joining 邻接法 MP,Maximum parsimony 最大简约法 ME,Min ...
- IQtree|构建进化树的软件
1.IQtree是什么 利用最大似然法构建系统发生树的软件. 2.IQtree下载安装 conda create -n iqtree conda activate iqtree conda ins ...
- 系统进化树构建及数据分析的简介
本文引用地址:http://blog.sciencenet.cn/blog-344534-352797.html 1.涉及基本概念 例如,"分子进化与生物进化是不是一个概念",&q ...
- 系统进化树的构建步骤和常用软件
系统进化树的构建步骤和常用软件 系统发生树(phylogenetic tree 或 evolutionary tree)又名分子进化树,被认为具有共同祖先的各物种间演化关系的树,它用来表示系统发生研究 ...
- linux系统发育树的构建步骤,手把手教你用 Fast Tree 快速构建序列进化树
常见的建树方法有:贝叶斯法(Bayesian),最大似然法(Maximum likelihood,ML),最大简约法(Maximum parsimony,MP),邻接法(Neighbor-Joinin ...
最新文章
- MySQL- SQL执行计划 统计SQL执行每阶段的耗时
- Algorithm之PGM之BNet:贝叶斯网络BNet的相关论文、过程原理、关键步骤等相关配图
- 刺猬文│从启动方式来看播客链的运行机制—设置验证者
- LOJ.6435.[PKUSC2018]星际穿越(倍增)
- ubuntu开机出现错误“Error found when loading /root/.profile”解决
- extjs中元数据_json – 如何配置ExtJS 4 Store(代理和阅读器)来读取元数据
- mysql ddl 失败_PHP加载mysql DDL会失败
- mysql time转换输出_MySQL将timediff输出转换为日,时,分,秒格式?
- dispatcherServlet 真正处理请求的源码解析
- [转贴] 电脑族请关爱自己的身体--远离“电脑病”完全实用手册
- 全网最全Linux命令总结!!(史上最全,建议收藏)
- 删除脚注分割线(word左下角的横线)
- uni-app 压缩图片
- 纽约市建筑能源得分预测代码分析
- 返利网的制作思路与原理
- 社群运营普遍存在的五个问题
- ps aux指令詳解 http://blog.csdn.net/hanner_cheung/article/details/6081440
- 评价的等级优良差_小学生期末评语-等级优良合格(最佳版本)
- Python从入门到精通 | Python作用与定义 | 三种注释的语法(三)
- 3D Object Detection 3D目标检测综述
热门文章
- 苹果27寸一体机拆机图解_苹果新iMac一体电脑拆机图解
- mro python_Python的mro
- python编程最大值_python求最大值最小值方法总结
- 今日头条、抖音推荐算法原理全文详解!
- 计算机软件毕业论文周记,毕业论文写作周记范文共12周
- 【MV-3D】----多视图的3D目标检测网络
- 互联网+创新创业大赛反思总结
- 海淀朝阳房山二手房价格反弹明显 短期内仍将持续
- linux软防火墙DDOS,Linux iptables防火墙详解 + 配置抗DDOS***策略实战
- u盘写保护+计算机管理,电脑如何去除u盘写保护?