MR——自动绘图

最新推荐文章于 2024-03-04 11:28:41 发布

帅锅不识字

最新推荐文章于 2024-03-04 11:28:41 发布

阅读量278

点赞数 1

分类专栏： MR 文章标签： mr

本文链接：https://blog.csdn.net/shuai19970911/article/details/131871391

版权

MR 专栏收录该内容

6 篇文章 3 订阅

订阅专栏

library(TwoSampleMR)
data = read.table('C:/Users/DELL/Desktop/lummr_result.csv',header = T,sep = ',')

for (m in c(147:152)) {
  i = data$id.exposure[m]
  exposure = extract_instruments(outcomes = i)
  outcome <- extract_outcome_data(snps = exposure$SNP, outcomes = 'ukb-b-2349')
  dat <- harmonise_data(exposure_dat = exposure, outcome_dat = outcome)
  filename <- paste0("outfig_", i, '.pdf')
  pdf(file = filename)
  a <- mr_scatter_plot(mr_results = mr(dat, method_list = c("mr_ivw", "mr_egger_regression", "mr_weighted_median")), dat)
  b <- mr_funnel_plot(singlesnp_results = mr_singlesnp(dat))
  c <- mr_leaveoneout_plot(leaveoneout_results = mr_leaveoneout(dat))
  print(a)
  print(b)
  print(c)
  dev.off()
  print(m)
}