Q-Q plot图的含义是什么

发布时间:2022-01-15 18:12:43 作者:柒染
来源:亿速云 阅读:178

小编今天带大家了解Q-Q plot图的含义是什么,文中知识点介绍的非常详细。觉得有帮助的朋友可以跟着小编一起浏览文章的内容,希望能够帮助更多想解决这个问题的朋友找到问题的答案,下面跟着小编一起深入学习“Q-Q plot图的含义是什么”的知识吧。

Q-Q plot是关联分析结果可视化的一种经典方案,这里的Q代表quantile, 分位数的意思,关联分析的Q-Q plot示意如下

Q-Q plot图的含义是什么

x轴代表期望p值,y轴代表实际p值。在解释这张图的含义之前,有必要先来了解下什么是分位数。

分位数,也称之为分位点,最常见的有中位数,四分位数等。以中位数为例,将数据集从小到大排列后,50%区域对应的点就是中位数。同理,四分位数分别对应25%, 50%, 75%, 依次称之为第一四分位数,第二四分位数,第三四分位数,其中第二四分位数就是中位数了,3个四分位数将数据划分为了4个区间,这也是其名字中四的由来。

四分位数只是分位数中的3个点,起始分位数可以是任何一个比例,比如10%的分位数,20%的分位数。给定一个数据集,在R中计算分位数的代码如下

Q-Q plot图的含义是什么

分位数可以很好的展示数据从最小值到最大值的跨度变化,在分位数点取值足够多的情况下, 可以用来代表整体数据。

Q-Q plot就是基于这样的原理,分别计算两个数据的分位数,然后绘制散点图。可以想象,如果两个总体完全一致,其Q-Q plot是一条y=x的直线,代码如下

Q-Q plot图的含义是什么

结果如下图所示

Q-Q plot图的含义是什么

在此基础上进一步推论,如果两个数据符合同一分布,则其分位数应该符合线性关系,验证如下

Q-Q plot图的含义是什么

输出结果如下

Q-Q plot图的含义是什么

上述代码中x和y由两个范围不同的均匀分布抽样产生,可以看到其Q-Q plot近似一条直线。通过这个例子可以看到,Q-Q plot的核心作用就是比较两个数据的分布是否一致。上述代码只是用于展示Q-Q plot的原理,分位数的求解过程本质是将数据从小到大排序之后计算百分比,在实际绘制Q-Q plot图示,直接将两个数据从小到大排序之后绘图就可以了,代码如下

Q-Q plot图的含义是什么

输出结果如下

Q-Q plot图的含义是什么

可以看到,趋势和用分位数来画是一样的。Q-Q plot适用范围广泛,可以比较任意两个数据集的分布。常见的用法有两种,第一种直接拿两个真实的数据集进行比较,查看分布是否一致,  第二种那一个真实数据,和一个基于理论分布计算的数据来比较,查看是否符合理论分布。

关联分析的Q-Q plot就是第二种用法,, 理论分布是均匀分布,拿实际关联分析的p值来和理论分布进行比较。为什么理论分布是均匀分布呢?

要判断一个数据的理论分布,最直观的方式是绘制实际数据的密度分布图,GWAS中p值分布的密度直方图如下

Q-Q plot图的含义是什么

Q-Q plot图的含义是什么

在上图中,每个bin内的密度基本是一样的,是典型的均匀分布的特征。将实际数据的密度分布图和各种理论分布的密度分布图进行比较,可以快速确定候选的理论分布。从峰型的比较来看,也只有均匀分分布和该数据的分布接近。用qqman中的数据集展示Q-Q plot绘制的方法,代码如下

Q-Q plot图的含义是什么

输出结果如下

Q-Q plot图的含义是什么

Q-Q pot绘制非常简单,三句代码就可以搞定。绘图时对p值做了-log10转换,所以图中越靠右的点代表的是p值显著的snp位点。可以看到左侧的点基本服从均匀分布,而显著的snp位点位于标准直线的上方,说明p值显著的点导致了实际的分位数小于了理论的分位数,示意图如下

Q-Q plot图的含义是什么

说明在较小的p值范围内,实际数据的分布更加紧密,当检测到关联信号时,关联位点的p值肯定较为紧密和偏小,数据产生上图的分布是正常现象。如果完全服从均匀分布,即基本和直线重合,说明p值的产生就是一个随机性的过程,gwas分析的可信度较低。

感谢大家的阅读,以上就是“Q-Q plot图的含义是什么”的全部内容了,学会的朋友赶紧操作起来吧。相信亿速云小编一定会给大家带来更优质的文章。谢谢大家对亿速云网站的支持!

推荐阅读:
  1. Q&A about Linux
  2. Q语言 我的命令库《通用找图》之后台找字函数源码

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

上一篇:MES修炼的五层境界是什么

下一篇:springboot整合quartz定时任务框架的方法是什么

相关阅读

您好,登录后才能下订单哦!

密码登录
登录注册
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》