Affy包中的MAplot函数的作用是什么？画出的MAplot图如何解读？为何要进行MAplot检查？ - 做专业的生信问答网站

登录

该问题已被锁定！

3: 关注

1734: 浏览

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

[attach]37[/attach] 如图所示，这个图该如何解读？

好问题 0 评论收藏举报

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轴的两侧。如何不满足上面的特点，就说明要对数据进行校正，使其满足我们的基本假设，才能够进行下一步的统计检验。

赞同 4 0评论

关于作者

: fcdslz 超级管理员
这家伙很懒，还没有设置简介

13: 回答

1: 文章

12: 问题

问题动态

发布时间: 2018-08-24 15:31

更新时间: 2018-08-24 23:09

关注人数: 3 人关注

相关问题

使用lapa进行APA分析: 1454 浏览 2 关注 1 回答 0 评论

怎么对基因启动子含有什么元件进行分析: 1393 浏览 2 关注 1 回答 0 评论

关于对3dDNA产生的hic文件进行纠错的问题: 1488 浏览 2 关注 1 回答 0 评论

如何对新建网站进行暴力测试？: 4304 浏览 8 关注 11 回答 0 评论

根据Barcode序列进行样本拆分？: 2065 浏览 2 关注 1 回答 0 评论

我用admixture进行群体结构的时候，为什么做出来到K=1到11中K=9最小，但是到K=12之后又更小了: 1992 浏览 2 关注 1 回答 0 评论

log2后的数据进行Wilcoxon秩和检验对结果存在什么影响？: 1746 浏览 2 关注 2 回答 0 评论

perl的如何进行读取和循环的: 1472 浏览 3 关注 2 回答 0 评论

在很多分析表达谱芯片的教程中，使用limma 包寻找差异表达基因，都是先找到差异表达探针，再把探针注释为基因id，为何不先进行注释再计算差异表达呢？这两种方法有什么区别呢？: 1341 浏览 2 关注 1 回答 0 评论

关于affy芯片Affymetrix Human Genome U133 Plus 2.0 Array的注释文件，为何GEO、affy官方与Bioconductor的注释数据有差别？: 1837 浏览 2 关注 2 回答 0 评论

推荐内容

问关于cufflinks中cuffdiff的使用问题。: 1248 浏览 2 关注 1 回答 0 评论

问请问芯片数据分析的一般流程及涉及的常用算法？: 1810 浏览 2 关注 3 回答 0 评论

问 ATAC-seq与RNA-seq的联合分析: 1046 浏览 1 关注 0 回答 0 评论

问 GEOquery下载GEO数据软件报错: 2111 浏览 4 关注 3 回答 0 评论

问表达谱芯片的类型及每种类型能够测定的基因数量如何查询？有没有网站进行过总结？: 1122 浏览 2 关注 1 回答 0 评论

问用R画RNA_seq中变化显著的heatmap聚类热图: 1797 浏览 4 关注 3 回答 0 评论

问转录组数据样本聚类结果不理想: 1591 浏览 2 关注 1 回答 0 评论

问求助，安装DESeq2遇到问题: 1108 浏览 2 关注 1 回答 0 评论

问多个处理组的RNAseq分析中关于counts表格合并的问题: 2084 浏览 2 关注 1 回答 0 评论

问使用MAS5.0处理过的芯片数据，能不能直接取对数后用limma包找差异表达基因，为什么不能？: 1457 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025