Next Article in Journal
A Deep Learning Neural Network Method Using Linear Eigenvalue Statistics for Schizophrenic EEG Data Classification
Next Article in Special Issue
A New Multi-Target Three-Way Threat Assessment Method with Heterogeneous Information and Attribute Relevance
Previous Article in Journal
Research on Time-Varying Meshing Stiffness of Marine Beveloid Gear System
Previous Article in Special Issue
New Hybrid EC-Promethee Method with Multiple Iterations of Random Weight Ranges: Applied to the Choice of Policing Strategies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Exploration of Prediction Performance Based on Projection Pursuit Regression in Conjunction with Data Envelopment Analysis: A Comparison with Artificial Neural Networks and Support Vector Regression

1
College of Humanities and Law, Shanghai Business School, Shanghai 200235, China
2
College of Economics and Management, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
3
School of Information Management, Shanghai Lixin University of Accounting and Finance, Shanghai 201209, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(23), 4775; https://doi.org/10.3390/math11234775
Submission received: 1 November 2023 / Revised: 17 November 2023 / Accepted: 22 November 2023 / Published: 26 November 2023

Abstract

:
Data envelopment analysis (DEA) is a leading approach in performance analysis and discovering newer benchmarks, and the traditional DEA models cannot forecast the future efficiency of decision-making units (DMUs). Machine learning, such as the artificial neural networks (ANNs), support vector machine/regression (SVM/SVR), projection pursuit regression (PPR), etc., have been viewed as beneficial for managers in predicting system behaviors. PPR is especially suitable for small and non-normal distribution samples, the usual cases in DEA analysis. This paper integrates DEA and PPR to cover the shortcomings we faced while using DEA and DEA-BPNN, DEA-SVR, etc. This study explores the advantages of combining these complementary methods into an integrated performance measurement and prediction model. Firstly, the DEA approach is used to evaluate and rank the efficiency of DMUs. Secondly, we establish two DEA-PPR combined models to describe the DEA efficiency scores (also called the production function) and the DEA-efficient frontier function. The first combined model’s input variables are input–output indicators in the DEA model, and the output variable is the DEA efficiency. In the second model, its input variables are input or output indicators in the DEA model, and the output variable is the optimal input indicator for input-oriented DEA or the output indicator for output-oriented DEA. We conducted positive research on two examples with actual data and virtual small, medium-sized, and large samples. Compared with the DEA-BPNN and DEA-SVR models, the results show that the DEA-PPR combined model has more vital global optimization ability, better convergence, higher accuracy, and a simple topology. The DEA-PPR model can obtain robust results for both small and large cases. The DEA-BPNN and DEA-SVR models cannot obtain robust results for small and medium-sized samples due to overfitting. For large samples, the DEA-PPR model outperforms DEA-BPNN, DEA-SVR, etc. The DEA-PPR combined model possesses better suitability, applicability, and reliability than the DEA-BPNN model, the DEA-SVR model, etc.

1. Introduction

The measurement of the efficiency or performance of homogeneous production or services, namely decision-making units (DMUs), has long been a hotspot or subject of interest for researchers in many fields. Since Charnes et al. [1] proposed the radial-based data envelopment analysis (called DEA, DEA-CCR, or CCR) under input-oriented cases with a constant return-to-scale (CRS), DEA is widely and successfully used to measure the performance or efficiency of DMUs [2,3]. DEA is a nonparametric basis for evaluating the multiple input–output efficiency of different DMUs. To better meet various situations and conditions, a variety of DEA models have been proposed successively [4]. Distinct calculation principles, a simple structure, fewer input–output indicators, low-cost data collection, and easy application characterize them.
Meanwhile, the DEA model also has some shortcomings. For example, DEA is very sensitive to data noise and cannot be used for prediction [5]. When new DMUs are added, they must be remodeled, and the efficiency of the original DMUs is changed. At the same time, the shortcomings of these DEA models happen to be the advantage of another nonparametric modeling approach—the artificial neural network (ANN) [6]. Furthermore, DEA models treat DMUs as a black box; inputs enter and outputs exit without considering the intervening steps [7]. Therefore, Athanassopoulos and Curram [8] proposed to combine the DEA and BPNN models for efficiency analysis and conducted empirical research on the operating efficiency of the Cobb–Douglas production function and bank; the results show that the DEA model has better performance in simulating production functions, and the BPNN model is equivalent to DEA in the efficiency ranking of DMUs; they conducted in-depth research on how to combine DEA and BPNN better to give full play to their respective advantages. Since then, extensive theoretical and applied research has been conducted on the combination of ANN and DEA [4,9,10,11,12,13,14,15,16]. Additionally, other machine learning methods (such as SVR and SVM) are also combined with DEA in applied and comparative studies [4,11,17,18,19,20]. Some scholars have optimal experimental research on the integrated model of RSM (response surface methodology) and ANN-DEA. So far, extensive and in-depth research has been conducted on combined models such as ANN-DEA, achieving good results. However, there are still areas to be improved and perfected. First, for small and medium-sized samples in DEA, the input–output data and the evaluation results (efficiency scores) usually do not obey normal distribution, and it is not easy to guarantee the generalization ability and practical value of ANN or SVR combined models. Although some studies divide the sample data into training and test subsets [11,13,21], all studies do not use validation subsets to monitor the training process to avoid over-training. In training the ANN-DEA combined model, we cannot judge whether over-training has occurred, so it is not easy to guarantee its generalization ability and practical value. Second, to establish SVR (SVM) and ANN, the principles of determining input and output variables in different studies are often inconsistent, resulting in unconvincing modeling results. Some scholars [11,13,21,22,23] took the input–output indicators of the DEA model as the input variables of the ANN model. The DEA model evaluated the efficiency scores as the output variables of the ANN for modeling, but did not discuss how to reasonably determine the number of hidden-layer neurons and avoid over-training. Farahmand et al. [24] and Zhu et al. [11] used the input–output indicators of the DEA model as the input variables of SVR and the efficiency value evaluated via DEA as the output variables of SVR. Through empirical research, Zhu et al. [11] believed that the performance of the ANN-DEA model was better than that of SVR-DEA. Zhong et al. [13] argued that the performance of the ANN-DEA model was better than that of SVR-DEA, random forest (RF-) DEA, gradient boosting regressor (GBR-) DEA, and other models. Bose [15] established an input-oriented CCR-DEA model, took values uniformly within the output-indicator range of DEA-efficient samples, generated a certain number of modeling virtual samples, and used the output indicators of DEA as the input variables of the ANN model and the input indicators of the DEA model as the output variables of the ANN model; without a test subset, the DEA-efficient frontier ANN model was established.
It can be seen from those mentioned above that although the academic community has carried out extensive research on the combined prediction model of ANN-DEA (including SVR-, RF-, etc.), there are still many problems to be further researched, especially on the efficiency of small and medium-sized samples that do not obey normal distribution; a better model should be established with a better generalization ability and robustness. On the other hand, projection pursuit regression (PPR) is particularly applied to modeling with nonlinear, small, and medium-sized samples that do not obey normal distribution [25,26,27,28]. Therefore, this paper introduces PPR into efficiency research to establish a DEA-PPR combined model to obtain more reliable and robust results. It expands not only the efficiency evaluation method but also the application field of the PPR model.
This study fulfills a practical need and improves benchmarking and decision-making processes by exploring an innovative performance measurement and prediction framework using DEA-PPR. The proposed combined model utilizes DEA as a preprocessor, and the subsequent PPR model conducts prediction tasks for the best performance output for each DMU. This paper provides insight into two examples to build the DEA-PPR combined models and compare the performance of different combined models such as DEA-BPNN, DEA-SVR, etc.
In summary, the primary purpose and motivations for this research paper are fivefold:
(1)
To present a new DEA-PPR combined model for performance measurement and prediction, thus bridging the research gap through methodological advancement;
(2)
To provide empirical support for the proposed model using two datasets through streamlining sequential processes of DEA measurement and DEA-PPR prediction;
(3)
To provide a way of thinking to improve managerial efficiency and enhance administrative flexibility in selecting actionable options from theoretical and practically feasible alternatives and potential progress monitoring via the DEA-PPR combined model;
(4)
To discuss the advantages and disadvantages of machine learning, such as BPNN, SVR, PPR, RF, etc.;
(5)
To put forward the basic principles and some matters needing attention for building machine learning, such as BPNN, PPR, SVR, etc.
The rest of this paper is organized as follows: Section 2 is a review of the literature on DEA combined models and PPR; Section 3 is about CCR-DEA and PPR modeling principles; Section 4 is on the empirical research on the production function of combined models of DEA and the DEA-efficient frontier function; Section 5 is the results and discussion; and Section 6 describes the limitations and future research fields.

2. Literature Review

2.1. DEA and Its Combined Models

