有没有办法只查看 tbl_summary 中某个因子变量的特定级别？

Question

我正在尝试比较并获取因子变量的 p 值，按 trt/control 分层。

为了可重复性，这是我如何加载数据的简化版本：

# DEMO
demo <- nhanes('DEMO')
demo_vars <- names(demo)
demo2 <- nhanesTranslate('DEMO', demo_vars, data = demo)

# PRESCRIPTION MEDICATIONS
rxq_rx <- nhanes('RXQ_RX')
rxq_rx_vars <- names(rxq_rx)
rxq_rx2 <- nhanesTranslate('RXQ_RX', rxq_rx_vars, data = rxq_rx)
rxq_rx2 <- rxq_rx2 %>% select("SEQN", "RXD240B") %>% filter(!is.na(RXD240B)) %>% group_by(SEQN) %>% dplyr::summarise(across(everything(), ~toString(na.omit(.))))

nhanesAnalysis = join_all(list(demo2, rxq_rx2), by = "SEQN", type = "full")

# Reconstructing survey weights for combining 1999-2018 - Combining ten survey cycles (twenty years) 
nhanesAnalysis$wtint20yr <- ifelse(nhanesAnalysis$SDDSRVYR %in% c(1,2), (2/10 * nhanesAnalysis$WTINT4YR), # for 1999-2002
                                          (1/10 * nhanesAnalysis$WTINT2YR)) # for 2003-2018


# sample weights 
nhanesDesign <- svydesign(id      = ~SDMVPSU,
                          strata  = ~SDMVSTRA,
                          weights = ~wtint20yr,
                          nest    = TRUE,
                          data    = nhanesAnalysis)

# subset
obese_adults  <- subset(nhanesDesign, (obesity == 1 & !is.na(BMXBMI) & RIDAGEYR >= 60))

我正在使用以下 tbl_summary() 函数来获取此表。

  ageDesign %>%
  tbl_strata(
    strata = obesity,
    .tbl_fun =
      ~ .x %>% tbl_svysummary(
        by = wlp_yn,
        percent = "row",
        include = c(phys_func2, standing_chair)
      )%>%
      add_overall() %>%
      add_p()
  )

我只想看因素“standing_chair”的“难度”级别。有什么办法可以做到这一点？谢谢！