Occupational burnout, flourishing and job satisfaction among HIV/AIDS healthcare workers in Western China: a network analysis

Background Healthcare workers suffered with high prevalence of occupational burnout, which might be related with their job satisfaction and well-being. This study aimed to provide evidence of complex interrelations among occupational burnout, flourishing, and job satisfaction, and identify key variables from the perspective of network structure among healthcare workers. Methods A cross-sectional study was conducted between July and October 2021, and 907 (the response rate was 98.4%) HIV/AIDS healthcare workers completed their sociodemographic characteristics, occupational burnout, flourishing and job satisfaction. Network analysis was conducted to investigate the interrelations of occupational burnout, flourishing, and job satisfaction communities, and identify central variables and bridges connecting different communities with different bridge strength thresholds in the network structure. The Network Comparison Test (NCT) was conducted to examine the gender differences in networks. Results In the network, feeling exhausted at work (strength: 1.42) and feeling frustrated at work (1.27) in occupational burnout community, and interested in daily activities (1.32) in flourishing community were central variables. Bridges in the network were job reward satisfaction (bridge strength: 0.31), satisfaction with job itself (0.25), and job environment satisfaction (0.19) in job satisfaction community, as well as interested in daily activities (0.29) and feeling respectable (0.18) in flourishing community, with bridges selected with top 20% bridge strengths. Feeling frustrated at work (0.14) in occupational burnout community and leading a purposeful and meaningful life (0.11) in flourishing community became bridges when using thresholds of top 25% and 30% bridge strengths, respectively. We also observed higher network densities in females (network density: 0.37) than that in males (0.34), and gender differences in the distribution of partial correlation coefficients (M = 0.27, P = 0.017). Conclusions In the network structure of occupational burnout-flourishing-job satisfaction, feeling frustrated at work in occupational burnout community and interested in daily activities in flourishing community were both central variables and bridges, which may be targeted variables to intervene to alleviate the overall level of symptoms in the network and therefore prevent poor health outcomes in healthcare workers. Supplementary Information The online version contains supplementary material available at 10.1186/s12888-023-04959-7.


