Study Design and Participants
This cross-sectional study was conducted among full employees and contractual workers of the Esfahan Steel Company (workforce of 16,000). The sample size was determined based on an epidemiological survey of chronic diseases in manufacturing employees by Roohafza et al. (ESCOME) (33) carried out to estimate the prevalence of psychological disorders (depression and anxiety) among the workforce. The sample size was estimated to be 3500 by considering 0.1 (19, 20), 0.05, and 0.01 as prevalence of psychological disorders, type one error rate, and sampling error rate, respectively.
Three thousand and sixty-three volunteers returned complete questionnaires (response rate: 0.87) and were included in statistical analysis. The inclusion criterion was work experience for at least one year and willing and agreeing to participate in the study. Volunteers who did not answer more than 10 % of the questionnaire were excluded from analysis.
We applied multi-stage cluster sampling, in which clusters were the main seven departments and their sections, stratified by job categories. Sample sizes in the clusters and strata were proportional to the size of respective departments. Due to a low number of women workers (n=800), we relied on convenience sampling to recruit women volunteers (n=260).
Demographic data were gathered through self-administered questionnaires at the company premises with the help of study coordinators, who also monitored questionnaire administration over the six months of data gathering. The data were quality checked for inclusion and exclusion criteria as they were entered in the computer for statistical analysis. The study design and its implementation has been presented elsewhere in more detail (33,34).
All participants were informed about the study protocol and signed informed consent. Medical research ethics committee of the Isfahan University of Medical Sciences approved the study protocol (projects numbers #87115 and #395482).
Study instruments and assessment of variables
Quality of life
A standardized and generic form of Euro QoL-five dimensions (EQ-5D-3L) (35,36) was used for assessing QoL. This self-report instrument comprises the two parts of self-classifier (descriptive system) and the visual analog scale (VAS).
The EQ-VAS was considered as a single value for assessing of overall health status. This is a visual analog scale (VAS), ranging from 0 (worst imaginable health) to 100 (best imaginable health). Respondents by marking the scale with a single vertical mark, rate their current health status. EQ-5D self-classifier describes health state of subjects in five domains: mobility, self-care care, normal activities, pain/discomfort, anxiety and depression. Each domain consists of three levels: no problem (1), some problem (2), severe problem (3). We combined the last two categories into a single category in latent class analysis (LCA) due to poor response rate in last category. Higher EQ-5D scores represent worse health status. EQ-5D has shown a good reliability and validity (37,38). Internal reliability of the questionnaire was assessed in a pilot sample of 300 participants in current study and Cronbach alpha was obtained to be 0.51.
Stressful Life Events
The number and intensity of experienced life stressors were measured by Stressful Life Events questionnaire (SLE) (39). Participants were asked about the occurrence of stressors within the past 6 months. The SLE is a 46-item self-administered scale which consists of eleven dimensions, including home life (measured with 7 items), financial problems (5 items), social relations (4 items), personal conflicts (5 items), job conflicts (4 items), educational concerns (4 items), job security (5 items), loss and separation (4 items), sexual life (4 items), daily life (2 items), and health concerns (2 items). The items are rated on a 6-point Likert scale (0: never, 1: very mild, 2: mild, 3: moderate, 4: severe, 5: very sever) and the higher score indicates higher stress level. Internal consistency of the SLE questionnaire was 0.92(39).
Assessment of other variables
Variables that were considered as potential confounders included demographics [age (years), gender (male/female), marital status (married/single), education (0–5 years / 6–12 years / over 12 years)], lifestyle variables [sleep duration (hours) and physical activity (hours per week) ,BMI (weight (kg)/ height (m2))], job-related variables [job stress (effort-reward imbalance), and second job (yes/no)]. Physical activity was evaluated with the International Physical Activity Questionnaire (IPAQ), which included 11 questions (40). The internal reliability of this questionnaire was reported good by Moghaddam et al. (41), based on Cronbach's alpha of 0.7 and Spearman Brown correlation coefficient of 0.9.
Statistical analysis
Two level latent class analysis (42–44) (employees as level 1 and job categories level 2 units) was employed to identify homogeneous latent classes of employee according to the their responses to 5 indicators of QoL. Multilevel latent class analysis is an extension of traditional latent class analysis (LCA) that handles situations where there is a multilevel construct (in our study, employees (level 1) are nested within job categories (level 2)).
The two-level latent class model not only classifies the employees but also the job categories based on the distribution of QoL of employees nested in job categories. The modeling process consists of the following steps: In the first step, with ignoring the nesting structure of data, LCA was used to classify individuals based on their response to the 5 items of the EQ-5D. To determine the appropriate number of classes at the employee-level, we started by fitting a one-class model and sequentially increased the number of classes until to yield the best fit. Bayesian Information Criterion (BIC) (45) was used to select the best-fitting model. Additionally, interpretability of the identified classes was also considered.
In the next step, two-level latent class model was estimated to take (the clustered nature of the data) into account the multilevel structure of our data. A two-level model with two classes at level one and two classes at the second level was selected as the best fitted model.
Finally, we used multilevel latent class regression (MLLCR) for evaluating the predictors of QoL. The proposed predictors of QoL in current study were stressors. For evaluating their association with QoL, at first, we used exploratory factor analysis on eleven domains of life stressors and two factors were extracted and labeled as “personal stressors” and “socioeconomic stressors” and these domains were used as latent predictors of QoL of employees.
We adjusted for the effects of potential confounding variables that were statistically significant at α=0.1 in univariate analyses.
Continuous and categorical variables were represented as mean (SD) and number (percentage) respectively. Independent samples t-test / Mann-Whitney U test and Chi-square test were used for comparing continuous and categorical variables between studied groups, respectively. MLLCR was fitted in Mplus 8 statistical software.