R 中卡方和 T 学生的置信区间

Question

所以我在 R 中有这段代码...

Tabla2 <- tbl_svysummary(data = Beta18.1Pond, by = CTBC, missing = "no", 
                         digits = list(all_continuous() ~ c(2,2)), 
                         include = c(SR, QS23, QS23R, QSSEXO, QS500, QS501U, QS501C, QS50328,
                                     QS505A, QS505B, QS505C, QS505D, QS506, QS2425N, QS25AA, 
                                     QS26, QS27R, QS102, QS109, HV270, Region, HV025, HV009, HV115, 
                                     CTBC)) %>% 
  add_p(test = list(all_continuous() ~ 
                      "svy.t.test", all_categorical() ~ 
                      "svy.adj.chisq.test"), include = everything()) %>%
  modify_header(label ~ "**Variable**") %>%
  modify_spanning_header(c("stat_1", "stat_2") ~ 
                           "**Conocimiento sobre tuberculosis pulmonar**")  %>%
  modify_caption("**Tabla 2. Análisis descriptivo**") 
Tabla2

我想添加一个命令来计算表中卡方的置信区间（在本例中为比例差异的置信区间）和 t 检验。有人知道我该怎么做吗？

感谢您的回复。

Answer 1

对于 t 检验，您可以简单地使用 add_ci()，因为它旨在处理连续数据。 对于分类数据，我认为您需要一些将卡方测试结果转换为允许置信区间计算的度量，例如将卡方测试转换为比例测试（如果适用）：

Tabla2 <- tbl_svysummary(data = Beta18.1Pond, by = CTBC, missing = "no", 
                         digits = list(all_continuous() ~ c(2,2)), 
                         include = c(SR, QS23, QS23R, QSSEXO, QS500, QS501U, QS501C, QS50328,
                                     QS505A, QS505B, QS505C, QS505D, QS506, QS2425N, QS25AA, 
                                     QS26, QS27R, QS102, QS109, HV270, Region, HV025, HV009, HV115, 
                                     CTBC)) %>% 
  add_p(test = list(all_continuous() ~ 
                      "svy.t.test", all_categorical() ~ 
                      "svy.adj.chisq.test"), include = everything()) %>%
  add_ci(conf.level = 0.95, test = list(all_continuous() ~ "svy.t.test")) %>% # Add CI for t-tests
  modify_header(label ~ "**Variable**") %>%
  modify_spanning_header(c("stat_1", "stat_2") ~ 
                           "**Conocimiento sobre tuberculosis pulmonar**")  %>%
  modify_caption("**Tabla 2. Análisis descriptivo**") 
Tabla2