Background
Healthcare workers, who are a significant part of the health system, are well-known to experience excessive workloads [1,2].Overloaded work, high work-related demand and unsupportive environments can induce healthcare workers' occupational burnout [3].Rotenstein and colleagues estimated a pooled prevalence of 21.3% on occupational burnout in practicing physicians [4].In Chinese healthcare workers, related studies also showed a high prevalence of occupational burnout ranging from 75.5% to 80.0% [5][6][7], and the prevalence in HIV/AIDS healthcare workers has reached 76.9% [8].
The ratio between healthcare workers and patients was imbalanced because of a serious shortage of HIV/AIDS healthcare workers in China [8].In addition, HIV/AIDS healthcare workers usually face occupational exposure risks for HIV infection, which may aggravate their occupational burnout [9].Occupational burnout is one of the critical reasons for poor quality of medical services, medical accidents, and low efficacy of job performance [10][11][12][13][14][15].It is urgent to understand occupational burnout and its influencing factors especially among HIV/AIDS healthcare workers.
Occupational burnout may be caused by negative psychological factors, such as anxiety, depression, and poor aptitude for interpersonal relationships [16,17].In contrast, positive psychological factors are beneficial in reducing occupational burnout [18].Flourishing, a state in which all aspects of one's experience are well, could reflect overall well-being comprehensively and has more advantages in explaining complex relations and connotations compared with other tools [19,20].Employees with flourishing would exhibit positive outcomes, such as higher job satisfaction [21][22][23][24][25].To date, evidence on flourishing and the association between flourishing and occupation burnout in HIV/AIDS healthcare workers is scarce.
Job satisfaction, defined as people's feelings about their jobs [26], is another factor that would greatly increase enthusiasm and reduce the absenteeism of workers, even in those with high occupational burnout [27][28][29][30][31][32].Job satisfaction was related to positive psychological factors [33,34].As described in the spillover model, higher job satisfaction leads to higher well-being [35].However, current studies estimated the relationships between these positive and negative psychological factors (i.e., occupational burnout, flourishing, and job satisfaction) based on the sum-scores of scales rather than individual items, and the subtle interrelationships among the dimensions of job satisfaction, flourishing and occupational burnout among HIV/AIDS healthcare workers remains to be elucidated clearly.
Network analysis has increasingly been used in psychological research to quantify the relationships among multi-dimensional psychological factors and identify the most interconnected dimensions at a time.Meanwhile, it is plausible that a group of nodes (variables) belong to a community (i.e., the same rating scale) where they are closely related, and network models can help identify central variables in these closely related variables [36].Moreover, network models describe an interacting web of symptoms and could help to understand the transdiagnostic symptom relationships that underlie the well-established association of occupational burnout, flourishing and job satisfaction, and can be a useful method to examine the bridges that are represented by one or more key variables bridging two or more communities [37].Deactivating bridges may constitute an effective strategy to block connections between items efficiently.For example, a previous study identified bridge symptoms (i.e., "irritability", "feeling afraid" and "sad mood") that could be targeted in specific treatment and preventive measures for comorbid depressive and anxiety symptoms [38].However, limited research, to our knowledge, used this method in occupational populations who have a high probability of occupational burnout and other negative emotions.
In light of the above gaps, we used network analysis to (1) estimate the associations of occupational burnout, flourishing, and job satisfaction in a constructed network ; (2) identify central variables in network structure of occupational burnout-flourishing-job satisfaction and (3) seek bridges linking occupational burnout, flourishing and job satisfaction.Previous studies found that there were differences between males and females in the prevalence of occupational burnout and mental health, and attitudes toward job satisfaction [39,40].We therefore test the differences in the network structure of occupational burnout-flourishing-job satisfaction in males and females.The findings of this study would help obtain a straightforward view of variables' interactions to clarify critical issues and guide more effective interventions to improve the health and professional happiness of healthcare workers in China.More broadly, our findings would help policy-making for improving the psychological health of other healthcare workers.

Study design and participants
This study is an observational and cross-sectional design in the Sichuan province of China, where is one of the most HIV-affected province in China [41,42].In such a scenario, HIV/AIDS healthcare workers in Sichuan province may undertake overloaded work and high workrelated demands.The current study used a stratified cluster sampling method to obtain a sample.First, we randomly selected three cities based on the prevalence of HIV as our study sites.Second, all related healthcare workers in selected cities were recruited for investigation (e.g., from the Centers for Disease Control and Prevention (CDC) and Designated Hospitals for HIV treatment) (Figure S1).Finally, 922 HIV/AIDS healthcare workers were enrolled to complete the electronic survey.Inclusion criteria of the study included at least 18 years old at the time of entry to this study, medical and management personnel involved in HIV/AIDS prevention and control work.Exclusion criteria included local HIV/AIDS health workers who were illiterate, and those who refused to participate in this survey.All the participants gave informed assent forms before data collection; it has been stated on the cover page that the participants cannot proceed with the online survey without consent.Ethics approval was obtained from the Ethics Committee West China School of Public Health and West China Fourth Hospital, Sichuan University (Gwll2021059).
The sample size was calculated based on n = . π = 76.9% (76.9% of HIV/AIDS healthcare workers met the criteria for burnout in a province of China) [8], δ = 0.05, α = 0.05, µ α/2 = 1.96, and the sample size was calculated to be 273.However, due to the large sampling error of the cluster, the sample size was usually 1.5 times of the sample size of the pure random sample, so the sample size of participants should be at least 409.

Data collection
Data for the analysis was collected between July 2021and October 2021 through an online questionnaire designed by an expert panel.The panel consisted of two epidemiologists, one health psychologist, and two HIV/AIDS healthcare workers.A pretest was conducted in 10 HIV/ AIDS healthcare workers to test the questionnaire's comprehension and the feasibility of investigation process.Their feedback was used to revise and finalize the questionnaire.Besides, the online questionnaire with attached instructions for completion was released by the HIV/ AIDS management department in the Health Commission of the three cities through WeChat and other social media platforms and forwarded to all the HIV/AIDS healthcare workers.The questionnaire took about 20 min to complete on average.

