Ecological Evaluation of Industrial Parks Using a Comprehensive DEA and Inverted-DEA Model

Data envelopment analysis (DEA) and inverted data envelopment analysis (inverted-DEA) are used so that the desirable and undesirable outputs of decision-making units (DMUs) exist simultaneously. We developed a new approach based on the concept of utilizing both DEA and inverted-DEA to enhance the discrimination power of DMUs with undesirable outputs. DMUs are ranked by the Z -score method and classiﬁed based on the eﬃciency scores of DEA and inverted-DEA. Then, the characteristics of the DMUs are analyzed based on the classiﬁcation result. This paper performs an eﬃciency evaluation of 21 industrial parks in China in 2017 using this new approach. The overall evaluation results indicate that the proposed new approach increases the discrimination ability in this empirical study.


Introduction
Data envelopment analysis (DEA) was proposed by three operational research experts, A. Charnes, W. Cooper, and E. Rhodes, in 1978 [1]. DEA is a new systematic analysis method used for evaluating the efficiency of similar decisionmaking units (DMUs) with multiple inputs and outputs.
is method does not need to restrict the production function and can avoid the subjective decision and objective factor dimensionality, and the impact of a unit's impact on evaluation results [2]. e classic DEA models first identify the production frontier for which the DMUs can be regarded as efficient. DMUs, which are outliers in terms of low inputs relative to the output level, map the efficient frontier. e DMUs on the frontier are efficient, and the DMUs not on the frontier are inefficient.
ose inefficient DMUs are compared with efficient DMUs to estimate their efficiency scores. DEA provides users with information about the efficiency scores and inefficiency scores and reference sets for inefficient units. One of the main features of DEA is to allow the DMUs to select their weights, which is a favorable approach for achieving maximum efficiency scores. However, this full flexibility may considerably reduce the discrimination power of DEA in the sense that too many DMUs often exist on the frontier, and some DMUs cannot be further ranked in standard DEA models.
To solve this problem, Allen et al. studied weight constraints and value judgments and concentrated on the implications of weight restrictions on the efficiency, targets, and comparisons of inefficient DMUs in DEA [3]. Additionally, Andersen and Petersen [4] employed a superefficiency DEA model to evaluate efficient DMUs; the model excludes the DMU being evaluated from the reference set.
is model was first used to identify outliers of observations by Banker and Gillord [5]. e model uses different reference sets to evaluate efficient and inefficient DMUs. Furthermore, Banker and Chang [6] reported that Andersen and Petersen's [4] procedure using superefficiency scores to rank efficiency observations yielded poor performance. Some scholars have tried to solve problems by using cross-efficiency methods. In such approaches, DMUs are evaluated based on their characteristics and those of other DMUs using a cross-efficiency matrix [7,8]. Although the crossefficiency method is often beneficial, the concept of crossefficiency score methods is considerabe from the basic principle of DEA.
Additionally, some scholars have attempted to construct new evaluation methods using DEA and inverted-DEA models. Amirteimoori [9] proposed the production possibility set and the quasi-production possibility set and constructed an alternative efficiency measure by using an efficiency frontier composed of the boundary points of the production possibility set and an antiefficient frontier composed of the boundary points of the quasi-production possibility set. However, this approach failed to solve the identification problem that occurs when many DMUs are efficient or inefficient. Zhang et al. [10] proposed an evaluation model of DMUs by using a good reference set and a bad reference set from the best and worst perspectives. To obtain the bad reference set, they simply treated the inputs and outputs of DMUs as undesirable. Cao et al. [11] used the evidential-reasoning (ER) approach to construct a performance indicator for combining the efficiency and antiefficiency values obtained by DEA and inverted-DEA models. Zhou et al. [12] used a DEA model without explicit inputs (see, e.g., Liu et al. [13], Liu et al. [14], and Yang et al. [15]) to combine efficient and antiefficient measures to rank DMUs. However, Shen et al. [16] easily verified that their approach could not significantly increase the discrimination power of DEA models and constructed three intuitive DEA performance indicators based on the distances to both the efficient and antiefficient frontiers. Although each performance indicator is useful for a specific sample size, no one approach can be referred to as a complete solution to all problems.
Furthermore, Entani et al. [17] employed both DEA and inverted-DEA models to obtain the upper and lower bounds of the interval efficiency of DMUs. ey argued that if the range of the interval efficiency is broad, then although the DMUs perform well from an optimistic perspective, they perform poorly from a pessimistic perspective. en, they used the interval efficiency to obtain a partial-order relation for DMUs. us, there has been a lack of clear distinctions among evaluations, and explicitly, certain information is difficult to obtain, providing little practical help to decisionmakers.
While the abovementioned techniques are useful in specific research areas, no one method provides a complete solution for all problems. People often have more than one reference perspective in assessing DMUs. e standard DEA models have employed the best-practice DMUs to construct the efficient frontier and have not fully taken advantage of the information implied in the data, especially for DMUs with undesirable outputs. In this paper, we explore a concept that involves enhancing the discrimination ability of DEA. In the case of the same inputs, the DEA model and the inverted-DEA model are used to enhance the discrimination ability of DMUs with desirable outputs and DMUs with undesirable outputs, respectively. e earliest work on the antiefficient frontier can be traced to the inverted-DEA model proposed by Yamada et al. [18]. Compared to the standard DEA models that evaluate DMUs from the perspective of optimism, the inverted-DEA model evaluates the performance from pessimism. e DEA method is not affected by the input and output dimensions of the problem and can comprehensively evaluate the data of different indicators. In particular, this approach has distinct advantages in dealing with multiple inputs and multiple outputs. However, most of the problems addressed with traditional DEA models have assumed that inputs and outputs are arbitrarily determined, the management activities of DMUs are controllable, and the output of DMUs should be optimized. However, in the actual production process, desirable and undesirable outputs may exist simultaneously. For example, a defective product is an undesirable output. In industrial production, economic benefits will also result in pollution, such as waste or smoke pollution, which is an undesirable output [19]. If the production processes that yield final products that generate wastes and pollutants are inefficient, the waste and pollutant outputs will be undesirable and should be reduced to improve performance. However, some wastes and pollutants are inevitably produced and cannot be reduced. us, the typical assumption of DEA is invalid. As a result, the undesirable and desirable outputs should be differently treated when DEA is used to evaluate the performance of DMUs. e most common approach is to consider only desirable outputs and ignore undesirable outputs [20]. However, such an evaluation method is overly simple and ignores essential information. erefore, performance evaluations based on such methods are not comprehensive.
Some researchers have suggested that some undesirable variables can be transformed [21]. For example, the ADD method proposed by Koopmans multiplies the undesirable output by −1 [19]. Similarly, it is also possible to add a translation vector based on negative transformations to keep the output data nonnegative [22]. Additionally, undesirable variables can be transformed nonlinearly, such as by multiplicative inverse operations [18]. However, these transformation methods may produce some unfavorable results [23]. Notably, transformations are often nonlinear and cannot retain convexity. In addition to applying transformation methods, undesirable outputs can also be regarded as inputs. If one treats the undesirable outputs as inputs, although the method is simple and easy to implement, the resulting DEA model does not reflect the actual production process. Additionally, the constraint function in linear programming can be adjusted, that is, the desirable and undesirable outputs classified, the output constraint function of the traditional DEA model is divided into two functions, and the desirable and undesirable outputs are constrained [24]. Furthermore, the distance measure can be adjusted to limit the range of undesirable outputs. e methods mentioned above for dealing with undesirable outputs in DEA mostly involve adjusting variables or models; notably, in the traditional model, undesirable outputs are processed through mathematical methods, but this approach does not reflect the actual situation. In this paper, we consider using an inverted-DEA to process undesirable outputs. Japanese scholar Yamada et al. proposed inverted-DEA in 1994 to evaluate the inefficiency score of DMUs, which is contrary to the concept of efficiency in DEA. When a DMU is inefficient, the model considers decreasing the output level and increasing the input level to improve efficiency, which reflects the objective of minimizing undesirable outputs under actual conditions. erefore, the use of an inverted-DEA model to address undesirable outputs can reflect the actual production process and is more straightforward and reasonable than other methods. e remainder of the paper is organized as follows. In Section 2, we present a comprehensive evaluation method for DMUs utilizing both DEA and inverted-DEA. Section 3 discussed a category analysis of DMUs. Section 4 applied the proposed method to an empirical data set consisting of 21 industrial parks in China in 2017. Finally, the discussion and conclusions are given in Section 5.