DEA, proposed in 1978 [1], provides an effective tool for measuring the efficiency of the homogeneous DMUs of non-profit (e.g., schools, hospitals, and local authorities) and for-profit (e.g., banks, public houses, corporates, listed companies) organizations or persons (e.g., scholars, employees). Since then, many new DEA models have been introduced, such as BCC, SBM, super-efficiency (SE-) DEA, SE-SBM, (SE-)SBM with undesirable output, dynamic DEA, and the window analysis technique [2,3,4]. DEA models with multiple input and output indicators have been widely used to assess DMUs’ efficiency.
On the other hand, the DEA models used for performance assessment have inherent disadvantages. First, DEA models are easily affected by data noise. Second, when new DMUs are added, we must re-establish the DEA, and the efficiency of all the original DMUs must be changed. It is not easy to compare the efficiency of different datasets. To overcome the above disadvantages, Athanassopoulos and Curram [8] thought that DEA and artificial neural networks (ANNs) were nonparametric methods, in that no assumption is made, and the inputs and outputs are used to describe an operational process. They started work to combine DEA with BPNN and studied the known Cobb–Douglas production function with two inputs and one output (simulated data). A series of 16 datasets were artificially generated and tested using the ANN-DEA combined models. The number of nodes in the input and output layers of the ANN model was equal to the DEA model. The results showed that the DEA performed more satisfactorily as a tool for estimating empirical production functions. Furthermore, they researched 250 commercial bank branches, using four inputs and three outputs. Take the network (topology) structure of the BPNN model as 4-10-3 (that is, the input layer, hidden layer, and output layer have four, ten, and three neurons, respectively, and the number of network connection weights is 5 × 10 + 11 = 61); the input-layer and output-layer neurons of the BPNN model are the same as the input–output indicators of the DEA model; that is, the BPNN model is used to simulate the operation of the bank branch. Randomly select a 20% (50) proportion of test (validation) subset data; it is believed that the more validation data, the better the model’s generalization ability is. The error changes in the validation data are used to monitor the training process. When the error reaches the minimum and does not decrease with further training, we take the network connection weights with the minimum error of the validation data (called the early-stopping method), eliminate the influence of over-training as much as possible, and improve the model’s generalization ability and practical value. The results show that the ANN is instead a tool for obtaining relative rankings of DMUs based on their predicted outcome. Since then, scholars have extensively researched other DEA combined models. (These previous studies added meaningful value to the existing literature, as summarized in Table 1).
It can be seen from Table 1 that there are more than 10 DEA combined models, including ANN (mainly BPNN and RBFNN models), SVR/SVM, RF, the bagging model, etc. There are three main functions (purposes) of establishing a DEA combined model. The first purpose is to reveal the function relationships between the input–output indicators of the DEA model; that is, the input and output indicators of the ANN model are exactly or partially the same as the DEA model [8,15]. The second purpose is to reveal the working principle of the DEA model [11,14,23,29]; that is, the efficiency score calculated via the DEA model is used as the output value of the machine learning model (ML) (e.g., BPNN, RF, SVR, etc.); input–output indicators of the DEA model are used as the inputs of the ML; to take all DMUs as modeling samples; and to divide the samples into training and test datasets (also called validation datasets in some studies, but they play the role of a test dataset, which is only used to evaluate the performance of the model instead of monitoring the training process to avoid over-training). This combined model is a subsystem of the DEA model, which can be used to calculate the value of input–output efficiency. Still, it cannot be used to simultaneously calculate input excesses and output shortages (or shortfalls). The third purpose is to construct the DEA-efficient frontier function of DEA obtained by establishing the DEA model, that is, to use the DEA-efficient (valid) DMUs as the modeling samples of the ML [9,15,16,23]. Since there are a relatively small number of valid DMUs, some scholars generate a certain number of virtual DMUs by uniformly taking values within the range of the indicator value of valid DMUs. Some scholars established a BPNN model for small samples. Therefore, building a combined DEA model is not a substitute for the DEA model but can augment rather than replace the DEA model. Na et al. [9] aimed to determine the total staff employed in government organizations in China; collected statistical data from 27 provinces, municipalities, and autonomous regions; set five input indicators (such as population) and two output indicators (administrative expenditure and the number of officially employed staff); and established the CCR and CGS-DEA models. There were 10 and 13 provinces that were valid for DEA. Then, according to the 10 and 13 valid DEA samples (all as training samples; no test and validation samples), the input layer had five neurons (that is, the indicators of the DEA model), and the output layer had two neurons (the output indicators of the DEA model); a DEA-efficient frontier BPNN model was established with a network structure of 5-3-2 and 26 connection weights, which is significantly larger than the number of training samples. It does not meet the requirements for BPNN modeling, as the number of training samples must be more than the number of network connection weights. Ren et al. [29] established a DEA model for a total of four DMUs (A, B, C, and D), with two input indicators and three output indicators; they generated uniformly distributed samples at five levels (which were all training samples; no test and validation samples), used a genetic algorithm to build a BPNN model with a 10-21-1 network topology, and constructed an efficient frontier of DEA with 253 network connection weights, which is 50 times the number of training samples; it does not meet the basic requirements for building a BPNN model. Na et al. [9] and Ren et al. [29] did not use validation samples to monitor the training process and could not judge whether over-training occurred in the training process. Because the training samples were less than the number of connection weights [9], and the number of training samples was even less than the number of hidden-layer neurons [29], over-training must occur during the training process. Therefore, it is impossible for BPNN models established in this way to have a generalization ability and practical value. In [15], the first case was to establish a DEA model based on the collected data of doctors and nurses (two input indicators) and the number of outpatients and inpatients (two output indicators) from 12 hospitals, and obtain three valid DEA samples (A, B, and D); then, they generated 125 virtual DEA-efficient DMUs using linear combinations of real DEA-efficient DMUs, A, B, and D. The researchers established a BPNN model to calculate the optimal numbers of doctors and nurses according to the number of inpatients and outpatients; that is, the output indicators of the DEA model were used as the input variables of the BPNN model, and the input indicator of the DEA-efficient DMUs was used as the output variable of BPNN. The second example studied in [15] is based on the borrowing data of 99 Indian microfinance institutions that had been collected; there were five input indicators (such as assets, etc.) and two output indicators (the number of borrowers and borrowings); according to the sample data of the six DEA-efficient DMUs, a total of 729 virtual DEA-valid samples were generated. The authors used the DEA model inputs as BPNN inputs and the DEA outputs as the BPNN outputs, established the BPNN model, and tried to obtain the optimal levels of inputs under a certain number of borrowers and borrowings—assets, offices, personnel, and expenses. For the two examples in Ref. [15], the authors did not explain whether the samples were divided into training, test, and validation subsets, nor how to determine a reasonable number of hidden-layer neurons and avoid over-training. Therefore, the validity and reliability of the results of Ref. [15] require further investigation. Tsolas et al. [30] studied the operational efficiency of 160 bank branches in Greece based on operational data; there were three inputs, namely personnel expenses, rents and depreciation, and operational expenses, and two outputs, namely net interest income and non-interest income (fee and trading income). The researchers established input-oriented DEA models, mainly including radial CCR, BCC, NIRS, NDRS, and non-radial FDH. Russell and others obtained the input–output efficiency score, and divided the efficiency score into four levels of 1–4; the authors took the input–output indicators of the DEA model as the inputs of the BPNN model (six in total) and the efficiency score level as the output (one), and randomly divided the 160 samples into training and test subsets at a ratio of 8:2. They established a BPNN model with certain results; however, since the test subset was not used to monitor the training process in real time, it was impossible to determine whether over-training occurred or not during the training process. Know et al. [23] employed ten years of longitudinal data from Japanese electronic manufacturing firms from 2003 to 2012, and each firm-year was treated as an individual DMU; the researchers eliminated non-normal data and obtained 1419 samples for modeling. There were three inputs for the DEA model (namely employees, total assets, and operating expenses) and two outputs (namely revenue and market value); they established output-oriented CCR and BCC models to produce input and output projections for each inefficient DMU to become DEA-efficient. Then, two BPNN models were established; the inputs of the BPNN1 model were the input–output indicators (five) of the DEA model, and the output variable was the efficiency score of the DEA model, revealing the production function of the DEA model. When building the BPNN1 model, the samples were randomly partitioned into training and test datasets at a 7:3 ratio, without a validation dataset; the BPNN1 model had a higher prediction accuracy, with only six DMUs beyond 10% error and with a maximum error of 17%. However, the input and output variables of the BPNN2 model were the input–output indicators of the DEA-valid DMUs, revealing the frontier function of the DEA model; the BPNN2 split 1419 DMUs into 772 training and 647 test samples, the established BPNN model demonstrated a high prediction accuracy of less than 10% error for 95% of the DMUs, but Ref. [23] had no validation samples to monitor the training process in real time, and it is impossible to determine whether over-training occurred or not during the training process. Similarly, Kwon (2014) also established two types of BPNN models that reveal the production function and frontier function of the DEA model based on the input–output indicators of the DEA model and their modeling results. The author split the data into training and test datasets at a ratio of 7:3; 54 samples were used to train the network to learn underlying patterns, and the remaining 24 samples were used to test the network. A total of five BPNN models were established, with network topologies of 8-11-1, 8-6-1, 8-4-1, 8-18-1, and 8-14-1, and connection weights of 111, 61, 41, 181, and 141, respectively; except for the third model, the numbers of network connection weights were larger than the number of training samples, and so they do not meet the requirements for building a BPNN model. A further review of the research on combining the DEA model with neural networks such as BPNN, RBFNN, PNN (probabilistic neural network), and SOM to establish DEA-efficient frontier functions or DEA production functions for efficiency prediction, efficiency classification, the screening of training data, and data processing refers to Kwon [16]. The combined DEA models with ANN are also used for risk management, feature selection, etc. In recent years, ML, such as SVR/SVM, RF, bagging, RSM, lightGBM, etc., has been introduced in combination with various DEA models, such as SBM, SESBM with undesirable output, and window analysis [4,11,22,31]. After removing outliers, Zhong et al. [13] obtained 710 sample data points from 95 rural commercial banks from 2011 to 2018 in Guangdong Province, China. The input indicators of the SESBM model were the number of employees, fixed assets, and intermediate business expenses, and the output indicator was the intermediate business income. The dataset was randomly split into training and test subsets with a proportion of 3:1. The authors used SESBM to measure the efficiency of DMUs in the training subset and then adjusted all the training datasets to the efficiency frontier based on projection values. The DEA model’s input and output indicators were the machine learning model’s input and output variables, respectively, and 15 machine learning models, such as BPNN and SVR, were established. Four performance metrics were used to evaluate the different ML algorithms, namely the mean square error (MSE), the mean absolute error (MAE), the root mean square error (RMSE), and the coefficient of correlation R2. After comparing 15 MLs, the results show that the top four MLs for efficiency evaluation performance are BPNN, the extra trees regressor (ETR), RF, and the gradient boosting regressor (GBR). SVR is the worst, and BPNN is the best. Zhong et al. [13] also established combined models between SESBM with undesirable output and SVR, BPNN, etc.; the BPNN model is better than the other models.
To sum up, many scholars have conducted in-depth research on DEA combined models and achieved good results, which have played an essential role in management and other fields. However, the following problems still exist:
(1)
The model’s generalization ability and prediction accuracy are difficult to guarantee without a validation dataset. To establish MLs such as BPNN, except for [8], no other authors divided the samples into training, validation, and test datasets with similar properties. No study discusses how, according to the error changes of the validation dataset, to use the early-stopping or regularization methods to prevent over-training, and use the trial-and-error approach to determine a reasonable number of hidden-layer neurons to ensure the generalization ability and practical value [6,39,40,41,42,43]. Although, most studies clearly state (and some studies do not) to randomly divide the samples into training and test datasets (in some studies, the validation dataset is a test dataset) with similar properties according to the ratios of 8:2, 7:3, or 3:1 (or according to some algorithm). There are also studies stating that over-training should be avoided. Still, since no validation dataset is used to monitor the training process, it is impossible to judge whether over-training has occurred. When the number of training samples is less than three times the number of connection weights due to the extensive network topology, over-training is easy to occur during training. In the case of over-training, even if the error of the training dataset is minimal, and even if the RMSE of the test dataset is casually small, the established model has no generalization ability and practical value. Scholars should pay more attention to this problem.
(2)
It is difficult to establish reliable and effective DEA-ML combined models for DEA modeling problems, usually with only small and medium-sized samples or frontier modeling with only a few DEA-valid samples. Most DEA efficiency modeling uses small or medium-sized samples. For modeling DEA-efficient frontier functions, small samples are usually used. To establish the DEA-BPNN combined model, most of the literature does not meet the basic requirement—the number of training samples must be more than the number of connection weights. Ref. [42] put forward a rule of thumb: you should aim to have at least five times as many cases (training samples) as connection weights in the network, and preferably ten times as many for establishing a reliable and effective BPNN model. According to the rule of thumb, you can determine the reasonable number of neurons in the hidden layer through trial-and-error. To establish DEA combined with SVR, RF, etc., you should avoid over-training and overfitting to carefully determine the model’s parameters with small and medium-sized samples.
(3)
It is not easy to judge which ML model is better than others. Among the ML models currently used for DEA combined modeling, some studies believe that BPNN performs better [10,12,30], while others believe that SVR and other models are better [13]. Which model has better performance is also a question worthy of study.
It can be seen that there are still several tricky problems that need to be improved in the establishment of reliable and effective DEA-combined MLs with a good generalization ability and production precision, particularly for frontier surface functions in the case of small and medium-sized samples.

2.2. Projection Pursuit Regression (PPR) Model

PPR is consistent with the BPNN model regarding the nonlinear approximation ability, especially for the nonlinear modeling of small and medium-sized samples that do not obey normal distribution [25,26,27,28,44,45]. The constraint condition of the PPR model is that the sum of the squares of the weights of all independent variables equals one. The modeling usually starts from linear to quadratic; cubic; polynomial ridge function (PRF); and one to two, three, etc., PRFs. So, there is a low possibility of over-training or overfitting during modeling [26,27,44,46]. It has been widely used in small and medium-sized samples such as agricultural engineering, water conservancy, demography, and earthquakes [26,44,46]. However, there are no studies on DEA-PPR combined modeling. Therefore, in this paper, the PPR model is introduced into efficiency evaluation research for the first time to construct the DEA-PPR combined model, hoping to improve the generalization ability, robustness, and reliability, mainly to solve the efficiency production function and frontier function of DEA combined modeling with small and medium-sized samples, and expand the DEA combined model and the new method of input–output efficiency research, as well as the application field of the PPR model.
The following are the innovations and contributions of this paper:
(1)
We are the first to establish the DEA-PPR combined model to effectively and reliably solve the problem of input–output efficiency of small and medium-sized samples that do not obey the normal distribution, to overcome the disadvantages of MLs, such as BPNN, SVR, etc., that can only be applied to large samples.
(2)
We proposed the modeling principles and steps for establishing a BPNN model with good generalization ability. Through empirical research, the reliable and effective DEA-BPNN, DEA-SVR, and DEA-PPR combined models were found to verify the effectiveness of each other. The DEA-PPR model has a relatively better generalization ability and prediction accuracy among them.
(3)
We established the DEA-PPR combined model to simulate the production function by setting the efficiency score to one, adopting the optimization technology to obtain the frontier function directly, and realizing the unification of the production and frontier surface functions.
(4)
We established the DEA-PPR combined model of the frontier function of the DEA by generating virtual samples according to the valid DEA samples. According to the input-oriented DEA-PPR combined model, the optimal input quantity can be obtained; on the contrary, if the output-oriented DEA-PPR combined model is used, the optimal output quantity can be obtained, providing a decision-making basis and technical paths for DMUs to organize production, strengthen management, improve efficiency, and reduce costs.

3. Methodology and Principle

In this paper, for the prediction research of input–output efficiency, we establish combined models of DEA with BPNN, PPR, and SVR to improve the generalization ability, robustness, and practical value.
This paper briefly introduces some typical models to evaluate input–output efficiency.

3.1. DEA-CCR Model

