RAxML-建立进化树
(2016-11-15 20:19:35)
标签:
bioinfo |
分类: Bioinformatics |
数据的输入RAxML
的数据位 PHYLIP 格式
-s sequenceFileName
-n outputFileName
-m substitutionModel
[-a weightFileName]
[-b
[-c numberOfCategories] 设定位点变化率的等级
[-d] -d 完全随机的搜索进化树,而不是从 maximum parsimony tree 开始。在 100 至 200 个分类单元间,该选可能会生成拓扑结构完全不同的局部最大似然树。
[-e likelihoodEpsilon]默认值为 0.1
[-E excludeFileName]
[-f a|b|c|d|e|g|h|i|j|m|n|o|p|s|t|w|x] 算法
-f a rapid Bootstrap
-f b draw the bipartitions using a bunch of topologies
-f c checks if
-f d rapid hill-climbing algorithm
-f e optimize the model parameters
-f g compute the per–site log Likelihoods for one ore more trees passed via -z.
-f h compute a log likelihood test (SH-test [21]) between a best tree passed via -t and a bunch of other trees passed via -z.
-f i performs a really thorough
standard
[-g groupingFileName]
[-h] program options
[-i initialRearrangementSett
[-j]
[-k] optimize branchlength and model parameters
on
[-l
sequenceSimilarityThresh
[-M]
-m GTRCAT: GTR approximation
-m GTRMIX: Search a good topology under GTRCAT
-m GTRGAMMA: General Time Reversible model of nucleotide subistution with the gamma model of rate heterogeneity.
-m GTRCAT_GAMMA:
-m GTRGAMMAI: Same as GTRGAMMA, but with estimate of proportion of invariable sites
-m GTRMIXI: Same as GTRMIX, but with estimate of proportion of invariable sites.
-m GTRCAT GAMMAI: Same as GTRCAT_GAMMA, but with estimate of proportion of invariable sites.
-n outputFileName
-o outgroupName(s) 设定外类群
[-p parsimonyRandomSeed] [-P proteinModel]
[-q multipleModelFileName]
-q multiple modelfile name
如将以下信息拷贝到另存为文件
DNA, rbcLa = 1-526
DNA, matK = 527-1472
调用方法
-m GTRGAMMA
[-r binaryConstraintTree]
-s sequenceFileName
[-t userStartingTree]
[-T numberOfThreads]
[-u multiBootstrapSearches] Specify the number of multiple BS
searches per replicate to obtain betterML
[-w workingDirectory] 将文件写入的工作目录
[-x
rapidBootstrapRandomNumb
[-y] -y 只输出简约树拓扑结构,之后推出,该树也可以用于 GARLI 等软件
[-z multipleTreesFile]
[-#|-N numberOfRun]
4 分析实例
若当前已经有比对好的序列,名为 test1.phy 文件
raxmlHPC -x 12345 -p 12345 -# 100 -m GTRGAMMA out1 -s test1.phy -f d
-q gennames -n TEST
将以上语句粘贴到记事本中,另存为 test1.bat 文件,保存到 raxmlHPC.exe 相同的文件夹,双击 test1.bat
即可运行。
-x 用快速方法进行 Bootstrap
-p 设定随机数
-# Bootstrap100 次
-m GTRGAMMA
-o out1 将 out1 序列设置为外类群
-s test1.phy
-n TEST
-q gennames
-f d

加载中…