TwoSampleMR 两样本孟德尔随机化 在使用twosampleMR这个包时,暴露文件是提取的本地文件,如何能筛选出P<5*10-8的SNP?
读取本地暴露文件:
exp_dat <- read_exposure_data(
filename = 'Firmicutes.txt',
sep= "\t",
snp_col = "rsID",
beta_col = "beta",
se_col = "SE",
effect_allele_col ="eff.allele",
other_allele_col = "ref.allele",
#eaf_col = "EAF",
pval_col = "P.weightedSumZ"
)
exp_dat$exposure <- "Firmicutes"
read_exposure_data()只能提取SNP,P值没办法筛选,和extract_instruments()函数不太像。
想问有可以筛选出本地暴露文件P值的函数吗?