Identifying the risk features for occupational stress in medical workers: a cross-sectional study.

OBJECTIVE: Occupational stress is considered a worldwide epidemic experienced by a large proportion of the working population. The identification of characteristics that place people at high risk for occupational stress is the basis of managing and intervening in this condition. In this study, we aimed to identify and validate the risk features for occupational stress among medical workers using a risk model and nomogram.
METHODS: This cross-sectional study included 1988 eligible participants from Henan Province in China. Occupational stress and worker-occupation fit were measured with the Depression, Anxiety and Stress Scales (DASS-21) and Worker-Occupation Fit Inventory (WOFI). The identification of risk features was achieved through constructing multiple logistic regression model, and the risk features were used to develop the risk model and nomogram. Receiver operating characteristic (ROC) curves and calibration plots were generated to assess the effectiveness and calibration of the risk model.
RESULTS: Among 1988 participants in our study, there were 42.5% (845/1988) medical workers experienced occupational stress. The risk features for occupational stress included poor work-occupation fit (WOF score < 25, expected risk: 77.3%), nurse population (expected risk: 63.1%), male sex (expected risk: 67.2%), work experience duration of 11-19 years (expected risk: 54.5%), experience of a traumatic event (expected risk: 65.3%) and the lack of a regular exercise habit (expected risk: 60.2%). For medical workers who have these risk features, the expected risk probability of occupational stress would be 90.2%.
CONCLUSION: The current data can be used to identify medical workers at risk of developing occupational stress. Identifying risk features for occupational stress and the work-occupation fit can support hierarchical stress management in hospitals.

Introduction

Occupational stress is generally acknowledged as a global phenomenon that occurs in both developed and developing countries and has significant health and economic consequences (Suleman et al. 2018). In 2019, the State Council of the People's Republic of China issued a public policy entitled Healthy China 2030 (2019–2030) with the aim of promoting mental health screening, and supervisory departments need to give more attention to students and relevant professions. Healthy China 2030 (2019–2030) requires organization managers to regularly assess and take steps to prevent occupation-related psychological disorders (General Office of the State Council of PRC 2019). Occupational stress is a typical occupation-related psychological disorder. At present, the features of groups at high risk for occupational stress have not been determined. It is worthwhile to determine the typical features of groups at high risk, as these risk features can serve as the basis of management in the workplace. Therefore, identification of the risk features and development of tools that can predict the risk of occupational stress are needed to improve targeted prevention and intervention measures.

Comparing the occupational stress level among different professions, scientific researchers (Kageyama et al. 2001), policemen (Queirós et al. 2020) and petroleum industry workers (Ning et al. 2020) would be considered as high-risk professions, especially medical workers (Lexén et al. 2020; Tsai and Liu 2012). Many studies have reported that medical workers have a high-to-moderate level of occupational stress (Chatzigianni et al. 2018; Kaburi et al. 2019; Nabirye et al. 2011; Tadesse et al. 2016). In the United States, the report of Exposure to Stress-Occupational Hazards in Hospitals by the National Institute for Occupational Safety and Health showed that the common stressors experienced by medical workers in hospital settings consist of inadequate staffing levels, long work hours, shift work, role ambiguity and exposure to infectious diseases (National Institute for Occupational Safety and Health 2008). In the United Kingdom, a report indicated that higher staff turnover, more absence due to illness, decreased performance, and more complaints and grievances were signs of stress (Health and Safety 2019). The above-mentioned report showed that medical workers have a high level of occupational stress, and many countries have begun to give more attention to the prevention of occupational stress in hospitals.

Previous studies on occupational stress have mainly referred to three aspects: the prevalence of stress in different occupations (Araújo et al. 2020; Choy and Wong 2017; Cordioli et al. 2019; van der Wal et al. 2018), the associations between stress and related diseases (Han et al. 2019; Sara et al. 2018; Yang et al. 2018) and stress interventions (Babanataj et al. 2019; Basu et al. 2017; Bresesti et al. 2020; Clough et al. 2017; Nowrouzi et al. 2015). These findings are insufficient to serve as the foundation for the development of targeted intervention strategies because we still lack effective methods of distinguishing subgroups with different levels of the risk of occupational stress, and management depends on understanding these subgroups. The demographic factors, occupational factors and health behaviour factors among medical workers who are at high risk for occupational stress need to be identified. Moreover, a predictive tool that can be used to evaluate the risk of occupational stress in the work force is needed. We believe that a stress nomogram would be a relatively better predictive tool for implementation in hospitals.

