gtsummary中的纵向数据表

Question

library(data.table)
library(gtsummary)    
mm<-structure(list(time = c(1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 
    2, 3, 1, 2, 3, 1, 2, 3), meas1 = c(59, 0, 0, 109, 6.8, 0, 98, 
    11, 2.7, 121, 9.9, 0, 136, 63, 22.9, 130, 0, 0, 86, 23, 0), ID = c("001", 
    "001", "001", "002", "002", "002", "003", "003", "003", "004", 
    "004", "004", "005", "005", "005", "006", "006", "006", "007", 
    "007", "007"), gr = c("A", "A", "A", "B", "B", "B", "A", "A", 
    "A", "C", "C", "C", "C", "C", "C", "B", "B", "B", "B", "B", "B"
    )), row.names = c(NA, -21L), class = c("data.table", 
    "data.frame"))

现在，我想生成一个表格，其中包含每个不同时间点的不同组 (gr) 级别的中位数和 IQR。

g <- function(x) {
  k <- as.list(quantile(x, (1 : 3) / 4, na.rm=TRUE))
  names(k) <- .q(Q1, Q2, Q3)
  k
}
mm2 <- mm[, g(meas1), by=.(gr, time)]
mm2[,iqr:=paste0(Q2,"(",Q1,"-",Q3,")")]
dcast(mm2,time~gr,value.var = "iqr")

在 gtsummary 中是否有一种简单的方法可以做到这一点？ 最好不要夸大 header 列中的 N 以及如果范围可以以中位数为中心的换行符的额外奖励？

Answer 1

使用tbl_strata()您可以获得两个或多个变量的摘要。

library(data.table)
library(gtsummary)  

mm<-structure(list(time = c(1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 
                            2, 3, 1, 2, 3, 1, 2, 3), 
                   meas1 = c(59, 0, 0, 109, 6.8, 0, 98, 
                             11, 2.7, 121, 9.9, 0, 136, 63, 22.9, 130, 0, 0, 86, 23, 0), 
                   ID = c("001", 
                          "001", "001", "002", "002", "002", "003", "003", "003", "004", 
                          "004", "004", "005", "005", "005", "006", "006", "006", "007", 
                          "007", "007"), gr = c("A", "A", "A", "B", "B", "B", "A", "A", 
                                                "A", "C", "C", "C", "C", "C", "C", "B", "B", "B", "B", "B", "B"
                          )), row.names = c(NA, -21L), class = c("data.table", 
                                                                 "data.frame"))

tbl <- 
  mm %>%
  mutate(gr = paste("Group", gr)) %>%
  tbl_strata(
    strata = gr,
    ~ .x %>%
      tbl_summary(
        by = time, 
        include = -ID,
        type = meas1 ~ "continuous2",
        statistic = meas1 ~ c("{median} ({p25}, {p75})", "{min}, {max}"),
        label = list(meas1 = "Measure 1")
      ) %>%
      modify_header(all_stat_cols() ~ "**Time {level}**")
  )

^{使用reprex v2.0.2创建于 2022-09-04}