Is it possible to include multiple random effects? about finalfit HOT 4 CLOSED

# Summary, univariable and multivariable analyses of the form:
# lme4::glmer(dependent ~ explanatory + (1 | random_effect), family="binomial")

explanatory = c("age.factor", "sex.factor", "obstruct.factor", "perfor.factor")
random_effect = "(1 | hospital) + (1 | rx.factor)"
dependent = "mort_5yr"
colon_s %>%
   finalfit(dependent, explanatory, random_effect = random_effect)

colon_s %>%
  glmmixed(dependent, explanatory, random_effect = random_effect) 

colon_s %>%
  finalfit(dependent, explanatory, random_effect = random_effect, keep_fit_id = TRUE) %>% 
    ff_merge(
      explanatory %>% 
        purrr::map_df(~ glmmixed(colon_s, dependent, .x, random_effect = random_effect) %>% 
          fit2df(estimate_suffix = " (univariable with RE)")
        ), 
      last_merge = TRUE
    ) %>% 
  dplyr::relocate(7, .before = 6) # reorder columns

  Dependent: Mortality 5 year          Alive       Died          OR (univariable)  OR (univariable with RE)           OR (multilevel)
                  Obstruction  No 408 (56.7) 312 (43.3)                         -                         -                         -
                              Yes  89 (51.1)  85 (48.9) 1.25 (0.90-1.74, p=0.189) 1.23 (0.83-1.83, p=0.310) 1.23 (0.82-1.83, p=0.313)
                  Perforation  No 497 (56.0) 391 (44.0)                         -                         -                         -
                              Yes  14 (51.9)  13 (48.1) 1.18 (0.54-2.55, p=0.672) 1.07 (0.44-2.57, p=0.888) 1.02 (0.42-2.46, p=0.969)