0

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

输出

在此处输入图像描述

4

1 回答 1

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=)参数传递另一个函数来舍入使用科学记数法的估计值，这样打印出来的效果并不疯狂。

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

快乐编程！

于 2020-12-28T23:13:59.993 回答