首页 / R语言实战教程 / 孟德尔随机化 R 路线图

孟德尔随机化(MR)R 路线图

从暴露/结局数据整理、工具变量筛选,到 TwoSampleMR 主分析、敏感性分析和结果报告,形成完整落地链路。

① 先修: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 报告清单。