95%CI输出问题-R中的gtsummary包

Question

有时，我会使用gtsummary包找到像这样的95% CI输出(见下图)。有时我会注意到这是因为它不喜欢变量名中的冒号、空格或撇号，或者因为在因子级别中有NA。当我更改名称时，CI是固定的。然而，这一次似乎不是这样。你知道这个输出的原因是什么吗？

谢谢

fit3 <- glmer(cclintr ~ SEX + Rg + agroup + death +
              psource_group + region + BMI + HOSPVENT + 
              PATMANICU + fibandflutter + MEDHISTO_03 + 
              MEDHISTO_35 + Cerebrovasc_group
              smo_vape + DOCUSYMP_12 + admonth + (1|FACILITY_DISPLAY_ID),
             data= filthosp,
             family = binomial,
             nAGQ = 1, 
             control=glmerControl(optimizer="bobyqa",optCtrl=list(maxfun=2e5)))

#Summary into table 
library(gtsummary)
f2model_tab <- fit3 %>%
  tbl_regression(exponentiate = TRUE ) %>%
  bold_labels() %>%
  bold_p(t= 0.05) 

​

​

​

​

Answer 1

从我们在这里看到的，gtsummary正在报告它应该报告的内容。问题是admonth变量的BM级别的置信度非常大。系数为-13.8，标准误差为74。如果我们自己计算置信上限，它将是：

-13.8 + 1.96 * 74 = 131.24

然后我们求幂，得到赔率比：

exp(131.24) = 9.92e+56    !!

默认情况下，tbl_regression()函数使用style_ratio()函数对边界进行格式化和四舍五入，该函数不使用科学记数法(因此会打印巨大的数字)。您可以使用tbl_regression(estimate_fun=)参数传递另一个函数来舍入使用科学记数法的估计值，这样打印结果就不会太夸张。

除此之外，我建议您尝试一下您的模型，看看您是否可以获得更合理的估计。

祝您编程愉快！