Charnes et al. [1] first proposed the DEA-CCR model. There are n DMUs, and each DMU has m inputs and q outputs. We set x i and y r to represent the inputs and outputs of the DMU, n × m is the input matrix, and n × q s is the output matrix. For each DMU, we can obtain the ratio r = 1 q u r y r k i = 1 m v i x i k of all its outputs to inputs (where u r is the output weight and v i the input weight); then, the problem of the CCR model is converted to the problem of selecting the best weights. We take the input orientation as an example and obtain the specific planning equation, as in Equation (1):
max θ k = r = 1 q u r y r k i = 1 m v i x i k s . t . r = 1 q u r y r j i = 1 m v i x i j 1   u 0 , v 0   i = 1,2 , . . . , m ; r = 1,2 , . . . , q ; k , j = 1,2 , . . . , n
The goal of the CCR model is to maximize the efficiency value of the rth DMU under the restriction that the other DMUs’ efficiency score is less than or equal to 1. We can solve the problem by transforming it into linear programming. Thus, by constructing an effective frontier, all DMUs either fall within the effective frontier (CCR-effective, CCR-efficient, efficient, and valid) or outside the effective frontier (CCR-ineffective, CCR-efficient, invalid, and inefficient).
In practice, we would like to replace Equation (1) with the corresponding linear programming pairwise model and obtain Equation (2).
min θ s . t . j = 1 n λ j x i j θ x i k j = 1 n λ j y r j = y r k i = 1,2 , . . . , m ; r = 1,2 , . . . , q ; k , j = 1,2 , . . . , n ; λ 0
In Equation (2), λ is the coefficient of the DMU, and the optimal solution (θ) of Equation (2) represents the efficiency score. The range of the θ is (0,1]. For a specific DMU, the DMU is efficient when and only when θ = 1 , and if θ < 1 , the DMU is inefficient.
The CCR model is developed based on the constant return-to-scale (CRS) assumption. The DEA-BCC (or BCC) model, on the other hand, is based on a variable return-to-scale (VRS). The BCC and CCR models differ only in that the former, but not the latter, includes the convexity condition of j = 1 n λ j = 1 , λ j 0 , j in its constraints. Thus, as might be expected, they share properties in common and exhibit differences. They also share properties with the corresponding additive models [2,3,47].
Using the CCR and BCC models helps to determine the DMU’s overall technical and scale efficiencies, and whether the data exhibit a VRS.
Furthermore, the CCR and BCC models fail to distinguish the DMUs with the highest efficiency scores of “1”; therefore, to overcome this problem, the super-efficiency CCR or BCC models were proposed in Ref. [2].
The CCR, BCC, and super-efficiency models are radial models and cannot fully consider the effect of slackness on efficiency. The slack-based-measured (SBM) DEA model (or SBM) and the super-efficiency SBM (called SE-SBM) were proposed to solve this problem [2]. Based on the slackness measurement, SBM is a non-radial method suitable for measuring efficiency when the input and output vary non-proportionally. SE-SBM is a model that combines super efficiency and the SBM.
In the super-efficiency evaluation method, the efficient DMUs are removed from the set, and the efficiency of the DMUs is re-evaluated; thus, the original non-efficient DMUs remain unchanged, and the original efficient DMUs can be greater than 1, then they can be compared [2,13].
According to Cooper et al. [3], the DEA models have a rule of thumb—the minimum number of DMUs should be equal to or greater than max q × m , 3 × q + m .

3.2. Machine Learning

Machine learning (ML) is a multi-disciplinary subject involving many disciplines, such as probability theory, statistics, approximation theory, convex analysis, and algorithm complexity theory. After decades of continuous development, ML is a well-known method that “uses algorithms to parse data, learn from it, and then make decisions or predictions about something unknown in the world” [10,40,42]. An important limitation of ML, such as BPNN, RF, SVR, etc., is that of over-training or overfitting. An ML trained to generalize well on new data (called the generalization ability) will produce a correct input–output mapping, even when the input differs slightly from the examples used to prepare the ML. However, when an ML model learns too many input–output examples, the ML model may memorize training examples that are not true of the underlying function that is to be modeled. Such a phenomenon is referred to as over-training or overfitting. During the training process or optimizing the parameters, the error of the training subset continuously decreases, even almost to zero. The error of the validation subset first decreases and then rises again. This is a sure sign that over-training or overfitting is occurring, and you should stop training once deterioration in the error of the validation subset is observed [40,42]. When an ML model is over-trained, it loses the ability to generalize between similar input–output patterns [6,40,42]. The most effective approach to ensure the generalization ability is to reserve data to cross-validate the ML’s performance. First, the available dataset is randomly partitioned into training, validation, and test subsets. Second, the training subset is used to train the ML model, the validation subset is used to validate the ML model, and the generalization performance of the established ML model is measured on the test subset. Third, if the validation and test errors are close, the trained ML model will likely generalize well; otherwise, we must re-train the ML model [40,42]. If there is no validation or test subset, judging whether over-training has occurred is challenging. BPNN and SVR are the most widely used MLs, and the PPR is not used to combine with the DEA model. We briefly introduce the principles of the BPNN and PPR models as follows.

3.2.1. BPNN Model

The BPNN model, proposed by Rumelhart et al. [6,39,40,42,43,46], is currently the most widely used. The BPNN model consists of an input layer, a hidden layer (usually 1, but 2–3 layers are optional), and an output layer. Each layer consists of multiple neurons. It simulates the human processing process of external input information.
Corresponding to the BPNN with one hidden layer, the output of the jth output node of the pth training sample is
O p j = f 2 h = 1 H ω j h f 1 i = 1 n ω h i x p i + θ h + θ j
In Equation (3), ω h i and ω j h are the connection weights between the input layer and the hidden layer and between the hidden layer and the output layer; θ j and θ h are the thresholds of the hidden layer and the input layer; and f 2 and f 1 are the transfer functions of the output and hidden layers, respectively. n and H are the number of neurons in the input and hidden layers, respectively.
According to the principle of least squares, the objective function is obtained as
Q ω = m i n p = 1 n E p 2 = m i n p = 1 n j = 1 m t p j O p j 2
In Equation (4), t p j is the expected output value; m is the number of neurons in the output layer.
According to the existence theorem proposed by Hornik et al. [48], as long as sufficient neurons are in the hidden layer, a BPNN can always make the error of the training samples as small as possible (even close to 0). However, if the hidden neurons are too many, the BPNN will likely remember the pattern of the training samples and fail to generalize the training samples. Therefore, the network topology must be as compact as possible (the number of hidden layers and neurons should be as small as possible). Second, the collected data must be split into training, validation, and test subsets with similar properties. The validation subset is used to monitor the training process. At the beginning of the training process, the error of training and validation subsets monotonically decrease. The error of the validation subset reduces to a minimum and then rises as the training continues, indicating that over-training occurred. The connection weights before over-training are taken. Refs. [39,42,43,46] proposed the basic principles and steps to be followed when building a BPNN model.

3.2.2. PPR Model

The PPR model was proposed by Friedman and Stuetzle [25]. Because its constraint is that the sum of the squares of the weights of each independent variable is equal to 1, over-training is less likely to occur. Therefore, the PPR model has better reliability and robustness and is especially suitable for small and medium-sized (SM) samples that do not obey the normal distribution. In addition, in comparison with other techniques such as Jackknife, Bootstrap, Monte Carlo, Lasso, and Robust Regression, the PPR model does not sample from the total dataset and can reveal the natural structure characteristics of the data; of course, the modeling process is also reliable and unbiased [46,49,50]. The other techniques do sample to model and some bias exists. Furthermore, the distribution of the efficiency from DEA and the economies is very non-normal. Methods such as Jackknife, Lasso et al. are mainly suitable for normal distribution data.
Assuming the normalized data of the independent variable x ( i , j ) , we can obtain the sample projection value z i = j = 1 m a ( j ) x ( i , j ) and establish the PPR model based on the cubic polynomial ridge function (PRF) of the sample projection value, as follows [25,26,46]:
f z i = c 0 + c 1 z i + c 2 z i 2 + c 3 z i 3   = c 0 + c 1 j = 1 m a ( j ) x ( i , j ) + c 2 j = 1 m a ( j ) x ( i , j ) 2 + c 3 j = 1 m a ( j ) x ( i , j ) 3 .
In Equation (5), c 0 ~ c 3 is the coefficient of the cubic PRF, and m is the number of independent variables. We use the least squares method and take the minimum of the sum of squares fitting errors as the objective function:
Q a , C = m i n i = 1 n y i f z i 2
Based on the parasitism–predation algorithm [51], we obtain the optimal global solution of Equation (6), the coefficient ( c 0 ~ c 3 ) of the cubic PRF of Equation (5), and the optimal projection vector coefficient a ( 1 ) , a ( 2 ) , , a ( m ) as well.
Due to the length limitation of this article, more detailed information on MLs can be found in Refs. [52,53].

4. Empirical Researches

Whether the number of DMUs is small or large, the number of DEA-efficient DMUs is small. To establish the combined model of the DEA-efficient frontier function, we must generate sufficient virtual DEA-efficient samples.

4.1. Empirical Illustration Using Hospital Data

We use the data for 12 hospitals as an example from Cooper et al. [3] that depict the CCR efficiencies in Refs. [3,6]. The samples’ data are shown in Table 2.
The primary purposes of establishing the two DEA-PPR combined models are as follows:
(1)
To analyze the relationship between the efficiency and the DEA input–output indicators, judge the importance of the inputs and outputs, and predict the efficiency of new hospital data;
(2)
To study how to improve the hospitals’ managerial efficiency, and to provide the lower bounds for the inputs (the numbers of doctors and nurses) for each inefficient DMU to produce or service its current level of outputs (the numbers of inpatients and outpatients).
According to Refs. [3,6], we also run a input-oriented CCR under CRS, applying the MaxDEA 9.0 software [2]. Table 2 shows the DEA efficiencies of the twelve DMUs. The efficiency scores of the three DEA-efficient DMUs (A, B, and D) equal “1”.
As aforementioned, we established the DEA-CCR model and obtained the efficiency of each DMU, determined whether one DMU is efficient or inefficient, constructed the efficient frontier of DEA, calculated the inputs’ excesses of each decision unit, etc. However, the optimal solution is not single for the DEA model, is used to establish discontinuously changing frontiers and production functions (models for calculating input–output efficiency), and cannot predict the efficiency values of other (new) DMUs and their optimal input quantities (for input-oriented and specific outputs, the output-oriented is the optimal input). So, we should establish combined models such as DEA-PPR to overcome the disadvantages of DEA. This study will select the DEA-ANN, DEA-PPR, and DEA-SVR combined models and compare their performance.

4.1.1. To Establish the DEA-PPR Combined Model of the DEA Production Function

According to the CCR model mentioned above, we obtain the efficiencies of 12 DMUs. Establishing a robust and reliable combined model is challenging due to the small samples. Therefore, we establish a PPR model with PRF by taking the efficiencies of 12 DMUs as the values of the output variable of PPR and the input–output indicators as the input variables of the PPR model. Because the CCR model describes a complex input–output relationship, the PPR model with a linear PRF cannot satisfy the accuracy. In this paper, we tried to establish a PPR model with a quadratic PRF. We input the above data into our compiled program based on the PPA [26,51] and obtain the optimal global solution, the first-dimension PPR model, and the DMUs’ efficiency scores (ES),
E S P P R - 1 = 1.00007 + 0.06054 z 1 0.36947 z 1 2
where, z 1 = 0.05135 I 1 0.71259 I 2 + 0.43926 O 1 + 0.54464 O 2 . The ES is shown in the E S P P R - 1 column of Table 2. Since the coefficient of the quadratic term of the Formula (7) is less than 0, Formula (7) has a maximum value. When z 1 = 0.081928, the maximum value of E S P P R - 1 is 1.00751.
It can be seen that, for the PPR model of the input-oriented DEA model, the two outputs are positive to the efficiency, and the two inputs are negative to the efficiency. If we take appropriate measures to reduce the input values of one DMU (or increase the output values), the z 1 will increase, and the ES rise; otherwise, the ES will decrease. The relationship between the ES and the input variables is entirely consistent with the production function of the DEA model. It can be seen from z 1 that the impact of I2 (nurse) on ES is more significant than that of I1 (doctor). Therefore, reducing the number of nurses is more conducive to improving the ES of these hospitals, while the two output indicators are not much different.
According to the PPR model (7), the mean absolute error (MAE) is calculated as MAE = 0.018, the mean absolute percentage error as MAPE = 2.2%, and the root mean squared error as RMSE = 0.024. That is to say, Equation (7) has revealed the production function of the DEA model well, but significant errors still exist in very few DMUs. For example, the ES of DEA-inefficient DMU K (1.002) is greater than DEA-efficient DMU B, which cannot accurately reflect the ES of the CCR. Then, we further build a PPR model with the second quadratic PRF to describe the ES:
E S P P R - 2 = 0.00310 + 0.07887 z 2 + 0.03985 z 2 2
where z 2 = 0.87938 I 1 + 0.31175 I 2 + 0.31072 O 1 + 0.18248 O 2 . The ES in Equation (8) is shown in the E S P P R - 2 column of Table 2. Since the coefficient of the quadratic term is greater than 0, Equation (8) has a minimum value. When z 2 = −0.98959, the minimum value E S P P R - 2 is −0.03592. Thus, the PPR model with two quadratic PRFs reveals the production function of the DEA model as follows:
E S P P R = E S P P R - 1 + E S P P R - 2   = 1.00381 + 0.06054 z 1 0.36947 z 1 2 + 0.07887 z 2 + 0.03985 z 2 2
The ES calculated according to Equation (9) is shown in the E S P P R column of Table 2. Therefore, for the DEA-PPR combined model of the original DEA model, if we take appropriate measures to reduce the input (or increase the output), the z 1 or z 2 increase, and the ES will rise; otherwise, the ES will decrease. Because the z 2 of all DMUs are greater than −0.98959, the larger the z 2 is, the higher the ES of the PPR model is. The relationship between the ES of the PPR model and the input variables is precisely the same as the production function of the original DEA model. Meanwhile, according to the combined DEA-PPR model in Equation (9) with two PRFs, the impact of I2 (nurse) on the hospital’s ES is reduced because its coefficient in z 2 is greater than 0. In contrast, the impact of I1 (doctor) is significantly enhanced because the coefficient of z 2 reaches 0.87938, which is significantly larger than that of I1. Therefore, the impacts of the number of nurses and doctors on the hospital’s ES are reduced, and the effects of the two output indicators are not much different, which is consistent with the actual operation of the hospital.
The linear relationship between the ES calculated with the DEA-PPR combined model and the ES of the original DEA model is E S P P R = 0.0303 + 0.968 E S D E A , and the correlation coefficient is R = 0.9737, indicating that the ES values calculated according to the two models are significantly correlated at the 0.01 level. For the DEA-PPR combined model, MAE = 0.014, MAPE = 1.5%, and RMSE = 0.017. The performance metrics, such as MAE, MAPE, etc., are better than that of Equation (7). The ES of DMU K from Equation (9) equals 0.984 and is still more significant than that of benchmark DMU L (0.972), but smaller than that of other DEA-efficient and benchmark DMUs. Therefore, the performance metrics and the distribution of the ES of DMUs show that the established DEA-PPR combined model has a high prediction accuracy, reveals the production function characteristics of the CCR, and suggests good practical value.
We may suppose that, through investigation or other methods, we have obtained the numbers of the outpatients, inpatients, doctors, and nurses of two more hospitals, as shown in Table 2 and Table 3, named M and N. If we rebuild the CCR, the ES of the 14 DMUs is shown in Table 2 (in the ESN column). Except for DMUs A, B, G, I, and L, the other DMUs’ ES has changed. Furthermore, DMU D changes from being DEA-efficient to inefficient. So, we cannot rank the DMUs objectively when adding new DMUs.
According to the above-established DEA-PPR combined model (9), we can quickly obtain the ES of DMU M and N to be 0.991 and 0.249, respectively.
For the above problem, since there are only 12 samples and four input–output indicators, even if we adopt the leave-one-out cross-validation method, the number of neurons in the hidden layer is one (there are seven connection weights); it does not meet the rule of thumb for building a BPNN model—the number of training samples should be more than five times the number of connection weights, so it is impossible to construct a DEA-BPNN, as well as DEA-SVR and DEA-RF models, with good generalization ability.

