孟德尔随机化中的无效工具变量检验

欢迎关注”生信修炼手册”!

两样本的孟德尔随机化研究只需要基于gwas summary数据,就可以研究暴露因素和结局变量之间的因果关系,是最广泛使用的研究手段之一。要保证MR研究结果的可靠性,需要在分析的各个环节进行有效的质控。

孟德尔随机化研究中以遗传变异为工具变量,而合格的工具变量需要符合以下几个条件

  1. 与暴露因素的强关联,称之为关联性假设

  2. 与混杂因素相独立,称之为独立性假设

  3. 只会通过暴露因素对结局变量造成影响,称之为排他性假设

众所周知,基因具有多效性(pleiotropy effect), 当把具有多效性的基因纳入工具变量的集合时,会对MR的分析结果造成影响。当多效性的基因可以通过混杂因素影响结局变量时,违背了独立性假设;当多效性的基因可以通过代谢通路等其他途径来影响结局变量时,违背了排他性假设。无法满足以上3点的遗传变异,称之为无效的工具变量。

为了确保分析结果的准确性,需要对无效工具变量进行检验, 有多种检验方式,其核心思想是异质性检验, 假设每个工具变量的因果效应估计值是相似的,当存在无效工具变量时,其计算得到的因果效应估计值与有效的工具变量存在较大差异,总体上看就会存在显著的异质性。

在进行无效工具变量的显著性检验时,常用的方法为Cochran’s Q检验,这一检验异质性的方法在meta分析中广泛使用。Q检验适用于处理如下数据

subject表示样本,其他列表示不同的条件。Q检验要求条件的个数要大于等于2,而且每个样本的数值为分二类的数值。其空假设为所有条件下数值的分布相同,其统计量称之为Q, 计算公式如下

其中

N表示样本总数,k表示条件的个数,i表示行,j表示列,Y表示i行j列对应的值。当样本量很大时,Q统计量服从自由度为k-1的卡方分布,对应p值的计算公式如下


当pvalue小于0.05 ,认为异质性不显著。除了Q检验外,还有一种在Q统计量的基础上发展而来的I square统计量,计算公式如下

通常情况下,I square的值在0-0.25之间,认为是较低异质性;0.25-0.5之间,中等程度的异质性,0.5到0.75之间,真实存在异质性,0.75以上,异质性很大。

在MR分析中,为了确保分析结果的可靠性,有必要进行无效工具变量的检验。

·end·

—如果喜欢,快分享给你的朋友们吧—

本公众号深耕耘生信领域多年,具有丰富的数据分析经验,致力于提供真正有价值的数据分析服务,擅长个性化分析,欢迎有需要的老师和同学前来咨询。

  更多精彩

  写在最后

转发本文至朋友圈,后台私信截图即可加入生信交流群,和小伙伴一起学习交流。

扫描下方二维码,关注我们,解锁更多精彩内容!

生物信息入门

只差这一个

公众号

点击下方“阅读原文”,我们一起进步!

  • 1
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
输入暴露和结局的数据,然后进行孟德尔随机分析,得出暴露和结局之间的因果关系。下面是使用TwoSampleMR包进行孟德尔随机分析的步骤: ```R # 安装TwoSampleMR包 install.packages("TwoSampleMR") 多变量孟德尔随机是# 加载TwoSampleMR包 library(TwoSampleMR) # 输入暴露和结局的数据 exposure_data <- read一种用于实验设计的方法,可以通过随机分配处理组来减少实验的.table("exposure_data.txt", header=TRUE) outcome_data <- read.table("outcome_data.txt", header=TRUE) # 进行孟混杂因素。下面是一个使用R语言进行多变量孟德尔随机的代码示例:\n\```R\# 安装并加载“randomizR”包\ins.packges(\randomizR\")\ibrary德尔随机分析 mr_results <- mr_egger(exposure_data, outcome_data) # 查看结果 summary(mr(randomizR)\n\# 创建一个包含3个处理组和2个因素的数据框\f <- .fram(\ group_results) ``` = (\A\", \B\", \C\"),\ factor1 = (1, 2, 3),\ factor2 = 需要注意的是,输入的数据必须符合一定的格式要求,具体可以参考TwoSampleMR包的文档。此外, (\X\", \Y\", \Z\")\)\n\# 进行多变量孟德尔随机\randomiz(df, str还有其他一些孟德尔随机的R包,例如MendelianRandomization和TwoSampleMRGUI等, = \factor1\", block = \factor2\", s = 123)\```可以根据具体需求选择使用。 \n\在这个示例,我们使用了R语言的“randomizR”包来进行多变量孟德尔随机。首先,我们创建了一个包含3个处理组和2个因素的数据框。然后,我们使用“randomiz”函数进行随机,其“str”参数指定了分层因素,
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值