双样本孟德尔随机研究随笔记录

随时学随时记录

rm(list = ls())
getwd()
setwd("D:\\科研——实验与课题\\孟德尔随机研究\\R目录")
install.packages("remotes")
library(remotes)
install_github("MRCIEU/TwoSampleMR")
library(TwoSampleMR)
#读取暴露相关的SNP
exposure_dat <- extract_instruments(outcomes = "ieu-a-2",
                                    clump = TRUE,r2=0.0001,kb=1000,
                                    access_token = NULL)
write.csv(exposure_dat,"exposure_dat.csv")
#结局
outcome_dat <- extract_outcome_data(snps=exposure_dat$SNP, 
                                    outcomes = "bbj-a-109")
write.csv(outcome_dat,"outcome_dat.csv")
#合并数据,action = 3去除回文序列
dat <- harmonise_data(exposure_dat, outcome_dat,action = 3)
res <- mr(dat)
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容