今天分享一下如何复现一篇孟德尔随机化的文章，这篇文章研究的是BMI与骨关节炎的因果关系。首先，我们在MR-Base platform数据库上找到合适的BMI数据集作为暴露，然后再找到合适的骨关节数据集作为结局。接下来就是点击几下鼠标或者复制代码的事情，下面我们用代码来演示：
#安装相关的R包install.packages("devtools")devtools::install_github("MRCIEU/TwoSampleMR")#读取数据与统计分析library(TwoSampleMR)ao <- available_outcomes()exposure_dat <- extract_instruments(c('ieu-a-2'))outcome_dat <- extract_outcome_data(exposure_dat$SNP, c('ukb-b-14486'), proxies =1, rsq =0.8, align_alleles =1, palindromes =1, maf_threshold =0.3)dat <- harmonise_data(exposure_dat, outcome_dat, action =2)mr_results <- mr(dat)
根据结果来撰写论文：
其实，这样的孟德尔随机化文章复现是非常容易的，也是非常快的，基本不用十分钟。复现操作是非常容易，没有什么难度，最难的就是选题，选题决定着你这篇文章的命运。BMI与骨关节炎这个选题没有什么好的选题依据，不是一个好的选题，所以只能发在低分的OA期刊上。
现在的孟德尔随机化就相当于10年前的meta分析、五年前的生信数据挖掘，从大体行情来说，肯定还是处理于好发的时期，再过几年就不敢说了。
