mr_heterogeneity(dat) # 从上述结果来看,这些IV之间存在很强的异质性(Q_pval远小于0.05),这时候我们需要剔除某些outcome的P值非常小的SNP,或者直接使用随机效应模型来估计MR效应量,代码如下: mr(dat,method_list=c('mr_ivw_mre')) #使用随机效应模型 # 从随机效应...
进一步采用IVW-MRE方法,表明超重对FA有因果效应(β=0.029,95%CI=0.013-0.046,P=0.0005)。3、超重和寿命 MR分析显示,超重对寿命(第90生存百分位数和第99生存百分位数)有明显的因果关系。超重与预期寿命减少有关(第90生存百分位数,β=-0.220,95%CI=-0.323至-0.118,P<0.0001;第99生存百分位数,β=-0.389,95%...
mr_penalised_weighted_median Penalised weighted median FALSE FALSE mr_ivw Inverse variance weighted TRUE TRUE mr_ivw_radial IVW radial FALSE TRUE mr_ivw_mre Inverse variance weighted (multiplicative random effects) FALSE FALSE mr_ivw_fe Inverse variance weighted (fixed effects) FALSE FALSE mr_si...
# 8 mr_ivw # 9 mr_ivw_radial # 10 mr_ivw_mre # 11 mr_ivw_fe # 12 mr_simple_mode # 13 mr_weighted_mode # 14 mr_weighted_mode_nome # 15 mr_simple_mode_nome # 16 mr_raps # 17 mr_sign # 18 mr_uwr # 省略 如果你要更改方法也是非常简单: mr(dat, method_list = c("mr...
从上述结果来看,这些IV之间存在很强的异质性(Q_pval远小于0.05),这时候我们需要剔除某些outcome的P值非常小的SNP,或者直接使用随机效应模型来估计MR效应量,代码如下: 代码语言:javascript 代码运行次数:0 运行 AI代码解释 mr(mydata,method_list=c('mr_ivw_mre'))#使用随机效应模型 ...
het<-mr_heterogeneity(mydata)hetmr(mydata,method_list=c('mr_ivw_mre'))#使用随机效应模型 pleio<-mr_pleiotropy_test(mydata)pleio single<-mr_leaveoneout(mydata)mr_leaveoneout_plot(single) 第一部分 绘制散点图 散点图的绘制其实非常简单,只需如下一行代码即可: ...
write.table(het,file="het.txt",col.names=T,row.names=F,sep="\t",quote=F)#mr(mydata,method_list=c('mr_ivw_mre')) #使用随机效应模型#多效性检验pleio<-mr_pleiotropy_test(mydata)pleio write.table(pleio,file="pleio.txt",col.names=T,row.names=F,sep="\t",quote=F)#逐个剔除检验si...
[1] 80 16mydata <- harmonise_data( exposure_dat=bmi_exp, outcome_dat=t2d_out, action= 2 )res <- mr(mydata)reshet <- mr_heterogeneity(mydata)hetmr(mydata,method_list=c('mr_ivw_mre')) #使用随机效应模型pleio <- mr_pleiotropy_test(mydata)pleiosingle <- mr_leaveoneout(mydata)...
从上述结果来看,这些IV之间存在很强的异质性(Q_pval远小于0.05),这时候我们需要剔除某些outcome的P值非常小的SNP,或者直接使用随机效应模型来估计MR效应量,代码如下: mr(mydata,method_list=c('mr_ivw_mre')) #使用随机效应模型 从随机效应模型的结果不难看出,BMI和二型糖尿病确实存在因果关系(pval < 0.05)...
BI7-EM18WD-AP6X/S90 NSD 编码器 MRE-32SP101LKR10-G HERZOG传感器8-6841-335700-1 SIEMENS 6SL3766-1FA00-0AA0 Parsytec Computer GmbH变送器Y9002-HM-D变送器Parsytec Computer GmbH SIEMENS进口型号3RK1402-3CE00-0AA2 SCHMERSAL限位开关T.250-11Z-UE-966限位开关seli GmbH DUNGS gas leak detection de...