元分析亚组比较

作者

[编辑] 郑虎;

[审核] .

注记

Hiplot 网站

本页面为 Hiplot Meta-Subgroup Analysis 插件的源码版本教程，您也可以使用 Hiplot 网站实现无代码绘图，更多信息请查看以下链接:

https://hiplot.cn/basic/metawho?lang=zh_cn

metawho 的目标是提供“使用元分析方法识别治疗获益最大的人群”的 R 实现。

metawho 由 R 包 metafor 提供支持，目前不支持包含个体的数据集。如果您更熟悉 Stata 代码，请使用 Stata 包 ipdmetan。

环境配置

  • 系统: Cross-platform (Linux/MacOS/Windows)

  • 编程语言: R

  • 依赖包: meta; cowplot

# 安装包
if (!requireNamespace("metawho", quietly = TRUE)) {
  install.packages("metawho")
}
if (!requireNamespace("cowplot", quietly = TRUE)) {
  install.packages("cowplot")
}

# 加载包
library(metawho)
library(cowplot)

数据准备

# 加载数据
data <- read.delim("files/Hiplot/121-metawho-data.txt", header = T)

# 整理数据格式
data = deft_prepare(data, conf_level = 1 - 0.95)
res = deft_do(data, group_level = unique(data$subgroup))

# 查看数据
head(data)
                 entry         trial subgroup   hr ci.lb ci.ub  ni   conf_q
1      Rizvi 2015-Male    Rizvi 2015     Male 0.30  0.09  1.00  16 1.959964
2    Rizvi 2015-Female    Rizvi 2015   Female 0.11  0.02  0.56  18 1.959964
3      Rizvi 2018-Male    Rizvi 2018     Male 1.25  0.82  1.90 118 1.959964
4    Rizvi 2018-Female    Rizvi 2018   Female 0.63  0.42  0.95 122 1.959964
5   Hellmann 2018-Male Hellmann 2018     Male 0.90  0.41  1.99  37 1.959964
6 Hellmann 2018-Female Hellmann 2018   Female 0.28  0.12  0.67  38 1.959964
          yi       sei
1 -1.2039728 0.6142831
2 -2.2072749 0.8500678
3  0.2231436 0.2143674
4 -0.4620355 0.2082200
5 -0.1053605 0.4030005
6 -1.2729657 0.4387290

可视化

# 元分析亚组比较
p1 <- deft_show(res, element = "all")
p2 <- deft_show(res, element = "subgroup")
p <- plot_grid(p1, p2, nrow = 2)

p
图 1: 元分析亚组比较