DEA.
DEA is an approach for analyzing the relative efficiency of peer DMUs that have multiple inputs and outputs. e evaluation of the DMU's efficiency is carried out by measuring the distance of this unit from the efficiency frontier created on the basis of the best units in the group, serving as benchmarks. In DEA, the maximum ratio of outputs is assumed to be efficiency, which is calculated from the optimistic perspective for each DMU. e efficiency for DMU 0 , which is analyzed as an object, is evaluated based on the efficiency values of the other DMUs. e following basic DEA model evaluates the efficiency of DMU 0 with s dimensional input vectors and m dimensional output vectors: where the decision variables are the weight vectors u and v; x j and y j are the input and output vectors for DMU j , respectively; and x 0 and y 0 are the input and output vectors for DMU 0 under evaluation, respectively. Each element of x j and y j is positive. To address many inputs and outputs, we consider the weighted sums of inputs and outputs as a hypothetical input and a hypothetical output. e efficiency is obtained by maximizing the ratio of weighted sum of outputs to that of inputs for DMU 0 under the condition that ratios for all DMUs are less than or equal to 1. is fractional programming problem is replaced with the following LP problem, which is called the C 2 R model, by limiting the denominator of the objective function to 1: For which the LP dual problem is When the optimal value of the objective function is equal to 1, DMU 0 is considered efficient, and otherwise, it is not deemed efficient. Specifically, "efficient" in this paper includes "weakly efficient." If we consider the slacks of inputs and outputs, we can then introduce the variables s − i and s + r and transform model (3) into the following model: where ε is a non-Archimedean infinitesimal and θ is the DEA efficiency score.  (4) are zero.