The structure of work is changing at a rapid pace and the mismatch between workers and their work environments poses a threat to their health. In this study, we propose a new factor affecting occupational stress, namely worker-occupation fit (WOF). In our early pilot study, we constructed a theoretical model of the effects of WOF on occupational stress and related disorders (Sun 2020a). In addition, we also reported the association between WOF and occupational stress: a lower level of WOF was associated with a higher level of occupational stress, indicating that occupational stress has a strong negative correlation with the level of WOF (Sun 2020b). WOF was defined in accordance with the theoretical concept of occupational stress as the match between a worker's characteristics, needs, and abilities and the culture, supplies and demands of the occupational environment. There are three types of WOF: characteristic fit, need-supply fit and demand-ability fit. A WOF misfit can have negative effects on society and on the mental health, physical health or behaviour of workers. As a result, low-level WOF can increase the risk of occupational stress incidence. In this study, we expanded the sample size to explore whether a low-level WOF is a risk feature for occupational stress. The findings are expected to provide new avenues for research on occupational stress and stress prevention in hospitals.

Nomograms, which are graphical illustrations of a mathematical model in which different factors are combined to predict a definite endpoint, have been utilized as convenient and reliable tools for the prediction of the survival and mortality of cardiovascular disease and cancer patients (Caulfield et al. 2018; Jalali et al. 2019; Zi et al. 2020). Researchers have constructed nomograms for the prediction of depression and anxiety in cancer outpatients (Lima et al. 2016). Guo et al. developed and validated a prognostic nomogram and evaluated its discrimination ability to improve the prediction of 30-day survival among critically ill myocardial infarction (MI) patients (Guo et al. 2020). In that study, the nomogram exhibited reasonably accurate discrimination (AUC: 0.765, 95% CI 0.716–0.814) and calibration (C-index: 0.758, 95% CI 0.712–0.804) in the validation cohort. In another study, Pan et al. constructed a predictive nomogram for mortality in patients with COVID-19, and the areas under the ROC curves for the nomogram were 0.988 (95% CI 0.972–1.000) and 0.956 (95% CI 0.874–1.000) in the primary and validation groups, respectively (Pan et al. 2020). Nomograms have relatively better performance for the prediction of outcomes in patients with cardiovascular disease and cancer, and it is worthwhile exploring whether a nomogram can perform well with regard to identifying groups of medical workers at high risk for occupational stress. Currently, the occupational stress scale is the major method used to assess whether workers are experiencing occupational stress. Workers who are identified as experiencing stress then directly enter the intervention stage. Because a nomogram could identify the high-risk group at an early stage, hospital administrators could implement control and intervention measures in a timely manner to prevent the occurrence of occupational stress. Thus, identifying the risk features for occupational stress in medical workers and being able to regularly screen workers could facilitate the early recognition of this problem, which would support the prevention and treatment of occupational stress and improve the quality of health care provided to patients.

Accordingly, this study aimed to identify and validate the risk features for occupational stress among medical workers based on a risk model and nomogram. The results may provide evidence and technical support that can be used when implementing a model of hierarchical management for occupational stress in hospitals.

Materials and methods

Study design and participants

This cross-sectional study was carried out from October 2019 to June 2020. A multicentre survey was conducted among medical workers from four targeted hospitals (1 general hospital and 3 specialized hospitals) in Henan, China.

The inclusion criteria were as follows: (1) individuals who had more than 1 year of work experience; (2) individuals who had no family history of and were not taking medications for a mental disorder; (3) individuals who were regular employees of the targeted hospitals; and (4) individuals who were willing to participate in this survey. The exclusion criteria were as follows: (1) participants who were absent during the survey and (2) participants whose questionnaires were ineligible or were < 80% completed.

There were 2050 medical workers in the surveyed hospitals. Nine participants did not agree to participate in the survey, and 32 participants were absent from the hospital during the survey. A total of 2009 questionnaires were distributed, and 21 participants returned incomplete (< 80%) questionnaires. Finally, 1988 questionnaires were considered valid, corresponding to a response rate of 97.0% (1988/2050).

From among the 1988 eligible participants, we randomly selected 80% of the subjects (n = 1568) as the training set to establish the risk model and nomogram, and the remaining 20% of the subjects (n = 420) as the validation set to verify the risk model. The study flowchart is presented in Fig. 1.

Fig. 1

Study flowchart

Data collection tools and measurements