Occupational burnout scale
The widely used scale to measure occupational burnout is the version of the Maslach Burnout Inventory-General Survey (MBI-GS) by Maslach and Jackson [43], which has good reliability and validity in China [44].Our study used a self-designed scale developed based on the MBI-GS to better adapt the local culture of HIV/AIDS healthcare workers, and to ensure the high-quality results, that is, within the endurance limits to be interviewed.The scale used 11 items to assess the occurrence of occupational burnout, including emotional exhaustion (3 items), work fatigue (1 item), emotional numbness (3 items), work frustration (2 items), and work pressure (2 items), and it was a five-point scale ranging from strongly disagree (1-point) to strongly agree (5-point).A sum score ranging from 11 to 55 was calculated, and higher scores indicated more severe occupational burnout (Table S1).The items in occupational burnout scale specifically described the individual's experiences and feelings about burnout.By using these items as analyzed variables, we could capture individual responses to different aspects of burnout, and further understand the network structure and characteristics of burnout.The occupational burnout scale demonstrated high internal consistency, with Cronbach's alpha of 0.895 in this study.The KMO coefficient was 0.878, and Bartlett's sphericity test was significant (P < 0.05), indicating good construct validity for the scale.

Flourishing scale
The flourishing scale (FS) was used to assess the flourishing of HIV/AIDS healthcare workers.The flourishing scale consists of 8 items describing essential aspects of human functioning and human needs, including the meaning of life, interpersonal relationship, life and work engagement, helpful to others, competence, self-esteem, optimism and respect from others [45].Each item was answered on a 7-item scale ranging from strong disagreement (1-point) to strong agreement (7-point).The total score of all items ranged from 8 (strong disagreement) to 56 (strong agreement).A higher total score indicated a greater sense of thriving, which showed a person had a positive attitude toward life and owed many psychological resources and advantages (Table S2).The Chinese version of flourishing scale has good validity and reliability [46,47].The flourishing scale had high internal consistency with Cronbach's alpha of 0.932 in this study.The KMO coefficient was 0.919, and Bartlett's sphericity test was significant (P < 0.05), indicating good construct validity for the scale.

Job satisfaction scale
Job satisfaction was measured by a self-developed scale and was revised and finalized after the pilot study.Three dimensions of job satisfaction, including satisfaction with job itself (10 items), job environment satisfaction (6 items) and job reward satisfaction (8 items), were recorded.All items had five degrees ranging from disagree (1-point) to agree (5-point).The total score of each dimension was calculated, and a higher score indicated higher job satisfaction in the corresponding dimension (Table S3).The job satisfaction scale consisted of 24 items of three dimensions.Since more variables may increase the difficulty of interpretation and lead to over-fitting and very unstable estimates in the network analysis [48], we therefore used dimensions to estimate the individual's job satisfaction of different aspects, increasing the visualization and interpretability of the results.The Cronbach's alpha coefficients for the three dimensions of job satisfaction were 0.853, 0.834 and 0.808 in this study, respectively.The KMO coefficient for the three dimensions were 0.858, 0.841 and 0.812, respectively, and Bartlett's sphericity test was significant for each dimension of the scale (P < 0.05).

Socio-demographic and work-related characteristics
The collected socio-demographic characteristics included age, sex, marital status, educational level, personal monthly income and living situation.Besides, work-related factor included types of institutions.All relevant descriptions and explanations were detailed in Table S4.

Statistical analysis Descriptive statistics
Data were summarized by descriptive statistics, with frequency and percentage for categorical variables, and median and interquartile range (IQR) for continuous variables (Table S5).R version 4.0.3 for Windows was used for data management and all statistical analyses, with P value < 0.05 indicating statistical significance.

