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

登录

该问题已被锁定！

3: 关注

2346: 浏览

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 人关注

相关问题

进行转录组数据分析时，进行cuffdiff后的输出文件gene_exp.diff中，一个基因出现了两个不同的表达量数据，应该如何处理？: 1744 浏览 2 关注 1 回答 0 评论

在riboseq分析中，如何对没有起始密码子和终止密码子的转录本进行三碱基准确性分析？: 1929 浏览 2 关注 1 回答 0 评论

二代测序得到测序数据组装完成后如何进行丰度计算: 2228 浏览 1 关注 1 回答 0 评论

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

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

R进行GO时出现No gene can be mapped....: 2485 浏览 3 关注 3 回答 0 评论

如何针对affy表达谱芯片计算结果分析不同样本之间是否表达差异: 1606 浏览 2 关注 1 回答 0 评论

请问一下，我的噬菌体基因组fasta文件还是config打头的，是不是需要进一步拼接成scaffold？还是挑选最大的config进行后续分析？: 1758 浏览 2 关注 1 回答 0 评论

生物信息学需要对哪些方面的数学知识进行深入研究: 2609 浏览 3 关注 2 回答 0 评论

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

推荐内容

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

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

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

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

问芯片数据去除批次效应一般用什么软件如何操作？: 1713 浏览 1 关注 1 回答 0 评论

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

问 R进行GO时出现No gene can be mapped....: 2485 浏览 3 关注 3 回答 0 评论

问求助 limma计算DEG最后结果中的调整p值没有小于0.01的怎么办。。。: 1730 浏览 2 关注 1 回答 0 评论

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

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025