0

以下代码生成 N 列“N”和“事件 N”作为单变量回归表的一部分。我有一个案例控制数据集,我希望让“案例”和“控制”列给出案例和控制的数量。

“案例”和“控制”由下面代码中的变量“响应”确定。例如 response(1) = "cases",而 response(0) = "controls"。

我怎样才能做到这一点?

谢谢,耐莉

tbl_uv_nevent_ex <-
trial[c("response", "trt", "age", "grade")] %>%
tbl_uvregression(
method = glm,
y = response,
method.args = list(family = binomial)
) %>%
add_nevent()
4

1 回答 1

2

您可以通过将控件的数量添加到.$table_body数据框中来做到这一点。我在下面提供了一个示例。目前有一个症结......在将控件数量添加到将要打印的数据框后,我们需要将新列添加到内部指令集以在 gtsummary 中打印。这一步现在很头疼,但我们正在研究一种解决方案,让用户可以访问它。这是该功能的初稿:http ://www.danieldsjoberg.com/gtsummary/dev/reference/modify_table_header.html

同时,您可以通过以下方式完成此操作:

library(gtsummary)

tbl <-
  trial[c("response", "trt", "age")] %>%
  tbl_uvregression(
    method = glm,
    y = response,
    method.args = list(family = binomial),
    exponentiate = TRUE
  ) %>%
  add_nevent()

# add the number of controls to table
tbl$table_body <-
  tbl$table_body %>%
  dplyr::mutate(
    n_nonevent = N - nevent
  ) %>%
  dplyr::relocate(n_nonevent, .after = nevent)

# updating internal info with new column (this part will not be required in the future)
tbl$table_header <- 
  gtsummary:::table_header_fill_missing(tbl$table_header, 
                                        tbl$table_body)

# print tbl with Case and Control Ns
tbl %>%
  modify_header(
    list(nevent ~ "**Case N**",
         n_nonevent ~ "**Control N**")
  )

在此处输入图像描述

于 2020-10-16T00:09:28.660 回答