孟德尔随机化代码 数据提取-在线_孟德尔clump 在线读取

setwd("D:/Software/R/projects/MR")
library(TwoSampleMR)
exposure <- extract_instruments(outcomes='ieu-b-40',
                               p1 = 5e-08,
                               clump = TRUE,
                               p2 = 5e-08,
                               r2 = 0.001,
                               kb = 10000)

outcome <- extract_outcome_data(snps=exposure$SNP, 
                                outcomes="ieu-b-2", 
                                proxies = TRUE,
                                maf_threshold = 0.01,#最小等位基因频率
                                access_token = NULL)#是否通过谷歌访问，国内选NULL

#之后进行harmonise_data()等等,同本地
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16