如果我想执行逻辑回归,我可以使用三种不同的方式:
- 手动:
fit1 <- (glm(outcome~A+B,family=binomial,data=database))
summary(fit1)
exp(confint(fit1))
exp(coefficients(fit1)
- 包
tbl_regression
中的函数gtsumary
fit1 <- (glm(outcome~A+B,family=binomial,data=database))
tbl_regression(fit1, exponentiate = TRUE)
- 包
logistic.display
中的函数epiDisplay
fit1 <- (glm(outcome~A+B,family=binomial,data=database))
logistic.display(fit1)
调整后的 OR 和 p 值与三种方法相同。
但我很惊讶地看到 95%CItbl_regression
和手动计算与 95%CI 不同logistic.display
...
但是,只有来自函数的 95%CIlogistic.display
与 p 值匹配(如果 p 值 > 0.05,则应该包括 1)。
有人知道为什么吗?