Inverted-DEA.
Since the inverted-DEA was first proposed in Japanese [18], we will illustrate the inverted-DEA method here. In contrast to DEA, which evaluates DMU 0 from the optimistic perspective, inverted-DEA is formulated as follows: Mathematical Problems in Engineering is fractional programming problem is also replaced with the following LP problem: By using the duality principle of linear programming, model (6) can be equivalently changed into the following linear programming model: where θ ′ represents the inverted-DEA inefficiency score of DMU 0 . We can then introduce variables s ′− i and s ′+ r and transform model (7) into the following model: We use one-dimensional input and two-dimensional output data to illustrate the difference between DEA and inverted-DEA. Figure 1 shows the processed data, the efficient frontier based on DEA, and the inefficient frontier based on inverted-DEA.

Comprehensive Evaluation.
Using the DEA model and inverted-DEA model, an efficiency analysis of DMUs is conducted with both desirable and undesirable outputs. Next, we will introduce four indicators (DEA analysis and inverted-DEA analysis for DMUs with inputs/desirable outputs and DEA analysis and inverted-DEA analysis for DMUs with inputs/undesirable outputs), which constitute a new efficiency measure, and rank DMUs according to this efficiency measure.
For n production units, or DMUs, we denote the input vector of DMUj as x j (j � 1, . . . , n), where x j ∈ E m , and the output vectors of DMUj with desirable outputs and undesirable outputs as y g j and y b j , where y g j , y b j ∈ E s . We denote the efficiency scores of the DMUs with desirable outputs and undesirable outputs as θ g and θ b , respectively. Similarly, we denote the inefficiency score of the DMUs with desirable outputs and the inefficiency score of the DMUs with undesirable outputs as θ ′ g and θ ′ b , respectively.
Using model (4), the DEA efficiency scores of each DMU with desirable and undesirable outputs are calculated as θ g and θ b , respectively. Similarly, by using model (8), the inverted-DEA inefficiency scores of each DMU with desirable and undesirable outputs are calculated as θ ′ g and θ ′ b , respectively. According to the properties of models (1) and (5), we In the sense of evaluating the efficiency of DMUs, larger values of θ g and θ ′ b are preferred, and small values of θ b and θ ′ g are ideal. us, the new model is as follows: It can be verified that e 0 is invariant to the units of the data. Furthermore, it holds that 0 ≤ e 0 ≤ 1.
Based on the information above, we can obtain the following definitions. Each index score of a DMU is marked on the corresponding coordinate axis, and these points are connected to form a planar geometric graph. e corresponding shaded area of the graph is the comprehensive evaluation score of the DMU.

