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突变所产生的遗传作用造成的

关于统计学中q-q图为什么正态分布是一条直线(R语言绘图说明)

q-q图的概念说明
  • rav009
  • rav009
  • 2014年10月13日 15:51
  • 9651

SPSS PP图和QQ图

P-P图和Q-Q图都是用来观察变量是否服从正太分布的。 选择菜单中的分析-描述统计-P-P图,则如下图 1变量选项栏:选择绘制pp图的变量,可多选 2检验分布选项栏:用于检测待检验分布的类型,如be...
  • myl1992
  • myl1992
  • 2015年05月13日 08:58
  • 4216

R语言自定义函数,利用QQ图和shaprio.test检验正态分布

代码 norm.test if(pic==TRUE){#画图形 dev.new() par(mfrow=c(2,1)) qqnorm(data,main="qq图") qqline(data) his...
  • SunCherryDream
  • SunCherryDream
  • 2016年04月26日 16:25
  • 4270

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

  • 2013年08月28日 14:32
  • 22.55MB
  • 下载

强化学习系列之九:Deep Q Network (DQN)

文章目录 [隐藏] 1. 强化学习和深度学习结合 2. Deep Q Network (DQN) 算法 3. 后续发展 3.1 Double DQN 3.2 Prioriti...
  • bbbeoy
  • bbbeoy
  • 2018年01月16日 10:45
  • 7

量化分析师的Python日记【Q Quant兵器谱 -之偏微分方程1】

从今天开始我们将进入一个系列 —— 偏微分方程。作为这一系列的开篇,我们以热传导方差为引子,引出: 如何提一个偏微分方程的初边值问题;利用差分格式将偏微分方程离散化;显示差分格式;显示差分格式的条件...
  • artemis_123
  • artemis_123
  • 2016年12月21日 10:06
  • 835

三菱Q系列PLC的modbusTCP通讯

BCNet-Q是一款新型的以太网通讯处理器,是为满足日益增多的工厂设备信息化需求(设备网络监控和生产管理)而设计,用于三菱Q0x/Q0xU/Q0xUD系列PLC的以太网数据采集,BCNet-Q支持工控...
  • wxbca
  • wxbca
  • 2017年01月11日 13:32
  • 3235

[IMX6Q]fastboot下载u-boot.bin失败提示太大原因

u-boot版本: v2009.08 在调试其他功能时增加了u-boot.bin容量到1M,使用fastboot download时提示: sending 'bootloader' (1191 K...
  • kris_fei
  • kris_fei
  • 2016年01月08日 14:12
  • 2392

imx6q内核定时器实现led闪烁

/************************************************************************* > File Name: led_drv.c ...
  • u011583798
  • u011583798
  • 2017年03月12日 23:04
  • 155

HDU 4500 小Q系列故事――�丝的逆袭

毕业于普通本科的小Q一直自称是资深�丝,不仅学校不知名,甚至他自己在这个普通学校也是默默无闻――直到临近毕业的时候,班里5朵金花中的2位甚至从没和他说过话!   谁又能想到,如此不起眼的小Q在历经重...
  • u011217342
  • u011217342
  • 2013年07月31日 21:53
  • 1265
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:Q-Q plot 在GWAS中的意思
举报原因:
原因补充:

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