Can gtsummary be used to predict an ordinal variable (several predictors of all kinds in one model), adjusted for confounding factors

Question

i am trying to bulid a prediction model for an oridinal variable. I now that MASS:polr() function tragets this issue but i want to present it in a more approchable way. i thought gtsummary package may be sutiable

my code-

reg_tb<-tbl_uvregression(
  reg_df,
  include = c(a,b,c,d), 
  method = polr,
  y = e,
  exponentiate = TRUE,
  pvalue_fun = ~style_pvalue(.x, digits = 2))

i now that tbl_uvregression() is a univariate model but under 'methods = ' i used the 'polr' option. i suspected polr can't be used in tbl_uvregression() to do an adjusted prediction model because after including 15 predictors, they all remained significant when runing the model (not reasonbele, several factors are strongly associated with each other).