4.1.2. To Establish the DEA-PPR Combined Model of the DEA-Efficient Frontier Function

Because there is usually a relatively small number of DEA-efficient DMUs, we cannot directly build a DEA-BPNN combined model based on the efficient DMUs only. Ref. [15] constructed 125 virtual efficient DMUs based on three DEA-efficient DMUs, obtained 128 samples, and established a DEA-BPNN combined model. However, since Ref. [5] did not specify whether to divide the samples into training, validation, and test subsets with similar properties, nor how to prevent over-training (or overfitting, over-learning), the generalization ability and robustness of the established DEA-BPNN combined model is questionable and not explicit.
Based on Equation (9), the DEA-PPR combined model of the DEA production function, if we set the ES to be 1, we can obtain the frontier function model of the DEA-efficient DMUs. If I2 (I1), O1, and O2 are assumed, I1 (I2) can be obtained, and the rest can be obtained similarly. Of course, with the virtual samples of Ref. [15], we can also establish the DEA-PPR combined model of the DEA-efficient frontier functions. Referring to [15], we divided the closed interval [0,1] into five equal subintervals, each of the five values (0.2, 0.4, 0.6, 0.8, 1.0) for DEA-efficient A, B, and D, respectively. We generated 5 3 = 125 combinations of DEA-efficient DMUs in all. We have a matrix of the above proportions, P, consisting of 125 rows and three columns, and matrix S of three DEA-efficient DMUs, as follows:
P = 0.2 0.4 0.6 0.8 1.0 0.2 0.4 . . . 1.0 1.0 1.0 1.0 0.2 0.2 0.2 0.2 0.2 0.4 0.4 . . . 1.0 1.0 1.0 1.0 0.2 0.2 0.2 0.2 0.2 0.2 0.2 . . . 0.4 0.6 0.8 1.0 T S = 20 151 100 90 19 131 150 50 27 168 180 72
Thus, we form the matrix V = P 125 × 3 S 3 × 4 and generate 125 ( 5 3 ) virtual linear combinations of the DEA-efficient DMUs. Along with the matrix representing the DEA-efficient DMUs taken singly, we thus have 125 + 3 = 128 linear combinations.
We take the above 128 samples as modeling samples and the nine benchmark DMUs as test samples (or solve the optimal input of these DMUs), input the data into the PPA-based program, and establish the DEA-PPR combined model with one quadratic PRF of the DEA-efficient frontier function. For the input-oriented CCR model, referring to [14], we, respectively, establish, under the given output conditions, the combined DEA-PPR models for the optimal number of doctors and nurses ( D P P R , N P P R ) (represented by the letters D and N, respectively), as follows:
D P P R = 38.6034 + 8.2998 z D 0.00461 z D 2
N P P R = 263.057 + 54.931 z N + 0.0048 z N 2
where z D = 0.85137 O 1 + 0.52457 O 2 , z N = 0.65370 O 1 + 0.75676 O 2 ; O 1 , O 2 has been centered on having a mean zero and scaled to have a standard deviation of 1. From the formulas of z D and z N , it can be seen that the number of outpatients significantly impacts the optimal number of doctors. In contrast, the number of inpatients significantly impacts the optimal number of nurses, but the difference between the two indicators is insignificant.

4.1.3. To Determine the Optimal Number of Doctors and Nurses

According to Equations (10) and (11), the frontier function of the DEA-PPR combined model, we can determine the optimal numbers of doctors and nurses. The optimal numbers of doctors and nurses for nine inefficient DMUs (verification samples) are shown in the D P P R and N P P R columns of Table 3.

4.1.4. To Compare the Performance of the Different Models

Based on the optimal numbers of doctors and nurses obtained from the DEA combined models, the performance metrics of nine verification samples, such as the MAE, MAPE, RMSE, correlation coefficient (R), maximum absolute error ( E A - m a x ), and maximum relative (percentage) error ( E R - m a x ) of the DEA-PPR combined model, were obtained, as shown in the DPPR and NPPR rows of Table 4. The performance metrics of the DEA-BPNN model established in Ref. [15] are shown in the DB and NB rows of Table 4.
From the performance metrics of different models shown in Table 4, it can be seen that the generalization ability of the DEA-PPR combined model with one quadratic PRF of the DEA-efficient frontier function is significantly better than that of the DEA-BPNN combined model in Ref. [15]. For the optimal number of doctors and nurses, the DEA-BPNN combined model has a relatively poor generalization ability, with the maximum relative errors reaching 29.6% and 49.0%, respectively, which shows that over-training occurs when the BPNN model is established (but Ref. [15] does not explicitly state whether they use the validation or test subsets to monitor the training process and how to avoid over-training).

4.1.5. To Analyze the Robustness and Reliability of the DEA-PPR Combined Models

To further verify the robustness and reliability of the PPR model for small samples, we once divided the closed interval [0,1] into three and two equal subintervals, respectively, and generated virtual DEA-efficient samples. The former three values are 0.33, 0.66, and 1, and the latter two values are 0.5 and 1. We can generate 33 = 27 and 23 = 8 virtual samples, a total of 30 and 11 samples. We take the nine DEA-inefficient DMUs as the test subset and establish the DEA-PPR combined model with one quadratic PRF of the DEA-efficient frontier function based on the 27 and 8 virtual samples (respectively, represented with the subscripts 27 and 8), as follows:
D P P R - 27 = 38.666 + 9.0173 z D - 27 0.0196 z D - 27 2
N P P R - 27 = 263.554 + 59.576 z N - 27 + 0.0203 z N - 27 2
where z D - 27 = 0.8515 O 1 + 0.5243 O 2 , z N - 27 = 0.6543 O 1 + 0.7562 O 2 ; O 1 , O 2 has been normalized to have a mean zero and scaled to have a variance of 1, the same as follows:
D P P R - 8 = 36.366 + 9.9389 z D - 8 0.0743 z D - 8 2
N P P R - 8 = 247.515 + 65.574 z N - 8 + 0.0489 z N - 8 2
where z D - 8 = 0.8497 O 1 + 0.5274 O 2 , z N - 8 = 0.6526 O 1 + 0.7577 O 2 .
Since the virtual samples and quantities used to build DEA-PPR combined models differ, the normalized values ( O 1 , O 2 ) are also slightly different.
It can be seen from Table 4 that the performance metrics, such as the MAE, MAPE, and RMSE, of models with different virtual samples have a good consistency. Even if only eight virtual samples are generated, the established DEA-PPR combined models have a high fitting accuracy and generalization ability, indicating that the DEA-PPR combined models have good robustness and reliability.

4.1.6. To Establish the DEA-SVR Combined Model of the DEA-Efficient Frontier Function

Theoretically, an SVR model can minimize the structural risk, so it has a good nonlinear approximation ability, and has also been widely used in nonlinear data modeling, such as constructing DEA frontier functions [11,54,55]. However, in practice, the optimal values of the penalty factor and the width coefficient of the radial basis kernel function must be obtained via the multi-fold cross-validation method [54,55]; otherwise, there is a high possibility of overfitting during optimization. The performance metrics of the SVR model are directly related to whether the values of parameters such as the penalty factor and width coefficient are reasonable, so the model does not have strong stability [11]. Due to the length limitation, we only built the DEA-SVR combined model for the case of 125 virtual samples (taking nine inefficient DMUs and two new DMUs, M and N, to be test samples). If the value ranges of the penalty factor and the width coefficient are both [10−2, 105], we obtained the optimal solutions of the two parameters, 0.019 and 105 (the width coefficient is equal to the boundary value, which is not reasonable). The mean square error of the training samples is already less than 10−7, but the error of the above nine benchmark verification samples is vast; the MAE of the optimal number of nurses is 21.01, the MAPE is 13.91%, the RMSE is 39.07, the E A - m a x is 112.65, and the E R - m a x is 74.63; clearly, overfitting occurred during the optimization process. Suppose the value ranges of the penalty factor and the width coefficient are adjusted to [10−2, 103]. In that case, we obtain the optimal solutions of two parameters, 586.7 and 103 (the width coefficient is still equal to the boundary value, which is not reasonable), and the mean square error of the training sample is 0.36; the MAE, MAPE, RMSE, E A - m a x , E R - m a x , and other data on the optimal number of nurses in the nine benchmark verification samples of the DEA-SVR combined model are shown in the N S V R row of Table 4. Similarly, when the DEA-SVR combined model of the optimal number of doctors is established when the value ranges of the penalty factor and the width coefficient are both [10−2, 105], the optimization process has overfitting; when the value ranges are adjusted to [10−2, 103], the optimal solutions are 484.0 and 103, respectively; the MAE, MAPE, RMSE, E A - m a x , E R - m a x and other data of the established DEA-SVR combined model are shown in the D S V R row of Table 4. It can be seen that the DEA-SVR combined model has poorer performance than the DEA-PPR combined model.
It can be seen from the above modeling process that if the value of the penalty factor is too small, there is a high possibility of overfitting. Unfortunately, there is no commercial software for the SVR model that can, while optimizing the parameters, monitor the error changes of the test samples to prevent overfitting. It brings great difficulties and challenges to establish a DEA-SVR combined model.

4.1.7. To Establish the DEA-BPNN Combined Model of the DEA-Efficient Frontier Function

