赞
踩
打开网址:IEU OpenGWAS project (mrcieu.ac.uk)
在框内输入你想查询的暴露因素,例如Mineral and other dietary supplements(矿物质及其他膳食补充剂)
接着选取自己想要研究的的暴露因素的【GWAS ID】,以第一个ukb-b-7043为例。
然后检索结局指标,【stroke】
选择数据集ebi-a-GCST90038613作为结局指标。
打开RStudio我分享的脚本文件,复制暴露因素的GWAS ID到下图的指定位置。
复制结局指标的GWAS ID到下图的指定位置
加载TwoSampleMR、ieugwasr两个包
加载访问令牌
下载暴露因素的数据
展示第一列前四行的数据
获取所有的列名
获取结局数据
统一效应等位与效应量(Harmonization),协调或整合暴露数据和结局数据,以确保它们具有相同的格式和结构
展示数据内容
默认使用五种方法开始MR分析
查看当前支持的统计方法
指定其他方法,以英文逗号为间隔。
beta 转换 OR 值
异质性检验
也可以像 mr() 函数那样,指定方法进行检验。
水平多效性检验
散点图
要是只想展示某几种方法,在进行 MR 分析 的步骤时,指定方法即可
森林图
和上面一样,想展示指定方法的结果,在函数中指定即可。
留一图
漏斗图
选择导出图像后,选择图像的格式
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。