R语言实现多变量孟德尔随机化分析（1）

最新推荐文章于 2024-11-07 18:36:05 发布

weixin_49320263

最新推荐文章于 2024-11-07 18:36:05 发布

阅读量7.1k

点赞数 1

分类专栏： GWAS 文章标签： r语言

本文链接：https://blog.youkuaiyun.com/weixin_49320263/article/details/134490073

版权

GWAS 专栏收录该内容

10 篇文章

订阅专栏

本文介绍了多变量孟德尔随机化（MVMR）方法，用于调整潜在混杂因素对遗传关联研究中暴露和结局关系的影响。通过实例展示如何使用`TwoSampleMR`和`MVMR`包来处理多个暴露变量和结局数据，以控制混杂偏倚并得出独立效应的结论。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

多变量孟德尔随机化分析调整了潜在混杂因素的影响。

1、调整哪些因素？参考以往文献。可以分别调整，也可以一起调整。

2、解决了什么问题？某个暴露相关的SNP，往往与某个或者某几个混杂因素相关。可以控制混杂偏倚。

3、如何解释结果？若该暴露的P值小于0.05，则可以说明该暴露独立于其他暴露对结局产生影响。否则是通过其他因素对结局产生影响。

#多变量孟德尔随机化（MVMR）
library(TwoSampleMR)
#提取多个暴露变量工具
#body mass index:ieu-b-40；
#hypertension:ebi-a-GCST90038604
#creatinine：ebi-a-GCST90025946
exposure_dat_mv<-mv_extract_exposures(c("ieu-b-40",
                                        "ebi-a-GCST90038604",
                                        "ebi-a-GCST90025946")) #Serum creatinine levels、Smoking initiation

#提取结局信息
outcome_dat_mv<-extract_outcome_data(exposure_dat_mv$SNP,"ebi-a-GCST90013862") #colorectal cancer

#整合数据
mvdat<-mv_harmonise_data(exposure_dat_mv,
                         outcome_dat_mv,
                         harmonise_strictness = 2)

#进行MVMR的分析
res <- mv_multiple(mvdat)

#提取结果
result<-res$result
#install package
# remotes::install_github("WSpiller/RMVMR",
#                         build_opts=c("--no-resave-data", "--no-manual"),
#                         build_vignettes = TRUE)
library(MVMR)
help(package="MVMR")
wer <- format_mvmr(BXGs = mvdat[["exposure_beta"]],
                      BYG = mvdat[["outcome_beta"]],
                      seBXGs = mvdat[["exposure_se"]],
                      seBYG = mvdat[["outcome_se"]],
                      RSID = rownames(mvdat[["exposure_beta"]]))
#IVW多变量孟德尔随机化结果
ivw_mvmr(wer)
#计算F值
Fz<- strength_mvmr(r_input = wer, gencov = 0)
#异质性检验
pres <- pleiotropy_mvmr(r_input = wer, gencov = 0)