Although Ref. [15] generated 125 virtual samples and established a DEA-BPNN combined model of the frontier function, the errors of the nine verification samples were large (see Table 4), among which, for the relative errors of the optimal numbers of doctors and nurses in benchmark K, are as high as 22.8% and 32.9%, respectively. The relative error of the optimal number of nurses in benchmark L is also 22.8%, indicating that over-training may have occurred in the training process. So, we try to rebuild the DEA-BPNN combined model. Although many scholars have used the BPNN model to construct the frontier function and production function of the CCR model, few follow the basic principles and steps of establishing a reliable and robust BPNN model. According to the existence theorem proposed by Hornik et al. [48], as long as the hidden layer has sufficient neurons, the error of the training subset of a BPNN can reach as small a value as possible (even close to 0). However, the established BPNN model may have no generalization ability. Therefore, the core of building a BPNN model is to prevent over-training during the training process. We establish a DEA-BPNN combined model with good reliability, robustness, and generalization ability, and follow the basic principles and steps as follows:
(1)
We randomly divide the samples into training, validation, and test subsets with similar properties according to the ratio of 2:1:1. (The number of the validation and test subsets should account for at least 15%, respectively.);
(2)
We use the trial-and-error method and make the BPNN topology as compact as possible (usually one hidden layer and the number of neurons in the hidden layer is as small as possible). The ratio of the number of training subsets to the number of connection weights must be greater than one and should be more than five, preferably ten, according to the rule of thumb;
(3)
We use the training subset to adjust the connection weights to reduce the sum of squares error (SSE) of the training subset and the validation subset to monitor the training process. Along with the training process, the SSE of the training subset gradually decreases, and the SSE of the validation subset first falls to a specific minimal value, and then begins to rise again, which is a sure sign that over-training is occurring. We stop training (called the early-stopping method). To take the network weights before the SSE begins to rise, we establish the BPNN model;
(4)
We use the test subset to measure the prediction ability and performance of the BPNN model. If the SSE of the test subset is reasonably close to or slightly larger (generally less than 1.3 times) than the SSE of the training and validation subsets, the established BPNN model has a good generalization ability, reliability, robustness, and prediction ability, as well as practical value. Otherwise, we should restart the process from (3) until the BPNN model has good generalization and prediction abilities.
The above principles and steps must be wholly followed; otherwise, the generalization ability and practical value of the established BPNN cannot be guaranteed.
We randomly selected 68, 30, and 30 training, validation, and test samples from the above 128 samples. The three subsets had the same properties (the mean and standard deviation were almost the same; if they were inconsistent, we should resample). According to [15], the input layer had two neurons, outpatient and inpatient, and the output layer had one output, the optimal number of nurses or doctors. We established two BPNN models to predict them, respectively. According to the rule of thumb, the number of training samples must reach more than five times the number of connection weights, and the number of hidden-layer neurons cannot exceed three. Using the STATISTICA Neural Networks (SNNs) software [42], we used the logistic transformation function for the hidden layer and output layer; considering that the model must have a certain degree of extrapolation needs, the input and output data were linearly transformed into the range of [0.2, 0.8]. We used the quasi-Newton optimization algorithm and monitored the training process with the validation subset. We compared and studied the performance of the models with 1, 2, and 3 hidden neurons. Under the premise of following the above principles and no over-training, the model with one neuron in the hidden layer had a high enough performance; if the number of hidden-layer neurons was two or three, there is a high possibility of over-training, and the performance was not improved significantly. We established the DEA-BPNN combined models for the optimal numbers of nurses and doctors, respectively, and the performance metrics of nine benchmark verification samples are shown in the N N N and D N N rows of Table 4.
The performance metrics, such as the MAE, MAPE, RMSE, R, E A - m a x , and E R - m a x of the training, validation, and test subsets of the DEA-BPNN combined models for the optimal numbers of doctors and nurses, are shown in Table 5.
It can be seen from Table 5 that the performance metrics of the training, validation, and test subsets are relatively close, indicating that the established BPNN models have a good generalization ability and reliability.
In cases where only 27 and 8 virtual samples are generated, even if we use the leave-one-out cross-validation method, it is impossible to establish DEA-BPNN and DEA-SVR combined models with generalization ability and reliability due to the small samples.

4.1.8. To Establish the Combined Models of the DEA-Efficient Frontier Function for Large Samples

To verify the suitability of the DEA-PPR combined model for large samples, we generated 729 (93) samples taking each of the nine values (0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1) for efficient DMUs A, B, and D, respectively. The minimal values of the samples were equal to that of 125 samples. We randomly divided the 731 samples into 371, 180, and 180 training, validation, and test datasets according to the ratio of 2:1:1. We established the DEA-PPR combined models for doctors and nurses, and the performance metrics of the nine benchmark verification samples are shown in the DPPR-729 and NPPR-729 rows of Table 4. By obeying the principles and steps, we established the DEA-BPNN combined models, and the performance metrics are shown in the DNN-729 and NNN-729 rows of Table 4. The performance metrics of the three subsets are shown in the DNN-729 and NNN-729 rows of Table 5. The optimal number of doctors and nurses is shown in the DNN-729 and NNN-729 columns of Table 3. We set the search range of parameters to be [10−2, 103], optimized the SVR models, and obtained the optimal values to be (7.762, 0.626) and (7.762, 3.206), respectively, which show that overfitting did not occur. The performance metrics of the nine benchmark verification samples are shown in the DSVR-729 and NSVR-729 rows of Table 4.
From the performance metrics of different models in Table 4, it can be seen that, for generating 125 and 729 virtual samples, respectively, the DEA-PPR, DEA-BPNN, and DEA-SVR combined models we established without over-training or overfitting have quite a high prediction accuracy and generalization ability. However, the DEA-PPR combined model slightly outperforms the DEA-BPNN and DEA-SVR for medium-sized samples, and the DEA-PPR combined model outperforms the DEA-BPNN and DEA-SVR for large samples. Furthermore, the DEA-BPNN combined model has a better generalization ability than the DEA-SVR combined models under medium-sized and large samples. More importantly, due to the strongly nonlinear approximation ability of the BPNN and SVR models, there is a high possibility of over-training and overfitting during training or optimizing the parameters within unreasonable search ranges. Cross-validation methods with validation subsets must be used, which require a cumbersome modeling process and cannot guarantee the generalization and prediction abilities of the model without following the principles and steps. Because the PPR model with a quadratic PRF is characterized by a relatively simple structure and moderate nonlinear approximation ability, it is especially suitable for modeling problems that are not very complex, such as the frontier surface, the production function, etc., of the DEA model. Whether for 729, 125, 27, or 8 virtual samples, the established DEA-PPR combined models have quite a good generalization ability and prediction ability. The PPR model is especially suitable for small and medium-sized samples that do not obey normal distribution. Furthermore, the DEA model is also ideal for analysis and modeling to research the input–output efficiency of small and medium-sized samples. Therefore, for the problem of small and medium-sized samples, the DEA-PPR combined model has unique advantages, such as good generalization, robustness, and reliability.

4.2. Empirical Illustration Using China’s Provincial Carbon Dioxide Emission Quotas

As the world’s top CO2 emitter (accounting for 26.5% of the world’s total emissions) and the second largest economy, China’s achievement of peak CO2 emissions and carbon neutrality (called the “dual carbon” target) will play a pivotal role in reducing the world’s CO2 emissions. Therefore, reasonably allocating China’s provincial carbon emission quotas (CPCEQs) plays a vital role in promoting high-quality development and striving to achieve a “dual carbon” target. In this study, firstly, it is necessary to judge and evaluate the advantages and disadvantages of different allocation methods. Secondly, it is required to study the carbon emission efficiency of different periods (e.g., comparing efficiency between 2019 and 2030), and so on. DEA is the most used model to measure the efficiency of carbon emissions. It is a pity we cannot use the DEA model to compare the efficiency change between different periods and the advantages and disadvantages of other allocating methods.
The primary purposes of establishing the two DEA-PPR combined models are as follows:
(1)
To build a relationship between the carbon emissions efficiency and the DEA inputs and outputs, judge the importance of the inputs and outputs, and predict the efficiency of new carbon-allocating methods and the quotas in 2030;
(2)
To study ways to improve China’s provincial carbon emissions efficiency and help to implement the “dual carbon” target, and provide lower bounds for the inputs (the carbon emissions) for each inefficient DMU to produce its current level of outputs (the GDP and population).

4.2.1. Data Resource