Category Analysis of DMUs
In the section, two critical steps are performed. First, all of the DMUs are categorized based on the ranks obtained with equation (9). Second, all of the DMUs are categorized by using θ g , θ b , θ ′ g , and θ ′ b .
To compare the production level of several production units under similar conditions, the same type of experiment for multiple production units must be performed to evaluate the production skills by comparing the results of the experiments. For such joint detection, an evaluation rule is required. A general evaluation rule is based on the following factors: (a) expert opinions; (b) rules confirmed by authoritative departments; and (c) minimizing the influence of extreme values.
In this paper, the results are compared by using the previous result set, which is the Z-score-based method of evaluation.
If the data of z j (j � 1, . . . , n) follow a normal distribution, the Z-score is calculated using the quartile method according to the following steps.
Step 3: calculate the Z-score of the measurement value z j according to where μ is the criterion value and σ is the standard deviation.
According to the relation of the standard deviation σ, quartile, and IQR (interquartile range), the following equation can be obtained: us, In equation (10), if μ � z N2 , the Z-score of z j is as follows: Regarding the production units evaluated, we obtain two evaluation criteria based on equation (13) and |Z| � 3: where Z and Z are the reference values of Z-score Z j corresponding to Z � −3 and Z � 3, respectively. According to Z and Z, the Z-score can be divided into three categories: "excellent," "ordinary," and "pending improvement." Next, we categorize all of the DMUs by using θ * and θ ′ * .
is category analysis includes analyzing DMUs with desirable outputs and undesirable outputs. e specific method involves finding the appropriate thresholds α and β and dividing the DMUs into four categories according to the efficiency scores of DEA and inefficiency score of inverted-DEA. e specific categories are shown in Table 1. For a DMU with desirable outputs and undesirable outputs, the category analysis based on DEA/inverted-DEA is shown in Figure 3.
According to the category analysis of DMUs based on DEA/inversed DEA, combined with the comprehensive score of DMUs obtained using e 0 , the category of the DMU can be determined.

Case Study
We use the DEA model and the inversed DEA model to calculate the efficiency and inefficiency scores of 21 industrial parks, all of which are the first batch of circular economy pilot parks in China in 2017. Relevant data of the study come from research reports and statistical yearbooks of the samples. e input and output variables are defined as follows. e input variables are divided into four classes: resources (X 1 ), energy (X 2 ), land (X 3 ), and water resources (X 4 ). e output variables are divided into two classes: desirable outputs (D) and undesirable outputs (U). e desirable outputs include the resource output rate (D 1 ) and total production value (D 2 ), which reflect the economic output per unit input of material resources and the value of all final products and services produced by the park economy during a certain period. e undesirable outputs include sulfur dioxide emissions (U 1 ), industrial solid waste emissions (U 2 ), and industrial wastewater discharge (U 3 ).  Table 2.
According to the data from 21 industrial parks in Table 2 and model (4) and model (8), the relative efficiency scores and inefficiency scores of DMUs can be obtained (see Table 3).
As shown in Table 3, for the desirable outputs, the relative efficiency scores of nine DMUs have reached 1 for model (4); that is, nine DMUs are efficient based on DEA. Similarly, for the undesirable outputs, we have two inefficient DMUs based on the inverted-DEA and model (8). As the fourth and fifth columns of the table indicate, for the undesirable outputs, ten DMUs lie on the efficient frontier, so the efficiency scores of these units are equal to 1. Conversely, five DMUs lie on the antiefficient frontier, so the inefficiency scores of these units are equal to 1. Among them, DMU 2 and DMU 10 are not only DEA efficient but also inverted-DEA inefficient, which means that these two DMUs have efficient economic outputs but also produce efficient waste discharge. Many DMUs exist on the frontier that cannot be further ranked in standard DEA models and the inverted-DEA model. Additionally, the DMUs on both frontiers will have the same performance scores, which cannot be differentiated in traditional DEA methods.
To solve the above problem, we use the performance indicators proposed in Table 3 to calculate the comprehensive evaluation scores of 21 industrial parks based on model (9). e computed results are shown in Table 4.
To use the quartile method to categorize the data in Table 4, the distribution of these data must be analyzed. e analysis results are shown in Figure 4.
Because the comprehensive evaluation scores tend to be normal distribution, we can use the quartile method to calculate the Z-score and categorize the overall evaluation scores according to the Z-score. e specific steps are as follows.
Step 1: rank the comprehensive evaluation scores z j (j � 1, . . . , n) by ascending sorting. e ranked results are shown in Table 5.
Step 2: the Z-score Z j of z j is calculated according to equation (13). From the comprehensive evaluation scores, find the values for a quarter of z 6 � 0.027, two quarters of z 11 � 0.165, and three quarters of z 16 � 0.339 of the population.
Step 3: the evaluation criteria (Z � −0.528 and Z � 0.859) are determined using equations (14) and (15). According to Z and Z, the Z-scores are divided into three categories, namely, "excellent" (Z j < Z) for DMUs 2, 4, 6, 10, and 20; "ordinary" (Z < Z j < Z) for DMUs 1, 3, 5, 7, 8, 9, 16, 17, 18, and 21; and "pending improvement" (Z < Z j ) for DMUs 11,12,13,14,15,and 19. e computed results are shown in Table 5. e next steps involve categorizing all of the DMUs using θ g , θ b , θ ′ g , and θ ′ b in Table 3. e results are shown in Figure 5. A comprehensive feature analysis of DMUs is performed according to the category assessment standard given in Table 1. Table 6 gives the location of each DMU. Table 6 and Figure 5 show that the "excellent" DMUs in Table 5 are mostly distributed in the "excellent" area of Table 6, especially DMU 10 , which is "exemplary" in both the DEA evaluation and inverted-DEA evaluation because the input/output scale is small compared to that of other DMUs. Additionally, the two undesirable outputs (sulfur dioxide emissions and industrial solid waste emissions) are small (the sulfur dioxide emissions are 0). e "pending improvement" DMUs in Table 5 are mostly distributed in the "pending improvement" area of Table 6. In particular, DMU 15 is "pending improvement" in terms of evaluation of the desirable outputs of DEA and undesirable outputs of inverted-DEA because the input indexes (energy and land) are high (in comparison with the average) and the desirable output indexes (resource output rate and the total production value) are relatively low; moreover, the undesirable output indexes (industrial solid waste emissions and industrial wastewater discharge) are relatively high. Additionally, DMU 3 is in the "ordinary" area in terms of the desirable output evaluation based on DEA and undesirable output evaluation based on inverted-DEA, and it is also in the middle position in the rankings in Table 5. For DMU 3 , in addition to the high resource input (compared to the  average), the three input indexes (energy, land, and water resources) are generally low, and the two desirable output indexes (resource output rate and total production value) are relatively low; moreover, the three undesirable output indexes are relatively high, except for sulfur dioxide emissions. e other two indexes (industrial solid waste emissions and industrial wastewater discharge) of DMU 3 are low, especially industrial solid waste emissions (0.0103), for which the value is much lower than the average (0.6148).
It is also worth noting that DMU 11 and DMU 19 are "pending improvement" DMUs in Table 5. In Table 6, DMU 19 is in the "excellent" area for desirable outputs based on DEA and the "pending improvement" area for undesirable outputs based on inverted-DEA. e scenario for DMU 11 is the opposite of that for DMU 19 . e four input indexes of DMU 11 are higher than the average, and the desirable output indexes are lower than the average, so the efficiency score of inverted-DEA is high (1.00), and the efficiency score of DEA is very low (0.030). For the undesirable output indexes, except for higher sulfur dioxide emissions (compared with the average), other industrial solid waste emissions and industrial wastewater discharge are low, so the efficiency score of DEA of DMU 11 is high   Table 4 indicate that the proposed models increase the discrimination ability in this empirical study.  (1.00), and the efficiency score of inverted-DEA is deficient (0.056). us, DMU 11 should reduce the related inputs according to these findings and increase the related desirable outputs; specifically, reducing the level of sulfur dioxide emissions should be prioritised. e four input indexes of DMU 19 are lower than the average values, and the output indexes, whether desirable or undesirable, are lower than the average index values, so the efficiency scores of DEA are all 1 and the efficiency scores of inverted-DEA are 0.081 and 0.056 (lower than the  average). Although the three undesirable output indexes of DMU 19 are lower than the average values, the values of the desirable output indexes are relatively high. For DMU 19 , while improving the inputs, efforts should be made to improve the desirable output indexes and reduce the undesirable outputs.