A three-section cross-sectional survey was administered by researchers and investigators to medical workers in the four targeted hospitals. The following survey instruments were used to collect data:

Basic information

A self-administered basic information questionnaire was divided into demographic and occupational characteristics for each participant. The demographic characteristics included sex (male or female), age, marital status (single, cohabiting, married, divorced and living alone, widowed and living alone), children’s grade level (no children, infant stage, primary school, junior-senior high school, graduate school or higher), education (diploma or less, bachelor’s degree or higher), monthly income (< 3000 yuan, 3000–3999 yuan, 4000–4999 yuan, 5000–5999 yuan, ≥ 6000 yuan) and whether the participant had a healthy lifestyle (e.g., exercise, smoking and drinking). The occupational characteristics included work experience, department, professional title, night shift frequency and work hours.

Occupational stress investigation

Occupational stress is defined as harmful physical and emotional responses that occur when the requirements of a job do not match the capabilities, resources or needs of the worker (National Institute for Occupational Safety and Health 1999). Occupational stress was measured with the Depression, Anxiety and Stress Scale (DASS-21), which was developed by Lovibond (Lovibond and Lovibond 1995) and has been proven to be a valid measure for the assessment of occupational stress in the Chinese population (Lu et al. 2018). Furthermore, we emphasized that each response to stress was related to occupational features in hospitals. In this study, we used the stress subscale (seven items: 1, 6, 8, 11, 12, 14 and 18) from the DASS-21. Each item was scored from 0 (it does not apply to me at all in the last week) to 3 (it applies to me perfectly in the last week). The DASS-21 score was multiplied by two to calculate the final subscale score, yielding a maximum of 42 points (Lovibond and Lovibond 1995). Participants with a total score from 0 to 14 were considered to have a normal level of stress, those with scores from 15 to 18 were considered to have mild stress, those with scores from 19 to 25 were considered to have moderate stress, and those with scores greater than 26 were considered to have severe stress (Lovibond and Lovibond 1995).

Worker-occupational fit inventory (WOFI)