Referring to Gomes et al. [56], we take the gross domestic production (GDP) and the population of provinces as the output indicators and the actual carbon emissions as the input indicator. The data of the output indicators are from the China Statistical Yearbook (2020) (http://www.stats.gov.cn/sj/ndsj (accessed on 8 July 2023)), and the input data are from Chinese carbon emission accounts & datasets (https://www.ceads.net (accessed on 8 July 2023)). The data are shown in Table 6.
We establish the input-oriented BCC model to measure the ES of the Chinese provincial carbon emissions in 2019, referring to Refs. [56,57] and the ES shown in Table 6. According to Refs. [56,57], we built the ZSG-DEA to allocate the provincial quotas and obtain the quotas, as shown in the “quotas” column of Table 6. If we rebuild the BCC model using the quotas and make the ES values all equal to 1, we cannot judge whether the ES will rise or decline. Meanwhile, we can allocate the provincial quotas by establishing an allocation system consisting of multiple variables [46] and obtain the quotas denoted in the “Quotas-y” column of Table 6. But how do we judge the ZSG-DEA’s quotas and the quotas of the allocation system? We can obtain the provincial quotas in 2030 according to the GDP and carbon emission intensity in the Action Plan for Carbon Peaking before 2030 (namely “Plan”), as well as the 14th Five-Year Plan for Economic and Social Development and the Outline of the Vision for 2035 (namely “Outline”) of provinces in China. We cannot compare and evaluate the ES of different allocation methods and periods solely through the DEA model because adding more DMUs changes the ES of the original DMUs. So, we must establish combined models to do so.

4.2.2. To Build the Combined Model Characterizing the DEA Production Function

(1)
To build a DEA-PPR combined model
Referring to Section 4.1, we take the ES as the output variable and the three input–output indicators as the input variables, and build the DEA-PPR combined model with one quadratic PRF, and the 30 training and ten verification samples as shown in Table 6 (noted “b”), as follows:
E S P P R - 1 = 0.7099 0.5571 z 1 + 0.1613 z 1 2
where z 1 = 0.8622 I 1 0.4605 O 1 0.2110 O 2 , I 1 denotes CO2 (Mt), O 1 denotes the population, and O 2 denotes the GDP, as shown in Table 6. The performance metrics of the training samples, such as the RMSE, AAE, MAPE, and EA-max, are 0.1129, 0.0834, 17.48%, and 0.3305, respectively. The performance metrics of the verification samples are 0.0737, 0.0615, 6.15%, and 0.1172, respectively. The performance of the DEA-PPR model is not good enough, and we rebuild the DEA-PPR model with the second quadratic PRF as follows:
E S P P R - 2 = 0.0557 + 0.0727 z 2 + 0.0876 z 2 2
where z 2 = 0.7195 I 1 + 0.4261 O 1 0.5484 O 2 . The performance metrics of the training samples, such as the RMSE, AAE, MAPE, and EA-max, are 0.0871, 0.0621, 14.83%, and 0.2448, respectively. The performance metrics of the verification sample are 0.0552, 0.0501, 5.00%, and 0.0819, respectively. Although we can build the third PRF to improve the performance, it is shown that the performance of the DEA-PPR model is good enough to describe the DEA production function. That is to say, the DEA-PPR combined model can approximately characterize the DEA production function as follows:
E S P P R = E S P P R - 1 + E S P P R - 2
The calculated ES of the provincial quotas, quotas-y, and quotas—2030 are all shown in the “ES-q”, “ES-y”, and “ES-2030” columns of Table 6.
According to the allocation system, the average ES of the quotas is greater than that of the actual carbon emissions, and shows that the allocated quotas are reasonable and improve the efficiency of carbon emissions. The average ES in 2030 is less than the actual carbon emissions in 2019 and shows that the carbon emissions efficiency in 2030 is lower than in 2019. To raise the carbon emissions efficiency in 2030, we must take measurements to reduce carbon emissions. If we can reduce carbon emissions by 5% in 2030, the average ES calculated with the DEA-PPR combined model is 0.6446, greater than the average ES in 2019.
The average ES of the quotas from the ZSG-DEA model is lower than the average ES of the actual carbon emissions, and shows that the ZSG-DEA model for carbon emission quotas can equalize the carbon emission efficiency of provinces, which does reduce carbon emission efficiency.
(2)
To build the DEA-SVR and DEA-BPNN combined models
There are too few samples to build DEA-BPNN and DEA-RF combined models. We build an SVR model according to 4.6. We set the parameters’ search range to be [0.01, 10] and obtain the optimal solution; the penalty factor equals 0.0258, and the width coefficient equals 10. The performance metrics of the training samples, such as the RMSE, AAE, and MAPE, are 0.0155, 0.0123, and 3.15, respectively. The performance metrics of the verification sample are 0.3205, 0.3116, and 31.16, respectively. Overfitting occurred during the optimization process.

4.2.3. To Build the Combined Models Characterizing the DEA Frontier Function

We cannot build the combined models of the DEA-efficient frontier function using the five DMUs only. Referring to Section 4.2, we also generate 25 or 35 virtual DEA-efficient DMUs by dividing the value interval [0, 1] into two (0.5, 1) or three (0.3333, 0.6667, 1) equal parts in a linear combination with the data of five DEA-efficient samples. We generated 32 and 243 virtual samples, respectively.
(1)
To build a DEA-PPR combined model
We take the DEA model’s input indicator as the PPR model’s output variable and the two output indicators as the input variables. We take the ten samples in Table 6, denoted as “b”, as the verification samples.
We build the PPR model with the first quadratic PRF as follows:
C P P R - 32 = 677.96 + 272.23 z 32 + 0.9597 z 32 2
C P P R - 243 = 765.99 + 172.26 z 243 + 0.4784 z 243 2
where z 32 = 0.8784 O 1 + 0.4779 O 2 and z 243 = 0.8792 O 1 + 0.4765 O 2 . The performance metrics of the training and validation subsets, such as the RMSE, AAE, MAPE, and EA-max, are shown in Table 7, respectively. The performance metrics show that the DEA-PPR combined model has good performance, generalization ability, and robustness for 32 or 243 virtual samples. The DEA-PPR combined model with one quadratic PRF already has a high enough prediction accuracy.
(2)
To build the DEA-SVR and DEA-BPNN combined models
We build an SVR model according to Section 4.1.6 and Section 4.1.7. When the parameters’ search range is [0.01, 8.7096], the RMSE equals the RMSE of the DEA-PPR combined model. The optimal solution is that the penalty factor equals 0.2041, and the width coefficient equals 8.7096; that is, the optimal width coefficient equals the upper bound value. The performance metrics of the training and validation subsets are shown in Table 7.
It can be seen that the DEA-SVR combined model has almost the same performance for the training subset as DEA-PPR but has inferior performance for the validation subset compared to the DEA-PPR combined model. Overfitting occurred during the optimization process.
Furthermore, the main disadvantage of the SVR model is that its performance is directly dependent on the search range of the parameters. If we change the upper boundary limit, the performance will change, too. If we set the upper boundary limit to 100, the SVR has a poorer generalization ability.
(3)
To build the DEA-BPNN and DEA-RF combined models
The 32 virtual samples are too few to build DEA-BPNN and DEA-RF combined models. The 243 virtual samples are barely suitable for building the DEA-BPNN and DEA-RF combined models.
We build a DEA-BPNN combined model. We randomly split the 248 (243 + 5) samples into training, validation (especially for the BPNN and RF model), and test subsets in the proportion of 6:2:2 [31,42,45] and obtain 148, 50, and 50 training, validation, and test samples. According to the rule of thumb, we should preferably have ten times as many training samples as the number of the connection weights. We determine the number of hidden neurons to be three; the connection weights are 13. The transfer functions in the hidden and output layers are Sigmoid applying the SNN [42]. We train the BPNN, monitor the error change in the validation subset, and build the BPNN model without over-training. The performance metrics of the DEA-BPNN combined model of the training and validation subset are shown in Table 7. The performance metrics of the validation and test subsets are 11.41, 1.61%, 14.64, 47.81, and 111.00, and 11.71, 3.78%, 13.83, 25.79, and 3.85, respectively. The performance metrics show that the built BPNN model has a good generation ability but is still poorer than the DEA-PPR combined model. According to the built model, we can conclude that GDP significantly impacts efficiency more than the population.
We build a DEA-RF combined model. We randomly split the 248 (243 + 5) samples into training and validation subsets in the proportion of 1:1 [46]. We take the default value of the parameters of the DPS software [46], the number of trees being 300, and build the DEA-RF combined model. The performance metrics of the DEA-RF combined model of the training and validation subsets are shown in Table 7. The validation subset’s performance metrics are greater than the training subset’s, indicating that overfitting has occurred in the optimization. The built DEA-RF combined model has a poor generalization ability.
The comparison of the performance metrics of different combined models shows that the DEA-PPR model has the best performance and generalization ability, robustness, and prediction accuracy, whether for large samples or small and medium-sized samples. The DEA-BPNN model without over-training has the second-best performance for medium-sized and large samples. The DEA-SVR combined model has a poorer generalization ability than the DEA-PPR and DEA-BPNN combined models, and overfitting occurs in optimizing the parameters under medium-sized samples. The DEA-RF combined model has the poorest generalization ability of the four combined models for large samples. We cannot build DEA-BPNN, DEA-SVR, and DEA-RF combined models for small samples.

5. Results and Concluding Remarks

5.1. The PPR and DEA Models Have Similarities in Frontier Morphology and Theoretical Consistency

According to the DEA modeling principle, the frontier surface of the input-oriented (or output-oriented) DEA is the top surface of a convex polyhedron that is convex (or concave) to the coordinate origin. The PPR model based on the quadratic PRFs comprises multiple (at least one) quadratic surfaces, which are convex polyhedra. It can be seen that the convex polyhedron of the PPR model and the top surface of the frontier surface of the DEA model have a similar morphology and consistent theory. Therefore, there is a low possibility of overfitting in building the DEA-PPR combined model, which thus has a good generalization ability and robustness. The widely used BPNN, RF, and SVR models with very complex curved surfaces have a low morphological similarity and theoretical consistency with the DEA model.

5.2. The Characterization Ability of the DEA-PPR Combined Model to the DEA Production Function

The main disadvantage of DEA is the static analysis of the existing DMUs. The production function constructed via DEA is discontinuous, so the newly added DMUs cannot be analyzed. If the DEA model is re-established with the newly added DMUs, the efficiency and ranking of the original DMUs will be changed, leading to a lack of comparability. For example, after adding new DMUs (M and N) and rebuilding the DEA model, the ES of the seven original DMUs changed, and DMU D changed from being DEA-efficient to inefficient.
To overcome its disadvantages, scholars have established combined models, such as DEA-BPNN, DEA-SVR, DEA-RF, etc., based on the input–output indicators and the ES of the DEA model that have played a good role. However, most combined models, such as BPNN, RF, etc., are mainly suitable for large samples and not suitable for small and medium-sized samples. Therefore, under the cases of small and medium-sized samples, we are the first to propose an innovative idea of establishing a DEA-PPR combined model and its modeling steps. For the two actual examples, we established reliable and effective combined models, describing the production function and the frontier surface function of the CCR, respectively. Especially in the cases of small and medium-sized samples that do not obey the normal distribution, it is challenging to establish a reliable DEA-BPNN, DEA-SVR, or DEA-RF combined model for the production or frontier surface functions of the CCR with generalization ability and prediction accuracy.
It can be seen from Table 2 and Table 6 that the established DEA-PPR combined model can characterize the production function of the DEA model, has sufficient fitting accuracy, and predicts the ES of new DMUs. We input the data of the new DMUs (M and N) into the DEA-PPR combined model (Equation (10)), and obtained their ES values of 0.991 and 0.249, respectively.
Therefore, the DEA-PPR combined model can exploit the advantages of both nonparametric models and overcome their shortcomings.

5.3. The Characterization Ability of the DEA-PPR Combined Model to DEA-Efficient Frontier Function

We establish the DEA model to obtain the DEA-efficient frontier function consisting of a series of line segments of the polyline (approximately). There are only a few DEA-efficient DMUs, and we must fill the gap between the broken line and the actual situation. Of course, if new DMUs are added, it is impossible to analyze whether it is on the frontier surface or to calculate its optimal input amount. Therefore, according to the DEA-efficient DMUs, some scholars have divided each indicator into equal parts to generate virtual DEA-efficient DMUs. According to different equivalent fractions, we can establish combined models such as DEA-PPR, DEA-BPNN, DEA-SVR, etc., to describe the DEA-efficient frontier function. To establish the BPNN, RF, and SVR combined models, we must randomly divide the samples into training, validation, and test subsets with similar properties. We use cross-validation to train or optimize the parameters to avoid over-training.
We usually establish combined models for large samples, such as DEA-PPR, DEA-BPNN, DEA-SVR, etc., with a good generalization ability; otherwise, it is difficult to guarantee the generalization ability. The DEA-PPR combined model has a better generalization ability than the DEA-BPNN model, the DEA-SVR model, etc.
The PPR model is especially suitable for small and medium-sized samples of nonlinear and non-normal distribution. Our empirical research generates medium-sized and medium-sized virtual DMUs with a linear combination of DEA-efficient DMUs. We established reliable DEA-PPR combined models with a good generalization ability. Differently to the BPNN and SVR combined models, according to the DEA-PPR combined model, we can easily judge the importance and ranking of the indicators according to the coefficients of different inputs (or outputs). The DEA-PPR combined model outperforms the DEA-BPNN and DEA-SVR combined models.
In short, the DEA-PPR combined model outperforms DEA-BPNN, DEA-SVR, etc., for small, medium-sized, and large samples, and should be recommended for use regarding generalization ability, robustness, and accuracy.
According to the DEA-PPR combined model, we can obtain the optimal inputs (outputs) for new DMUs and take measures to improve efficiency and management.

6. Limitations and Future Research

This study mainly focuses on overcoming the weaknesses of the DEA model and the shortcomings of the DEA-BPNN and DEA-SVR combined models not being well applied to small and medium-sized samples. We introduce the PPR model and first establish a DEA-PPR combined model with a good generalization ability. However, there are some limitations. Future research should use more examples to verify the generalization ability and accuracy between DEA-BPNN, DEA-SVR, and DEA-PPR for establishing production and frontier functions. Second, we take the DEA-CCR model as an example, and in the next step, we should study other DEA models to verify the applicability of the DEA-PPR combined model. This paper is the first to focus on the DEA-PPR combined model, which requires more empirical research to confirm its universality and reliability.

Author Contributions

Conceptualization, W.L. and X.Y.; methodology, X.Y.; software, W.L.; validation, X.Y.; formal analysis, W.L.; investigation, X.Y.; resources, X.Y.; data curation, X.Y.; writing—original draft preparation, W.L. and X.Y.; visualization, X.Y.; supervision, W.L.; project administration, X.Y.; funding acquisition, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by Shanghai Business School Plateau Discipline Project of Business Administration, No. SHGDD-GSGL-2023-02-02.

Data Availability Statement

The data are available on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Charnes, A.; Cooper, W.W.; Rhodes, E. Measuring the efficiency of decision making units. Eur. J. Oper. Res. 1978, 2, 429–444. [Google Scholar] [CrossRef]
  2. Cheng, G. Data Envelopment Analysis: Methods and MaxDEA Software; Intellectual Property Press: Beijing, China, 2014; Available online: http://www.maxdea.cn/ (accessed on 12 June 2023).
  3. Cooper, W.; Seiford, L.; Tone, K. Data Envelopment Analysis—A Comprehensive Text with Models, Applications, References and DEA-Solver Software; Klumer Academic Publishers: Boston, MA, USA, 2007. [Google Scholar]
  4. Panwar, A.; Olfati, M.; Pant, M.; Snasel, V. A review on the 40 years of existence of data envelopment analysis models: Historic development and current trends. Arch. Comput. Methods Eng. 2022, 29, 5397–5426. [Google Scholar] [CrossRef] [PubMed]
  5. Wu, D.; Yang, Z.; Liang, L. Using dea-neural network approach to evaluate branch efficiency of a large Canadian bank. Expert Syst. Appl. 2006, 31, 108–115. [Google Scholar] [CrossRef]
  6. Bishop, C. Neural Networks for Pattern Recognition; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  7. Lewis, H.F.; Sexton, T.R. Network DEA: Efficiency analysis of organizations with complex internal structure. Comput. Oper. Res. 2004, 31, 1365–1410. [Google Scholar] [CrossRef]
  8. Athanassopoulos, A.D.; Curram, S.P. A comparison of data envelopment analysis and artificial neural networks as tool for assessing the efficiency of decision making units. J. Oper. Res. Soc. 1996, 47, 1000–1016. [Google Scholar] [CrossRef]
  9. Na, J.; Zhou, Z.; Zhou, H. A decision-analysis approach to determine the total staff employ in government organization—Establishment of local government organizations. Chin. J. Manag. Sci. 1997, 5, 7–17. [Google Scholar]
  10. Ma, C.; Wu, C.; Zhang, S.; Yang, Y.; Li, H.; Han, Y. Decision making method for variable-rate fertilization based on data envelopment analysis and artificial neural network. Trans. CSAE 2004, 20, 152–155. [Google Scholar]
  11. Zhu, N.; Zhu, C.; Emrouznejad, A. A combined machine learning algorithms and DEA method for measuring and predicting the efficiency of Chinese manufacturing listed companies. J. Manag. Sci. Eng. 2021, 6, 435–448. [Google Scholar] [CrossRef]
  12. Olanrewaju, O. Integrated index decomposition analysis-artificial neural network-data envelopment analysis (IDA-ANN-DEA)—Implementation guide. Energy Effic. 2021, 14, 71–78. [Google Scholar] [CrossRef]
  13. Zhong, K.; Wang, Y.; Pei, J. Super efficiency SBM-DEA and neural network for performance evaluation. Inf. Process. Manag. 2021, 58, 102728. [Google Scholar] [CrossRef]
  14. Anouze, A.L.M.; Bou-Hamad, I. Data envelopment analysis and data mining to efficiency estimation and evaluation. Int. J. Islam. Middle East. Financ. Manag. 2019, 12, 169–190. [Google Scholar] [CrossRef]
  15. Bose, A.; Patel, G.N. “NeuralDEA”—A framework using neural network to re-evaluate DEA benchmarks. OPSearch 2015, 52, 18–41. [Google Scholar] [CrossRef]
  16. Kwon, H. Exploring the predictive potential of artificial neural networks in conjunction with DEA in railroad performance modeling. Int. J. Prod. Econ. 2017, 183, 159–170. [Google Scholar] [CrossRef]
  17. Hong, H.K.; Ha, S.H.; Shin, C.K.; Park, S.C.; Kim, S.H. Evaluating the efficiency of system integration projects using data envelopment analysis (DEA) and machine learning. Expert Syst. Appl. 1999, 16, 283–296. [Google Scholar] [CrossRef]
  18. Tsaples, G.; Papathanasiou, J.; Georgiou, A.C. An exploratory DEA and machine learning framework for the evaluation and analysis of sustainability composite indicators in the EU. Mathematics 2022, 10, 2277. [Google Scholar] [CrossRef]
  19. Mirmozaffari, M.; Shadkam, E.; Khalili, S.M.; Kabirifar, K.; Yazdani, R.; Gashteroodkhani, T.A. A novel artificial intelligent approach: Comparison of machine learning tools and algorithms based on optimization DEA Malmquist productivity index for eco-efficiency evaluation. Int. J. Energy Sect. Manag. 2021, 25, 523–550. [Google Scholar] [CrossRef]
  20. Yang, X.; Dimitrov, S. Data envelopment analysis may obfuscate corporate financial data: Using support vector machine and data envelopment analysis to predict corporate failure for nonmanufacturing firms. INFOR Inf. Syst. Oper. Res. 2017, 55, 295–311. [Google Scholar] [CrossRef]
  21. Emrouznejad, A.; Shale, E. A combined neural network and DEA for measuring efficiency of large scale datasets. Comput. Ind. Eng. 2009, 56, 249–254. [Google Scholar] [CrossRef]
  22. Saeidi, S.; Jouybanpour, P.; Mirvakilli, A.; Iranshahi, D.; Klemeš, J.J. A comparative study between modified data envelopment analysis and response surface methodology for optimisation of heterogeneous biodiesel production from waste cooking palm oil. J. Clean. Prod. 2016, 136, 23–30. [Google Scholar] [CrossRef]
  23. Kwon, H.; Lee, J.; Roh, J. Best performance modeling using complementary DEA-ANN approach—Application to Japanese electronics manufacturing firms. Benchmarking Int. J. 2016, 23, 704–721. [Google Scholar] [CrossRef]
  24. Farahmand, M.; Desa, M.; Nilashi, M. A combined data envelopment analysis and support vector regression for efficiency evaluation of large decision making units. Int. J. Eng. Technol. 2014, 6, 2310–2321. [Google Scholar]
  25. Friedman, J.; Stuetzle, W. Projection pursuit regression. J. Am. Stat. Assoc. 1981, 76, 817–823. [Google Scholar] [CrossRef]
  26. Lou, W. The Projection Pursuit Theory Based on Swarm Intelligence Optimization Algorithms—New Developments, Applications, and Software; Fudan University Press: Shanghai, China, 2021. [Google Scholar]
  27. Hwang, T.; Lay, S.; Maechler, M. Regression modeling in back-propagation and projection pursuit learning. IEEE Trans. Neural Netw. 1994, 5, 342–353. [Google Scholar] [CrossRef] [PubMed]
  28. Zhan, H.R.; Zhang, M.K.; Xia, Y.C. Ensemble projection pursuit for general nonparametric. arXiv 2022, arXiv:2210.14467. [Google Scholar]
  29. Ren, H.; Ma, X.R.; Li, H.B. Improvement of input evaluation for giant projects based on GA-BP neural network. Syst. Eng.—Theory Pract. 2015, 35, 1474–1481. [Google Scholar]
  30. Tsolas, I.E.; Charles, V.; Gherman, T. Supporting better practice benchmarking: A DEA-ANN approach to bank branch performance assessment. Expert Syst. Appl. 2020, 160, 113599. [Google Scholar] [CrossRef]
  31. Zhang, Z.; Xiao, Y.; Niu, H. DEA and Machine Learning for Performance Prediction. Mathematics 2022, 10, 1776. [Google Scholar] [CrossRef]
  32. Fallahpour, A.; Kazemi, N.; Molani, M.; Nayyeri, S.; Ehsani, M. An Intelligence-Based Model for Supplier Selection Integrating Data Envelopment Analysis and Support Vector Machine. Iran. J. Manag. Stud. 2018, 11, 209–241. [Google Scholar]
  33. Yazdanparast, R.; Tavakkoli-Moghaddam, R.; Heidari, R.; Aliabadi, L. A hybrid Z-number data envelopment analysis and neural network for assessment of supply chain resilience: A case study. Cent. Eur. J. Oper. Res. 2021, 29, 611–631. [Google Scholar] [CrossRef]
  34. Sreekumar, S.; Mahapatra, S. Performance modeling of Indian business schools: A DEA-neural network approach. Benchmarking Int. J. 2011, 18, 221–239. [Google Scholar] [CrossRef]
  35. Kao, H.; Huang, C.; Chen, J. Classification using DEA and SVM approaches: The empirical study of higher education. Information 2013, 16, 7801–7810. [Google Scholar]
  36. Barros, C.; Wanke, P. Insurance companies in Mozambique: A two-stage DEA and neural networks on efficiency and capacity slacks. Appl. Econ. 2014, 46, 3591–3600. [Google Scholar]
  37. Sanei, R.; Hosseinzadeh lotfi, F.; Fallah, M.; Sobhani, F. An estimation of an acceptable efficiency frontier having an optimum resource management approach, with a combination of the DEA-ANN-GA technique (A case study of branches of an insurance company). Mathematics 2022, 10, 4503. [Google Scholar] [CrossRef]
  38. Liu, Q.; Shang, J.; Wang, J.; Niu, W.; Qiao, W. Evaluation and prediction of the safety management efficiency of coal enterprises based on a DEA-BP neural network. Resour. Policy 2023, 83, 103611. [Google Scholar] [CrossRef]
  39. Lou, W.G. Evaluation and prediction of soil quality based on artificial neural network in the Sanjiang Plain. Chin. J. Manag. Sci. 2002, 10, 79–83. [Google Scholar]
  40. Haykin, S. Neural Networks and Learning Machines; China Machine Press: Beijing, China, 2009. [Google Scholar]
  41. Zhang, G.; Patuwo, E.; Hu, M. Forecasting with artificial neural networks: The state of the art. Int. J. Forecast. 1998, 14, 35–62. [Google Scholar] [CrossRef]
  42. StatSoft, Inc. STATISTICA Neural Networks; StatSoft, Inc.: Tulsa, OK, USA, 2011. [Google Scholar]
  43. Zhong, C.; Lou, W.G.; Wang, C. Neural Network-Based Modeling for Risk Evaluation and Early Warning for Large-Scale Sports Events. Mathematics 2022, 10, 3228. [Google Scholar] [CrossRef]
  44. Singhee, A. SiLVR: Projection Pursuit for Response Surface Modeling. In Machine Learning in VLSI Computer-Aided Design; Elfadel, I., Boning, D., Li, X., Eds.; Springer Nature: Cham, Switzerland, 2019. [Google Scholar]
  45. Hall, P. On projection pursuit regression. Ann. Stat. 1989, 17, 573–588. [Google Scholar] [CrossRef]
  46. Yu, X.H.; Xu, H.Y.; Lou, W.G.; Xu, X.; Shi, V. Examining energy eco-efficiency in China’s logistics industry. Int. J. Prod. Econ. 2023, 258, 108797. [Google Scholar] [CrossRef]
  47. Banker, R.D.; Charnes, A.; Cooper, W.W. Some models for estimating technical and scale inefficiencies in data envelopment analysis. Manag. Sci. 1984, 30, 1078–1092. [Google Scholar] [CrossRef]
  48. Hornik, K.; Stinchcombe, M.; White, H. Multilayer feedforward networks are universal approximators. Neural Netw. 1989, 2, 359–366. [Google Scholar] [CrossRef]
  49. Chen, C.; Tuo, R. Projection Pursuit Gaussian Process Regression. IISE Trans. 2023, 55, 901–911. [Google Scholar] [CrossRef]
  50. Tao, H.; Tao, J.; Li, Q.; Aihemaiti, M.; Jiang, Y.; Yang, W.; Wei, J. Average relative flow of single-wing labyrinth drip irrigation tape based on projection pursuit regression. Sci. Rep. 2022, 12, 8543. [Google Scholar] [CrossRef] [PubMed]
  51. Mohamed, A.A.A.; Hassan, S.A.; Hemeida, A.M.; Alkhalaf, S.; Mahmoud, M.; Eldin, A.M.B. Parasitism–Predation algorithm (PPA): A novel approach for feature selection. Ain Shams Eng. J. 2020, 11, 293–308. [Google Scholar] [CrossRef]
  52. Athey, S.; Imbens, G. Machine learning methods that economists should know about. Annu. Rev. Econ. 2019, 11, 685–725. [Google Scholar] [CrossRef]
  53. Marsland, S. Machine Learning: An Algorithmic Perspective; CRC Press: Boca Raton, FL, USA; Taylor & Francis Group: Abingdon, UK, 2015. [Google Scholar]
  54. Valero-Carreras, D.; Aparicio, J.; Guerrero, N.G. Support vector frontiers: A new approach for estimating production functions through support vector machines. Omega 2021, 104, 102490. [Google Scholar] [CrossRef]
  55. Gomes, E.G.; Lins, M.P.E. Modelling undesirable outputs with zero sum gains data envelopment analysis models. J. Oper. Res. Soc. 2008, 59, 616–623. [Google Scholar] [CrossRef]
  56. Zhou, X.; Niu, A.Y.; Lin, C.X. Optimizing carbon emission forecast for modelling China’s 2030 provincial carbon emission quota allocation. J. Environ. Manag. 2023, 325, 116523. [Google Scholar] [CrossRef]
  57. Tang, Q.; Zhang, C. Data processing system (DPS) software with experimental design, statistical analysis and data mining developed for use in entomological research. Insect Sci. 2013, 20, 254–260. [Google Scholar] [CrossRef]
Table 1. The detailed information on the articles establishing a combined DEA model with machine learning such as BPNN, RF, SVR, etc.
Table 1. The detailed information on the articles establishing a combined DEA model with machine learning such as BPNN, RF, SVR, etc.
ReferencesDEA ModelTopology Employed in a Combined ModelCombined ModelNumber of Samples (Validation, Test)To Obey the Rule of ThumbModel for Efficiency Score (ES) and Efficient Frontier (EF)
Wu et al. [5]CCR (2, 3)5-10-1, 5-4-1two BPNNs142 (84)noES
Athanassopoulos et al. [8]CCR (2, 1) *, DEA (4, 3)2-3-1 **, 4-10-3BPNN250 (50)yesES, DEA > BPNN; efficiency rank, BPNN > DEA
Na et al. [9]CCR (5, 2)5-3-2BPNN13N, NnoES
Ma et al. [10]CCR (6, 1)4-4-3BPNN38 (5)noES
Zhu et al. [11]CCR (2, 2)4-?-1BPNN, GANN, SVR, ISVR948 (48)yesES, GANN > BPNN > ISVM > SVM
Olanrewaju [12]CCR (1, 1)1-5-1BPNN8noES
Zhong et al. [13]SESBM (3, 1)3-?-115 MLs such as CART, CIT, Bagging, RF, BPNN, etc.710 (3:1)noEF, BPNN > ET > RF > GBR
Anouze et al. [14]CCR (5, 4), and 15 environmental variables9-5-115 MLs such as CART, CIT, Bagging, RF, BPNN, etc.151 (2:1)noThe bagging and RF are better than BPNN, CIT, etc.
Bose et al. [15]CCR (2, 2)2-3-2-2, 5-4-3-2Two BPNNs12, 99noEF
Kwon [16]two DEAs, CCRs (4, 1)5-7-1, 5-8-1, 5-8-1, 5-7-1; 4-3-1, 4-2-1, 4-4-1, 4-5-1eight BPNNs56 (17)noES, EF
Hong et al. [17]CCR (4,4)/SOMc50//
Yang et al. [20]SBM (5, 5)trial-and-errorSVMccreating at most 500 instances, 10-fold CV/ES
Saeidi et al. [22]CCR (4, 2)6-?-1BPNN26noES and EF
Kwon et al. [23]CCR (3, 3), CCR (3, 1)15-9-1, 6-30-1, 15-3-1, 15-3-1four BPNNs181 (37, 36) #noES, EF
Ren et al. [29]CCR (5, 5)10-21-1BPNN5N, NnoES
Tsolas et al. [30]eight DEAs, CCR (3, 2)3-3-1, 6-?-1two BPNNs160 (4:1)noES, EF
Zhang et al. [31]SESBM (3, 2)5-10-1; 5-10-20-1; 5-10-10-10-1; 5-10-10-20-1; 5-10-20-30-111 MLs such as BPNN, SVR, etc.420 (30), 5-fold CVnoES, BPNN is the optimal
Fallahpour et al. [32]CCR (3, 3)trial-and-errorSVM48 (12)/ES
Yazdanparast et al. [33]Z-DEA (1, 17)17-?-1BPNN150 (45)noEF
Sreekumar et al. [34]CCR (3, 8), BCC (3, 8)11-?-1GRNN49/ES
Kao et al. [35]CCR (10, 2)12-?-1two SVMsc91/ES
Barros et al. [36]PCA-CCR (2, 2)4-20-1BPNN50noES
Sanei et al. [37]SBM (3, 3)5-6-1three BPNNs155 (46)noEF
Liu et al. [38]BBC (4, 4)7-8-4BPNN120 (20)noES
Notes: *: CCR(2,1) denotes the used CCR-DEA model with two inputs and one output. **: the “2-3-1” BPNN model denotes the used BPNN model with two inputs, three neurons in the hidden layer, and one output; the “?” means that the article does not specify the number of neurons in the hidden layer. c: the “c” denotes the model used for classification. #: 181 (37, 36) denotes that the numbers of the total, validation, and test subsets are 181, 37, and 36, respectively; 250 (50) denotes the numbers of the total and test subsets are 250, and 50, respectively; 160 (4:1) denotes that the total data are divided into training and test subsets according to a 4:1 ratio.
Table 2. The DMUs’ data and their DEA efficiency, as well as the predicted results.
Table 2. The DMUs’ data and their DEA efficiency, as well as the predicted results.
DMUsOriginal DataESEfficient/
Inefficient
PPR Model
I1 *I2O1O2PPR-1PPR-2PPRESN
A20151100901Efficient0.9960.0181.0141
B19131150501Efficient0.9740.0170.9921
C25160160550.883Inefficient0.902−0.0120.8910.874
D27168180721Efficient0.995−0.0080.9870.941
E2215894660.764Inefficient0.787−0.0110.7770.748
F55255230900.835Inefficient0.8070.0230.8310.791
G33235220880.902Inefficient0.930−0.0020.9280.902
H31206152800.796Inefficient0.806−0.0240.7820.752
I302441901000.960Inefficient0.8970.0220.9190.960
J502682501000.871Inefficient0.905−0.0280.8770.819
K533062601470.955Inefficient1.002−0.0350.9680.873
L382842501200.958Inefficient0.9790.0080.9840.958
B-C **20.9141.2160551Efficient0.9800.0120.9931(M) #
B-E16.8120.694661Efficient0.9820.0150.9970.963(N)
B-F33.8212.9230901Efficient1.002−0.0120.990
B-G29.8208.6220881Efficient1.0010.0131.013
B-H24.7164152801Efficient0.9960.0010.997
B-I28.8207.11901001Efficient1.0020.0141.016
B-J37.5233.32501001Efficient1.002−0.0180.985
B-K43.3292.32601471Efficient0.984−0.0120.972
B-L36.4259.52501201Efficient1.0000.0111.011
Notes: *: I1 and I2 denote two inputs, the doctor and nurse, respectively; O1 and O2 denote two outpatient and inpatient outputs, respectively; ES denotes the efficiency score. **: B-C, B-E, …, and K-L represent the benchmark of DMU C, E, …, and K, respectively. #: the letters M and N in the bracket denote the new DMUs, and the column of ESN is the efficiency scores after adding new DMUs.
Table 3. Comparison of the optimal number of doctors and nurses based on PPR and BPNN combined models.
Table 3. Comparison of the optimal number of doctors and nurses based on PPR and BPNN combined models.
DMUsI1I2O1O2DPPRNPPRDBNBDPPR-27NPPR-27DPPR-8NPPR-8DNN-729NNN-729
A201511009020.04150.9520.00150.1920.02150.9820.00151.0020.26149.17
B191311505020.08129.9119.01131.3120.07129.9220.04129.9520.19129.94
D271681807225.57169.4526.93167.1725.55169.4725.52169.5025.29168.20
C251601605521.63140.4421.06140.1621.62140.4521.59140.4821.47139.30
E22158946616.61120.8214.68122.4716.60120.8316.58120.8617.97122.48
F552552309032.51214.1833.84194.5132.48214.2032.45214.2332.56214.94
G332352208831.31207.0432.97191.4331.29207.0731.26207.1031.33207.61
H312061528023.84164.9024.86163.2623.82164.9223.79164.9523.58163.41
I3024419010029.86206.0630.86186.8929.84206.0929.81206.1229.85206.51
J5026825010035.61235.2535.42200.6435.59235.2735.55235.3135.69236.11
K5330626014742.15293.4233.42196.1542.12293.4642.09293.4842.05292.72
L3828425012037.99257.9334.96200.3837.96257.9637.93257.9938.03258.24
M251501707925.45272.52//22.43172.5425.40172.5725.18171.38
N303009013023.85191.45//23.81191.4923.79191.5123.73191.02
Table 4. Comparison of performance between different models with different virtual combinations.
Table 4. Comparison of performance between different models with different virtual combinations.
MAEMAPE (%)RMSER E A - m a x E R - m a x (%)
D B 2.3468.03.6950.8989.88029.6
D P P R 1.0685.31.1900.9891.5897.9
D S V R 1.0835.41.1930.9571.8279.1
D N N 1.1634.21.2810.9832.2626.8
D P P R - 27 1.0683.71.1880.9891.8825.6
D P P R - 8 1.0653.71.1820.9891.8205.3
D P P R - 729 1.0703.661.1920.9891.8945.8
D N N - 729 1.2094.391.2830.9881.8077.0
D S V R - 729 1.6216.471.7360.9833.16317.0
N B 27.4814.140.620.89596.1549.0
N P P R 1.0810.71.20211.9461.3
N S V R 1.0810.71.2120.9991.9991.3
N N N 1.1110.61.2530.9992.3791.0
N P P R - 27 1.0780.61.19911.9390.9
N P P R - 8 1.0750.61.19311.6720.9
N P P R - 729 1.0830.561.20411.9520.9
N N N - 729 1.2980.751.50412.8071.6
N S V R - 729 2.6481.653.2300.9996.6294.4
Table 5. Comparison of performance metrics of three subsets of DEA-BPNN combined models.
Table 5. Comparison of performance metrics of three subsets of DEA-BPNN combined models.
ModelTraining SubsetValidation SubsetTest Subset
MAEMAPEEA-maxER-maxMAEMAPEEA-maxER-maxMAEMAPEEA-maxER-max
D N N 0.4961.521.59711.680.5131.491.5154.350.5271.331.8226.75
N N N 2.1702.063.2458.602.4172.176.139.472.6492.6610.0610.29
D N N - 729 0.4291.132.25014.000.4291.313.2824.90.4301.212.0213.26
N N N - 729 0.8330.3715.7515.400.7870.309.322.130.8800.398.556.37
Table 6. The input–output data and ES of the DEA model and the predicted ES with the DEA-PPR combined model.
Table 6. The input–output data and ES of the DEA model and the predicted ES with the DEA-PPR combined model.
ProvincesCO2 (Mt)Population
(104 Persons)
GDP (CNY 100 Million)ESQuotasES-qQuotas-yES-yQuotas-2030ES-2030SVR
Beijing89.2215435,371189.20.6906426.50.651684.30.95780.9877
Tianjin158.5156214,1040.41766.10.7636179.20.6503157.70.60940.4351
Hebei914.2759235,1050.3133286.40.3407453.00.8803910.00.45060.3410
Shanxi564.9372917,0270.2562144.70.5744207.50.8229727.00.32220.2874
Inner Mongolia794.3254017,2130.1280101.70.6701547.60.2707712.20.34840.1726
Liaoning533.4435224,9090.3146167.80.5214361.30.6240504.00.24910.3401
Jilin203.7269111,7270.5245106.80.6702218.60.5558213.50.52360.5353
Heilongjiang278.2375113,6130.5231145.50.5779290.60.5776262.90.45800.5302
Shanghai192.9242838,1550.5581107.70.6301374.80.7978173.00.75420.5641
Jiangsu804.6807099,6320.6410515.81.2084787.10.9723760.20.86520.6389
Zhejiang381.4585062,3520.7201274.70.4306546.30.8834379.70.85580.7108
Anhui408.1636637,1140.5923241.70.4090366.20.9428427.80.45550.5950
Fujian278.1397342,3950.5922164.70.5287352.60.9392285.60.66150.5945
Jiangxi242.3466624,7580.7393179.10.4996276.80.8492267.50.64070.7402
Shandong937.110,07071,0680.4823452.00.4691824.40.7480885.40.58380.4947
Henan460.6964054,2590.9059417.30.3018528.11458.50.82360.8966
Hubei354.8592745,8280.6636235.40.4297417.20.8935371.90.57740.6600
Hunan310.6691839,7520.8430261.90.3874385.30.9731325.70.75270.8359
Guangdong569.111,521107,6711569.11.2367826.41510.31.46830.9877
Guangxi246.7496021,2370.7689189.70.4799232.61258.70.71770.7617
Hainan43.19455309143.10.854274.50.690067.80.82210.9860
Chongqing156.3312423,6060.7891123.30.6135238.10.8199155.50.85150.7736
Sichuan315.2837546,6161315.20.3491572.20.7915313.70.99020.9877
Guizhou261.1362316,7690.5394140.90.5831165.61288.30.50350.5414
Yunnan186.0485823,2241186.00.4869367.00.6421216.10.79920.9881
Shaanxi296.3387625,7930.5085150.70.5553308.70.7174294.90.51610.5184
Gansu164.5264787180.6397105.20.6820108.81172.40.66700.6323
Qinghai51.860829660.832243.10.8170118.20.321348.90.83530.8346
Ningxia212.469537480.202843.10.825738.01222.70.39590.2404
Xinjiang455.3252313,5970.2214100.80.6800271.30.4695408.20.24310.2574
Tianjin b 156214,104160.70.9783 0.6401
Hebei b 759242,2581285.70.9669 0.7387
Shanxi b 372920,7561140.30.9747 0.6924
Inner Mongolia b 254017,213196.70.9757 0.6219
Liaoning b 435224,9091164.00.9673 0.8048
Jilin b 269114,9781101.30.9869 0.6233
Heilongjiang b 375120,8781141.20.9745 0.6931
Shanghai b 251538,1551107.70.9107 0.8277
Jiangsu b 10,47999,6321515.80.9770 0.622
Zhejiang b 585062,3521274.70.8479 0.6197
Average///0.6239/0.6089/0.7828/0.6567/
Notes: b denotes the benchmark of the provinces of DEA-efficient DMUs. The quotas are obtained according to Refs. [56,57] using the ZSG-DEA model, and its ES is ES-q; we obtain the quotas (denoted in “Quotas-y” column of Table 6) according to an allocation system consisting of multiple variables. The quotas—2030 are obtained according to the GDP and carbon emission intensity of the “Outline” and “Plan”.
Table 7. Performance comparison of the PPR, BPNN, and SVR models in different virtual combinations.
Table 7. Performance comparison of the PPR, BPNN, and SVR models in different virtual combinations.
Training SubsetValidation Subset
ModelMAEMAPE (%)RMSE E A - m a x E R - m a x (%)MAEMAPE (%)RMSE E A - m a x E R - m a x (%)
C P P R - 32 12.303.2915.5343.2036.418.375.0514.6635.0027.37
C S V R - 32 10.925.0315.0045.78106.342.5427.4353.41132.249.06
C P P R - 243 12.081.8914.6945.3341.819.525.7314.8637.5122.82
C S V R - 243 11.402.1014.6751.29119.136.4524.3344.44106.244.77
C B P N N - 243 11.741.6014.5152.789.2722.7714.0330.8369.6552.94
C S V R - 243 7.391.8810.5471.50166.034.3528.8037.1359.4497.98
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, X.; Lou, W. An Exploration of Prediction Performance Based on Projection Pursuit Regression in Conjunction with Data Envelopment Analysis: A Comparison with Artificial Neural Networks and Support Vector Regression. Mathematics 2023, 11, 4775. https://doi.org/10.3390/math11234775

AMA Style

Yu X, Lou W. An Exploration of Prediction Performance Based on Projection Pursuit Regression in Conjunction with Data Envelopment Analysis: A Comparison with Artificial Neural Networks and Support Vector Regression. Mathematics. 2023; 11(23):4775. https://doi.org/10.3390/math11234775

Chicago/Turabian Style

Yu, Xiaohong, and Wengao Lou. 2023. "An Exploration of Prediction Performance Based on Projection Pursuit Regression in Conjunction with Data Envelopment Analysis: A Comparison with Artificial Neural Networks and Support Vector Regression" Mathematics 11, no. 23: 4775. https://doi.org/10.3390/math11234775

APA Style

Yu, X., & Lou, W. (2023). An Exploration of Prediction Performance Based on Projection Pursuit Regression in Conjunction with Data Envelopment Analysis: A Comparison with Artificial Neural Networks and Support Vector Regression. Mathematics, 11(23), 4775. https://doi.org/10.3390/math11234775

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop