tb2 <- glm(formula_current, data = train, family = "binomial") |>
tbl_regression(exponentiate = TRUE, pvalue_fun = ~style_sigfig(., digits = 3))
tb1 <- train[, c("phando2", all_feats)] |>
tbl_uvregression(
method = glm,
y = phando2,
method.args = list(family = binomial),
exponentiate = TRUE,
hide_n = TRUE
)