Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Psychol., 06 November 2025

Sec. Organizational Psychology

Volume 16 - 2025 | https://doi.org/10.3389/fpsyg.2025.1670195

Predicting teacher turnover in private universities: a machine learning approach based on 10 years of data and satisfaction factors

Wang Jingwen,
Wang Jingwen1,2*Liu YiLiu Yi1Yang XiaohongYang Xiaohong2
  • 1General Office of the Party and Government, Xi'an Fanyi University, Xi'an, China
  • 2College of Educational Technology, Northwest Normal University, Lanzhou, China

Background: Teacher turnover poses a significant challenge to the sustainable development of private universities in China. While machine learning (ML) has been increasingly applied to turnover prediction, existing studies often overlook psychological factors and lack longitudinal analysis.

Methods: This study integrates a 10-year longitudinal dataset with satisfaction surveys from a private university in Western China. Exploratory Factor Analysis (EFA) was employed to extract key dimensions influencing turnover. Three ML models—K-Nearest Neighbors (KNN), Naive Bayes (NB), and Backpropagation Neural Network (BPNN)—were constructed and evaluated using accuracy, F1-score, and AUC.

Results: The KNN model achieved the highest predictive performance (accuracy = 83.64%, F1 = 84.16%, AUC = 0.901). The “Compensation, Benefits, and Development” dimension was identified as the most influential factor, accounting for 25.41% of the variance.

Conclusion: This study proposes an “EFA + ML” hybrid approach that enhances feature interpretability and prediction robustness, offering practical insights for human resource management in private higher education institutions.

1 Introduction

