我从多灾多难的2020年的11月开始正式学习构建系统发育树,拖拖拉拉地学到12月。星座运势告诉我:“你想掌握地大部分技能都可以通过持续地坚持来获得。”我心想:信了你的邪。
下面进入正题。
最大似然法的原理在之前的文章有过介绍。这里只讲建树方法。最大似然法的计算强度大,非常耗时,用自己的电脑跑太耗时耗力,因此在这里给出几个线建树的方法。
IQ-TREE有本地地软件(下载地址:http://www.iqtree.org/),也有在线地建树平台IQ-TREE web server (http://iqtree.cibiv.univie.ac.at/),可以通过阅读它提供的教程(http://www.iqtree.org/doc/Web-Server-Tutorial)学习如何使用,下面我也说一下:
在页面上方有三个主要的选项:Tree Inference、Model Selection和Analysis Results。这里我们主要用的是Tree Inference。
Model Selection是当你不想建树,只要选择最佳的核苷酸替代模型。什么是最佳的核苷酸替代模型?为什么要选择最佳的核苷酸替代模型?请参考本人之前对最大似然法介绍的文章https://www.jianshu.com/p/23198baef227和https://www.jianshu.com/p/8a26a5c24232。

在Input data选项框中,点击Alignment file后的Browse…选择要进行构建系统发育树的序列文件(注意,是已经完成多重序列比对的序列文件)。
Sequence type可以选择也可以不选,默认是自动识别。
在Substitution Model Options选项框中的Substitution model一栏中提供了核苷酸替代模型选项,默认是Auto,就是它会自动选择最佳的核苷酸替代模型。
其他的选项默认就好,其实到头来发现,我们除了选择要分析的序列外,别的都是默认就好,还是很方便的。
在最下面的Email中给出自己的邮箱,任务结束时会发邮件提醒。点击邮件中的链接,进入下方的网页:

点击Full Result栏,仔细阅读一下里面的内容,会告诉你为你的数据选择了哪个最优核苷酸替代模型等信息,还有用到哪些相关程序,在论文的参考文献部分要怎么引用。
点击左下角的DOWNOAD SELECTED JOBS下载压缩文件,解压后,用查看系统发育树的软件打开其中的.treefile文件进行查看。
我发现我的树上有两个支持率的值,阅读Full Result栏的解释后发现是SH-aLRT support (%) / ultrafast bootstrap support (%)
IQ-TREE只能进行最大似然法的建树,但是CIPRES能进行多个建树方法的建树。
这个需要注册,然后登陆,点击首页蓝色背景的Use the CRIPRESScience Gateway

先Create New Folder,在新建文件夹下会有Data和Tasks两个文件

点击左边的Data文件,点击Upload Data上传自己的序列文件(已经完成多重序列比对):

点击左边Tasks文件,点击Create New Task,在Description一栏中,可以写入对这个任务的描述:

在Input一栏中,选择要分析的序列文件。在Tool一栏中,选择建树的方法,弹出以下页面:

一般选择RAxML-HPC2 on XSEDE。
然后是Ser Parameters,设置参数:

需要改动的参数有
Maximum hours to run:128
Set a name for output files:给输出结果一个名字
Enter the number of patterns in your dataset:这里给出的是进行多序列比对之后的总体序列长度
Please select the data type:选择数据类型,这里默认是核苷酸序列
Estimate proportion of invariable sites:默认是No,选择Yes。
最后选择save and run task。
任务完成之后邮箱会收到邮件提醒。
然后打开对应的任务,出现如下界面,点击右面的View Output:

之后会出现一堆output文件,我们要下载的,包含支持率的树文件是:RAxML_bipartitions.result
点击右面的Download,就可以下载了。
不过好像CIPRES最近开始收费了,唉。
https://www.hiv.lanl.gov/content/sequence/PHYML/interface.html
https://www.hiv.lanl.gov/content/index
这个网站主要是收集HIV遗传数据的,但是也有一些在线的分析软件
构建最大似然树的软件包括:
PhyML:https://www.hiv.lanl.gov/content/sequence/PHYML/interface.html
IQ-tree:https://www.hiv.lanl.gov/content/sequence/IQTREE/iqtree.html
但是数据太大的话,可能会受到限制,不能跑。
https://raxml-ng.vital-it.ch/#/

这个是用RaxML构建最大似然树的在线平台,但同样的,如果数据太大的话,可能会受到限制,不能跑。
在Ruby中可以使用哪些替代方法来ping一个ip地址?标准库“ping”库的功能似乎非常有限。我对在这里滚动我自己的代码不感兴趣。有没有好的gem?我应该接受它并忍受它吗?(我在Linux上使用Ruby1.8.6编写代码) 最佳答案 net-ping值得一看。它允许TCPping(如标准rubyping),但也允许UDP、HTTP和ICMPping。ICMPping需要root权限,但其他则不需要。 关于ruby-Pingruby网站?,我们在StackOverflow上找到一个类
在编写Ruby(客户端脚本)时,我看到了三种构建更长字符串的方法,包括行尾,所有这些对我来说“闻起来”有点难看。有没有更干净、更好的方法?变量递增。ifrender_quote?quote="NowthatthereistheTec-9,acrappyspraygunfromSouthMiami."quote+="ThisgunisadvertisedasthemostpopularguninAmericancrime.Doyoubelievethatshit?"quote+="Itactuallysaysthatinthelittlebookthatcomeswithit:themo
电脑0x0000001A蓝屏错误怎么U盘重装系统教学分享。有用户电脑开机之后遇到了系统蓝屏的情况。系统蓝屏问题很多时候都是系统bug,只有通过重装系统来进行解决。那么蓝屏问题如何通过U盘重装新系统来解决呢?来看看以下的详细操作方法教学吧。 准备工作: 1、U盘一个(尽量使用8G以上的U盘)。 2、一台正常联网可使用的电脑。 3、ghost或ISO系统镜像文件(Win10系统下载_Win10专业版_windows10正式版下载-系统之家)。 4、在本页面下载U盘启动盘制作工具:系统之家U盘启动工具。 U盘启动盘制作步骤: 注意:制作期间,U盘会被格式化,因此U盘中的重要文件请注
在应用开发中,有时候我们需要获取系统的设备信息,用于数据上报和行为分析。那在鸿蒙系统中,我们应该怎么去获取设备的系统信息呢,比如说获取手机的系统版本号、手机的制造商、手机型号等数据。1、获取方式这里分为两种情况,一种是设备信息的获取,一种是系统信息的获取。1.1、获取设备信息获取设备信息,鸿蒙的SDK包为我们提供了DeviceInfo类,通过该类的一些静态方法,可以获取设备信息,DeviceInfo类的包路径为:ohos.system.DeviceInfo.具体的方法如下:ModifierandTypeMethodDescriptionstatic StringgetAbiList()Obt
需求:要创建虚拟机,就需要给他提供一个虚拟的磁盘,我们就在/opt目录下创建一个10G大小的raw格式的虚拟磁盘CentOS-7-x86_64.raw命令格式:qemu-imgcreate-f磁盘格式磁盘名称磁盘大小qemu-imgcreate-f磁盘格式-o?1.创建磁盘qemu-imgcreate-fraw/opt/CentOS-7-x86_64.raw10G执行效果#ls/opt/CentOS-7-x86_64.raw2.安装虚拟机使用virt-install命令,基于我们提供的系统镜像和虚拟磁盘来创建一个虚拟机,另外在创建虚拟机之前,提前打开vnc客户端,在创建虚拟机的时候,通过vnc
我正在尝试在配备ARMv7处理器的SynologyDS215j上安装ruby2.2.4或2.3.0。我用了optware-ng安装gcc、make、openssl、openssl-dev和zlib。我根据README中的说明安装了rbenv(版本1.0.0-19-g29b4da7)和ruby-build插件。.这些是随optware-ng安装的软件包及其版本binutils-2.25.1-1gcc-5.3.0-6gconv-modules-2.21-3glibc-opt-2.21-4libc-dev-2.21-1libgmp-6.0.0a-1libmpc-1.0.2-1libm
因为我现在正在做一些时间测量,我想知道是否可以在不使用Benchmark类或命令行实用程序time的情况下测量用户时间或系统时间。使用Time类只显示挂钟时间,而不显示系统和用户时间,但是我正在寻找具有相同灵active的解决方案,例如time=TimeUtility.now#somecodeuser,system,real=TimeUtility.now-time原因是我有点不喜欢Benchmark,因为它不能只返回数字(编辑:我错了-它可以。请参阅下面的答案。)。当然,我可以解析输出,但感觉不对。*NIX系统的time实用程序也应该可以解决我的问题,但我想知道是否已经在Ruby中实
我需要用任何语言编写一个算法,根据3个因素对数组进行排序。我以度假村为例(如Hipmunk)。假设我想去度假。我想要最便宜的地方、最好的评论和最多的景点。但是,显然我找不到在所有3个中都排名第一的方法。Example(assumingthereare20importantattractions):ResortA:$150/night...98/100infavorablereviews...18of20attractionsResortB:$99/night...85/100infavorablereviews...12of20attractionsResortC:$120/night
在Ruby中,以毫秒为单位获取自纪元(1970)以来的当前系统时间的正确方法是什么?我试过了Time.now.to_i,好像不是我想要的结果。我需要结果显示毫秒并且使用long类型,而不是float或double。 最佳答案 (Time.now.to_f*1000).to_iTime.now.to_f显示包含十进制数字的时间。要获得毫秒数,只需将时间乘以1000。 关于ruby-以毫秒为单位获取当前系统时间,我们在StackOverflow上找到一个类似的问题:
关闭。这个问题需要更多focused.它目前不接受答案。想改进这个问题吗?更新问题,使其只关注一个问题editingthispost.关闭8年前。Improvethisquestion我们有以下(以及更多)系统,我们将数据从一个应用推送/拉取到另一个:托管CRM(InsideSales.com)Asterisk电话系统(内部)横幅广告系统(openx,我们托管)潜在客户生成系统(自行开发)电子商务商店(spree,我们托管)工作板(本土)一些工作网站抓取+入站工作提要电子邮件传送系统(如Mailchimp,自主开发)事件管理系统(如eventbrite,自主开发)仪表板系统(大量图表和