Scale-construction process (construct validity and reliability)
For the self-designed scales, psychometric analyses were performed in several steps: the first step consisted of removing items with a rate of missing values (missing value and "I don't know" choice) > 20%, or a floor or ceiling effect > 50% [49], and no item was removed.Second, Kaiser-Meyer-Olkin (KMO) coefficients and Cronbach's alpha coefficients were computed for evaluating the validity and internal consistency of scales, respectively.Third, an exploratory principal component analysis (PCA) using a varimax rotation for the remaining items of scales such as the job satisfaction scale was used in the study [49].The number of dimensions was determined using a scree plot, and factor loading more than 0.60 was attributed to as a dimension [50].

Network estimation
We assessed the associations among occupational burnout, flourishing and job satisfaction based on the network of Graphical Gaussian Model (GGM).Partial correlation analysis was conducted to indicate the association of each pairwise variable and form networks; nodes in the network represented variables, and edges represented partial correlation coefficients between two variables.Stronger correlations were shown in thicker and more saturated edges.Positive and negative correlations were shown in green and red, respectively.In this network estimation, an Extended Bayesian Information Criterion (EBIC) model with the least absolute shrinkage and selection operator (LASSO) were used to get a sparse and intelligible network [51,52].We also estimated the predictability, the upper bound of variance (measured in R 2 ) of a given variable explained by all the other variables in the network [53].

Centrality and bridge symptoms
For the constructed networks, we calculated strength centrality (i.e., the sum of the absolute value of all partial correlation coefficients for a given variable) to identify the most central variables, using the centralityplot function in the qgraph package in R [54].In addition, we conducted a stability and reliability analysis of our results using a commonly applied bootstrapping procedure in R with bootnet package, which showed whether the networks remained stable when dropping 75% of the sample, and whether the results of bootstrap 95% confidence interval (CI) or edges were narrow, indicating the trustworthy of edges [55].
We used the term community to designate a group of items or dimensions, which were supposed to be related according to scale classification, independent of the actual network structure [36].We also estimated variables that acted as bridges connecting communities, which allowed us to identify which items were the most interconnected across occupational burnout, flourishing and job satisfaction communities.The bridges were estimated with different scoring variable thresholds (i.e., variables that are most strongly connected to all the variables of different communities) using the bridge function in the networktools package in R [56].

Network comparison test
Network Comparison Test (NCT) is a permutation-based hypothesis test for invariance of network structure (i.e., how the connections between variables within a network differ across samples) and global strength of connections (i.e., how the density of the network differs across samples -the sum of all edge strengths).To examine sex differences in the structure of networks, we ran the NCT using the Network Comparison Test package in R [57].We compared the distribution of partial correlation coefficients in each network to characterize the network structure.Then, we compared differences in strength for each edge of networks between males and females, after controlling for multiple tests using a Holm-Bonferroni correction.Additionally, the network densities (the actual number of connections between variables in the network / the theoretical number of connections between variables in the network [58]) between males and females were also compared.

Sociodemographic and work characteristics
A total of 907 (the response rate was 98.4%) HIV/AIDS healthcare workers were included in the final analysis, with a mean age of 38.0 ± 9.4 years.Among them, 69.1% (n = 627) were females, 49.5% (n = 449) graduated from junior colleges and 27.9% (n = 253) had a personal income below 3,000 Renminbi (RMB) per month.Further details were provided in Table 1.

Network structure and centrality measures
The network structure showed a high degree of interrelations between variables from occupational burnout, flourishing and job satisfaction communities, and the variables in the same communities tended to cluster together.Most of the associations between variables of different communities were negative (Fig. 1).The predictabilities of variables were shown as ring shaped pie charts, and variables with the highest centrality were E2 (feeling exhausted at work in occupational burnout community; strength: 1.42), H3 (interested in daily activities in flourishing community; strength: 1.32) and E10 (feeling frustrated at work in occupational burnout community; strength: 1.27).We also found the stability of the network remained stable (i.e., case-dropping coefficient = 0.75) even dropping large proportions of the sample (Fig. 2).The numerical interactions were shown in Table S6 using weighted adjacency matrix that represented the strength of associations between variables.The results of the reliability and stability analyses were presented in Figures S2,  S3 and S4.

