Testing for diversifying selection for two clades with a background clade
在利用branch-site检测趋同进化的时候 .可以将各个趋同进化枝分别进行检测,分析的时候不去除某一趋同枝系
.在分析的时候,需要去除其他趋同枝系的影响
I have sequences of a bacteria gene each from one strain. Among them, there are five strains belong to the same sequence type (ie highly similar) and the other four strains belong to another sequence type. The rest of the strains are from different sequence types such that they are pretty divergent among themselves and also the aforementioned strains I classify the first five as clade , the next four as clade and the rest as clade which serves as a background clade by labeling the tree as instructed in the manual. So now I have a Clade Model C with three clades. I am mostly interested in whether there is a diversifying selection going on between clade and clade . What is the right way to do this? I came across this paper saying I can use a null model that combines clade and clade into one clade such that we can simulate the situation of w3 == w4. Then by using an alternative model that has three clades, I can run LRT with df= to test for diversifying selection
http://www.biomedcentral.com/1471-2148/12/206 Does this make sense? Is this the right way for me to test for diversifying selection between two clades? Assume you have a simple tree with three labeled clades, as follows: ((A1,A2)$ , (B1,B2)$ , (C1,C2)$); Clade A is labeled with $, clade B with $, and clade C with $. This last labelling doesn't need to be specified—codeml will automatically label the unlabelled branches/clades with ''. First, to test for significant variation among clades you can compare the fit of Clade model C (using the labeled input tree shown above) versus M2a_rel. M2a_rel assumes that $, $, $, etc, are all evolving under the same selection pressures. This test should have degrees of freedom. Second, to test for significant variation between clades A and B while simultaneously allowing clade C to be different, you can compare the fit of CmC when run using the tree provided above versus CmC when run using a simpler tree. In this case, the simpler tree would assign clades A and B to the same group, like so: ((A1,A2)$ , (B1,B2)$ , (C1,C2)$); . This test should have degree of freedom.
So with a three clades model C against M2a_rel, you can test if each of the three clades undergoing diversifying selection themselves to adapt to the environment. Not exactly. You can't say that each of the clades is divergent. If you have three clades, you have three freely-estimated 'site class 2' parameters under CmC: w2, w3, and w4. M2a_rel assumes that w2 = w3 = w4. Comparing these two models tests whether the M2a_rel assumption holds, and if the test is significant then you can conclude that the assumption doesn't hold. However, there are several ways to violate this M2a_rel assumption. For example:
w2 and w3 and w4 might all be different from each other
w2 might equal w3, with w4 being different
w2 might equal w4, with w3 being different
etc...
By analogy, think of an ANOVA where you're comparing the mean value in three groups. If you have a significant test result, you can then conclude that there is significant variation among the mean values for the groups. However, you can't know for sure which groups are significantly different unless you restructure your test by designing a more appropriate null model or unless you conduct pairwise comparisons between groups. If we test the three clades model C against a model with two of the clades merged, then we are testing whether there are diversifying selection between the two clades that are merged.
Testing for diversifying selection for two clades with a background clade的更多相关文章
- semantic框架 - selection dropdown自定义样式
semantic和Bootstrap是我最常用的前端框架. 这次需求,是要修改selection dropdown. 直接贴代码,不解释. .ui.form input:not([type]):foc ...
- js编写规范
JavaScript编码规范 Bug----33条 1. 不要使用’==’和’!=’,使用’===’和’!==’替代 等级:Major 原因:==和!=在判断值相等前会判断类型是否相等.这容易因为类型 ...
- MSDS596 Homework 9
MSDS596 Homework 9 (Due in class on Nov 14) Fall 2017Notes. The lowest grade among all twelve homewo ...
- 树莓派配置文档 config.txt 说明(转)
原文连接:http://elinux.org/RPi_config.txt 由于树莓派并没有传统意义上的BIOS, 所以现在各种系统配置参数通常被存在"config.txt"这个文 ...
- 树莓派进阶之路 (012) - 树莓派配置文档 config.txt 说明
原文连接:http://elinux.org/RPi_config.txt 由于树莓派并没有传统意义上的BIOS, 所以现在各种系统配置参数通常被存在”config.txt”这个文本文件中. 树莓派的 ...
- Raspberry Config.txt 介绍
原文连接:http://elinux.org/RPi_config.txt Config.txt 由于树莓派并没有传统意义上的BIOS, 所以现在各种系统配置参数通常被存在"config.t ...
- wingide 显示中文 及 配色方案
http://lihuipeng.blog.51cto.com/3064864/923231 网上收集的方法: 显示中文: 任意文本编辑器打开:x:\Wing IDE\bin\gtk-bin\etc\ ...
- PHP- 如何在终端输出带颜色的字体?
转自: http://www.neatstudio.com/show-2568-1.shtml 终端显示颜色,在以前的想法当中,都是因为有了.profile的配色方案.而我一般也都是 采用默认的(sn ...
- 你应该知道的25个非常有用的CSS技巧
在我们的前端CSS编码当中,经常要设置特殊的字体效果,边框圆角等等,还要考虑兼容性的问题, CSS网页布局,说难,其实很简单. 说它容易,往往有很多问题困扰着新手,在中介绍了非常多的技巧,这些小技巧与 ...
随机推荐
- Linux基础入门 - 3
第四节 Linux 目录结构及文件基本操作 4-1.Linux目录结构 Linux 的目录与 Windows 的目录的实现机制是完全不同的.一种不同是体现在目录与存储介质(磁盘,内存,DVD 等)的关 ...
- Eclipse PHPEclipse 配置
最近偶来兴致趁着有些时间,看了看php的书. 说到php就不得不提php的开发环境了,一般的都是采用apache做服务器.mysql做数据库,再加上php组合成一个完备的运行环境,但是好像没有写代码的 ...
- list 去重并保持原来排序
public <T> List<T> removeDuplicateKeepOrder(List<T> list){ /* Set set = new HashSe ...
- Static和static block(静态块)的用法
一.用法:是一个修饰符,用于修饰成员(成员变量 成员函数)被动态所共享 当成员被静态修饰后,就多了一种调用方式,除了可以被对象调用外,还可以直接被类名调用. 类名.静态成员 二.static特点: 1 ...
- django admin后台插件:django-suit入门
去年9月底开始用django来做公司内部项目,开始对django有了一些了解,感觉django真的蛮强大的(也有很多人推荐flask,将来有空的话我会试试).今天的话只是介绍一个小东西,django管 ...
- Spring framework 反序列化的漏洞
理解这个漏洞需要先看freebuff上的jdni的小例子. jndi注入在jdk8u121绕过参考这俩篇文章: https://bl4ck.in/tricks/2019/01/04/JNDI-Inje ...
- burpsuit常用功能
1.生成GET数据包:复制url -> 打开burp -> repeater -> 右键paste url as request 2.生成POST数据包:生成一个GET数据包 -&g ...
- 10-RabbitMQ-整合SpringBoot
RabbitMQ整个SpringBoot SpringBoot因其配置简单.快速开发,已经成为热门的开发之一 消息中间件的工作过程可以用生产者消费者模型来表示.即,生产者不断的向消息队列发送信息 而消 ...
- javascript第一课
JavaScript 一個完整的JavaScript實現由以下三個不同的部分組成 ECMAScript(ECMA歐洲計算機製造商協會):規定了這門語言的語法,類型,語句,關鍵字,保留字,操作符,對象等 ...
- MAX7456和DP7456的对比
DP7456 是一款集成了 EEPROM 的单通道.单色随屏显示发生器,集成了视频驱动器.同步分离器.视频分离开关以及 EEPROM,提高了系统的集成度,有效降低了系统成本. DP7456 采用符合 ...