该问题已被锁定!
3
关注
2916
浏览

Affy包中的MAplot函数的作用是什么?画出的MAplot图如何解读?为何要进行MAplot检查?

为什么被折叠? 0 个回复被折叠
孟浩巍 超级管理员 用户来自于: 北京市
2018-08-24 23:09
MAplot是一个非常好的质控图,无论是在芯片数据分析中还是在RNA-Seq数据分析中都是经常使用的。   假设有ctrl和treat两组数据,gene K的表达量分别是 ctrl[K] 和 treat[K],   MAplot的X轴就是每一个gene的几何平均数的log2值, M[K] = 0.5 * log2(ctrl[K] * treat[K]); 这个几何平均数的log2值就可以衡量这个gene的大概表达量,是高表达gene还是低表达gene.   MAplot的Y轴就是每个gene的log2 FoldChange, A[K] = log2(treat[K] / ctrl[K]), 这个量表达的是变化关系,0表示没变化,大于0表示treat中升高,小于0表示降低。   那么MAplot的意义是什么呢?这个主要和我们的基本假设有关系,我们一般在做芯片或者是RNA-Seq分析的时候,都默认绝大多数的gene不发生变化,高表达的gene不发生变化。因此,在MAplot中,我们应该看到:   1. 绝大多数的点,Y轴都在0附近; 2. X轴很靠右的点,几乎不怎么发生变化; 3. 发生显著性变化的点,比较均匀地分布在X,Y轴的两侧。   如何不满足上面的特点,就说明要对数据进行校正,使其满足我们的基本假设,才能够进行下一步的统计检验。

关于作者

问题动态

发布时间
2018-08-24 15:31
更新时间
2018-08-24 23:09
关注人数
3 人关注

推荐内容

安捷伦(Agilent)的数据处理软件Feature Extraction Software的默认的数据标准化的算法是什么?最后得到的值是否做了log2处理?是否是RMA或者MASS的一种?
请问芯片数据分析的一般流程及涉及的常用算法?
RNA-seq logFC与pvalue
转录组数据样本聚类结果不理想
表达谱芯片的类型及每种类型能够测定的基因数量如何查询?有没有网站进行过总结?
WGCNA中的TOM热图绘制
RNA-seq count值和fpkm值不对应
小麦660K芯片里,参考碱基那一列的A/C,哪个才是参考基因组的碱基?
关于cufflinks中cuffdiff的使用问题。
Rstudio怎么打开project
All Rights Reserved Powered BY WeCenter V4.1.0 © 2026