Bridges and bridge centrality measures
Figure 3 showed the bridges and bridge centrality indices.The top 20% scoring variables on bridge strength were W1 (satisfaction with job itself; bridge strength: 0.25), W2 (job environment satisfaction; bridge strength: 0.19) and W3 (job reward satisfaction; bridge strength: 0.31) in job satisfaction community; H3 (interested in daily activities; bridge strength: 0.29) and H8 (feeling respectable; bridge strength: 0.18) in flourishing community (Fig. 3A).When identifying bridges with the top 25% bridge strengths, we observed that E10 (feeling frustrated at work; bridge strength: 0.14) in the occupational burnout community was additionally added as a bridge linking other communities (Fig. 3B).Furthermore, H1 (leading a purposeful and meaningful life of flourishing community; bridge strength: 0.11) was added as a bridge when using the top 30% bridge strengths (Fig. 3C).Bridge strength was reported in Fig. 3D.

Network comparisons between males and females
The network structures and network centrality indices of males (n = 607) and females (n = 280) were showed in Fig. 4, and the network structures were both stable (Figures S5, S6, S7 and S8).The number of edges in female network were 86, and in male network were 79.The network density of males and females were 0.37 and 0.34, respectively.The results of the NCT indicated significant sex differences in the distribution of partial correlation coefficients (M = 0.27, P = 0.017), and no significant  S9).We also calculated bridge centrality in males and females' networks, separately (Fig. 5).No bridge was found with the top 20% and 25% bridge strengths, but seven bridges were found in both males and females' networks with the top 30% bridge strength.However, there was a significant sex difference.For example, E10 (feeling frustrated at work of occupational burnout community; bridge strength: 0.14) was a bridge in the network of females, but no bridge was found in occupational burnout community in the network of males.Besides, H2 (having supportive and rewarding social relationships; bridge strength: 0.12) was a bridge in flourishing community in females' network, however, H1 (leading a purposeful and meaningful life; bridge strength: 0.09) and H4 (contributing to the happiness and well-being of others; bridge strength: 0.08) were bridges in males' network.

