欢迎光临散文网 会员登陆 & 注册

个人学习记录:建树流程2——Bayes建树前模型选择(Mrmodeltest2的使用)

2022-05-02 21:44 作者:猫腻需要更多的学习  | 我要投稿

2023.12.8加笔置顶

替代模型是一样的,不同软件可能基于的标准不同,大多都是AIC标准,但是MEGA好像就不是。所以最佳模型是根据序列得到的,不是说换个软件就不同了,没有说贝叶斯要用的模型就必须是MrModeltest得到的。

核算替代模型和氨基酸替代模型都推荐使用 Modelfinder,操作简单不用几个软件反复捣腾,好用的软件层出不穷,说不定还会有更简单的一键式软件出现。


下面的基本上不用看了,都已经是历史了,用处不大。

  1. 用MEGA将fas文件转换为PAUP的分析文件:

分析好的序列点进去

点保存,保存为PAUP  4.0的格式,保存到本地就行了。

2.下载MrModeltest2 (https://github.com/nylander/MrModeltest2)

直接打包下载到本地,解压使用,找到他的执行文件。

使用前先看看使用手册肯定没有害处

作者让你改名字,按照他说的来就好,我懒,改成了mt2.exe。(改成什么,之后运行的时候输入的命令就是什么,相当于重命名命令的名字)

doc文件夹里有个MrModelblock文件,看一眼,右面要用到。

2.下载使用PAUP完成scores文件的输出

下载安装PAUP(window版本:http://phylosolutions.com/paup-test/paup4-setup.msi)

打开PAUP,把刚刚保存的  .nexus文件拖入打开,把上面那个MrModelblock文件也拖进去,在Output Display的框下输入命令载入。

输好了点运行或者回车,运行完后里面会有如图里提示的completed

输入exe MrModelblock ,回车运行。运行结果输出文件mrmodel.scores

出现如下报错,可以先检查一下序列数据是否有成功导入,有的时候序列命名不规范可能导致数据导入失败。


3.使用Mrmodeltest2来分析模型的选择。

将mrmodel.scores复制到Mrmodeltest2执行文件所在目录

打开命令提示符,cd  到该目录下,然后输入执行命令行mrmodeltest2 < mrmodel.scores > out

这个命令的第一部分是 mrmodeltest2_32bit.exe 这个执行文件我起名为 mt2.exe 所以使用mt2调用,< mrmodel.scores >表示分析对象是之前我们用paup得到的那个文件,后面跟输出文件名,建议起名为.txt 格式文件方便查看。

可以打开查看一下

22.8.6 补充奇怪的bug示例

PAUP应该是对序列根据各模型的计算打分,这里需要注意的一点是,生成的mrmodel.scores文件里不同的版本可能对模型的描述略有偏差,较新的MrModeltest能够识别,但是稍微旧一些的版本可能无法识别,导致输入指令后并没有显示已经完成,并输出一个空的文件。不过这个不能识别的解释只是推测,总之使用最新的版本就不会遇到这种情况(大概)。

如果mrmodel.scores文件为空或数据不全,输出的文件里会提示你不全或有问题。

这两种情况在命令框里都会直接换行,正常情况会返回program is done


个人学习记录:建树流程2——Bayes建树前模型选择(Mrmodeltest2的使用)的评论 (共 条)

分享到微博请遵守国家法律