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

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

查看全部 1 个回答

孟浩巍 超级管理员 用户来自于: 北京市
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 人关注

推荐内容

关于ggplot2画条图的问题
用R画RNA_seq中变化显著的heatmap聚类热图
ATAC-seq与RNA-seq的联合分析
RNA-seq reads和表达量在转座子上的分布meta图
RNA-seq差异分析
请问芯片数据分析的一般流程及涉及的常用算法?
求助 limma计算DEG最后结果中的调整p值没有小于0.01的怎么办。。。
RNA-seq count值和fpkm值不对应
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024