Início / coding / Perguntas / 79553847
Accepted
BlueKnight2025
Asked: 2025-04-04 03:36:35 +0800 CST

Criando tabela Resumo usando GTSUMMARY

  • 772

Estou tentando obter uma tabela como a abaixo usando GTSUMMARY. Observe que há 2 camadas de recuo necessárias.

insira a descrição da imagem aqui

Eu tentei usar tbl_hierarchicalpara conseguir isso, pois não consigo encontrar uma maneira de tbl_summarycriar várias camadas de indentação. Abaixo está meu código:

# Denominator
adsl.d <- pharmaverseadam::adsl %>%
  select(c(USUBJID, TRT01A))

#Table body input data
adsl1 <- pharmaverseadam::adsl %>%
  select(c(USUBJID, TRT01A, RFSTDTC, EOSSTT, DTHFL)) %>%
  mutate(
    TOTPART = 'TOTAL NUMBER OF PARTICIPANTS',
    RANDFL = 'RANDOMIZED',
    ITTFL = 'INTENT-TO-TREAT',
    SAFFL = if_else(!is.na(RFSTDTC) & is.na(DTHFL), 'SAFETY', NA_character_),
    COMPLFL = case_when(
      EOSSTT == 'COMPLETED' ~ 'Y',
      EOSSTT == 'DISCONTINUED' ~ 'N'),
    DSREAS = case_when(
      grepl('Screen', TRT01A) ~ TRT01A,
      DTHFL == 'Y' ~ 'Death',
      EOSSTT == 'DISCONTINUED' & substr(RFSTDTC, 1, 4) == '2012' ~ 'Adverse event',
      EOSSTT == 'DISCONTINUED' & substr(RFSTDTC, 1, 4) == '2013' ~ 'Withdrawal by subject',
      EOSSTT == 'DISCONTINUED' & substr(RFSTDTC, 1, 4) == '2014' ~ 'Lost to Follow-up')
  ) %>%
  set_variable_labels(DSREAS = "PRIMARY REASON FOR DISCONTINUATION")

#Manipulate data to create ITT & SAF counts
adsl2 <- adsl1 %>%
  select(c(USUBJID, TRT01A, TOTPART, ITTFL, SAFFL)) %>%
  pivot_longer(
    cols = c("ITTFL", "SAFFL"),
    names_to = "EOSSTT",
    values_to = "RANDFL")

adsl <- bind_rows(adsl1, adsl2) %>%
  filter(!grepl('Screen', TRT01A)) %>%
  set_variable_labels(DSREAS = "PRIMARY REASON FOR DISCONTINUATION")

# Generate hierarchical summary
tbl.hrch <- adsl %>% tbl_hierarchical(
  variables = c(TOTPART, RANDFL, EOSSTT),
  id = USUBJID,
  denominator = adsl.d,
  by = TRT01A,
  include = c(TOTPART, RANDFL, EOSSTT),
  statistic = c(TOTPART, RANDFL, EOSSTT) ~ c("{n} ({p})"),
  digits = list(RANDFL = c(1, 2),
                EOSSTT = c(1, 1))
) %>% add_overall(
  last = TRUE,
  col_label = "**Overall**  \nN = {style_number(N)}",
)

tbl.hrch$table_body <- tbl.hrch$table_body %>%
  mutate(across(starts_with('stat_'), ~ if_else(variable == 'TOTPART', NA, .))) %>%
  filter(!(variable == 'EOSSTT' & label %in% c('ITTFL', 'SAFFL')))

O tbl.hrchme fornece o corpo da tabela pretendido. Tenho duas perguntas.

  1. Existe uma maneira de conseguir isso sem manipular os tbl.hrch$table_bodymetadados?
  2. Podemos personalizar o tipo de variáveis ​​em vez de alfabético? Eu obtenho a linha ITT primeiro, o que não é a intenção do shell.

insira a descrição da imagem aqui

Obrigado por qualquer ajuda!

1 respostas

  1. Best Answer

    Eu usaria a tbl_summary()função básica para esta tabela. No exemplo abaixo, tbl_summary()nos dá a estrutura básica da tabela que precisamos. Então chamamos algumas funções para estilizá-la como a especificação, por exemplo, removendo um cabeçalho de linha desnecessário, adicionando recuo extra, etc.

    Boa programação!

    set.seed(123)
library(gtsummary)
library(dplyr)
packageVersion("gtsummary")
theme_gtsummary_compact()

tbl <- 
  pharmaverseadam::adsl |> 
  # remove sreen failures
  filter(TRT01A != "Screen Failure") |> 
  # create fake data for table
  mutate(
    RANDFL = sample(c(TRUE, FALSE), size = n(), replace = TRUE, prob = c(0.95, 0.5)),
    ITTFL = sample(c(TRUE, FALSE), size = n(), replace = TRUE, prob = c(0.95, 0.5)) %>%
      {ifelse(RANDFL == FALSE, NA, .)},
    COMPLFL = 
      sample(c(TRUE, FALSE), size = n(), replace = TRUE, prob = c(0.95, 0.5)) %>%
      {ifelse(RANDFL == FALSE, NA, .)} |> 
      factor(levels = c(TRUE, FALSE), labels = c("COMPLETED", "DISCONTINUED"))
  ) |> 
  # create a basic summary table
  tbl_summary(
    by = TRT01A,
    include = c(RANDFL, COMPLFL, ITTFL, SAFFL),
    value = list(SAFFL = "Y"),
    label = list(RANDFL = "RANDOMIZED",
                 ITTFL = "INTENT-TO-TREAT",
                 SAFFL = "SAFETY"),
    missing = "no"
  ) |> 
  # add column for all subjects
  add_overall(last = TRUE) |> 
  # update header to match spec
  modify_header(
    all_stat_cols() ~ "{level}  \n(N={n})  \n n (%)",
    label = ""
  ) |> 
  # remove header row for COMPLFL
  remove_row_type(variables = "COMPLFL", type = "header") |> 
  # group summaries together under one heading
  add_variable_group_header(
    header = "TOTAL NUMBER OF PARTICIPANTS",
    variables = everything()
  ) |> 
  # indent the COMPLFL rows under RANDFL
  modify_column_indent(
    columns = "label",
    rows = variable == "COMPLFL",
    indent = 8
  ) |> 
  # remove default footnote
  remove_footnote_header(columns = all_stat_cols())

    insira a descrição da imagem aqui Criado em 2025-04-03 com reprex v2.1.1

    • 1

relate perguntas