create make_table_33_gtsum

Author: aassuied-ps

NAMESPACE

@@ -31,6 +31,7 @@ export(make_table_22)
 export(make_table_32)
 export(make_table_32_gtsum)
 export(make_table_33)
+export(make_table_33_gtsum)
 export(make_table_34)
 export(make_table_35)
 export(make_table_36)

R/fda-table_33.R

@@ -70,3 +70,81 @@ make_table_33 <- function(advs,
 
   tbl
 }
+
+#' @describeIn make_table_33 Create FDA table 33 using functions from `gtsummary`.
+#'
+#' @return
+#' * `make_table_33_gtsum` returns a `gt` object
+#'
+#' @examples
+#' tbl <- make_table_33_gtsum(advs = advs)
+#' tbl
+#'
+#' @export
+make_table_33_gtsum <- function(advs,
+                                adsl,
+                                id_var = "USUBJID",
+                                arm_var = "ARM",
+                                saffl_var = "SAFFL",
+                                lbl_overall = NULL) {
+  checkmate::assert_subset(c(
+    saffl_var, "AVISITN", "PARAMCD", "AVAL", "AVALU", arm_var, id_var
+  ), names(advs))
+  assert_flag_variables(advs, saffl_var)
+
+  advs_all <- advs %>%
+    filter(
+      .data[[saffl_var]] == "Y",
+      AVISITN >= 1,
+      PARAMCD %in% c("DIABP", "SYSBP")
+    ) %>%
+    group_by(.data[[id_var]], PARAMCD) %>%
+    mutate(
+      MAX_DIABP = if_else(PARAMCD == "DIABP", max(AVAL), NA_real_),
+      MAX_SYSBP = if_else(PARAMCD == "SYSBP", max(AVAL), NA_real_)
+    ) %>%
+    ungroup() %>%
+    mutate(
+      SBP90 = formatters::with_label(PARAMCD == "SYSBP" & MAX_SYSBP < 90, "SBP <90"),
+      DBP60 = formatters::with_label(PARAMCD == "DIABP" & MAX_DIABP < 60, "DBP <60")
+    )
+
+  advs_diabp <- advs_all %>%
+    filter(PARAMCD == "DIABP")  %>%
+    distinct(.data[[id_var]], .keep_all = TRUE) %>%
+    select(all_of(id_var), DBP60, ARM, AVALU)
+
+  advs_sysbp <- advs_all %>%
+    filter(PARAMCD == "SYSBP")  %>%
+    distinct(.data[[id_var]], .keep_all = TRUE) %>%
+    select(all_of(id_var), SBP90)
+
+  advs_combined <-
+    full_join(advs_diabp, advs_sysbp, by = "USUBJID") %>%
+    select(SBP90, DBP60, ARM, AVALU)
+
+  avalu <- unique(advs_combined$AVALU)[1]
+  advs_combined <- advs_combined %>% select(-AVALU)
+
+  tbl <- advs_combined %>%
+    tbl_summary(
+      by = arm_var,
+      statistic = list(all_categorical() ~ "{n} ({p}%)"),
+      digits = everything() ~ c(0, 1)
+    ) %>%
+    modify_header(label ~ paste0("**Blood Pressure (", avalu, ")**")) %>%
+    modify_header(all_stat_cols() ~ "**{level}**  \nN = {n}") %>%
+    gtsummary::modify_column_alignment(columns = all_stat_cols(), align = "right")
+
+  if (!is.null(lbl_overall)) {
+    tbl <- tbl %>%
+      add_overall(last = TRUE, col_label = paste0("**", lbl_overall, "**  \n N = {n}"))
+  }
+
+  tbl <- tbl %>% modify_footnote(update = everything() ~ NA)
+
+  gtsummary::with_gtsummary_theme(
+    x = gtsummary::theme_gtsummary_compact(),
+    expr = as_gt(tbl)
+  )
+}