① 先修:R统计基础
② 主线:MR 专题模块
③ 常用流程骨架
library(TwoSampleMR) exp <- read_exposure_data("exposure.txt", sep = "\t", snp_col="SNP", beta_col="beta", se_col="se", effect_allele_col="ea", other_allele_col="oa", eaf_col="eaf", pval_col="p") out <- read_outcome_data("outcome.txt", snps = exp$SNP, sep = "\t", snp_col="SNP", beta_col="beta", se_col="se", effect_allele_col="ea", other_allele_col="oa", eaf_col="eaf", pval_col="p") dat <- harmonise_data(exp, out)
res <- mr(dat)
pleio <- mr_pleiotropy_test(dat)
het <- mr_heterogeneity(dat) 建议在完整结果中补充:leave-one-out、MR-PRESSO、森林图、漏斗图、STROBE-MR 报告清单。