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

翻译 2012年03月21日 14:15:14

GIF值通常用在GWA过程中,用来控制该过程质量;鉴别出那些低质量的基因标记;例如在统计学中,GIF值大于1.0表示结果有低质量数据的出现;

那么我们如何计算gwa中的GIF呢;在matlab中,我们使用qqplot()来生成QQplot图,但是这并不意味这我们能够很容易的计算并得到GIF值,我们生成QQplot以后才能从头计算GIF。

第一步就是,我们首先得到正态分布的期望顺序统计量;正态分布的第i次顺序统计量服从参数(i,n-i+1)的贝塔分布,期望值是:

E \left\{p_{(i)} \right\} = i / (n+1) \; \; \; \; (i=1,\ldots,n)

这个是比较容易计算的。下面是个简单的MATLAB程序来显示GWA分析的QQplot并计算GIF:

function gif=compute_gif(pvals)
 
% number of p-vals
n=length(pvals);
 
% expected order statistics
es=(1:n)' ./ (n+1);
 
% x-axis
x = -log10(es);
y = -log10(sort(pvals(:)));
 
% compute GIF
gif =(x'*y)/(x'*x);
 
% QQ-plot
figure;
hold;
grid;
maxh=ceil(log10(n));
xlim([0 maxh]);
ylim([0 maxh]);
plot(x,y,'bx');
plot(x,x,'r-');
xlabel('-log10(Expected Order Statistics)');                                      
ylabel('-log10(Observed Order Statistics)');
title('QQ Plot');
 
% done
return;
函数compute_gif() 只有一个参数,这是p值的list,对p值进行排序并使用该排序生成QQplot计算GIF。

Q-Q plot 在GWAS中的意思

Q-Q plot 即Quantile-Quantile Plot。它在各类研究中经常用到,主要是直观的表示观测值与预测值之间的差异。 在SPSS中很容做,Analysis - Descriptive...
  • likelet
  • likelet
  • 2012年03月21日 13:32
  • 11621

4th Oct 2014:R语言中的factor类型

Factor在回归分析中的应用今天做回归的时候要用到这个函数,被help文档以及不专业的翻译搞昏了factor这个函数到底在干嘛,真是佩服将其翻译成因子。 实际上很简单,factor(C(…))是用...
  • u011022115
  • u011022115
  • 2015年05月28日 21:03
  • 2178

ROS naviagtion analysis: costmap_2d--InflationLayer

这是costmap_2d的最后一个重要的类了。InflationLayer 本身不存储map数据,所以所谓的地图层的概念,仅仅指的是一种对地图的操作,并不是数据层面的地图。 首先来看一个数据结构,知...
  • u013158492
  • u013158492
  • 2016年01月11日 10:21
  • 1182

XP/Win7主机 + Oracle VM VirualBox + MAC OSX Lion

更多精彩内容,请见:http://www.16boke.com 安装过程经历了一些错误,不过结局是好的一切都是值得的。 之所以遇到一些错误是对MAC OS X不了解,找了无数帖子,按照别人一个一个的...
  • one_in_one
  • one_in_one
  • 2013年01月30日 10:28
  • 5058

Score Inflation_usaco3.1.2_dp

DESCRIPTIONThe more points students score in our contests, the happier we here at the USACO are. We ...
  • jpwang8
  • jpwang8
  • 2016年08月28日 23:08
  • 246

OpenCV双边滤波测试

Opencv双边滤波测试
  • bxinquan
  • bxinquan
  • 2017年09月26日 11:40
  • 139

R语言中的factor

对于初学者来说,R语言中的factor有些难以理解。如果直译factor为“因子”,使得其更加难以理解。我倾向于不要翻译,就称其为factor,然后从几个例子中理解: data ...
  • yeweiyangblog
  • yeweiyangblog
  • 2016年10月17日 14:35
  • 15687

R语言的factor是什么含义(1)

首先,c( )是R的一个函数,表示将括号里面的内容连接起来成为一个向量。 例如 先输入 >  y 再输入 > y 得到的结果是 [1] "郭靖 "黄蓉" "老顽童" 有数...
  • weblwh
  • weblwh
  • 2016年01月02日 11:50
  • 373

R语言从基础入门到提高(五)factor 因子

第1程序 What's a factor and why would you use it? 100xp In this chapter you dive into the wo...
  • qq_26816591
  • qq_26816591
  • 2016年08月31日 20:05
  • 2310

R语言入门--R中的factor

如果学过Java或其他高级语言,应该有枚举类型这个概念,个人理解,R中的factor类似于这个枚举类型,但又有不同,毕竟Java中的enum本质上是整数,可以当成整数使用,但R中的factor却不可以...
  • pc_amoon
  • pc_amoon
  • 2017年07月19日 10:12
  • 606
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:GIF(genomic inflation factor)值的计算方法
举报原因:
原因补充:

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