基于未加权结果编辑含多变量的加权gt表格

要求将未加权样本量n<2的单元格转换为NA。已有针对仅含单一变量行情况的解决方案。

数据 #

df <- data.frame(a= c("m","f","f","m", "m", "f", "f"),b= c("g1","g2","g1","g2", "g2", "g2", "g2"),
                 c=c("c", "c", "c", "d", "c", "e", "e"),
                 d =c("t", "s", "t", "t", "s", "s", "t"),
                 weight = c(1.1, 0.8, 2.2, 4, 3, 0.3, 1.9))

加权表格代码 #

library(gtsummary)
library(srvyr)
library(dplyr)

tbl_summary_object <- df %>%
  as_survey_design(1, weight = weight) %>%
  gtsummary::tbl_svysummary(
    missing="no",
    by=a,
    include=-weight,
    label=list(
      b ~ "group",
      c ~ "group2",
      d ~ "group3"
    ),
    percent="row"
  )

tbl_summary_object

未加权表格代码及结果 #

tbl1 <- df %>%
  gather(key, value, b,c, d) %>%
  group_by(a, value) %>%
  summarise(n=n()) %>%
  spread(a, n) %>%
  dplyr::ungroup() %>%
  arrange(match(value, c("g1", "g2", "c", "d", "e", "s", "t"))) %>%
  rename(stat_1=f, stat_2=m) %>%
  add_row( value=NA,stat_1=NA, stat_2=NA, .before=3) %>%
  add_row( value=NA,stat_1=NA, stat_2=NA, .before=7)

> tbl1
# A tibble: 9 × 3
  value stat_1 stat_2
  <chr>  <int>  <int>
1 g1         1      1
2 g2         3      2
3 NA        NA     NA
4 c          2      2
5 d         NA      1
6 e          2     NA
7 NA        NA     NA
8 s          2      1
9 t          2      2 

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