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

【问题标题】：Problem with 95% CI output - gtsummary package in R95% CI 输出的问题 - R 中的 gtsummary 包
【发布时间】：2020-12-28 21:26:02
【问题描述】：

有时，我使用 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)

【问题讨论】：

有没有办法可以发布一个可重现的例子？你也可以分享broom::tidy(fit3)的输出吗？我们最近在内部进行了巨大的更改，这将更好地识别具有奇怪名称的变量。请查看开发版remotes::install_github("ddsjoberg/gtsummary")是否仍然存在问题。我们将在下周 CRAN 假期归来时发布新版本。
@DanielD.Sjoberg，感谢您的快速回复。我添加了 broom::tidy(fit3) 的屏幕截图。处理可重现的示例，但不确定我是否可以上传，因为我在锁定平台上工作。
我了解您是否无法上传此特定示例，但类似的示例会有所帮助（至少是相同的模型类型，希望有相同的问题）。模型类型是什么？开发版的 gtsummary 问题是否仍然存在？
您能否更新broom::tidy() 输出以包含列标题？
感谢您添加更多信息。我将在下面发布答案。

标签： r confidence-interval gtsummary

【解决方案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=) 参数传递另一个函数来对使用科学记数法的估计值进行四舍五入，这样打印出来的效果不会很疯狂。

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

编程愉快！

【讨论】：