如何在gtsummary中按分层分类变量展示对应连续变量统计值

需求说明 #

基于trial数据集生成表格：

• 按grade分层展示
• 按trt分组对比
• 每行对应stage的一个水平（T1-T4），展示该水平下marker变量的中位数(IQR)
• 移除原表格中第一行的marker整体统计行

修改后代码 #

library(gtsummary)

trial %>%
  select(trt, marker, stage, grade) %>%
  tbl_strata(
    strata = grade,
    ~ .x %>%
      tbl_continuous(
        variable = marker,
        by = trt,
        group = stage,
        statistic = all_continuous() ~ "{median} ({p25}, {p75})",
        digits = all_continuous() ~ 2,
        label = marker ~ "Marker Level (ng/mL)"
      ) %>%
      modify_header(all_stat_cols() ~ "**{level}**")
  )

关键调整说明 #

• 改用tbl_continuous函数：该函数专为展示分类变量各水平下的连续变量统计设计，完美匹配需求场景
• group = stage：指定以stage的不同水平作为行维度，每行对应一个T1-T4分组
• statistic参数：明确指定输出格式为中位数(IQR)，替代默认的观测数占比
• 自动移除了原代码中的marker整体统计行，无需额外处理

内容的提问来源于stack exchange，提问作者Juicebox