With the continuous development of China’s higher education system, private universities have become an indispensable part of this landscape, playing a crucial role in promoting educational equity and diversity. According to the “2023 National Education Development Statistical Bulletin” released by the Ministry of Education, there are 789 private universities in China, accounting for 25.67% of the total number of universities. Private universities also enroll 26.34% of the total student population across ordinary, vocational, and junior colleges, and employ 20.13% of all full-time teachers in Chinese universities (Ministry of Education of the People's Republic of China, 2024).

Despite substantial growth in both number and student enrollment, private universities continue to experience difficulties in recruiting and retaining high-quality faculty. This challenge has emerged as a key constraint on the sustainable development of these institutions. Owing to their “non-enterprise” and “extra-system” status, faculty in private universities often experience systemic inequities compared to their public-university counterparts in areas such as political standing, social recognition, and welfare benefits. These disparities not only undermine job satisfaction and loyalty among teachers, but are also associated with higher turnover rates. Studies show that 33.5% of teachers in private universities exhibit turnover intention (Qi, 2023), with actual annual turnover rates reaching as high as 10% in some cases (Bai, 2018). High faculty attrition increases recruitment and training costs, and adversely affects teaching quality, academic development, and long-term institutional planning.

In response to the issue of teacher turnover, researchers in China and abroad have adopted a variety of approaches, including both traditional statistical methods and machine learning techniques. Conventional methods such as descriptive statistics, regression analysis, factor analysis, and mediation analysis have played significant roles in predicting employee turnover. However, their applicability and predictive power are often limited by strict assumptions, such as normality and linearity. In contrast, machine learning methods—such as decision trees, random forests, K-Nearest Neighbors (KNN), and Naive Bayes (NB), Backpropagation Neural Network (BPNN)—can more effectively capture non-linear relationships and complex decision boundaries, thereby enhancing predictive accuracy and efficiency. Nonetheless, existing machine learning studies on teacher turnover have primarily focused on corporate employees and relied heavily on objective data such as salary levels, while paying insufficient attention to psychological traits such as job satisfaction and professional efficacy. Moreover, the absence of feature importance analysis in many studies limits their utility in formulating targeted teacher retention strategies.

Based on the aforementioned research gap, this study proposes two research questions: (1) What are the key dimensions influencing faculty turnover in private universities, and what is the relative importance (predictive weight) of each dimension in predicting turnover? (2) Among the three machine learning algorithms—KNN, NB, and BPNN Neural Network—which model demonstrates superior accuracy and stability in predicting faculty turnover in private universities?

Furthermore, grounded in classical theoretical frameworks such as Herzberg’s Two-Factor Theory and Maslow’s Hierarchy of Needs, this study puts forward three main research hypotheses: H1: Job satisfaction is negatively correlated with turnover intention. H2: The “Compensation, Benefits, and Development” dimension (including hygiene factors such as salary, security, and promotion) has a stronger predictive power on turnover than other dimensions. H3: The KNN model outperforms both the NB and BPNN models in terms of accuracy, F1-score, and AUC value in predicting faculty turnover.

The introduction of a longitudinal tracking design is a key methodological feature of this study. Multiple studies have demonstrated that, compared to cross-sectional data, longitudinal data can more effectively capture the dynamic evolution of employee attitudes and behaviors, thereby enhancing the timeliness and robustness of turnover predictions (Rubenstein et al., 2017; Chen et al., 2011). Particularly for highly mobile occupational groups, long-term tracking data facilitates the identification of both cumulative and stage-specific factors leading to turnover (Hom et al., 2017). Specifically, the operational objectives of this study are as follows: (1) to revise a questionnaire based on established scales and systematically collect data on teacher satisfaction and turnover over a 10-year period; (2) to employ Exploratory Factor Analysis (EFA) to extract core influencing factors from the multidimensional questionnaire and clarify the predictive weight of each dimension on turnover; (3) to construct three machine learning prediction models and compare their performance differences; and (4) to propose actionable teacher retention intervention strategies by integrating the model results with qualitative analysis.

By integrating theoretical frameworks from psychology and organizational management, this study not only addresses the limitations of existing turnover prediction models but also offers practical strategies for human resource management in private universities. The principal innovations of this work are threefold: First, it introduces, for the first time, a 10-year longitudinal dataset focused on faculty in private universities, thereby enhancing the timeliness and robustness of the findings. Second, it proposes a hybrid feature extraction method based on factor analysis, overcoming the feature selection constraints inherent in traditional models, and further validates the effectiveness of the EFA framework through qualitative analysis. Third, it innovatively integrates teachers’ psychological characteristics with organizational management factors, and via a comparative analysis of multiple machine learning models, demonstrates the superiority of the KNN model in predicting faculty turnover. Moving beyond the reliance of conventional models solely on objective data, this research provides valuable guidance for developing teacher retention strategies in private universities, bearing significant theoretical and practical implications.

2 Literature review

2.1 Theoretical foundations of teacher turnover research

Job satisfaction among teachers is a multidimensional concept influenced by various psychological, organizational, and environmental factors. Herzberg’s Two-Factor Theory posits that job satisfaction is driven by two types of factors: hygiene factors (such as salary and working conditions) and motivators (such as recognition and career development) (Herzberg, 1966). Hygiene factors prevent employee dissatisfaction, while motivators actively enhance satisfaction and loyalty. In the context of higher education, studies suggest that compared to base salary levels, faculty satisfaction in private universities is more sensitive to motivational factors such as research autonomy and career development opportunities (Ali and Botha, 2021). This theory supports the inclusion of factors like “career development” and “administrative support” in satisfaction analyses. Maslow’s Hierarchy of Needs further emphasizes the hierarchical nature of human motivation, indicating that unmet lower-level needs (e.g., job security) can hinder the pursuit of higher-level needs (e.g., self-actualization through teaching or research) (Maslow, 1943). For faculty in private institutions, job stability (a security need) often conflicts with academic freedom (a self-actualization need). When universities fail to balance these priorities, unique tensions arise that may lead to faculty turnover (Kim and Lee, 2020). The Price-Mueller model integrates environmental, individual, and organizational variables to predict turnover behavior (Price and Mueller, 1986). This model identifies job satisfaction and organizational commitment as key mediating factors between external opportunities (e.g., alternative employment options) and turnover decisions. Research shows that even in an unfavorable external job market, dissatisfaction with workload or administrative policies can significantly increase teachers’ turnover intention (Chen et al., 2022).

2.2 Applications and research advances of machine learning in turnover prediction

In recent years, machine learning (ML) methods have been widely applied in employee turnover prediction research due to their powerful data-driven and pattern recognition capabilities (Akasheh et al., 2024). As a core branch of artificial intelligence, machine learning can automatically learn patterns from historical data and build predictive models, demonstrating particular strength in handling high-dimensional features and non-linear relationships (Hastie et al., 2009). Based on the availability of data labels and learning objectives, machine learning can be categorized into four main types: supervised learning, unsupervised learning, semi-supervised learning, and reinforcement learning. In the field of employee turnover prediction, supervised learning algorithms such as K-Nearest Neighbors (KNN), Naive Bayes (NB), and Neural Networks (NN) have been shown to be highly effective. For instance, Lazzari et al. (2022) employed seven classifiers, including k-nearest neighbors and decision trees, to analyze multinational survey data and construct a model of employee turnover intention. Their findings indicated that logistic regression and LightGBM performed particularly well in predicting employees’ turnover intention. Park et al. (2024) utilized survey data on career mobility among university graduates from the Korean Employment Information Service, applying logistic regression (LR), K-nearest neighbors (KNN), and Extreme Gradient Boosting (XGB) classifiers for modeling. Their study identified job security as the most significant predictor. Notably, Yiğit and Shourabizadeh (2017) compared the performance of multiple classifiers in predicting employee churn and found that Support Vector Machine (SVM) achieved the best results in terms of accuracy, precision, and F1-score. Fallucchi et al. (2020) validated the advantage of Gaussian Naive Bayes in recall using the IBM dataset, while Al-Darraji et al. (2021) adopted a deep neural network model, achieving an accuracy of 89.11% in employee turnover prediction.

Specifically focusing on the three algorithms examined in this study:

(1) K-Nearest Neighbors (KNN) is an instance-based lazy learning algorithm that classifies data by measuring the distance between different feature values. It is suitable for small-scale datasets and is insensitive to outliers (Cover and Hart, 1967).

(2) Naive Bayes (NB) is based on Bayes’ theorem and assumes feature independence. It is simple, efficient, and well-suited for processing high-dimensional data (Jadhav and Channe, 2016).

(3) Backpropagation Neural Network (BPNN) is a multi-layer feedforward neural network that adjusts weights through the error backpropagation algorithm, exhibiting strong non-linear fitting capability (Rumelhart et al., 1986).

These algorithms have been proven effective in multiple turnover prediction studies.

2.3 Research review

The organizational environment is a critical factor influencing faculty turnover. A positive school climate, adequate resources, and supportive administrative management can significantly reduce teacher attrition rates. Conversely, chaotic management, insufficient resources, and stringent accountability systems contribute to teacher dissatisfaction and turnover (Jackson, 2012). Research within the Chinese context indicates that urban living cost pressures, blocked promotion pathways, and a lack of organizational cultural identity constitute three major triggers of faculty attrition in private institutions (Yang and Lu, 2024). Based on survey data from 1,161 teachers in private universities, Chen et al. (2023) found that nearly one-third of the respondents reported turnover intention, with doctoral-degree holders exhibiting the strongest inclination to leave.

However, existing studies utilizing machine learning for predicting faculty turnover have predominantly focused on corporate employees and relied heavily on objective data such as salary levels. These studies often overlook psychological characteristics specific to teachers—such as job satisfaction and professional efficacy—and fail to provide feature importance rankings. As a result, they offer limited practical guidance for developing targeted teacher retention strategies in educational settings (Kim and Lee, 2020). Furthermore, few existing models integrate Exploratory Factor Analysis (EFA) with machine learning methods, thereby inadequately addressing the multicollinearity issues inherent in psychological characteristics (Holton et al., 2022).

To address these limitations, this study proposes a hybrid predictive framework that combines EFA and machine learning. First, EFA is employed to extract latent influencing factors from multidimensional satisfaction questionnaires, effectively mitigating feature redundancy. Subsequently, comparative models are constructed using K-Nearest Neighbors (KNN), Naive Bayes (NB), and Backpropagation Neural Networks to validate the advantages of nonlinear algorithms in predicting faculty turnover. This approach not only enhances the theoretical depth of turnover prediction research but also promotes practical innovation in the field.

3 Research object and methods

3.1 Research objects

This study focuses on a private university in Western China, referred to as University A, as the empirical research site. University A is one of the first private universities approved by the Ministry of Education. As of December 2023, the university enrolls approximately 20,000 full-time undergraduate students and holds complete academic qualifications. It enjoys a strong social reputation, with no significant public controversies or regulatory penalties during its operation.

The study adopts a longitudinal tracking design, with data collection spanning from January 2015 to December 2024. As shown in Table 1, over the past decade, University A has recruited 497 young and middle-aged faculty members aged between 26 and 45. Of these, 291 teachers (58.55%) remain, while 206 teachers (41.45%) have left. Personnel records were exported from the university’s human resources management system, and data cleaning was performed using a double-blind entry method. Ethical approval for the study was obtained from the Academic Committee of University A. All participants were informed about the survey, provided informed consent, and were assured of the confidentiality and anonymity of their data. In addition, this study commits to providing feedback on the overall research findings to participants after the completion of the study. Specific methods include disseminating results via university email systems and organizing dedicated results-briefing sessions, ensuring that no personal data is involved during the feedback process. These measures further protect participants’ right to information and privacy.

Table 1
www.frontiersin.org

Table 1. Data on new hires and departures from university a in the past decade.

3.2 Research methods

The study designed the “Teacher Satisfaction Survey for Young and Middle-Aged Faculty in Private Universities,” proposing a feature selection method based on factor analysis and identifying four key feature variables. Subsequently, three machine learning models were used to fit the teacher turnover prediction. The models’ prediction performance was evaluated and compared using metrics such as accuracy, F1 score, and AUC value.

3.2.1 Questionnaire survey method

This study was based on the core dimensions of the short-form Minnesota Satisfaction Questionnaire (MSQ) (Weiss et al., 1967), which contains 20 items, and the Job Satisfaction Survey (JSS) (Spector, 1985). The questionnaire was adapted to align with the organizational characteristics of private universities. The resulting “Teacher Satisfaction Survey for Young and Middle-Aged Faculty in Private Universities” was developed (see Table 2). The survey used a Likert 5-point scale and consisted of 21 items. The Cronbach’s α coefficient for the questionnaire was 0.92, indicating strong internal consistency, with a content validity of 0.91. Additionally, the item discrimination was found to be satisfactory.

Table 2
www.frontiersin.org

Table 2. TEACHER satisfaction survey for young and middle-aged faculty in private universities.

3.2.2 Questionnaire distribution and recovery

This study utilized the university’s human resource management system to conduct an annual tracking survey of faculty satisfaction each December. To ensure a fair comparison of satisfaction levels between teachers who left and those who remained, the analysis uniformly used the most recent questionnaire results obtained during each teacher’s period of employment. Specifically: (1) for teachers who had left, data from their most recent annual survey before departure were used; (2) for teachers still employed, data from the latest survey (December 2024) were used. This approach ensured that all satisfaction data were sourced from the most recent available point in each teacher’s career, effectively capturing their latest state either prior to leaving or while currently employed. This method effectively mitigates measurement errors that could arise from differences in career stages or recent events. Following the above screening and data cleaning procedures, a total of 428 valid questionnaires were included in the final analysis, comprising 227 from currently employed teachers and 201 from those who had left. The distribution of the surveyed teachers by gender, age, educational background, and professional title is presented in Table 3.

Table 3
www.frontiersin.org

Table 3. BASIC information of the survey.

In the “Satisfaction Survey for Young and Middle-aged Faculty in Private Universities” employed in this study, Item 11 (“I find my current work tedious and busy”) and Item 13 (“My compensation and benefits fail to meet expectations”) were designed as reverse-scored items. This design aimed to control for response bias and avoid habitual responding from participants. The scoring method for the scale was as follows: As the study focuses on turnover prediction, positively worded items were scored as: “Strongly Disagree” = 5 points, “Disagree” = 4 points, “Neutral” = 3 points, “Agree” = 1 point, “Strongly Agree” = 1 point. Conversely, reverse-scored items were scored as: “Strongly Disagree” = 1 point, “Disagree” = 2 points, “Neutral” = 3 points, “Agree” = 4 points, “Strongly Agree” = 5 points. This scoring scheme ensures that the direction of scoring is consistent with the conceptual meaning of the measured dimensions across all items.

3.2.3 Feature extraction method

To address the high-dimensional characteristics and potential information redundancy of the “Teacher Stability Survey for Young and Middle-Aged Faculty in Private Universities,” this study proposes a hybrid feature optimization framework based on factor analysis (Wang et al., 2024). Using EFA for dimensionality reduction, key information from the data can be preserved while reducing its dimensionality. This approach decreases computational complexity and helps classifiers perform classification and prediction more accurately, thereby assessing employee turnover.

Suppose there are n samples and p indicators. Let the random vector X=(X1,X2,,Xp)T, and the common factors to be found be f=(f1,f2,,fm)T. Then the calculation model is shown in Equation 1.

{x1=u1+a11f1+a12f2++a1mfm+ε1x2=u2+a21f1+a22f2++a2mfm+ε2xp=up+ap1f1+ap2f2++apmfm+εp    (1)

It can also be expressed as x = u + Af + ε, where f is the common factor vector, ε is the specific factor vector, and aij is the factor loading matrix.

The degree to which the original variable xi depends on the common factors is expressed as Equations 2 and 3:

hi2=j=1maij2    (2)
Var(xi)=hi2+εi2(i=1,2,,p)    (3)

hi2 represents the variance contribution of the common factors to xi, also known as the communality. And εi2 represents the variance contribution of εi2 to xi, also known as the uniqueness. Moreover, hi2+εi2=1.

The factors obtained from the principal component analysis were used as the extracted features, which were used as independent variables in the subsequent model fitting for turnover prediction.

3.2.4 Prediction models

3.2.4.1 K-Nearest neighbors (KNN)

KNN is a lazy machine learning method that supports incremental learning. It can be used for non-linear classification and can handle complex decision spaces. When the data volume is small but representative, it has a good classification effect. KNN classifies a sample to be classified by calculating the distances between it and the training samples and selecting the KNN (Everitt et al., 2011). When calculating the sample distances, distance calculation methods such as the Manhattan distance, Euclidean distance, and Chebyshev distance can be used, and the Euclidean distance is generally used. The disadvantages of the KNN are the high computational cost and poor interpretability. In the prediction of teachers’ turnover, the KNN algorithm can intuitively reflect the similarity between samples. Its main calculation process is as follows:

For two n-dimensional samples, and y=(y1,y2,,yn) the Euclidean distance between them is shown in the following Equation 4:

d(x,y)=i=1n(xiyi)2    (4)
3.2.4.2 Naive Bayes (NB)

Naive Bayes is a classification method based on Bayes’ theorem, assuming that all features are independent of each other (Panda and Patra, 2008). In the prediction of teachers’ turnover, the Naive Bayes algorithm is simple and efficient and is suitable for processing large-scale data sets. However, this algorithm has a strong assumption of feature independence; therefore, the correlation between features needs to be carefully evaluated in practical applications. The calculation process of the Bayes’ Equation is shown in the following Equation 5.

d(x,y)=i=1n(xiyi)2    (5)

P(yx) Where is the posterior probability, is the likelihood probability, is the prior probability, and P(x) is the evidence factor.

3.2.5 Backpropagation neural network (BPNN)

A neural network is a computational model that mimics the connection structure of human brain neurons and has a powerful non-linear mapping ability. In principle, first, the feature items are input, that is, the independent variable items are introduced. When using a neural network model, “pseudo features” can be constructed from the feature items. For example, if the inputs are gender, age, salary, years of work, etc., some “pseudo feature items” (i.e., feature items that do not actually exist and are completely constructed by the model and cannot be interpreted) can be constructed in combination with the activation function. Moreover, there can be multiple levels for constructing “pseudo feature items” (i.e., there can be multiple layers of “hidden layer neurons”), and each layer can have multiple neurons. Finally, the prediction item is obtained through calculation using a mathematical optimization algorithm. In the prediction of teachers’ turnover, the neural network can automatically extract complex patterns in the data, thus achieving high-precision prediction (Liu et al., 2015).

The Equation of the neuron model (taking a single neuron as an example) is See Equation 6:

y=f(i=1nwixi+b)    (6)

Where xi is the input, wi is the corresponding weight, b is the bias, f(·)is the activation function, and y is the output of the neuron.

The forward propagation calculation Equations of a multilayer neural network are as follows (where x is the input layer, h is the hidden layer, and y is the output layer):

The Equation for the hidden layer can be seen in Equation 7.

hj=f(i=1nwij(1)xi+bj(1)),j=1,2,,m    (7)

where m is the number of hidden layer neurons.

The Equation for the Output layer can be seen in Equation 8.

yk=f(j=1mwjk(2)hj+bk(2)),k=1,2,,l    (8)

Where l is the number of output layer neurons.

3.2.6 Evaluation metrics

Commonly used evaluation indicators for machine—learning classification models include accuracy (Accuracy, ACC), ROC curve (Receiver Operating Characteristic Curve), F1 score, etc. (Ramakrishnan et al., 2024).

3.2.6.1 Accuracy

According to Equation 9, accuracy refers to the proportion of the number of correctly classified samples to the total number of samples. It intuitively reflects the correctness of the model’s overall classification of the samples.

Accuracy=TP+TNTP+TN+FP+FN    (9)

where TP (True Positive) represents the number of true positives, that is, the number of samples that are actually positive and are predicted as positive by the model; TN (True Negative) represents the number of true negatives, that is, the number of samples that are actually negative and are predicted as negative by the model; FP (False Positive) represents the number of false positives, that is, the number of samples that are actually negative but are predicted as positive by the model; FN (False Negative) represents the number of false negatives, that is, the number of samples that are actually positive but are predicted as negative by the model.

3.2.6.2 F1 Score

The F1 score is an indicator to measure the performance of a classification model, the details are presented in Equation 10. It is the harmonic mean of precision and recall. The maximum value is 1, and the minimum value is 0. It comprehensively considers the accuracy and completeness of the model, and a higher value indicates better overall performance. Among them, precision represents the proportion of actually positive samples among the samples predicted as positive, and recall represents the proportion of samples that are actually positive and are correctly predicted as positive.

F1=2×Precision×RecallPrecision+Recall    (10)
3.2.6.3 ROC curve (receiver operating characteristic curve)

The ROC curve is plotted with the False Positive Rate (FPR) on the x-axis and the True Positive Rate (TPR) on the y-axis. This shows the model’s ability to distinguish between positive and negative examples under different classification thresholds. The True Positive Rate (TPR) represents the proportion of actually positive samples that are correctly predicted as positive. The Equation for TPR is refer to Equation 11:

TPR=TPTP+FN    (11)

The False Positive Rate (FPR) represents the proportion of actually negative samples that are incorrectly predicted as positive. The Equation for FPR is refer to Equation 12:

FPR=FPFP+TN    (12)

The AUC (Area Under the Curve) is the area under the ROC curve, and its value ranges from 0.5 to 1. The larger the AUC, the better the classification performance of the model.

4 Results and discussion

4.1 Feature extraction

4.1.1 Factor analysis

This section examines the information compression of the questionnaire data for young and middle-aged faculty in private universities. First, the suitability of the factor analysis for the standardized scores of each sample was tested, including the KMO test and Bartlett’s sphericity test. The results are as follows: (1) The KMO sampling adequacy statistic value is 0.93, indicating a high degree of information overlap and strong correlations between the variables. (2) The Bartlett test yielded a chi-square value of 7855.43, with a significant p-value less than 0.05, indicating that the correlations between the variables are significant. Both the KMO and Bartlett’s sphericity tests confirm that the sample data are suitable for factor analysis.

Using the maximum variance rotation method, the variance explained by the first four factors is presented in Table 4. After rotation, the first four factors had eigenvalues greater than 1, and the cumulative variance explained reached 75.59%. Thus, the first four factors were extracted, which successfully condensed the information from the stability survey for young and middle-aged faculty. The variance explained by factor 1 after rotation was 25.41%, which was higher than the variance explained by the other factors, indicating that factor 1 is the most important factor influencing teacher turnover.

Table 4
www.frontiersin.org

Table 4. Variance explanation rates of the factors.

Table 5 presents the rotated factor loading coefficients. The results indicate that the communality value for each item exceeded 0.4, signifying a strong association between the items and the factors, and confirming that the factors effectively extract meaningful information. Based on the rotated solution, Factor 1, which primarily loads on items related to salary, social security and housing fund, post-appointment training, and academic credential advancement, was labelled “Compensation, Benefits, and Development.” Factor 2, characterized by high loadings on items concerning sense of achievement, job challenge, and utilization of professional expertise, was named “Professional Efficacy.” Factor 3, which is predominantly associated with school administration, performance evaluation systems, and teaching schedule arrangements, was designated the “School Management” factor. Finally, Factor 4, mainly relating to office facilities and interpersonal environment, was labelled “Work Environment.”

Table 5
www.frontiersin.org

Table 5. Factor loading coefficients after rotation.

Through dimensionality reduction via factor analysis, redundant features were effectively eliminated, allowing the model to focus on the most discriminative factors (F1: Compensation, Benefits, and Development; F2: Professional Efficacy; F3: School Management; F4: Work Environment). This outcome aligns with Herzberg’s Two-Factor Theory, wherein both hygiene factors (e.g., salary, work environment) and motivators (e.g., career development, sense of efficacy) collectively influence faculty turnover.

4.1.2 Textual analysis

Based on the dimensional framework derived from the EFA, this study conducted a systematic coding analysis of the textual data from the open-ended question “Other suggestions for enhancing the job satisfaction of young and middle-aged teachers” in the questionnaire. Two independent coders were selected to perform the coding task. Both coders had backgrounds in higher education research (one doctoral candidate and one lecturer). They received unified training prior to coding. Inter-rater reliability was assessed using Cohen’s Kappa coefficient (κ = 0.82, p < 0.001), indicating high reliability of the coding results. The final coding results are presented in Table 6.

Table 6
www.frontiersin.org

Table 6. Text coding multiple response table.

As shown in Table 6, the dimension of welfare benefits and career development (F1) received the highest level of attention (n = 56, 48.7%). Core appeals focused on optimizing the compensation system (frequency 28), improving social security benefits (frequency 18), and reforming the establishment system (frequency 10), among other aspects. Representative excerpts included:

The school management dimension (F3) ranked second (n = 29, 25.2%). Respondents primarily suggested improvements regarding the rigidity of the quantitative assessment mechanism (frequency 15), the pressure from the “up-or-out” policy (frequency 9), and the excessive burden of administrative tasks (frequency 8). Representative excerpts included:

The occupational efficacy dimension (F2) garnered 14 suggestions (12.2%), emphasizing investment in disciplinary construction and the allocation of research resources. A representative excerpt included:

"The basic salary leaves little remains after deducting rent and living expenses. We hope the school can raise the salary standard and adjust the social security and housing fund contribution base to the actual income level. The gap in benefits compared to public institutions is simply too large."

"We, teachers in private institutions, lack the status of public institutional establishment, and there is no security for post-retirement benefits. If the school could promote reform of the establishment system or supplement commercial pension insurance, it would certainly help retain personnel."

The school management dimension (F3) ranked second (n = 29, 25.2%). Respondents primarily suggested improvements regarding the rigidity of the quantitative assessment mechanism (frequency 15), the pressure from the "up-or-out" policy (frequency 9), and the excessive burden of administrative tasks (frequency 8). Representative excerpts included:

"The school's evaluation overemphasizes teaching hours and paper counts. Even with excellent teaching outcomes, promotion is impossible without sufficient research points. This 'quantification-only' standard is unsuitable for humanities teachers."

"Filling out five or six administrative forms weekly and attending various irrelevant meetings leaves less than half the time for actual lesson preparation and research. We hope to streamline administrative procedures and return energy to teaching and research."

The occupational efficacy dimension (F2) garnered 14 suggestions (12.2%), emphasizing investment in disciplinary construction and the allocation of research resources. A representative excerpt included:

“The school’s laboratory equipment is outdated, and the process for applying research funding is cumbersome. There are simply no conditions for conducting innovative research. It would be helpful to increase seed funding for research and equipment investment.”

The work environment dimension (F4) received 16 suggestions (13.9%), primarily focusing on practical needs such as housing security (frequency 5) and commuting convenience (frequency 3). A representative excerpt included:

“The school is located too far from the city center, surrounding rentals are expensive and inconvenient. We hope the school can provide faculty transitional housing or housing subsidies.”

Qualitative analysis further revealed that the F1-F4 dimensions covered 92.4% (n = 115) of all suggestions, and their internal item factor loadings were all >0.65, confirming the effectiveness of the EFA model. The results also provide empirical evidence for private university administrators to build a multidimensional retention mechanism.

4.2 Model prediction results

The factors obtained from the factor analysis (F1 Welfare and Development, F2 Career Effectiveness, F3 School Management, F4 Work Environment) were used as independent variables, with employee turnover (whether the employee left or not) as the dependent variable. The training set proportion was set to 0.8, and three models—K-Nearest Neighbors (KNN), Naive Bayes (NB), and Backpropagation (BPNN) neural network—were used for modeling. The main parameter settings for each machine learning model are shown in Table 7.

Table 7
www.frontiersin.org

Table 7. MAIN parameter settings of the models.

4.2.1 Model performance comparison

The prediction results of the models on the training and test sets are shown in Table 8. On the training set, both KNN and the BPNN demonstrated the highest classification capability, with accuracies of 84.21 and 83.04%, respectively. Their recall rates both exceeded 85%, indicating that these two types of models can effectively identify positive samples within the target category. On the test set, although NB had the lowest accuracy on the training set (77.19%), it exhibited generalization ability comparable to the best model on the independent test set, with an accuracy of 83.72% and an F1-score of 84.09%.

Table 8
www.frontiersin.org

Table 8. Model prediction results for the training and test sets.

Table 9 provides the prediction results of various models for the overall training and testing datasets. All models showed accuracy rates ranging from 79.05 to 83.64% on the overall dataset, with KNN achieving the best performance, with an overall accuracy rate of 83.64% and an F1 score of 84.16%.

Table 9
www.frontiersin.org

Table 9. Prediction results of multiple models for the overall data of the training and test sets.

4.2.2 ROC curve validation

ROC curves were constructed for KNN_Prediction_probability, NB_Prediction_probability, and BP_Prediction_probability, with the three models showing high diagnostic value for turnover prediction. As shown in Table 10, all models demonstrated strong discriminative ability, with the Area Under the Curve (AUC) values all exceeding 0.80, indicating high prediction accuracy. Specifically, the KNN model had the highest AUC value of 0.901 (95% confidence interval: 0.872 ~ 0.930), followed by the BPNN (AUC = 0.888; 95% confidence interval: 0.855 ~ 0.921).

Table 10
www.frontiersin.org

Table 10. ROC results AUC summary.

Figure 1 shows the ROC curves for the three models. The KNN model’s curve is positioned at the topmost, closest to the ideal model point (where sensitivity = 1 and specificity = 1), indicating that it performed the best in distinguishing whether teachers left or not.

Figure 1
Receiver Operating Characteristic (ROC) curve comparing three prediction models: KNN, BP, and NB. The x-axis represents 1-Specificity, and the y-axis represents Sensitivity. KNN is shown with a solid line, BP with a dashed line, and NB with a dash-dotted line. The KNN model shows a step-like pattern, while BP and NB have smoother curves.

Figure 1. Roc curve.

5 Discussion

This study reveals a significant association between the four dimensions extracted by EFA (welfare benefits and career development, occupational efficacy, work environment, and school management) and teacher turnover behavior, further supporting the applicability of classical organizational behavior theories in the context of private education. For instance, the “welfare benefits and career development” factor, as the dimension with the highest weight, strongly aligns with the fundamental role of “hygiene factors” in Herzberg’s Two-Factor Theory (Herzberg, 1966). This finding is also consistent with conclusions from numerous recent machine learning-based employee turnover prediction studies, where multiple reviews and empirical analyses indicate that salary is one of the most important and stable predictors (Akasheh et al., 2024; Najafi-Zangeneh et al., 2021). Simultaneously, the importance of the “occupational efficacy” factor echoes the persistent influence of higher-level self-actualization needs on career choices within Maslow’s Hierarchy of Needs theory (Maslow, 1943). This result aligns with the findings of Chen et al. (2023) regarding private university teachers, which identified professional development opportunities and work autonomy as key motivational factors affecting turnover intention. Furthermore, Ma’s (2013) empirical research on occupational stress and turnover intention among teachers in private higher vocational colleges indicated that school management, along with four other stress factors—teacher compensation, role responsibilities, personal development, student factors, and interpersonal relationships—significantly influenced turnover intention. This highly overlaps with the four dimensions identified in this study, further confirming the composite nature of turnover drivers among private university teachers. These drivers encompass basic hygiene factors such as compensation and security, as well as motivational factors involving management climate and career growth, fitting the multi-causal and hierarchical patterns revealed by general employee turnover prediction research (Akasheh et al., 2024; Rombaut and Guerry, 2017). This suggests a certain degree of stability and generalizability in the factors influencing turnover among private university teachers across different periods and research contexts.

Regarding model prediction performance: the KNN algorithm demonstrated relatively outstanding results in the overall tests. Particularly on the test set, all four key classification metrics reached relatively high values, with accuracy and F1-score reaching 83.64 and 84.16%, respectively. This indicates that the algorithm effectively captured the data characteristics within this study’s dataset, resulting in favorable predictive performance for teacher turnover. This finding aligns with conclusions from existing studies utilizing machine learning for employee turnover prediction, which suggest that KNN exhibits stable performance when handling small-to-medium-scale classification problems with distinct features (Yiğit and Shourabizadeh, 2017; Lazzari et al., 2022). The neural network model demonstrated strong generalization ability between the training and test sets, highlighting its capability to handle complex nonlinear relationships. By leveraging their powerful capacity for feature combination and mapping, neural networks have shown high discriminative power in several employee retention prediction studies (Al-Darraji et al., 2021; Srivastava and Eacchempati, 2021). However, its high computational complexity and reliance on large-scale data limited its practicality in the context of this study’s limited sample size. In contrast, the Naïve Bayes algorithm achieved the highest recall rate on the test set (92.50%), indicating its unique value in reducing the “risk of missed detection.” This makes it particularly suitable for institutions that prioritize the identification of high-risk teachers. Similarly, Fallucchi et al. (2020) also noted that Naïve Bayes performs excellently in terms of recall rate when processing highly imbalanced human resources data, making it especially apt for the preliminary screening of high-risk populations.

6 Conclusion

This study addresses the issue of teacher turnover prediction in private universities in China, using a questionnaire survey method and applying EFA for feature extraction. Three machine learning models were compared based on their ability to predict teacher turnover. The main conclusions are as follows:

The EFA framework helps institutions identify key drivers of turnover. For example, “Welfare and Development” (variance contribution rate: 25.41%) and “Career Effectiveness” (variance contribution rate: 17.80%) are significantly associated with negative evaluations that influence turnover decisions. The “Welfare and Development” factor highlights the importance of material factors, such as salary, training, and promotion opportunities, in teachers’ turnover decisions. The “Career Effectiveness” factor emphasizes the role of psychological factors, such as job satisfaction and work challenges, To a certain extent, it provides a clear direction for colleges and universities to develop targeted teacher retention strategies. Text coding reveals that 48.7% of improvement suggestions focus on issues related to the F1 factor, such as optimizing the salary system and enhancing social security benefits. Qualitative analysis further validates the effectiveness of the EFA framework.

In terms of model prediction, the K-Nearest Neighbors (KNN) algorithm is relatively optimal, with an accuracy of 83.64%, an F1 score of 84.16%, and an ROC curve AUC value of 0.901. A comparison of multiple machine learning models revealed that Naive Bayes (NB), with its high recall rate, is well-suited for initial risk screening, while KNN’s stability makes it ideal for in-depth analysis of key personnel. On the premise that data quality and feature structure are stable, it provides a quantitative analysis tool for human resource management in private colleges and universities.

The “EFA + Machine Learning” hybrid framework constructed in this study provides a reference for conducting turnover prediction in similar institutions. By employing EFA to eliminate redundant features and focus on influential factors, it addresses the issue of multicollinearity in psychological characteristic data. Combining this with multi-model comparative validation overcomes the reliance on a single algorithm common in traditional research. Furthermore, the non-intrusive design of the satisfaction questionnaire reduces teachers’ psychological defensiveness and enhances data authenticity, making this framework more practical in sensitive contexts like private universities where personnel data is highly guarded.

6.1 Research limitations

1. Limited sample representativeness: The research data were sourced from only one private undergraduate university in Western China. Consequently, the applicability of the conclusions is primarily limited to private universities similar to the research object in terms of scale and faculty structure. The framework’s adaptability to public universities or private universities in different regions and at different levels requires further validation.

2. Potential bias in data processing: In the data collection process, data from teachers who left was their last survey response before departure, while data for retained teachers came from the most recent survey (December 2024). The time difference and status disparity between these data points may introduce bias.

3. Limitations in variable dimensions: The questionnaire design primarily relied on subjective satisfaction feedback and did not incorporate key objective performance indicators (such as salary records, research output) or external factors like the job market environment.

4. Lack of analysis on interacting factors: The study identified influencing factors within dimensions like “Welfare Benefits and Career Development” and “Occupational Efficacy.” However, constrained by the article’s length, it did not further explore the interaction effects between demographic characteristics (such as gender, professional title, age) and these core influencing factors.

6.2 Future research directions

Based on the above limitations, future research could be deepened in the following directions:

(1) Expanding samples and data sources: Future studies could extend to multi-regional institutional comparisons, integrate multi-source data (e.g., salary, research achievements) to build dynamic prediction systems, and explore the application of machine learning models for cross-institutional generalization.

(2) Optimizing data processing and model design: Subsequent research should consider introducing time-series analysis methods to construct dynamic turnover prediction models capable of capturing the evolution of teachers’ turnover intentions.

(3) Deepening the analysis of interacting factors: Research should investigate the interaction mechanisms between demographic dimensions (such as gender, age, education level, professional title) and core factors like “Welfare Benefits” and “Occupational Efficacy.”

(4) Integrating policy and environmental dynamics into the framework: External factors, such as macro-policy changes (e.g., reforms of the faculty establishment system in private universities) and regional disparities in educational resource allocation, should be incorporated into the analysis. Exploring the moderating effects of the institutional environment on teachers’ turnover decisions will make the research conclusions more aligned with the practical development context of private universities.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by The Research Department of Xi’an FanYi University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

WJ: Writing – original draft, Writing – review & editing, Conceptualization, Investigation, Software. LY: Data curation, Writing – review & editing. XY: Methodology, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the Xi’an Fanyi University Horizontal Project (Grant No.: 25XYH113) and the Xi’an Fanyi University Teaching Reform Project (Grant No.: J24B48).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Akasheh, M. A., Malik, S. F., Hujran, O. A., and Zaki, N. (2024). A decade of research on machine learning techniques for predicting employee turnover: a systematic literature review. Expert Syst. Appl. 238:121794. doi: 10.1016/j.eswa.2023.121794

Crossref Full Text | Google Scholar

Al-Darraji, S., Honi, D. G., Fallucchi, F., Abdulsada, A. I., Giuliano, R., and Abdulmalik, H. A. (2021). Employee attrition prediction using deep neural networks. Computer 10:141. doi: 10.3390/computers10110141

Crossref Full Text | Google Scholar

Ali, H., and Botha, N. A. (2021). Motivators and hygiene factors in higher education: a comparative study of public and private universities. J. Educ. Adm. 59, 321–335. doi: 10.1108/Jea-08-2020-0175

Crossref Full Text | Google Scholar

Bai, W. H. (2018). The lack of career attraction and its improvement for teachers in private universities. Heilongjiang Res. Higher Educ. 36, 37–41.

Google Scholar

Chen, G., Ployhart, R. E., Thomas, H. C., Anderson, N., and Bliese, P. D. (2011). The power of momentum: a new model of dynamic relationships between job satisfaction change and turnover intentions. Acad. Manag. J. 54, 159–181. doi: 10.5465/amj.2011.59215089

Crossref Full Text | Google Scholar

Chen, J. P., Cheng, H. Y., and Chen, Y. (2023). The relationship between teachers' working environment, subjective well-being and turnover intention: based on the analysis of survey data of teachers in private universities. J. Chongqing Univ. Sci. Technol. (Soc. Sci. Ed.) 2, 111–118. doi: 10.19406/j.cnki.cqkjxyxbskb.2023.02.004

Crossref Full Text | Google Scholar

Chen, X., Li, J., and Liu, Y. (2022). Predicting faculty turnover in Chinese private universities: an extension of the Price-Mueller model. High. Educ. Policy 35, 489–507. doi: 10.1057/s41307-021-00254-w

Crossref Full Text | Google Scholar

Cover, T., and Hart, P. (1967). Nearest neighbor pattern classification. IEEE Trans. Inf. Theory 13, 21–27. doi: 10.1109/TIT.1967.1053964

Crossref Full Text | Google Scholar

Everitt, B. S., Landau, S., Leese, M., and Stahl, D. (2011). Cluster analysis. 5th Edn. Chichester: John Wiley & Sons, Ltd.

Google Scholar

Fallucchi, F., Coladangelo, M., Giuliano, R., and William De Luca, E. (2020). Predicting employee attrition using machine learning techniques. Computer 9:86. doi: 10.3390/computers9040086

Crossref Full Text | Google Scholar

Hastie, T., Tibshirani, R., and Friedman, J. (2009). The elements of statistical learning: Data mining, inference, and prediction. 2nd Edn. New York, NY: Springer-Verlag.

Google Scholar

Herzberg, F. (1966). Work and the nature of man. Cleveland, Oh: World Publishing.

Google Scholar

Holton, M. K., Barry, A. E., and Chaney, J. D. (2022). Mitigating social desirability bias in organizational surveys: a meta-analysis of questionnaire design effects. Organ. Res. Methods 25, 369–392. doi: 10.1177/1094428120987426

Crossref Full Text | Google Scholar

Hom, P. W., Lee, T. W., Shaw, J. D., and Hausknecht, J. P. (2017). One hundred years of employee turnover theory and research. J. Appl. Psychol. 102, 530–545. doi: 10.1037/apl0000103

PubMed Abstract | Crossref Full Text | Google Scholar

Jackson, K. M. (2012). Influence matters: the link between principal and teacher influence over school policy and teacher turnover. J. Sch. Leadersh. 22, 875–901. doi: 10.1177/105268461202200503

Crossref Full Text | Google Scholar

Jadhav, S. D., and Channe, H. P. (2016). Comparative study of K-Nn, naive bayes and decision tree classification techniques. Int. J. Sci. Res. 5, 1842–1845. doi: 10.21275/NOV153131

Crossref Full Text | Google Scholar

Kim, S., and Lee, J. (2020). Maslow’s hierarchy in academia: balancing security and self-actualization among adjunct faculty. J. Career Dev. 47, 456–470. doi: 10.1177/0894845319833164

Crossref Full Text | Google Scholar

Lazzari, M., Alvarez, J. M., and Ruggieri, S. (2022). Predicting and explaining employee turnover intention. Int. J. Data Sci. Anal. 14, 279–292. doi: 10.1007/s41060-022-00329-w

Crossref Full Text | Google Scholar

Liu, J., Wang, H., Sun, Y., Fu, C., and Guo, J. (2015). Real-coded quantum-inspired genetic algorithm-based bp neural network algorithm. Math. Probl. Eng. 2015, 1–10.

Google Scholar

Maslow, A. H. (1943). A theory of human motivation. Psychol. Rev. 50, 370–396. doi: 10.1037/h0054346

Crossref Full Text | Google Scholar

Ma, X. Y. (2013). An empirical study on occupational stress and turnover intention of teachers in private higher vocational colleges. J. Jinhua Polytech. 13, 7–10.

Google Scholar

Ministry of Education of the People's Republic of China (2024). Statistical Bulletin of the Development of National Education in 2023. Available online at: http://www.moe.gov.cn/jyb_sjzl/sjzl_fztjgb/202410/t20241024_1159002.html (Accessed October 24, 2024).

Google Scholar

Najafi-Zangeneh, S., Shams-Gharneh, N., Arjomandi-Nezhad, A., and Hashemkhani Zolfani, S. (2021). An improved machine learning-based employees attrition prediction framework with emphasis on feature selection. Mathematics 9:1226. doi: 10.3390/math9111226

Crossref Full Text | Google Scholar

Panda, M., and Patra, M. R. (2008). “A comparative study of data mining algorithms for network intrusion detection,” in 2008 First International Conference on Emerging Trends in Engineering and Technology, Nagpur, India, 16–18 July 2008 (IEEE). pp. 504–507.

Google Scholar

Park, J., Feng, Y., and Jeong, S. P. (2024). Developing an advanced prediction model for new employee turnover intention utilizing machine learning techniques. Sci. Rep. 14, 1221. doi: 10.1038/s41598-023-50593-4

Crossref Full Text | Google Scholar

Price, J. L., and Mueller, C. W. (1986). Absenteeism and turnover among hospital employees. Greenwich, CT: Jai Press.

Google Scholar

Qi, Y. K. (2023). Research on the logic and management strategies of teacher mobility in private universities based on rational choice theory. Heilongjiang Res. Higher Educ. 41, 42–47. doi: 10.19934/j.cnki.shxyj.2023.06.005

Crossref Full Text | Google Scholar

Ramakrishnan, J., John, R., Mavaluru, D., Ravali, R. S., and Srinivasan, K. (2024). Transforming groundwater sustainability, management and development through deep learning. Groundw. Sustain. Dev. 27:101366. doi: 10.1016/j.gsd.2024.101366

Crossref Full Text | Google Scholar

Rombaut, E., and Guerry, M. A. (2017). Predicting voluntary turnover through human resources database analysis. Manag. Res. Rev. 41, 96–112. doi: 10.1108/Mrr-04-2017-0098

Crossref Full Text | Google Scholar

Rubenstein, A. L., Eberly, M. B., Lee, T. W., and Mitchell, T. R. (2017). Surveying the forest: a meta-analysis, moderator investigation, and future-oriented discussion of the antecedents of voluntary employee turnover. Pers. Psychol. 71, 23–65. doi: 10.1111/peps.12226

Crossref Full Text | Google Scholar

Rumelhart, D. E., Hinton, G. E., and Williams, R. J. (1986). Learning representations by back-propagating errors. Nature 323, 533–536. doi: 10.1038/323533a0

Crossref Full Text | Google Scholar

Spector, P. E. (1985). Measurement of human service staff satisfaction: development of the job satisfaction survey. Am. J. Community Psychol. 13, 693–713. doi: 10.1007/BF00929796

PubMed Abstract | Crossref Full Text | Google Scholar

Srivastava, P. R., and Eacchempati, P. (2021). Intelligent employee retention system for attrition rate analysis and churn prediction: an ensemble machine learning and multi-criteria decision-making approach. J. Glob. Inf. Manag. 29, 1–29. doi: 10.4018/JGIM.20211101.oa23

Crossref Full Text | Google Scholar

Wang, J., Yang, X., and Liu, D. (2024). A comparison of the validity of online courses quality evaluation models. Int. J. Web-Based Learn. Teach. Technol. 19, 1–17. doi: 10.4018/Ijwltt.349984

Crossref Full Text | Google Scholar

Weiss, D. J., Dawis, R. V., England, G. W., and Lofquist, L. H. (1967). Manual for the Minnesota satisfaction questionnaire. Minneapolis, Mn: University of Minnesota, Industrial Relations Center.

Google Scholar

Yang, X. Q., and Lu, Y. H. (2024). Causes and countermeasures for the brain drain of young teachers in private universities: a case study of a private university in Shanghai. J. Wuhan Univ. Technol. (Soc. Sci. Ed.) 27, 1103–1106.

Google Scholar

Yiğit, I. O., and Shourabizadeh, H. (2017). “An approach for predicting employee churn by using data mining,” in 2017 International Artificial Intelligence and Data Processing Symposium (Idap), Malatya, Turkey, pp. 16–17 September 2017 (IEEE). pp. 1–4.

Google Scholar

Keywords: teacher turnover, private universities, machine learning, exploratory factor analysis, job satisfaction

Citation: Jingwen W, Yi L and Xiaohong Y (2025) Predicting teacher turnover in private universities: a machine learning approach based on 10 years of data and satisfaction factors. Front. Psychol. 16:1670195. doi: 10.3389/fpsyg.2025.1670195

Received: 21 July 2025; Accepted: 14 October 2025;
Published: 06 November 2025.

Edited by:

Jun Yang, University of North Carolina at Greensboro, United States

Reviewed by:

Arlette J. Ngoubene-Atioky, Goucher College, United States
Enrique H. Riquelme, Temuco Catholic University, Chile

Copyright © 2025 Jingwen, Yi and Xiaohong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wang Jingwen, MTA4NDMxMzc4NkBxcS5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.