WOF is an important factor to consider in the assessment of occupational stress, but it is still neglected in occupational stress research. The Worker-Occupation Fit Inventory (WOFI) was based on the questionnaire by Cable and DeRue (Cable and DeRue 2002), with some items and descriptions adapted for the Chinese culture and way of thinking to facilitate its suitability for Chinese participants. The WOFI consists of three sections with a total of 9 items: (1) characteristic fit, for example, “Do you think your job style fits your job?”; (2) need-supply fit, for example, “Do you think the job provides what you need?”; and (3) demand-ability fit, for example, “Do you think that your educational background meets the job demands?”. Each item is rated on a 5-point Likert-like scale (1 = extremely poor fit, 2 = poor fit, 3 = general fit, 4 = good fit, 5 = extremely good fit), and the total score is the sum of the scores for the 9 items, ranging from 9 to 45 points. The higher the score is, the better the WOF. The WOFI showed satisfactory internal consistency (Cronbach's α = 0.888). The items used in this section of the survey are described in Table S1.

Data analysis

Descriptive statistics were used to examine variable characteristics. Categorical variables are expressed as frequencies (%) and were compared via chi-square analysis. In this study, we constructed a risk model for occupational stress using a training set, with a binary outcome (“0” indicated an occupational stress score ≤ 14, and “1” indicated an occupational stress score ≥ 15) included in the multiple logistic regression model, and those variables (including sex, age, marital status, children’s grade level, education, monthly income, exercise, smoking, drinking, work experience, department, professional title, night shift frequency and work hours) with statistical significance in the univariate analysis were ultimately included in the multiple logistic regression model. The incidence risk of occupational stress was expressed as odds ratios (ORs) and 95% confidence intervals (95% CIs). Based on the risk model, a nomogram was constructed using R software with the rms package. The discriminatory ability of the model was quantified using the area under the ROC curve (AUC) for internal validation and external validation. Accuracy and the ROC curve were used to reflect the performance of the selected predictors. The calibration of the nomogram was assessed by plotting the observed outcome probabilities and the probabilities predicted by the logistic model in the training set and validation set. Based on the risk model for occupational stress, we calculated the prevalence and used calibration plots to assess the effectiveness of the models at different levels of risk. A two-tailed P value < 0.05 was considered statistically significant. All the above processes were performed using R software (version 4.0.2).

Results

Participant characteristics

The respondents’ demographics and occupational characteristics are shown in Table 1. Among the entire sample of 1988 participants, the mean age was 32.7 (SD = 7.8) years, with 214 (7.7%) males and 1774 (92.3%) females. The mean WOF score was 34.3 (SD = 6.1). A total of 26.4% (n = 525) of the respondents were physicians, and 73.6% were nurses. The majority of the participants had a bachelor’s degree (n = 1519, 76.4%) or higher. A total of 10.8% (n = 214) of the participants reported working a ≥ 60-h work week, and a monthly night shift frequency ≥ 6 (n = 928, 46.7%) was the most common.

Table 1

Demographic and occupational characteristics of the 1988 participants

Variable	Groups	n	%
Sex	Male	214	7.7
	Female	1774	92.3
Occupational categories	Physician	525	26.4
	Nurse	1463	73.6
Marital status	Single	511	25.7
	Cohabiting	19	1.0
	Married	1423	71.6
	Divorced and living alone	31	1.5
	Widowed and living alone	4	0.2
Grade level of children	None	668	33.6
	Infant stage	455	22.9
	Primary school	573	28.8
	Junior-senior high school	124	6.2
	Graduate school or higher	168	8.5
Education	Diploma or less	469	23.6
	Bachelor’s degree or higher	1519	76.4
Monthly income	< 3000 yuan	334	16.8
	3000–3999 yuan	686	34.5
	4000–4999 yuan	450	22.6
	5000–5999 yuan	315	15.9
	≥ 6000 yuan	203	10.2
Exercise	Yes	818	41.2
	No	1170	58.8
Smoking	Yes	81	4.1
	No	1907	95.9
Drinking	Yes	159	8.0
	No	1829	92.0
Traumatic events	Yes	570	28.7
	No	1418	71.3
Event types	No	1427	71.8
	Family event	331	16.7
	Occupation event	114	5.7
	Occupation-family event	116	5.8
Work experience*	< 1 year	150	7.6
	1–3 years	297	14.8
	4–10 years	757	38.1
	11–19 years	512	25.8
	≥ 20 years	272	13.7
Department of	Internal Medicine	355	17.9
	Surgery	325	16.3
	Obstetrics and Gynaecology	237	11.9
	Paediatric Surgery	83	4.2
	Psychiatry	248	12.5
	Infectious Diseases	111	5.6
	Emergency	140	7.0
	ICU	136	6.8
	Outpatient Clinic	79	4.0
	Others	274	13.8
Professional title	Primary	1125	56.6
	Intermediate	756	38.0
	Senior	107	5.4
Night shift frequency (times per month)	None	691	34.8
	≥ 2	194	9.7
	≥ 4	175	8.8
	≥ 6	928	46.7
Work hours (hours per week)	< 40	534	26.8
	40–49	940	47.3
	50–59	300	15.1
	≥ 60	214	10.8
WOFI score	< 25	352	17.7
	25–29	507	25.5
	30–34	381	19.2
	35–39	434	21.8
	≥ 40	314	15.8

WOF worker-occupation fit

*Classification of the number of years of work experience based on the hierarchy of training for medical workers

Proportion of participants with occupational stress

Of the 1988 medical workers, 42.5% (n = 845) experienced mild to severe levels of occupational stress. A total of 17.8% (n = 355) of the medical workers had mild stress, 15.1% (n = 300) had a moderate level of stress, and 9.6% (n = 190) had a severe level of stress (Fig. 2).

Fig. 2

Proportions of participants in each occupational stress category (%)

Incidence of occupational stress in medical workers with different characteristics

Table 2 shows the incidence of occupational stress in subgroups with different characteristics in the training set (n = 1568) and the validation set (n = 420). In the training set, occupational stress significantly differed among subgroups with different characteristics (P < 0.05). Higher detection rates of occupational stress were observed in participants who were males (53.8%, n = 91), were nurses (47.1%, n = 549), were aged 30–34 years (52.6%, n = 247), were widowed and living alone (87.5%, n = 7), had a bachelor’s degree or higher (44.7%, n = 529), did not exercise (50.5%, n = 462), smoked (63.6%, n = 49), drank (55.6%, n = 74), had a monthly income greater than 6000 yuan (46.9%, n = 76), had worked for 11–19 years (49.1%, n = 200), worked in a psychiatry department (49.5%, n = 103), had an intermediate professional title (47.9%, n = 286), worked more than 6 night shifts per month (47.0%, n = 340), worked more than 60 h per week (55.3%, n = 94) and had a WOFI score of 30–34 points (60.1%, n = 179). All differences were significant (chi-square test, P < 0.05). In the validation set, occupational stress did not differ among subgroups stratified by sex, age, marital status, education, income, smoking, drinking, professional title, night shifts and work hours (chi-square test, P > 0.05).

Table 2

Comparison of occupational stress among 1988 medical workers with different characteristics in the training and validation sets

Characteristic	Training set (n = 1568)			P*	Validation set (n = 420)			P*
	n	Positive	%		n	Positive	%
Sex				0.021				0.431
Male	169	91	53.8		47	23	48.9
Female	1399	585	41.8		373	160	42.9
Occupational categories				< 0.001				< 0.001
Physician	403	127	31.5		112	32	28.6
Nurse	1165	549	47.1		308	151	49
Age group (years)				< 0.001				0.052
> 25	272	89	32.7		74	26	35.1
25–29	302	120	39.7		73	32	43.8
30–34	470	247	52.6		124	66	53.2
35–39	276	126	45.7		89	39	43.8
≥ 40	248	94	37.9		60	20	33.3
Marital status				0.012				0.257
Single	409	150	36.7		105	42	40
Cohabiting	19	12	63.2		4	0	0
Married	1105	496	44.9		304	138	45.4
Divorced and living alone	27	11	40.7		7	3	42.9
Widowed and living alone	8	7	87.5		0	0	0
Education				0.019				0.093
Diploma or less	385	147	38.2		93	32	34.4
Bachelor’s degree or higher	1183	529	44.7		327	151	101.2
Monthly income				0.024				0.260
< 3000 yuan	260	89	34.2		78	26	33.3
3000–3999 yuan	548	245	44.7		133	61	45.9
4000–4999 yuan	343	148	43.1		106	48	45.3
5000–5999 yuan	255	118	46.3		63	32	50.8
≥ 6000 yuan	162	76	46.9		40	16	40
Exercise				< 0.001				< 0.001
Yes	654	214	32.7		169	55	32.5
No	914	462	50.5		251	128	51
Smoking				0.007				0.956
Yes	77	49	63.6		14	6	42.9
No	1491	627	42.1		406	177	43.6
Drinking				0.024				0.185
Yes	133	74	55.6		33	18	54.5
No	1435	602	42.0		387	165	42.6
Traumatic events				< 0.001				< 0.001
Yes	434	250	57.6		140	79	56.4
No	1134	426	37.6		280	104	37.1
Work experience				< 0.001				0.046
< 1 year	121	32	26.4		33	11	33.3
1–3 years	238	87	36.6		62	20	32.3
4–10 years	582	273	46.9		173	87	50.3
11–19 years	407	200	49.1		101	46	45.5
≥ 20 years	220	84	38.2		51	19	37.3
Department				0.021				0.003
Internal medicine	279	125	44.8		70	37	52.9
Surgery	258	115	44.6		69	33	47.8
Obstetrics and gynaecology	184	68	37.0		54	22	40.7
Paediatric surgery	63	31	49.2		22	10	45.5
Psychiatry	208	103	49.5		42	24	57.1
Infectious diseases	87	30	34.5		19	6	31.6
Emergency	105	49	46.7		33	19	57.6
ICU	109	54	49.5		29	12	41.4
Outpatient clinic	64	20	31.3		17	1	5.9
Others	211	81	38.4		65	19	29.2
Professional title				0.008				0.567
Primary	885	354	40.0		239	99	41.4
Intermediate	597	286	47.9		158	74	46.8
Senior	86	36	41.9		23	10	43.5
Night shift (times per month)				0.007				0.651
No	556	208	37.4		137	55	40.1
≥ 2	151	68	45.0		46	20	43.5
≥ 4	137	60	43.8		37	19	51.4
≥ 6	724	340	47.0		200	89	44.5
Work hours (hours per week)				< 0.001				0.247
< 40	412	141	34.2		126	46	36.5
40–49	760	335	44.1		180	81	45
50–59	226	106	46.9		71	35	49.3
≥ 60	170	94	55.3		43	21	48.8
WOF score				< 0.001				< 0.001
< 25	280	119	42.5		75	32	42.7
25–29	394	201	51.0		111	56	50.5
30–34	298	179	60.1		83	51	61.4
35–39	337	125	37.1		95	39	41.1
≥ 40	259	52	20.1		56	5	8.9

WOF worker-occupation fit

The number indicates the number of individuals with occupational stress; the % indicates the incidence rate of occupational stress. *P value was analysed by the chi-square test, with significance defined at < 0.05

Risk model and nomogram

In the Table 3, multiple logistic regression was used to construct the risk model of occupational stress. Six variables were included in the final risk model for occupational stress (Table 3). A high WOF was associated with a decreased risk of occupational stress (P < 0.01). The variables were male sex (OR = 2.10, 95% CI 1.431–3.091), nurse population (OR = 6.04, 95% CI 4.043–9.031), lack of exercise (OR = 2.02, 95% CI 1.594–2.549), traumatic event experience (OR = 2.10, 95% CI 1.623–2.730), work experience duration of 1–3 years (OR = 1.71, 95% CI 0.958 − 3.032), work experience duration of 4–10 years (OR = 2.97, 95% CI 1.707–4.892), work experience duration of 11–19 years (OR = 3.61, 95% CI 2.104–6.322) and work experience duration of more than 20 years (OR = 3.45, 95% CI 1.932–6.732), which were all associated with an increased risk of occupational stress (all P < 0.05).

Table 3

The risk model for occupational stress by multiple logistic model in training set

Variable	B	OR	95% CI		P	SE	Z
			Lower	Upper
WOF (Ref: < 25)
25–29	– 0.660	0.52	0.352	0.759	0.001	0.10	– 3.37
30–24	– 1.077	0.34	0.205	0.567	< 0.001	0.09	– 4.14
35–39	– 1.872	0.15	0.093	0.253	< 0.001	0.04	– 7.36
≥ 40	– 2.760	0.06	0.037	0.107	< 0.001	0.02	– 10.24
Sex (Ref: Female)
Male	0.743	2.10	1.431	3.091	< 0.001	0.41	3.78
Occupational categories (Ref: Physician)
Nurse	1.799	6.04	4.043	9.031	< 0.001	1.24	8.77
Exercise (Ref: Yes)
No	0.744	2.02	1.594	2.549	< 0.001	0.24	5.86
Traumatic events (Ref: No)
Yes	0.701	2.10	1.623	2.730	< 0.001	0.06	– 5.61
Work experience (Ref: < 1 year)
1–3 years	0.566	1.71	0.958	3.032	0.030	0.46	2.17
4–10 years	0.819	2.97	1.707	4.892	0.001	0.54	3.43
11–19 years	0.973	3.61	2.104	6.322	< 0.001	0.65	3.93
≥ 20 years	1.104	3.45	1.932	6.732	< 0.001	0.81	4.13

OR odds ratio, 95% CI 95% confidence interval, SE standard error, WOF worker-occupation fit

Table 4 shows the risk features for occupational stress in the risk model. Based on the risk model, the subgroup with the highest probability of occupational stress would be the high-risk features. The features included: (1) a WOF score less than 25 points, which was associated with an expected risk probability of occupational stress of 77.3%; (2) nursing as an occupation, which was associated with an expected risk probability of occupational stress of 63.1%, (3) male sex, which was associated with an expected probability of occupational stress of 67.2%; (4) work experience duration from 11 to 19 years, which was associated with an expected risk probability of occupational stress of 54.5%, (5) the experience of traumatic events within the previous year, which was associated with an expected risk probability of occupational stress of 65.3%, and (6) no habit of regular exercise, which was associated with an expected risk probability of occupational stress of 60.2%. For medical workers with all these risk features, the expected risk probability of occupational stress would be 90.2%.

Table 4

Risk features for occupational stress in the training set

Variable	High-risk features	Expected risk prob. (%)
WOF	< 25	77.3
Sex	Male	67.2
Experience of traumatic events	Yes	65.3
Occupation category	Nurse	63.1
Exercise	No	60.2
Work experience	11–19 years	54.5
All features		90.2

WOF worker-occupation fit; expected risk prob.: incidence probability

Based on the risk model for occupational stress, a nomogram was developed (Fig. 3). The included variables were the six identified variables (sex, occupation, exercise, experience of traumatic events, work experience and WOF). Internal validation and external validation showed that the final risk model had excellent discriminatory ability, with AUCs of 0.747 (95% CI 0.723–0.772) and 0.799 (95% CI 0.756–0.841), respectively (Fig. 4).

Fig. 3

Nomogram used to predict the probability of occupational stress among medical workers. WOF worker-occupation fit. The seventh row (points) indicates the points that are assigned to each variable in rows 1–6, which are the variables that are included in the risk model. The assigned points for all variables are then summed, and the total value is shown as the total score. Once the total score is located, a vertical line is drawn down to the bottom line to obtain the predicted probability of occupational stress

Fig. 4

Receiver operating characteristic (ROC) curves showing the performance of the risk prediction model for the identification of occupational stress in the training set (AUC = 0.747) and validation set (AUC = 0.799)

Figure 5 shows the calibration curves indicating the performance of the model in the validation cohort. The model had excellent predictive performance throughout the range of predicted risks and was accurate through a range of predicted probabilities of occupational stress from 25% to approximately 70% in the training set. The calibration curve in the validation set showed a greater deviation than that in the training cohort because of the small sample size.

Fig. 5

Calibration curve for predicted versus observed risk of occupational stress in the training and validation cohorts. The risk model estimated probability is plotted on the X-axis, and the fraction corresponding to the positive probability is plotted on the Y-axis

To ensure the practical applicability of the risk model, a cut-off point of 56.4% was established based on the maximum Youden index and was used to stratify medical workers into high- and low-risk groups. The sensitivity at the cut-off point was 0.63, and the specificity at the cut-off point was 0.74.

Discussion

In our study, among the 1988 medical workers surveyed, the prevalence of occupational stress was 42.5% (n = 845), and 9.6% (n = 190) of the participants had a severe level of occupational stress. These results indicate that medical workers have substantial responsibilities and experience high levels of demands and overwork, predisposing them to have a high risk of occupational stress (Xu and Hu 2020; Jerg-Bretzke et al. 2020; Li et al. 2019).

Main results

In this study, we established a risk model and nomogram for occupational stress. We used the training set for internal validation and the validation set for external validation, and both showed good discriminative ability of the model, with AUCs of 0.747 and 0.799, respectively. The high-risk features for occupational stress were identified. A total WOF score less than 25 points was associated with an expected risk probability of individuals with occupational stress of 77.3%. The expected risk probability of nurses with occupational stress was 61.3%, which was higher than the risk of physicians. Among medical workers who had worked for 11–19 years, the risk with occupational stress was 54.5%, which was higher than the risks among groups with other levels of work experience. The traumatic events experienced was associated with an expected risk probability of occupational stress of 65.3%. The risk of medical workers with occupational stress was higher in the subgroup that did not regularly exercise (expected risk: 60.2%) than in the other groups. In a group of workers with all the identified risk features, the expected risk probability with occupational stress incidence would be 90.2%.

WOF level as a new feature of occupational stress

We found that the lower the WOF level was, the higher the risk of occupational stress in the study. WOF is an ignored risk factor for occupational stress, which is rarely mentioned in the study of occupational stress. Our results indicated that when the WOF score was less than 25 points was associated with an expected risk probability of individuals with occupational stress of 77.3%. The results suggested that a strong fit between workers and their occupational environment is a protective factor against occupational stress. The WOF score was less than 25 points can be used as an indicator to identify medical workers who are at high risk of occupational stress. Before the allocation of medical workers to certain departments, hospital administrators could consider whether the workers’ characteristics, abilities and needs are a good fit with the occupational environments’ demands and supplies. Using the strength of the fit to predict the risk of occupational stress incidence would not only protect the occupational health of each worker but also prevent a reduction in the quality of medical services provided. Furthermore, increasing the level of WOF when medical workers experience a high level of occupational stress would be a good way to resolve the problem.

Demographic and occupation features of the group at high risk for occupational stress

The demographic and occupation features included male sex, being a nurse, and work experience. Regarding to Male sex, the risk of occupational stress among male medical workers (expected risk: 67.2%) was higher than that among female medical workers. Previous studies reported that males were likely to suffer from occupational stress (Davey et al. 2014; Faraji et al. 2019). This could be due to males are likely to assume more social or family responsibilities in traditional Asian cultures, and the managers of hospitals tend to assign physical work to male medical workers. Regarding to be nurses, nurses always reported the highest level of occupational stress due to their heavy workload (Sandrin et al. 2019), reduced social support (Ma et al. 2020), staff shortages (Baye et al. 2020), and decreased amount of time to perform their work duties (Salilih and Abajobir 2014). As a result, the level of occupational stress was higher in nurses than in physicians, and the level of empowerment was lower in nurses than in physicians. Regarding to work experience, the work experience and occupation categories were contributors to the risk of occupational stress. In the risk model, medical workers with 4–9 years of work experience had a higher risk probability of occupational stress than the other groups. This group of medical workers, namely, primary medical workers and intermediate medical workers, is eager to seek promotions or increase their monthly income (Ta et al. 2018). Medical workers who have worked more than 20 years, however, have a reduced physical ability to perform healthcare services, and work fatigue may lead to occupational stress.

Health behaviours of the group at high risk for occupational stress

Exercising, drinking and smoking habits reflect whether people engage in healthy behaviours. This study showed that the risk of occupational stress incidence among medical workers who did not regular exercise (expected risk: 60.2%) and therefore had an unhealthy lifestyle was higher than that among medical workers who regular exercised. There is evidence that exercise is beneficial for mental health and can reduce the risk of mental illness (Rosenbaum et al. 2014). Chou et al. suggested that physical inactivity is the only significant factor correlated with occupational stress, and regular exercise could help medical workers resolve occupational stress (Chou et al. 2016). Soteriades and colleagues showed that every hour per week of physical activity performed by participants was associated with a 16% lower risk of occupational stress (Soteriades et al. 2019). This conclusion support for our results that obtained in this study. Exercise may be an effective approach to coping with occupational stress (Callaghan 2004). Consuelo Arbona et al. found that drinking was a strong risk factor for occupational stress in black firefighters. But in our study, drinking and smoking were not significantly associated with the risk of occupational stress. The reasons may include the following: (1) Firefighters are more often male, while physicians and nurses are more often female in China, and the proportion of those who smoke and drink is relatively small. (2) Medical workers are awareness of the harmful effects of smoking and drinking, so they consciously avoid such unhealthy lifestyle behaviours.

The results also indicated that exposure to traumatic event as a risk feature of occupational stress. Exposure to traumatic events in the workplace is common, and traumatic events may result in psychosomatic disorders (Lee et al. 2020). In our study, the results showed that the medical workers who had experienced traumatic events the expected risk of developed occupational stress was 65.3%. Studies on traumatic events have more often focused on firefighters (LaRocca et al. 2020). In South Korea, a nationwide, population-based survey of all employed firefighters found that the rate of PTSD was estimated to be 5.4% (Kim et al. 2018). Additionally, a cross-sectional survey in Australia indicated that the rates of PTSD and depression in current and retired firefighters were 8% and 5%, respectively (Jahnke et al. 2016). Both firefighters and medical workers have obvious occupational characteristics and are professionals at high risk for occupational stress. Firefighters and medical workers both experience higher demands at work than other professionals, and they may have experienced stressors such as continuously being at risk of dying. These results and the findings of studies on firefighters indicate that it would be beneficial for managers to look for signs that a nurse or physician may be experiencing occupational stress. We suggest that hospital managers assess occupational stress in workers who have experienced a traumatic event on a regular basis. In addition, psychological training to improve situational awareness and professional adaptability is needed for medical workers, as are interventions to help them overcome psychological crises or negative emotions.

Application of the nomogram for occupational stress

Nomograms have been widely used in clinical practice, especially for cancer, and have been found to have excellent predictive performance (Jiang et al. 2017; Kazushige et al. 2015). Currently, the measurement of occupational stress is based on a questionnaire survey that estimates the level of stress in an occupational group. Conducting such a survey requires the commitment of time, labour and material resources by professional survey administrators or researchers. The proposed nomogram for occupational stress is convenient for use by managers in hospitals; alternatively, it could be used for self-assessment. Our results showed that the cut-off point for the identification of a high risk of occupational stress is the expected risk was 56.4% by risk model. In a future study, we can explore the risk probability stratified by different characteristics in other professions and collect data on medical workers’ risk probability using the nomogram, with the probability considered a variable related to occupational stress.

Study limitations

Some limitations of this study should be noted. The cross-sectional design of this study limits its ability to confirm causal relationships among the variables that were investigated. Thus, the present study results should be interpreted cautiously, and future research using longitudinal and experimental designs is needed to clarify the direction of any causal relationships among the observed variables. Despite the limitations of our study, these findings contribute significantly to application-oriented research and provide practical support for the management and prevention of occupational stress in medical workers.

Conclusions

In this study, 42.5% (845/1988) of the medical workers experienced occupational stress. Having a poor WOF (WOFI score < 25), being male, being a nurse, having more work experience (11–19 years), experiencing traumatic events and without regular exercise habit were risk features for occupational stress. The identification of risk features and the development of a nomogram are valuable for enhancing application-oriented research and exploring ways to provide practice support to enable the targeted prevention of occupational stress in hospitals. Moreover, self-assessment with nomograms can also improve awareness of occupational stress in those performing routine work duties.

Below is the link to the electronic supplementary material.

Supplementary file 1 DOCX 15 KB)