Additional file 1: Table S1. Descriptive statistics for the sociodemographic and occupational variables (N = 40,986). Table S2. Multinomial logistic regressions exploring the sociodemographic and occupational associations with the identified classes of health (risk) behaviours for men (N = 25,788). Percentages are weighted with conditional probability weights. Unadjusted multinomial odds ratios (MOR) with 95% confidence intervals (CIs) are shown. Table S3. Multinomial logistic regressions exploring the sociodemographic and occupational associations with the identified classes of health (risk) behaviours for women (N = 15,198). Percentages are weighted with conditional probability weights. Unadjusted multinomial odds ratios (MOR) with 95% confidence intervals (CIs) are shown. Fig S1. Plot showing AIC, BIC, and SSABIC model fit criteria for each additional class.