trythis

library(gtsummary)
library(survey)

design.ACS2023 <- svydesign(
  strata = ~STRATA,
  id = ~SERIAL,
  weights = ~PERWT,
  nest = TRUE,
  data = data_ACS2023
)

design.ACS2023.NPs <- subset(design.ACS2023, OCC == 3258)

table1 <- tbl_svysummary(
  design = design.ACS2023.NPs,
  include = c(gender_f, AGE, race_f, hispanic_f, empstat_f, UHRSWORK, INCWAGE),
  missing = "no"
)