Conclusion
Based on the concept of using inverted-DEA to process DMUs with undesirable outputs, this paper proposes a method for the comprehensive evaluation of DMUs with both desired and undesired outputs using a DEA model and an inverted-DEA model. is approach enhances the recognition capabilities of DEA. e DMUs are sorted and categorized by the quartile method. Based on this category analysis, combined with the efficiency scores of DEA and inverted-DEA, the overall rank of DMUs is characterized. Finally, this paper uses the proposed method to evaluate the efficiency of 21 industrial parks in China in 2017, and the characteristics of some individual industrial parks are evaluated based on the input and output indicators. e findings highlight the advantages of the proposed method, which provides an improved discrimination evaluation approach compared to existing methods by utilizing DEA and inverted-DEA.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Additional Points
Highlights. (i) A new performance evaluation model based on both DEA and inverted-DEA is developed to improve the recognition capability of DEA. (ii) An effective method is proposed to address DMUs with undesirable outputs based on inverted-DEA. (iii) e quartile method is used to categorize the ranked DMUs. (iv) e DMUs are categorized using the efficiency score of DEA and the inefficiency score of inverted-DEA. (v) e proposed method is validated based on an empirical analysis of 21 industrial parks in China.

Conflicts of Interest
e authors declare no conflicts of interest.