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网页布局,说难,其实很简单. 说它容易,往往有很多问题困扰着新手,在中介绍了非常多的技巧,这些小技巧与 ...
随机推荐
- Python实例---简单的选课系统
要求 思路: 构造方法传递过去学校名称,同时利用UUID创建一个随机字符串,用这个字符串来作为要写入的文件名 利用类的__str__方法来实现类对象接口返回学校名称的操作 利用pickle的dumps ...
- 关于Oracle伪列rownum
rownum列跟rowid列不一样,虽然同为表的伪列.但是rowid列是列的一个固定属性,而rownum列是结果集的一个排序.所以像如下查询是不会返回结果的:select rowid,rownum,e ...
- MVC与EF结合:Contoso大学
中文教程 1.通过 MVC 5 使用 Entity Framework 6 Code First 入门 https://docs.microsoft.com/zh-cn/aspnet/mvc/over ...
- 异常处理与MiniDump详解(1) C++异常(转)
异常处理与MiniDump详解(1) C++异常 write by 九天雁翎(JTianLing) -- blog.csdn.net/vagrxie 讨论新闻组及文件 一. 综述 我很少敢为自己写 ...
- video文件格式说明(笔记)
video标签兼容IE8可使用html5media.js,具体demo可以下载文件中的压缩包 移动端兼容参考: http://www.xyhtml5.com/3252.html
- 关于使用 CALayer 中 mask 的一些技巧
CALayer 拥有 mask 属性,Apple 的官方解释如下: An optional layer whose alpha channel is used to mask the layer’s ...
- webview综述
nWebView 是webkit最核心的一个view,WebView管理WebFrameView和WebFrame之间的交互,一个WebView对象绑定一个window,并且要求MainFrame加载 ...
- 3171. [TJOI2013]循环格【费用流】
Description 一个循环格就是一个矩阵,其中所有元素为箭头,指向相邻四个格子.每个元素有一个坐标(行,列),其中左上角元素坐标为(0,0).给定一个起始位置(r,c) ,你可以沿着箭头防线在格 ...
- TensorFlow函数(四)tf.trainable_variable() 和 tf.all_variable()
tf.trainable_variable() 此函数返回的是需要训练的变量列表 tf.all_variable() 此函数返回的是所有变量列表 v = tf.Variable(tf.constant ...
- 自定义ViewPagerIndicator-视图指示器
ViewPagerIndicator.java public class ViewPagerIndicator extends LinearLayout { private Paint mPaint; ...