Discussion
This study estimated the network structure of occupational burnout, flourishing and job satisfaction, and firstly attempted to assess their interactions at the item level in HIV/AIDS healthcare workers.Our results indicated that variables feeling exhausted at work and feeling frustrated at work in occupational burnout community, and interested in daily activities in flourishing community were the most central variables.Job reward satisfaction, satisfaction with job itself and job environment satisfaction in job satisfaction community, as well as interested in daily activities and feeling respectable in flourishing community were bridges when identifying bridges with the top 20% bridge strengths.Feeling frustrated at work in the occupational burnout community and leading a purposeful and meaningful life in the flourishing community became bridges connecting other communities with the top 25% and 30% bridge strengths, respectively.Gender was associated with different distributions of partial correlation coefficients between their networks, and network densities in females were higher than that in males.
Central variables may give insights into the connectedness or importance of items within a network and may play a major role in causing the onset of and/or maintaining a syndrome [59].Bridges mediate the transition among different syndromes and may increase the risk of contagion to other disorders [60].In this study, feeling exhausted at work was the most central variable with high centrality in the network, followed by interested in daily activities and feeling frustrated at work.Meanwhile, we found that interested in daily activities and feeling frustrated at work were also bridges connecting other communities.Thus, there may have some overlaps between central variables and bridges, and these overlaps deserve special attention.A potential explanation is that interested in daily activities and feeling frustrated at work are two key mental statuses of coexistence that may affect other mental statuses.For HIV/AIDS healthcare workers in charge of disease prevention and control and patients' treatment, they are overburdened with heavy work demands, and the mismatch between effort and income may make them frustrated [61].Besides, previous studies showed that some healthcare workers in special departments experienced more psychological pressure, that is, much of the pressure mainly comes from daily activities [62,63].As Sichuan Province is one of the provinces with high HIV epidemic in China, the more attention AIDS prevention and treatment receives, the greater the responsibility HIV/AIDS healthcare workers assume.In addition, if the occupational population is interested in their daily work, they may have high job satisfaction and be less prone to occupational burnout since they would think the work being an emotional reward, and may be more motivated to complete [64].These central variables and bridges not only played an important role in understanding the structure of the network model, but also can serve as explicit intervention points to relieve all symptoms.It should be noted that bridges changed in different percentile cut-offs of bridge strength (80%, 75%, and 70%), which were different from Jones's study that only focused on one cut-off [36].The number of bridges would increase when the percentile cut-off point decreased, and therefore we identified at least one bridge from each community, which may guide us to carry out targeted interventions in the bridges corresponding to each community.In this study, the items belonging to the same communities tended to cluster together, and the relationships between occupational burnout, flourishing and job satisfaction were more subtle compared to previous studies' findings that only considered the total score or subscale total scores of the rating scales [21,65,66].Meanwhile, our study revealed that there were intrinsic interactions among variables of different communities.For example, leading a purposeful and meaningful life in flourishing community and feeling frustrated at work in occupational burnout community were related, which was consistent with a previous study [67].That is, personal value and a sense of self-fulfilled could help individuals find meaning and value in their work by increasing their satisfaction and sense of accomplishment at work, thereby reducing occupational burnout.Therefore, network analysis offers an additional tool to capture connecting patterns of occupational burnout, flourishing and job satisfaction, and intrinsic connections among variables of different communities suggesting that interventions targeting crucial variables in each community may be beneficial to the whole network.
Furthermore, our analysis documented both similarities and differences in the network structure of occupational burnout-flourishing-job satisfaction in males and females.The central variables were almost similar in both males and females.However, the bridges had a large difference in occupational burnout and flourishing community, which may be caused by the career expectations and attitudes towards problems between males and females The magnitude of the edge depicts the magnitude of the correlation, with thicker and saturated edges showing stronger correlations.Positive and negative correlations were shown in green and red, respectively.Variables identified as job satisfaction are colored in red circles; Variables identified as flourishing are colored in orange circles; Variables identified as occupational burnout are colored in green circles; Variables identified as bridges are drawn with rectangles.The top 30% scoring variables on a given statistic and selected as predicted bridges in the males and females [68].We also examined sex differences in the structure of networks and found network densities between males and females were different.Network density describes the portion of the potential connections in a network that are actual connections and refers to the mean strength of the connections between variables in networks [58].In the study, we found higher network densities in females than that in males, which showed that symptoms may transmit more quickly in females.The results were similar to previous studies which have found that males and females responded to conditions in different ways and females tend to develop internal symptoms [39,69].The differences may be derived from particular psychological characteristics, unique perspectives and experiences on work-related aspects in females which were different from males.
This study has several strengths.First, the network model was designed to consider individual symptoms as variables in the study.In addition, we compared bridges under different thresholds, which can help to find more important bridges in different communities and help design subtle interventions.Finally, NCT analysis explicitly tested gender as a sociodemographic influence on the network structure of occupational burnout-flourishing-job satisfaction of health workers, and indicated the heterogeneity in central and bridges for targeted intervention in two networks.
Several limitations should also be mentioned.First, the cross-sectional data is limited to infer causal relations.To fully conceptualize how the variables interact with one another, future experimental and longitudinal studies are needed.However, related research has shown that network structures among cross-sectional and longitudinal studies did not differ [70].Second, as the study relied on the self-report online questionnaire, the results may not rule out other possible factors, social desirability or recall biases, etc.There may be also understanding deviations of different respondents and the reliability and accuracy of the survey may be affected to some extent.Third, we only examined the difference in network structure between gender.However, other sociodemographic factors are also needed to examine in future studies.Finally, generalizations should be taken with caution since the sample comprised only a small proportion of healthcare workers in China.Therefore, researches on geographical and cultural diversity are warranted in the future.

