是否可以使用 tbl_svysummary() 创建分层 table (tbl_strata)？

Question

我对调查数据和 gtsumarry 包还很陌生。我正在尝试使用以下代码从调查数据中创建分层 table，但出现错误“错误：mutate() 输入 tbl 出现问题”。

# Reading the subset of the data

fileUrl <- "https://raw.github.com/Shadi-Sadie/Paper-1-Cancer-Screening-and-Immigrants/master/Cleaned%20Data/subset.csv"
SData<-read.csv( fileUrl , header = TRUE, sep ="," )

# Setting the weight 

options( "survey.replicates.mse" = TRUE)
svy <- svrepdesign(repweights = "PWGTP[0-9]+",
                   weights = ~PWGTP,
                   combined.weights = TRUE,
                   type = "JK1",
                   scale = 4/80, rscales = rep(1, 80),
                   data = SData)

# creating the table

SData %>%
        select(CITG, HICOV, ESRG , EXPANSION) %>%
        tbl_strata(
                strata = CITG,
                .tbl_fun =
                        ~ .x %>% tbl_svysummary(
                                by = EXPANSION,
                                include = c(CITG, HICOV, ESRG , EXPANSION),
                                label = list(CITG ~ "Nativity",
                                             HICOV~ "Any health insurance",
                                             ESRG~ "Employment",
                                             EXPANSION ~ "Expansion" )
                        )
        )

如果可以将 tbl_svysummary() 与 tbl_strata() 一起使用，谁能告诉我哪里做错了？

Answer 1

感谢您使用可重现的更新 post。我做了以下更改：

您正在将数据框传递给 tbl_strata()，它需要更新为调查设计对象。
分层变量不应列在 tbl_summary(include=) 参数中。

编程愉快！

library(gtsummary)
packageVersion("gtsummary")
#> [1] '1.6.0'

fileUrl <- "https://raw.github.com/Shadi-Sadie/Paper-1-Cancer-Screening-and-Immigrants/master/Cleaned%20Data/subset.csv"
SData <- read.csv(fileUrl, header = TRUE, sep = ",")

# Setting the weight
options("survey.replicates.mse" = TRUE)
svy <- survey::svrepdesign(
  repweights = "PWGTP[0-9]+",
  weights = ~PWGTP,
  combined.weights = TRUE,
  type = "JK1",
  scale = 4 / 80, rscales = rep(1, 80),
  data = SData
)

# creating the table
tbl <-
  svy %>%
  tbl_strata(
    strata = CITG,
    .tbl_fun =
      ~ .x %>% tbl_svysummary(
        by = EXPANSION,
        include = c(HICOV, ESRG, EXPANSION),
        label = list(
          HICOV = "Any health insurance",
          ESRG = "Employment",
          EXPANSION = "Expansion"
        )
      )
  )

^{由 reprex package (v2.0.1)}

创建于 2022-05-01

是否可以使用 tbl_svysummary() 创建分层 table (tbl_strata)？

Is it possible to create a stratified table (tbl_strata) using tbl_svysummary()?

r

gtsummary