R语言中的随机效应模型更改代码如下:fit <- summary(lm(b_out ~ -1 +b_exp, weights = 1/se_out^2))代码里的b_out表示结局的beta值,b_exp表示暴露的beta值,se_out就是结局的标准误,se_out^2就代表结局beta值的方差,而模型中的-1表示的就是去除截距项。R语言里lm()函数表示拟合线性模型(linear model),summary()函数是用来汇总回归模型拟合的结果。关于R语言的相关信息,请参考往期推送R语言入门系列和R语言进阶系列。那么这次回归得出来的beta,se和P值就是MR分析的结果。