Q-Q plot 在GWAS中的意思

转载 2012年03月21日 13:32:15

Q-Q plot 即Quantile-Quantile Plot。它在各类研究中经常用到,主要是直观的表示观测值与预测值之间的差异。

在SPSS中很容做,Analysis - Descriptive statistics - Q-Qplot。

Q-Q plot主要是用来估计数量性状观测值与预测值之间的差异。一般我们所取得的数量性状数据都为正态分布数据。在GWAS研究中Q-Q plot的X和Y轴主要是代表各个SNP的-lg P values。预测的线是一条从原点发出的45°角的虚线。实际观测值则是标的实心点。

Q-Q plot主要要点:

预测的虚线为什么是45°出来的呢?因为预测的线实际是通过在QQ图中第一象限作图得出。理论上一个点A在该图上的位置应该是A预测值=A实际值,转化为坐标就是A(x,y)x=y。所以预测的线是一条从原点发出的45°线。

观测值的点的坐标是怎么得出来的。同样设点A的坐标是(x,y)x为预测值,y为实际观测值。查了一下R 中qq plot的算法是这样的

pvals <- read.table("DGI_chr3_pvals.txt", header=T)

observed <- sort(pvals$PVAL)
lobs <- -(log10(observed))

expected <- c(1:length(observed))
lexp <- -(log10(expected / (length(expected)+1)))


具体解释是这样的,先把P值从小到大排序。lobs代表纵坐标,lexp代表横坐标,纵坐标就是观测P值的-log10,而横坐标则根据P值数目而定。比如,当只有3个P值 P1=0.0001 P2=0.001 P3=0.01,那么在这个P值组中,length(observed)=3,对于P1=0.0001 expected=1 lexp=-log10(1/3+1),对于P2=0.001 expected=2 lexp=-log10(2/3+1), P3=0.01 expected=3 lexp=-log10(3/3+1)。。。。。依此类推。

如果出现了偏离的情况说明实际值跟预测值有偏差,在GWAS研究中,那个SNP点出现了较大的偏离,则认为这个SNP位点的观测值的偏离是由这个SNP突变所产生的遗传作用造成的

相关文章推荐

利用Python的科学计算包scipy画Quantile-Quantile图

scipy包的画图功能非常强大,画出来的图形能够存为各种各样的格式,方便用户导入到各种文档中去。本文主要的讲解的是采用什么命令来画Quantile-Quantile图,以及画出来的图形如何保存为图片格...

Q-Q plot

Q-Q plot 即Quantile-Quantile Plot。它在各类研究中经常用到,主要是直观的表示观测值与预测值之间的差异。 在SPSS中很容做,Analysis - Descripti...

GIF(genomic inflation factor)值的计算方法

GIF值通常用在GWA过程中,用来控制该过程质量;鉴别出那些低质量的基因标记;例如在统计学中,GIF值大于1.0表示结果有低质量数据的出现; 那么我们如何计算gwa中的GIF呢;在matlab中,我...
  • likelet
  • likelet
  • 2012年03月21日 14:15
  • 3879

QQPlot/Quantile-Quantile Plot

QQPlot用于直观验证一组数据是否来自某个分布,或者验证某两组数据是否来自同一(族)分布。在教学和软件中常用的是检验数据是否来自于正态分布。 详细信息参考: http://onlinestatb...

扩增子图表解读4曼哈顿图:差异分类级别Taxonomy

作者: 刘永鑫 日期:2017-6-30 阅读时长:10min 背景介绍(Introduction)宏基因组学宏基因组学目前的主要研究方法包括:16S/ITS/18S扩增子、宏基因组、宏转录组和...

全基因组关联分析(GWAS)

全基因组关联分析是一种在人类或动植物全基因组中寻找变异序列的方法,全英文名为Genome-wide association study,缩写名为GWAS。 2005年,Science杂志报道了第一篇G...

quantile-quantile plot (qqplot) of the p-values

The QQ plot shows the expected distribution of association test statistics (X-axis) across the mil...

多重假设检验与Bonferroni校正、FDR校正

总结起来就三句话: (1)当同一个数据集有n次(n>=2)假设检验时,要做多重假设检验校正 (2)对于Bonferroni校正,是将p-value的cutoff除以n做校正,这样差异基因筛选的p-va...

Bootstrap-打造自己的轮播器

在上篇博客中介绍了Bootstrap框架的一些基本情况。也通过一个小例子说明了如何具体在项目中使用这个框架。相信大家对这套框架已经有了一定认识。这篇博客就介绍一下Bootstrap实际应用在项目中时所...

del /f /s /q 是什么意思_DOS命令教程

  • 2013年08月28日 14:32
  • 22.55MB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Q-Q plot 在GWAS中的意思
举报原因:
原因补充:

(最多只允许输入30个字)