Conclusions
Through network analysis, we estimated the complex association between occupational burnout, flourishing and job satisfaction among Chinese healthcare workers.Feeling frustrated at work of occupational burnout community and interested in daily activities of flourishing community were both the central variables and bridges.Specific strategies and measures targeting to these variables should be given such as improving levels of flourishing and job satisfaction and reducing occupational burnout.Additionally, government and healthcare agencies should address occupational burnout and improve flourishing, especially take initiatives to improve job satisfaction with effective interventions in healthcare workers.
• fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year

•
At BMC, research is always in progress.

Learn more biomedcentral.com/submissions
Ready to submit your research Ready to submit your research ?Choose BMC and benefit from: ? Choose BMC and benefit from:

Fig. 1
Fig. 1 Networks and strength centrality displaying the relationship between occupational burnout, flourishing and job satisfaction.A The regularized network between occupational burnout, flourishing and job satisfaction.B Standardized variable strength centrality for the network (z-scores).Variables are represented by nodes, numbers represent item numbers in the scale, partial correlations between the variables are represented by edges.The color of the edge indicates the direction of the correlation (red = negative, green = positive).The magnitude of the edge depicts the magnitude of the correlation, with thicker and saturated edges showing stronger correlations.Positive and negative correlations were shown in green and red, respectively.Variables identified as job satisfaction are colored in red circles; Variables identified as flourishing are colored in orange circles; Variables identified as occupational burnout are colored in green circles.The area in the rings around the variables indicates predictability (the upper bound of the variance of a given variable explained by the remaining variables in the network)

Fig. 2 Fig. 3
Fig. 2 Stability of centrality indices by case dropping bootstrap for the network shown in Fig. 1 (The bootstrap was repeated 1000 times).The x-axis represents the percentage of cases of original sample used at each step.The y-axis represents the average of correlations between the centrality indices from the original network and the centrality indices from the networks that were re-estimated after dropping increasing percentages of cases.Each line indicates the correlations of strength, while areas indicate 95%CI

Fig. 4
Fig. 4 Networks and strength centrality displaying the relationship between occupational burnout, flourishing and job satisfaction in males and females.A Graphical LASSO model of variables with occupational burnout, flourishing and job satisfaction in females (n = 627).B Graphical LASSO model of variables with occupational burnout, flourishing and job satisfaction in males (n = 280).C Comparison of network centrality indices between males and females.Variables are represented by nodes, numbers represent item numbers in the scale, partial correlations between the variables are represented by edges.The color of the edge indicates the direction of the correlation (red = negative, green = positive).The magnitude of the edge depicts the magnitude of the correlation, with thicker and saturated edges showing stronger correlations.Positive and negative correlations were shown in green and red, respectively.Variables identified as job satisfaction are colored in red circles; Variables identified as flourishing are colored in orange circles; Variables identified as occupational burnout are colored in green circles

Fig. 5
Fig. 5 Networks and bridge strength centrality displaying the relationship between occupational burnout, flourishing and job satisfaction in males and females.A Graphical LASSO model of variables with occupational burnout, flourishing and job satisfaction in females of the sample (n = 687).B Graphical LASSO model of variables with job satisfaction, flourishing and occupational burnout in males of the sample (n = 280).C Bridge centrality estimates for each variable in the network of females, ordered by highest value.D Bridge centrality estimates for each variable in the network of males, ordered by highest value.Variables are represented by nodes, numbers represent item numbers in the scale, and partial correlations between the variables are represented by edges.The color of the edge indicates the direction of the correlation (red = negative, green = positive).The magnitude of the edge depicts the magnitude of the correlation, with thicker and saturated edges showing stronger correlations.Positive and negative correlations were shown in green and red, respectively.Variables identified as job satisfaction are colored in red circles; Variables identified as flourishing are colored in orange circles; Variables identified as occupational burnout are colored in green circles; Variables identified as bridges are drawn with rectangles.The top 30% scoring variables on a given statistic and selected as predicted bridges in the males and females

Table 1
Baseline characteristics of the participantsAbbreviations: CDC Centers for Disease Control and Prevention; RMB renminbi a Others is a self-selected designation that indicates the institution is not listed