关于GWAS分析中PVE的计算方法:


我查了一下,大体计算PVE的方法有三种:第一种回归分析或者方差分析的方法,计算R方(GLM模型),第二种是根据effect,se,maf计算PVE,第三种是根据LMM的矩阵构建计算PVE。

汇总如下:

所以准备研究一下。

先看一个常见问题:GWAS分析中,SNP解释百分比(PVE)之和为何大于1?

问题来源:TASSEL的google group:

https://groups.google.com/g/tassel/c/v1aOPCYFyJE/m/HromqNnNIyMJ

问题描述:

Dearall, AgainIhaveaquestion. InassociationmappingresultsinTasselbyMLMmodelgivesMarkerR2 values.IfMarkerR2reportsproportionofphenotypicvariation (between0to1)thatisexplainedbycorrespondingmaker,whysumof markerR2valuesforallofmarkersismorethan1?

为何GLM或者MLM中的Marker R2之和会大于1,如果R2是解释的百分比,那应该是在0~1之间呀。

张志武老师回复:

The association studies in TASSEL (MLM or GLM) are performed on markers one at a time. Therefore, the sum of the R square (R2) of markers could be bigger than 100%. One of the reasons is due to linkage disequilibrium (LD) between markers. For example, if a marker has R2 of 20% and the marker is in complete LD with other five markers, then the five markers will have R2 sum to 120%.

在TASSEL中GLM或者MLM模型中,是单标记扫描,之所以SNP的R2(R square)之和会大于1,因为标记间存在LD,比如一个标记关联的基因能解释20%的变异,这个位点附近有6个标记都存在LD状态,那么这6个标记的解释百分比之和就会是120%。

下一节介绍一下GAPIT中GLM的PVE的计算方法,并用R语言实现。

GWAS分析中SNP解释百分比PVE | 第一篇,SNP解释百分比之和为何大于1?相关推荐

  1. GWAS分析中SNP解释百分比PVE | 第二篇,GLM模型中如何计算PVE?

    上一篇,介绍了一下显著性的SNP,他们的解释表型变异百分比(PVE)之和,为何可能大于1. https://yijiaobani.blog.csdn.net/article/details/12209 ...

  2. GWAS分析中SNP解释百分比PVE | 第三篇,MLM模型中如何计算PVE?

    之前,想研究一下GWAS分析汇中PVE(表型方差解释百分比)的计算方法,写了两篇: GWAS分析中SNP解释百分比PVE | 第一篇,SNP解释百分比之和为何大于1? GWAS分析中SNP解释百分比P ...

  3. GWAS分析中SNP解释百分比PVE | 第四篇,MLM模型中如何手动计算PVE?

    系列部分: GWAS分析中SNP解释百分比PVE | 第一篇,SNP解释百分比之和为何大于1? GWAS分析中SNP解释百分比PVE | 第二篇,GLM模型中如何计算PVE? GWAS分析中SNP解释 ...

  4. 笔记 GWAS 操作流程6-2:手动计算GWAS分析中的GLM和Logistic模型

    1. 名词解释 GWAS 全基因组关联分析 手动计算 使用R语言编程GLM模型和Logistic模型,提取Effect和Pvalue GLM 一般线性模型 Logistic 主要分析广义线性模型,Y变 ...

  5. GWAS分析中协变量的区分(性别?PCA?初生重?)

    1. 电子书领取 前几天发了一篇GWAS电子书分享,异常火爆,阅读量8000+,很多人评价比较基础.这本电子书主要特点是比较基础,GLM模型用软件和R语言进行比较,如何添加数字协变量.因子协变量.PC ...

  6. b2c项目基础架构分析(二)前端框架 以及补漏的第一篇名词解释

    b2c项目基础架构分析(二)前端框架 以及补漏的第一篇名词解释 继续上篇,上篇里忘记了也很重要的前端部分,今天的网站基本上是以一个启示页,然后少量的整页切换,大量的浏览器后台调用web服务局部.动态更 ...

  7. GWAS分析中0-1-2的编码问题

    昨天一个小伙伴给我写了封信,问我了几个GWAS的问题,我就回信说,答案可以写一下公众号了. 这不,就来了. 邓老师,您好,看了好多博文,学到好多,最近看一些GWAS原理相关的内容,有些疑惑.想请教一下 ...

  8. GWAS分析中使用PCA校正群体分层

    欢迎关注"生信修炼手册"! GWAS通过分析case/control组之间的差异来寻找与疾病关联的SNP位点,然而case和control两组之间,可能本身就存在一定的差异,会影响 ...

  9. python绘制基因结构图_使用Python绘制GWAS分析中的曼哈顿图和QQ图

    [前言]其实这篇文章是为了简单介绍一下geneview的用法,它是一个Python高级库,建立在matplotlib的基础之上,专门用于基因组数据的可视化,目的是为了使创建高大上(精致)的基因组数据图 ...

最新文章

  1. oracle--with as
  2. sha1withRSA md5withRSA分析
  3. Console-算法[for]-穷举法:百钱买百鸡
  4. 获取 metadata 过程详解 - 每天5分钟玩转 OpenStack(167)
  5. 解决Unity3D导出apk失败:Failed to re-package resources
  6. matlab 导入元胞,MATLAB导入xls文件以及cell的使用方法
  7. Paper Pal:一个中英文论文及其代码大数据搜索平台
  8. Java 算法 矩阵乘法
  9. TrustToken向Curve上tfTUSD贷款池新投入2400万美元资金
  10. java读二进制bin文件内容_利用Java读取二进制文件示例详细解说
  11. 拿什么拯救你,我的头发丝--PS抠图神器
  12. 非root用户安装或升级CUDA和CUDNN版本
  13. 第三阶段应用层——1.11 数码相册—setting_page设置页面的显存管理、页面规划、输入控制
  14. 第22章 声音与音乐
  15. 医学院交换去计算机学院,2017年北京协和医学院阜外心血管病医院823计算机原理考研强化模拟题...
  16. JTextField:单行文本框组件
  17. google搜索特殊字符
  18. 2017web前端校招面试总结 - 陌陌技术保障组
  19. Xcode-关于切换窗口的快捷键
  20. [C++]typedef typename什么意思?

热门文章

  1. AI in Finance<量化交易人工智能金融投资>(上)
  2. AH8652小封装220伏如何做到固定电压5伏呢?
  3. 虚拟机linux挂载虚拟光驱
  4. 第53届全球生态旅游大使世界总决赛落幕 俄罗斯佳丽夺冠
  5. python多因子量化选股模型_GitHub - zhangjinzhi/Wind_Python: 量化开发 多因子选股模型...
  6. python 报表框架_现在学Python还不晚!阿里巴巴推荐的Python学习教程,分享给大家!...
  7. 查找附近的人-geohash使用方法
  8. 【安全问道】系列:3分钟掌握一个互联网安全知识
  9. 数据结构大作业:纯C双链表实现贪吃蛇
  10. C常见问题之RAND_MAX到底是多少