代码里的b_out表示结局的beta值,b_exp表示暴露的beta值,se_out就是结局的标准误,se_out^2就代表结局beta值的方差,而模型中的-1表示的就是去除截距项。R语言里lm()函数表示拟合线性模型(linear model),summary()函数是用来汇总回归模型拟合的结果。关于R语言的相关信息,请参考往期推送R语言入门系列和R语言进阶...
clump = TRUE: 这个参数告诉函数是否进行 SNP clumping。SNP clumping 是一个常用的步骤,用于减少 SNP...
IVW是最早使用的方法,也是最常用的,它需要SNP完全符合MR研究三原则才能得到正确的因果估计;MR-Egger多加了截距项,其主要目的是判断水平多效性的有无;Weighted Median是利用大部分SNP(majority of genetic variants)来判断因果关系的有无。 接下来是重点部分: 首先推荐使用“MR-PRESSO”包(后面我会简单介绍它的使用...