全部案例
生信分析生信模式遗传流行病学真实运行

GWAS 汇总统计的孟德尔随机化分析

用公开的 GWAS 汇总统计,以 BMI 为暴露、2 型糖尿病为结局做两样本孟德尔随机化:筛选工具变量,跑 IVW、MR-Egger 和加权中位数,做异质性与多效性检验,画散点图和森林图。

智能体怎么做

  1. 下载真实数据暴露用 GIANT 2015 的 BMI 汇总统计(欧洲人群、不含 UK Biobank),结局用 FinnGen R12 的 2 型糖尿病(82,878 例 / 403,489 对照),LD 参考用千人基因组 EUR;网络慢就改为多线程分段下载,并按索引只取需要的位点。
  2. 筛选工具变量p<5×10⁻⁸ 的 1,860 个 SNP 按 r²<0.001 去连锁得到 68 个;结局缺失的用代理 SNP 补上;逐个对齐效应等位基因,按频率判断回文 SNP,剔除无法判断的 FTO 位点,最终 67 个,F 统计量最小 29。
  3. 三种方法估计BMI 每升高 1 个标准差,2 型糖尿病风险:IVW OR 1.95(1.43–2.66),加权中位数 2.17(1.93–2.44),MR-Egger 2.32(0.92–5.86),方向一致。
  4. 异质性与敏感性异质性很强(Q=1795,I²=96%),但 Egger 截距 p=0.70,没有定向多效性证据;MR-PRESSO 剔除 11 个离群 SNP、Steiger 方向过滤、留一法后 OR 都在 1.9–2.2 之间;局限里写明两份数据可能有少量芬兰队列重叠。

运行截图

PaperOK科研助手中这次对话的实际界面。

GWAS 汇总统计的孟德尔随机化分析 运行截图:先检查各数据源能否访问:S3 地址连不上,就换其他下载地址继续
先检查各数据源能否访问:S3 地址连不上,就换其他下载地址继续点击查看大图
GWAS 汇总统计的孟德尔随机化分析 运行截图:在右侧直接预览散点图 PDF,左侧是本轮生成的结果表和图
在右侧直接预览散点图 PDF,左侧是本轮生成的结果表和图点击查看大图

运行结果

以下文件由智能体在应用中实际生成,未做修改。

GWAS 汇总统计的孟德尔随机化分析:67 个工具变量的 SNP 效应散点图:IVW OR 1.95,Egger 截距接近 0(p=0.70);最偏离的点是 TCF7L2 位点
67 个工具变量的 SNP 效应散点图:IVW OR 1.95,Egger 截距接近 0(p=0.70);最偏离的点是 TCF7L2 位点 · 下载脚本、结果表与报告
GWAS 汇总统计的孟德尔随机化分析:单个 SNP 的估计与三种方法的合并结果:IVW 1.95、MR-Egger 2.32、加权中位数 2.17
单个 SNP 的估计与三种方法的合并结果:IVW 1.95、MR-Egger 2.32、加权中位数 2.17
GWAS 汇总统计的孟德尔随机化分析:剔除 11 个 MR-PRESSO 离群 SNP 后三种方法高度一致:IVW OR 2.08(1.87–2.30)
剔除 11 个 MR-PRESSO 离群 SNP 后三种方法高度一致:IVW OR 2.08(1.87–2.30)

把这个案例导入到软件里

这次运行的完整对话:打开后能看到每一步做了什么、打开生成的文件,还可以接着提问修改。

已经装好最新版 PaperOK科研助手:点「在软件中打开」,软件会自动下载并打开这个案例;也可以在软件左侧的「案例库」里找到它。点了没反应,就按下面的步骤手动导入:

  1. 安装 PaperOK科研助手需要 1.0.1 及以上版本,可以在下载页获取。
  2. 下载案例文件 paperok-case-gwas-mr.zip(125.5 MB),下载后不用解压。
  3. 导入打开软件,在左侧栏「对话」标题右边点上传图标(导入对话),选择刚下载的 zip。
  4. 查看或继续导入完成后,这个案例会出现在对话列表最上方:展开「工作过程」看每一步,右侧文件区可以打开生成的文件,也可以直接在输入框里接着提要求。

为控制体积,超过 100 MB 的原始数据和运行环境(如 Python 虚拟环境)没有打包进去,需要重新运行时模型会自动下载或安装。对话里出现的本机路径已做匿名处理。