Next Article in Journal
Development and Validation of a Safety Attitude Scale for Coal Miners in China
Next Article in Special Issue
A Novel Grey Wave Method for Predicting Total Chinese Trade Volume
Previous Article in Journal
A Holistic Analysis Approach to Social, Technical, and Socio-Technical Aspect of E-Government Development
Previous Article in Special Issue
Going Public and Industrial Upgrading of Traditional Clusters in Developing Countries: Rethinking the Dynamics of the ‘Jinjiang Model’ in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Interval Data Envelopment Efficiency Model Considering Different Distribution Characteristics—Based on Environmental Performance Evaluation of the Manufacturing Industry

Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, School of Economics and Management, China Institute for Manufacture Developing, Nanjing University of Information Science and Technology, Nanjing 210044, China
*
Author to whom correspondence should be addressed.
Sustainability 2017, 9(12), 2080; https://doi.org/10.3390/su9122080
Submission received: 7 October 2017 / Revised: 31 October 2017 / Accepted: 6 November 2017 / Published: 27 November 2017
(This article belongs to the Special Issue Transition from China-Made to China-Innovation )

Abstract

:
This study utilizes the Data Envelopment Efficiency (DEA) model to assess input–output efficiency from two perspectives. First, not considering the distribution of interval data, we introduce an adjusted parameter to transform interval data to determination data. Second, by contrast, we take into account the distribution characteristics of interval data and test the DEA model with interval data based on linear uniform distribution and normal distribution with uncertainty. Based on the normal distribution DEA evaluation model, this paper aims to evaluate the input–output performance of the manufacturing industry with the constraint of environmental pollution in the Yangtze River Delta (YRD) region, China. Research has shown that the optimal solution of the normal distribution model is better than that of linear distribution. Therefore, it is imperative to adopt an appropriate method to evaluate the energy and environmental efficiency of this region.

1. Introduction

To maintain sustainable economic development, China is now subject to a number of crucial constraints, including but not limited to: the outstanding supply–demand discrepancies in energy, severe environmental pollution, and regional discrepancies in energy efficiency. Enhancing energy efficiency and environmental efficiency is an important way to confront energy crisis and to prevent environmental pollution, and, ultimately, to boost China’s sustainable development capability and to promote economic restructuring. The Yangtze River Delta agglomeration is one of the regions with the highest level of economic development in China. The heightened levels of industrialization, urbanization and energy consumption have led to the dramatic decline in energy supply and the deterioration of environmental quality. Drawing upon an investigation on the energy efficiency with the environmental constraints of the Yangtze River Delta region, we explore the causes of low energy efficiency and seek practical recommendations to improve energy efficiency in order to alleviate some of the pressure on energy consumption. Many authors think it is useful to explain the productivity efficiency of Chinese manufacturing sectors using DEA [1,2,3,4]. Hence, the efficiency evaluation method of data envelopment analysis (DEA) is used to explore the energy efficiency differences of 16 cities in the Yangtze River Delta region from 2005 to 2014 under the environmental constraints.
Data envelopment analysis (DEA) was originally proposed by Charnes [5] and is a system analysis method based on “relative efficiency”. As a non-parametric evaluation method of the DEA model, the project efficiency evaluation problem with single-input and single-output gradually expanded into a complex system evaluation problem with multi-inputs and multi-outputs, which is widely used in many fields, such as mathematical economics and operations research. If the input–output data of decision-making units (DMUs) in the DEA model are accurate, the DEA method belongs to the definite type of DEA. However, in practice, due to measurement errors and data noise, amongst other reasons, the evaluation indicator values of each DMU often cannot be accurately determined, so the CCR [6,7,8], BCC [9,10,11], SBM [12], SBI DEA [13,14] and TurboDEA [15] models were derived.
Furthermore, in practice, problems such as policy strength assessment and investment projects evaluation arise and some or even all of the input–output indicators of DMUs can only be uncertain values due to the qualitative policy strength and uncertainty of economic phenomena. To overcome the limitations of the traditional DEA method, many authors have expanded to the random DEA [16,17], interval DEA [18,19,20,21] and normal DEA models [22,23,24], amongst others, which are mainly used to solve the problem relating to the DEA model with uncertain values.
(1) When the input and output variables are denoted as stochastic, the corresponding DEA model is known as the stochastic DEA model. The distribution of the input–output variables includes logarithmic normal, exponential and normal distributions. Cooper et al. [25] supposed that the variable obeys logarithmic normal distribution and as there was no abnormal variable, and the original addition model was converted to the multiplication model for efficiency solution. The study of the stochastic DEA model with logarithmic normal distribution also includes [26,27,28,29,30,31]. Zhang [32] further proposed a chance constraint to the stochastic DEA model based on exponential distribution, since input–output indicators obeyed the same distribution, and analyzed the effectiveness of corresponding DMUs in theory. The stochastic DEA model of exponential distribution also includes [33,34,35,36]. Jeang [37] assumed that the stochastic variables in the chance constraint model obeyed normal distribution. Consequently, the nonlinear stochastic DEA model was applied to solve the maximum output scale by using computer software (such as computer-aided engineering and DEA software).
(2) Under normal conditions, the stochastic DEA model is frequently solved by the stochastic simulation method. The concept of the satisfaction degree was introduced into the stochastic DEA problem and the consensus DEA model was proposed, e.g., by [38,39,40,41,42,43]. Furthermore, the confidence region or level was introduced into the stochastic DEA model, e.g., by [44,45,46,47,48]. The Monte Carlo simulation algorithm was used to solve the problem of stochastic input and output variables and the study of stochastic DEA is one of the directions of current research, such as studies by [49,50,51,52,53,54]. The stochastic DEA model including the Bayesian method, such as in studies by [55,56,57,58,59,60], was used to deduce the efficiency value of a finite sample.
(3) For special stochastic DEA models, the isometric transformation method is an effective method for solving the stochastic DEA model. Khodabakhshi and Asgharian [61] constructed a general deterministic stochastic DEA model that can replace the nonlinear model under normal conditions. Eslami et al. [62] dealt with a real decision problem with fuzzy constraints and uncertain information (stochastic data). The fuzzy probability constraint in the uncertain constraint DEA model was transformed into the equivalent form of certainty. Zhou et al. [63] took advantage of chance-constrained programming, the key value reduction method and equivalent transformation to solve stochastic uncertain problems and to measure the effectiveness, efficiency and productivity in uncertain environments and for different confidence levels. Yu et al. [64] combined time series and regression methods to construct a new stochastic chance constrained model and used it to evaluate and analyze the efficiency of 15 commercial banks from 2005 to 2011.
This study utilizes DEA model to assess input–output efficiency from two perspectives. First, not considering the distribution of interval data, we introduce an adjusted parameter to transform interval data to determination data. Second, by contrast, we take into account the distribution characteristics of interval data and test the DEA model with interval data based on linear uniform distribution and normal distribution with uncertainty. Research has shown that the optimal solution of the normal distribution model is better than that of linear distribution. Based on the normal distribution DEA evaluation model, this paper aims to evaluate the input–output performance of the manufacturing industry with the constraint of environmental pollution in the YRD region, China.
In practical efficiency evaluation, interval input–output variables often exist. In conventional interval DEA model processing, the efficiency solution of the DMUs is restricted by the interval endpoint values (upper and lower limits) if only the intervals are considered. In this study, we investigate the relative effectiveness of DMUs with different distributions, based on two cases of intervals obeying linear and normal distribution, respectively.

2. Preliminaries and Methodology

2.1. Standard DEA Model with Crisp Data

Supposing that there are n DMUs ( D M U j ( j = 1 , 2 , , n ) ) , m types of different resources invested and s kinds of different outputs yielded, input–output vectors of D M U j would be denoted separately as follows:
x i j = ( x 1 j , x 2 j , , x m j ) T , i = 1 , 2 , , m ,
y r j = ( y 1 j , y 2 j , , y s j ) T , r = 1 , 2 , , s .
For D M U j 0 , compared to other DMUs in efficiency solving, the traditional CCR model is constructed as follows [5]:
Z j 0 1 = M a x . { ( r = 1 s μ r y r j 0 ) / ( i = 1 m ω i x i j 0 ) , } s . t . ( r = 1 s μ r y r j ) / ( i = 1 m ω i x i j ) 1 , j = 1 , 2 , , n , ( 1 - 1 ) , μ r 0 , r = 1 , 2 , , s , ( 1 - 2 ) , ω i 0 , i = 1 , 2 , , m , ( 1 - 3 ) ,
where μ r and ω i are the weight coefficients of outputs and inputs, respectively, and their values are not negative. The objective function is the maximum efficiency for D M U j 0 , and D M U j 0 is evaluated DMU in this model. D M U j 0 is reference target for non-efficient DMU [65].
In Model (1), to keep the input unchanged and maximize the output, we construct a linear programming Model (2) as follows:
Z j 0 2 = M a x . { r = 1 s μ r j 0 y r j 0 } , s . t . i = 1 m ω i x i j 0 = 1 , ( 2 - 1 ) , r = 1 s μ r y r j i = 1 m ω i x i j 0 , j = 1 , 2 , , n , ( 2 - 2 ) , μ r 0 , r = 1 , 2 , , s , ( 2 - 3 ) , ω i 0 , i = 1 , 2 , , m , ( 2 - 4 ) .
Constraint (2-1) represents the unchanged input of the DMUs and the objective function is to maximize the output of the D M U j 0 , D M U j 0 is evaluated DMU in this model.
The duality of linear programming can be denoted by the real variable and the non-negative vector λ = ( λ 1 , λ 2 , , λ n ) [66]. To avoid the weakly efficient problem with a non-zero slack variable, we add the constraint of the sum of the maximum slack variables in the objective function. In this paper, Model (2) is transformed into Model (3) as follows:
Z j 0 3 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . j = 1 n λ j y r j S r + = y r j 0 , r = 1 , 2 , , s , ( 3 - 1 ) , j = 1 n λ j x i j + S i = θ x i j 0 , i = 1 , 2 , , m , ( 3 - 2 ) , λ j 0 , S i 0 , S r + 0 , j = 1 , 2 , , n , ( 3 - 3 ) .
D M U j 0 is evaluated DMU in this model. S r + and S i are introduced as slack variables in the model and the selection of slack variables do not affect the value of the optimal solution θ * . With the same output, the objective function of the dual model is to find the difference between the minimum variable elasticity and the slack variable. If θ * = 1 and all the slack variables are zero, i.e., S r + * = 0 and S i * = 0 , then it is absolutely efficient. If θ * = 1 and all the slack variables are not zero, i.e., S r + * 0 or S i * 0 , then the DMUs are weakly efficient.
In Model (3), the constraint j = 1 n λ j = 1 is added and the model could subsequently be transformed into the following BCC model:
Z j 0 4 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . j = 1 n λ j y r j S r + = y r j 0 , r = 1 , 2 , , s , ( 4 - 1 ) , j = 1 n λ j x i j + S i = θ x i j 0 , i = 1 , 2 , , m , ( 4 - 2 ) , j = 1 n λ j = 1 , j = 1 , 2 , , n , ( 4 - 3 ) , λ j 0 , S i 0 , S r + 0 , j = 1 , 2 , , n , ( 4 - 4 ) .
D M U j 0 is evaluated DMU in this model. Constraint (4-3) represents the variable scale returns, i.e., the ratio of the output increment of the DMUs is not equal to that of the increment of the input factor.

2.2. Interval DEA model

In Model (4), the input and output indicators are all accurate numbers. In practice, the input and output indicators may be intervals, and DEA Model (5) based on the interval indicators is constructed as follows:
Z j 0 5 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . j = 1 n λ j y r j S r + = y r j 0 , r = 1 , 2 , , s , ( 5 - 1 ) , j = 1 n λ j x i j + S i = θ x i j 0 , i = 1 , 2 , , m , ( 5 - 2 ) , y ˜ r j = [ y r j L , y r j U ] , ( 5 - 3 ) , x ˜ i j = [ x i j L , x i j U ] , ( 5 - 4 ) , j = 1 n λ j = 1 , j = 1 , 2 , , n , ( 5 - 5 ) , λ j 0 , S i 0 , S r + 0 , j = 1 , 2 , , n , ( 5 - 6 ) .
In Model (5), the input–output indicators, y ˜ r j = [ y r j L , y r j U ] ; x ˜ i j = [ x i j L , x i j U ] , are interval data, D M U j 0 is evaluated DMU in this model.
From the perspective of optimization, we first transform the input and output interval into an accurate number to obtain the optimal solution for Model (5). That is, let y ˜ r j = ( α r j y r j U + ( 1 α r j ) y r j L ) , 0 α r j 1 , x ˜ i j = ( α ¯ i j x i j U + ( 1 α ¯ i j ) x i j L ) , and 0 α ¯ i j 1 , respectively. The interval DEA model could be transformed into the conventional accurate DEA Model (6) as follows:
Z j 0 6 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . j = 1 n λ j y r j S r + = y r j 0 , r = 1 , 2 , , s , ( 6 - 1 ) , j = 1 n λ j x i j + S i = θ x i j 0 , i = 1 , 2 , , m , ( 6 - 2 ) , y ˜ r j = α r j y r j U + ( 1 α r j ) y r j L , ( 6 - 3 ) , x ˜ i j = α ¯ i j x i j U + ( 1 α ¯ i j ) x i j L , ( 6 - 4 ) , j = 1 n λ j = 1 , j = 1 , 2 , , n , ( 6 - 5 ) , λ j 0 , S i 0 , S r + 0 , j = 1 , 2 , , n , ( 6 - 6 ) .
D M U j 0 is evaluated DMU in this model.

2.2.1. Examples of Interval DEA Model

There were six DMUs and the input indicator x 1 , and output indicators y 1 and y 2 are intervals (Table 1).
Taking the input–output data in Table 1 as the example data, which are similar to those in the traditional Model (2), the interval DEA nonlinear model is constructed as follows:
Z j 0 7 = M a x . { u 1 ( 3 + 2 a 11 ) + u 2 ( 7 + a 21 ) } , s . t . w 1 ( 1 + a ¯ 11 ) = 1 , ( 7 - 1 ) , u 1 ( 3 + 2 a 11 ) + u 2 ( 7 + a 21 ) w 1 ( 1 + a ¯ 11 ) 0 , ( 7 - 2 ) , u 1 ( 5 + 2 a 12 ) + u 2 ( 6 + 2 a 22 ) w 1 ( 2 + a ¯ 12 ) 0 , ( 7 - 3 ) , u 1 ( 30 + 4 a 13 ) + u 2 ( 32 + 3 a 23 ) w 1 ( 4 + 2 a ¯ 13 ) 0 , ( 7 - 4 ) , u 1 ( 10 + 3 a 14 ) + u 2 ( 4 + a 24 ) w 1 ( 3 + a ¯ 14 ) 0 , ( 7 - 5 ) , u 1 ( 27 + 2 a 15 ) + u 2 ( 20 + a 25 ) w 1 ( 1 + 2 a ¯ 15 ) 0 , ( 7 - 6 ) , u 1 ( 25 + a 16 ) + u 2 ( 85 + 5 a 26 ) w 1 ( 3 + a ¯ 16 ) 0 , ( 7 - 7 ) , 0 a ¯ 11 1 , 0 a ¯ 12 1 , 0 a ¯ 13 1 , ( 7 - 8 ) , 0 a ¯ 14 1 , 0 a ¯ 15 1 , 0 a ¯ 16 1 , ( 7 - 9 ) , 0 a 11 1 , 0 a 12 1 , 0 a 13 1 , ( 7 - 10 ) , 0 a 14 1 , 0 a 15 1 , 0 a 16 1 , ( 7 - 11 ) , 0 a 21 1 , 0 a 22 1 , 0 a 23 1 , ( 7 - 12 ) , 0 a 24 1 , 0 a 25 1 , 0 a 26 1 , ( 7 - 13 ) .
D M U 1 is evaluated DMU in this model. Let A r j = μ r α r j , for 0 α r j 1 , so 0 A r j μ r ; and let B i j = ω i α ¯ i j , for 0 α ¯ i j 1 , so 0 B i j ω i .
The nonlinear interval DEA Model (7) is transformed into a linear Model (8) as follows:
Z j 0 8 = M a x . { 3 u 1 + 2 A 11 + 7 u 2 + A 21 } , s . t . w 1 + B 11 = 1 , ( 8 - 1 ) , ( 3 u 1 + 2 A 11 ) + ( 7 u 2 + A 21 ) ( w 1 + B 11 ) 0 , ( 8 - 2 ) , ( 5 u 1 + 2 A 12 ) + ( 6 u 2 + 2 A 22 ) ( 2 w 1 + B 12 ) 0 , ( 8 - 3 ) ( 30 u 1 + 4 A 13 ) + ( 32 u 2 + 3 A 23 ) ( 4 w 1 + 2 B 13 ) 0 , ( 8 - 4 , ) ( 10 u 1 + 3 A 14 ) + ( 4 u 2 + A 24 ) ( 3 w 1 + B 14 ) 0 , ( 8 - 5 ) , ( 27 u 1 + 2 A 15 ) + ( 20 u 2 + A 25 ) ( w 1 + 2 B 15 ) 0 , ( 8 - 6 ) , ( 25 u 1 + A 16 ) + ( 85 u 2 + 5 A 26 ) ( 3 w 1 + B 16 ) 0 , ( 8 - 7 ) , A 11 - u 1 0 ; A 12 u 1 0 ; A 13 u 1 0 , ( 8 - 8 ) , A 14 u 1 0 ; A 15 u 1 0 ; A 16 u 1 0 , ( 8 - 9 ) , A 21 u 2 0 ; A 22 u 2 0 ; A 23 u 2 0 , ( 8 - 10 ) , A 24 u 2 0 ; A 25 u 2 0 ; A 26 u 2 0 , ( 8 - 11 ) , B 11 - w 1 0 ; B 12 w 1 0 ; B 13 w 1 0 , ( 8 - 12 ) , B 14 w 1 0 ; B 15 w 1 0 ; B 16 w 1 0 , ( 8 - 13 ) .
D M U 1 is evaluated DMU in this model.
According to the dual equivalence principle, the linear interval DEA Model (8) is transformed into a dual linear Model (9) as follows:
Z j 0 9 = M i n . { y 1 } , s . t . y 1 y 2 2 y 3 4 y 4 3 y 5 y 6 3 y 7 y 20 y 21 y 22 y 23 y 24 y 25 0 , ( 9 - 1 ) , 3 y 2 + 5 y 3 + 30 y 4 + 10 y 5 + 27 y 6 + 25 y 7 y 8 y 9 y 10 y 11 y 12 y 13 3 , ( 9 - 2 ) , 7 y 2 + 6 y 3 + 32 y 4 + 4 y 5 + 20 y 6 + 85 y 7 y 14 y 15 y 16 y 17 y 18 y 19 7 , ( 9 - 3 ) , 2 y 2 + y 8 2 ; 2 y 3 + y 9 0 ; 4 y 4 + y 10 0 , ( 9 - 4 ) , 3 y 5 + y 11 0 ; 2 y 6 + y 12 0 ; y 7 + y 13 0 , ( 9 - 5 ) , y 2 + y 14 1 ; 2 y 3 + y 14 0 ; 3 y 4 + y 15 0 , ( 9 - 6 ) , y 5 + y 16 0 ; y 6 + y 17 0 ; 5 y 7 + y 18 0 , ( 9 - 7 ) , y 1 y 2 + y 20 0 ; y 3 + y 21 0 ; 2 y 4 + y 22 0 , ( 9 - 8 ) , y 5 + y 23 0 ; 2 y 6 + y 24 0 ; y 7 + y 25 0 , ( 9 - 9 ) .
Yagi et al. [67] regard the variables Y = { y 1 , y 2 , , y 25 } of the above model as shadow price, D M U j is efficient when the profit of shadow price equals zero. D M U 1 is evaluated DMU in this model. We use Lingo 11 software (LINDO SYSTEMS Inc., Chicago, IL, USA) to solve the aforementioned model and the optimal solution is shown in Table 2.
From Table 2, D M U 1 - D M U 4 are not efficient, whereas D M U 5 - D M U 6 are efficient; the efficiency from D M U 1 to D M U 6 therefore increases. In the next section, we discuss the improvement of the non-efficient DMUs in the DEA model.

2.2.2. The Improvement of the Non-Efficient DMU

Definition 1.
The DEA’s relative efficiency (efficient production frontier) is π : ω 0 T X 0 μ 0 T Y 0 = 0 . If the D M U K 0 is DEA efficient, and the optimal solution ω 0 T , μ 0 T of linear programming (P) satisfies the condition
V p = μ 0 T Y 0 = 1 , ω 0 0 , μ 0 0 .
Furthermore, ω 0 T X 0 = 1 , so ω 0 T X 0 = μ 0 T Y 0 . Therefore, the point ( X 0 , Y 0 ) is on the hyper plane π and the other points ( X , Y ) on the hyper plane π are also DEA efficient. Therefore, we use the method of “projection” on the relative efficient surface to improve the DMUs that were non-efficient for DEA.
Definition 2.
Assume that λ 0 , s 0 , s 0 + , θ 0 are the optimal solution to the duality of linear programming, and let
X ^ 0 = θ 0 X 0 s 0 , Y ^ 0 = Y 0 + s 0 +
We call ( X ^ 0 , Y ^ 0 ) of D M U K 0 the projection of the corresponding ( X 0 , Y 0 ) on the DEA relative efficient surface π.
Non-efficient D M U 1 could be constructed into a new D M U 7 , as seen in Table 3.
The dual nonlinear Model (10) of a new interval DEA is constructed as follows:
Z j 0 10 = M i n . { u 1 ( 5 + 2 a 17 ) + u 2 ( 8 + a 27 ) } s . t . w 1 ( 0.32 + 0.63 a ¯ 11 ) = 1 , ( 10 - 1 ) , u 1 ( 3 + 2 a 11 ) + u 2 ( 7 + a 21 ) w 1 ( 1 + a ¯ 11 ) 0 , ( 10 - 2 ) , u 1 ( 5 + 2 a 12 ) + u 2 ( 6 + 2 a 22 ) w 1 ( 2 + a ¯ 12 ) 0 , ( 10 - 3 ) , u 1 ( 30 + 4 a 13 ) + u 2 ( 32 + 3 a 23 ) w 1 ( 4 + 2 a ¯ 13 ) 0 , ( 10 - 4 ) , u 1 ( 10 + 3 a 14 ) + u 2 ( 4 + a 24 ) w 1 ( 3 + a ¯ 14 ) 0 , ( 10 - 5 ) , u 1 ( 27 + 2 a 15 ) + u 2 ( 20 + a 25 ) w 1 ( 1 + 2 a ¯ 15 ) 0 , ( 10 - 6 ) , u 1 ( 25 + a 16 ) + u 2 ( 85 + 5 a 26 ) w 1 ( 3 + a ¯ 16 ) 0 , ( 10 - 7 ) , u 1 ( 5 + 2 a 17 ) + u 2 ( 8 + a 27 ) w 1 ( 0.32 + 0.63 a ¯ 17 ) 0 , ( 10 - 8 ) , 0 a ¯ 11 1 , 0 a ¯ 12 1 , 0 a ¯ 13 1 , ( 10 - 9 ) , 0 a ¯ 14 1 , 0 a ¯ 15 1 , 0 a ¯ 16 1 , ( 10 - 10 ) , 0 a ¯ 17 1 , 0 a 11 1 , 0 a 12 1 , ( 10 - 11 ) , 0 a 13 1 , 0 a 14 1 , 0 a 15 1 , ( 10 - 12 ) , 0 a 16 1 , 0 a 17 1 , 0 a 21 1 , ( 10 - 13 ) , 0 a 22 1 , 0 a 23 1 , 0 a 24 1 , ( 10 - 14 ) , 0 a 25 1 , 0 a 26 1 , 0 a 27 1 , ( 10 - 15 ) ,
D M U 1 is evaluated DMU in this model. D M U 1 is evaluated DMU in this model. If the efficiency of D M U 7 is 1, then the new D M U 7 is efficient and D M U 1 could be constructed in accordance with the data of new DMU. Similarly, other non-efficient DMUs could be changed into new DMUs, as shown in Table 4.

2.3. Uncertain DEA Model

There are disordered, probable, vague and approximate attributes, i.e., uncertainty in the process of contact and development of objective things. Therefore, Liu [68] introduced the uncertainty theory mainly to solve decision-making optimization under the uncertain information environment. It is a mathematical theory concerning the quantitative characteristics of various uncertain phenomena.
In this section, the input–output objectives y ˜ r j = [ y r j L , y r j U ] ; x ˜ i j = [ x i j L , x i j U ] , are assumed separately to obey the uncertain linear distribution (as seen in Section 2.1) and the uncertainty of the normal distribution (as seen in Theorem 1 in Supplementary Materials).
Based on the dual Model (3), we construct the unchanged inputs, maximum outputs and opportunity constraint DEA Model (11) as follows:
Z j 0 11 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . M { j = 1 n λ j y ˜ r j y ˜ r j 0 + S r + } 1 α r j , r = 1 , 2 , , s , ( 11 1 ) , M { j = 1 n λ j x ˜ i j θ x ˜ i j 0 S i } 1 α ¯ i j , i = 1 , 2 , , m , ( 11 2 ) , j = 1 n λ j = 1 , λ j 0 , j = 1 , 2 , , n , ( 11 3 ) , 0 α r j 1 , 0 α ¯ i j 1 , ( 11 4 ) .
D M U j 0 is evaluated DMU in this model, where ε is a non-Archimedes infinitesimally small constant, namely 10 6 , and M is the uncertainty measure. Model (11-1) shows that the probability of occurrence of expected outputs of the DMUs being larger than the actual outputs is 1 α r j , which improves the expected outputs. Model (11-2) shows that when D M U j 0 ’s outputs y ˜ r j 0 remain unchanged, we should assure that the input values x ˜ i j 0 decrease as the same ratio θ , i.e., the probability of occurrence of expected inputs being less than the actual investment is 1 α ¯ i j and the corresponding adjustments are made to the DMUs that are deficient compared to DMUs that are efficient for DEA.
The objective function can be seen as having two parts, one of which is to minimize the real number θ of the linear programming dual model and the other is to maximize the value of the total slack variables in the input–output constraint. Assuming that the optimal problem is given by Model (11), we obtained the optimal solution that λ 0 , S i j 0 , S r j 0 + , θ 0 . Furthermore, there is θ 0 = 1 and S i j 0 = 0 , S r j 0 + = 0 . We subsequently call the evaluation DMU is efficient compared to other DMUs and satisfied both the technical and scale effectiveness. If θ 0 = 1 , but S i j 0 or S r j 0 + are not equal to the zero vector at the same time, then we call that the DMUs are weakly efficient and the DMUs are not technically and scale efficient at the same time. If θ 0 < 1 , then the DMUs are evaluated as DEA non-efficient.

2.3.1. Uncertain DEA Model Based on Linear Distribution

Definition 3.
Assume that an uncertain variable ξ obeys the linear uncertainty distribution as follows:
ϕ ( x ) = 0 , x a , ( x a ) / ( b a ) , a x b , 1 x b .
Consequently, ϕ ( x ) obeys the linear probability distribution, where a and b are definite numbers and a < b . Then, we call ξ obeys the linearly uniform distribution of a , b , expressed as ξ ( a , b ) .
Definition 4.
Let ξ be the uncertain variable obeying the linear uncertainty distribution ϕ ( x ) (Figure 1). The inverse function ϕ 1 ( x ) is called the inverse uncertain distribution function of an uncertain variable ϕ ( x ) . The inverse uncertain distribution function of linear uncertain distribution [ a , b ] is as follows:
ϕ 1 ( α ) = ( 1 α ) a + α b .
Based on the certain input and the largest output background, the uncertain DEA model based on linear distribution is constructed as follows:
Z j 0 12 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . M { j = 1 n λ j y ˜ r j y ˜ r j 0 S r + } 1 α r j , r = 1 , 2 , s , ( 12 - 1 ) , M { j = 1 n λ j x ˜ i j + S i θ x ˜ i j 0 } 1 α ¯ i j , i = 1 , 2 , , m , ( 12 - 2 ) , y ˜ r j [ y r j L , y r j U ] , x ˜ i j [ x i j L , x i j U ] , ( 12 - 3 ) , j = 1 n λ j = 1 , λ j 0 , j = 1 , 2 , , n , ( 12 - 4 ) , 0 α r j 1 , 0 α ¯ i j 1 , ( 12 - 5 ) , S i 0 , i = 1 , 2 , , m , ( 12 - 6 ) , S r + 0 , r = 1 , 2 , , s , ( 12 - 7 ) .
D M U j 0 is evaluated DMU in this model. Theorem 2 and Proof 0.2 (see Supplementary Materials).
According to Definition 5, Model (12) is transformed into a deterministic DEA Model (13).
Definition 5.
Uncertain linear distribution of DEA Model (12) with chance constraints is equivalent to the determined DEA Model (13):
Z j 0 13 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . j = 1 n λ j ( α r j y r j U + ( 1 α r j ) y r j L ) ( S r + ) ( ( 1 α r j ) y r j 0 U + α r j y r j 0 L ) , r = 1 , , s , ( 13 - 1 ) , j = 1 n λ j ( α ¯ i j x i j L + ( 1 α ¯ i j ) x i j U ) + ( S i ) θ ( ( 1 α ¯ ) i j x i j 0 L + α ¯ i j x i j 0 U ) , i = 1 , , m , ( 13 - 2 ) , j = 1 n λ j = 1 , λ j 0 , j = 1 , 2 , , n , ( 13 - 3 ) , 0 α r j 1 , 0 α ¯ i j 1 , ( 13 - 4 ) , S i 0 , r = 1 , 2 , , s , ( 13 - 5 ) , S r + 0 , i = 1 , 2 , , m , ( 13 - 6 ) .
D M U j 0 is evaluated DMU in this model. As α ¯ i j and α r j increases, the objective function φ is an increasing function of the probability α ¯ i j and α r j , and the objective function φ becomes larger. Similarly, if α ¯ i j increases, the objective function φ also becomes larger.

2.3.2. A Numerical Example of Uncertain Linear Distribution DEA Model Based on Opportunity Constraint

Assuming that there are still six DMUs for the interval model, the input indicators x 1 , and the output indicators y 1 and y 2 are the interval data in Table 1, which are subjected to linear uncertain distribution. The data can be seen in Table 5.

2.4. A Numerical Example of Equivalent DEA Nonlinear Transformation Model for Uncertainty Distribution

The input and output interval data (Table 5) indicators of the DMUs subjected to linear distribution are taken as the numerical example. We construct the equivalent nonlinear DEA Model (14) from the linear distribution model, which is similar to the certain linear distribution DEA Model (13):
Z j 0 14 = M i n . { x 7 0.000001 ( s 1 + s 2 + s 3 ) } , s . t . u 1 ( 3 + 2 a 11 ) u 2 ( 5 + 2 a 12 ) u 3 ( 30 + 4 a 13 ) u 4 ( 10 + 3 a 14 ) u 5 ( 27 + 2 a 15 ) u 6 ( 25 + a 16 ) + s 1 2 a 11 - 5 , ( 14 - 1 ) , u 1 ( 7 + a 21 ) u 2 ( 6 + 2 a 22 ) u 3 ( 32 + 3 a 23 ) u 4 ( 4 + a 24 ) u 5 ( 20 + a 25 ) u 6 ( 85 + 5 a 26 ) + s 2 a 21 8 , ( 14 - 2 ) , w 1 ( 2 a ¯ 11 ) + w 2 ( 3 a ¯ 12 ) + w 3 ( 6 2 a ¯ 13 ) + w 4 ( 4 a ¯ 14 ) , + w 5 ( 3 2 a ¯ 15 ) + w 6 ( 4 a ¯ 16 ) + s 3 x 7 ( 1 + a ¯ 11 ) , ( 14 - 3 ) , 0 a ¯ 11 1 , 0 a ¯ 12 1 , 0 a ¯ 13 1 , ( 14 - 4 ) , 0 a ¯ 14 1 , 0 a ¯ 15 1 , 0 a ¯ 16 1 , ( 14 - 5 ) , 0 a 11 1 , 0 a 12 1 , 0 a 13 1 , ( 14 - 6 ) , 0 a 14 1 , 0 a 15 1 , 0 a 16 1 , ( 14 - 7 ) , 0 a 21 1 , 0 a 22 1 , 0 a 23 1 , ( 14 - 8 ) , 0 a 24 1 , 0 a 25 1 , 0 a 26 1 , ( 14 - 9 ) .
D M U 1 is evaluated DMU in this model.The parameter “ x 7 ”is the representation of shadow price.

2.4.1. A Linear Equivalent Model for Nonlinear DEA Transformation Model

Let A r j = μ r α r j , for 0 α r j 1 , so 0 A r j μ r ; and let B i j = ω i α ¯ i j , for 0 α ¯ i j 1 , so 0 B i j ω i .
The linear equivalent Model (15) of the nonlinear DEA transformation model is as follows:
Z j 0 15 = M i n . { x 7 0.000001 ( s 1 + s 2 + s 3 ) } , s . t . 3 x 1 2 A 11 5 x 2 2 A 12 30 x 3 4 A 13 x 4 10 3 A 14 x 5 27 2 A 15 x 6 25 A 16 + s 1 2 a 11 5 , ( 15 - 1 ) , x 1 7 A 21 x 2 6 2 A 22 x 3 32 3 A 23 x 4 4 A 24 , x 5 20 A 25 x 6 85 5 A 26 + s 2 a 21 8 ( 15 - 2 ) , x 1 2 B 11 + x 2 3 B 12 + x 3 6 2 B 13 + x 4 4 B 14 + x 5 3 2 B 15 + x 6 4 B 16 + s 3 x 7 1 + B 11 , ( 15 - 3 ) , A 11 u 1 0 ; A 12 u 2 0 ; A 13 u 3 0 , ( 15 - 4 ) , A 14 u 4 0 ; A 15 u 5 0 ; A 16 u 6 0 , ( 15 - 5 ) , A 21 u 1 0 ; A 22 u 2 0 ; A 23 u 3 0 , ( 15 - 6 ) , A 24 u 4 0 ; A 25 u 5 0 ; A 26 u 6 0 , ( 15 - 7 ) , B 11 w 1 0 ; B 12 w 2 0 ; B 13 w 3 0 , ( 15 - 8 ) , B 14 w 4 0 ; B 15 w 5 0 ; B 16 w 6 0 , ( 15 - 9 ) .
D M U 1 is evaluated DMU in this model. Model (15) is solved by Lingo 11 software and the optimal solution is shown in Table 6.
As seen from Table 6, D M U 1 - D M U 6 are not efficient and the efficiency is lower than that of the original interval model (compared with those in Table 2). The efficiencies of D M U 1 to D M U 6 present an increasing trend.

2.4.2. The Concept of Normal Uncertainty Distribution

Definition 6.
Assume an uncertain variable ξ obeys normal uncertain distribution [68]:
ϕ ( x ) = ( 1 + e x p ( ( π ( e x ) ) / ( 3 σ ) ) ) 1 , x .
Consequently, ϕ ( x ) obeys normal distribution, where e denotes the mean value, and σ denotes the variance (which differ from the expression in probability). Then, we call ξ obeys the normal distribution of e , σ as mean and variance, denoted by ξ N ( e , σ ) .
The figure of normal uncertainty distribution function can be seen in Figure 2.
Definition 7.
The inverse uncertainty distribution of normal uncertain variable N ( e , σ ) is as follows:
ϕ 1 ( α ) = e + ( 3 σ ) / ( π ) l n ( ( α ) / ( 1 α ) ) .

2.5. Correlation between Interval Distribution and Normal Uncertainty Distribution

Assume that an uncertain variable ξ obeys the normal distribution, i.e., ξ N ( e , σ ) , and the distribution function is
ϕ ( x ) = ( 1 + e x p ( π ( e x ) ) / ( 3 σ ) ) ) 1 , x , e , σ R , σ > 0 .
The inverse function is as follows:
ϕ 1 ( α ) = e + ( 3 σ ) / ( π ) l n ( ( α ) / ( 1 α ) ) , e , σ R , σ > 0 .
If the function obeys the standard normal distribution, the probability of the normal uncertainty distribution variable x falling within [ μ 4 σ , μ + 4 σ ] (we call it “ 4 σ ”criterion) is as follows:
ϕ ( e + 4 σ ) ϕ ( e 4 σ ) = ( 1 + e x p ( ( π ( e ( e + 4 σ ) ) ) / ( 3 σ ) ) 1 ( 1 + e x p ( ( π ( e ( e 4 σ ) ) ) / ( 3 σ ) ) 1 .
This is reorganized as follows:
ϕ ( e + 4 σ ) ϕ ( e 4 σ ) = ( 1 + e x p ( 4 π ) / ( 3 ) ) 1 ( 1 + e x p ( ( 4 π ) / ( 3 ) ) 1 ) ) = 99.86 % ,
i.e.,
M { μ 4 σ x μ + 4 σ } = ϕ ( 4 ) ϕ ( 4 ) = 99.86 % .
According to the “small probability event “principle, the probability of x falling outside [ μ 4 σ , μ + 4 σ ] is less than 0.16 % . In reality, we would consider the corresponding event happening as being almost impossible.
Since the intervals themselves originate from the results of fuzzy judgments or random sampling, we only know the upper and lower limits of the interval. It is difficult to determine the specific value, so the interval can be regarded as fuzzy numbers and can also be seen as uncertain variables. It is possible to replace the intervals with uncertain variables, taking into account the good nature of normal distribution, for an uncertain variable ξ that follows the normal distribution N ( e , σ ) , according to the 4 σ principle, its probability of falling within [ μ 4 σ , μ + 4 σ ] is 99.86 % . Therefore, we can use the constant variable ξ obeying normal distribution N ( e , σ ) to approximately replace the interval [ a , a + ] :
[ a , a + ] = [ μ 4 σ , μ + 4 σ ] , μ = ( a + a + ) / 2 , σ = ( a + a ) / 8 .

2.5.1. Uncertain DEA Model Based on Normal Distribution

Definition 8.
The equivalent uncertain normal distribution of the interval y ˜ r j [ y r j L , y r j U ] , x ˜ i j [ x i j L , x i j U ] is as follows:
x ˜ i j N o r m a l ( e i j , σ i j ) ,   y ˜ r j N o r m a l ( e ¯ r j , σ ¯ r j ) ,
where
e r j = ( y r j L + y r j U ) / 2 , σ r j = ( y r j U y r j L ) / 8
and
e ¯ i j = ( x i j L + x i j U ) / 2 , σ ¯ i j = ( x i j U x i j L ) / 8 .
In the DEA model with interval input and output variables, where y ˜ r j [ y r j L , y r j U ] , x ˜ i j [ x i j L , x i j U ] , respectively, we construct an equivalent uncertain DEA Model (16) with a certain input and the largest output as follows:
Z j 0 16 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . M { j = 1 n λ j y ˜ r j y ˜ r j 0 S r + } 1 α r j , r = 1 , 2 , s , ( 16 - 1 ) , M { j = 1 n λ j x ˜ i j + S i θ x ˜ i j 0 } 1 ( α ¯ i j ) , i = 1 , 2 , , m , ( 16 - 2 ) , x ˜ i j N o r m a l ( e ¯ i j , σ ¯ i j ) , y ˜ r j N o r m a l ( e r j , σ r j ) , ( 16 - 3 ) , j = 1 n λ j = 1 , λ j 0 , j = 1 , 2 , , n , ( 16 - 4 ) , 0 α r j 1 , 0 α ¯ i j 1 , ( 16 - 5 ) , S i 0 , r = 1 , 2 , , s , ( 16 - 6 ) , S r + 0 , i = 1 , 2 , , m , ( 16 - 7 ) .
D M U j 0 is evaluated DMU in this model. The equivalent deterministic form of the Constraint (16-1) in the model is as follows:
j = 1 n λ j ( e y r j + 3 σ y r j / π l n ( ( α r j ) / ( 1 α r j ) ) ) ( S r + ) ( e y r j 0 3 σ y r j 0 / π l n ( ( α r j ) / ( 1 α r j ) ) )
(Theorem 3 and Proof 0.3, see Supplementary Materials)
The Model (16)’s equivalent deterministic DEA Model (17) is as follows:
Z j 0 17 = M i n . { θ ( ε r = 1 s S r + + ε i = 1 m S i ) } , s . t . j = 1 n λ j ( e y r j + ( 3 σ y r j ) / π l n ( ( α r j ) / ( 1 α r j ) ) ) , ( S r + ) ( e y r j 0 ( 3 σ y r j 0 ) / π l n ( ( α r j ) / ( 1 α r j ) ) ) , ( 17 - 1 ) , j = 1 n λ j ( e ¯ i j ( 3 σ ¯ i j ) / π l n ( ( α ¯ i j ) / ( 1 α ¯ i j ) ) ) + S i θ ( e ¯ x i j 0 ( 3 σ ¯ x i j 0 ) / π l n ( ( α ¯ i j ) / ( 1 α ¯ i j ) ) ) , ( 17 - 2 ) , j = 1 n λ j = 1 , λ j 0 , j = 1 , 2 , , n , ( 17 - 3 ) , 0 α r j 1 , 0 α ¯ i j , 1 , ( 17 - 4 ) , S i 0 , r = 1 , 2 , , s , ( 17 - 5 ) , S r + 0 , i = 1 , 2 , , m , ( 17 - 6 ) .
D M U j 0 is evaluated DMU in this model.

2.5.2. A Numerical Example of DEA Model Based on Normal Distribution

Assuming that there are still six DMUs of the interval model, the input indicators x 1 , and the output indicators y 1 and y 2 are interval data in Table 1, which are subjected to normal distribution. The data are shown in Table 7.

2.5.3. A Numerical Example of Equivalent Non-Linear DEA Model for Normal Distribution

Taking input and output indicators of the DMUs obeying the normal distribution of the interval data (Table 7) as the numerical example, similar to the normal distribution of the DEA Model (17), nonlinear DEA Model (18) for linear distribution is constructed as follows:
Z j 0 18 = M i n . { x 7 0.000001 ( s 1 + s 2 + s 3 ) } , s . t . u 1 ( 4 + 0.034 B 11 ) u 2 ( 6 + 0.034 B 12 ) u 3 ( 32 + 0.138 B 13 ) u 4 ( 11.5 + 0.078 B 14 ) u 5 ( 28 + 0.034 B 15 ) u 6 ( 25.5 + 0.086 B 16 ) + s 1 0.034 B 11 4 , ( 18 - 1 ) , u 1 ( 7.5 + 0.0086 B 21 ) u 2 ( 7 + 0.034 B 22 ) u 3 ( 33.5 + 0.078 B 23 ) u 4 ( 4.5 + 0.0086 B 24 ) u 5 ( 20.5 + 0.0086 B 25 ) u 6 ( 87.5 + 0.2154 B 26 ) + s 2 0.0086 B 21 7 . 5 , ( 18 - 2 ) , w 1 ( 1 . 5 0.0086 B ¯ 11 ) + w 2 ( 2.5 0.0086 B ¯ 12 ) + w 3 ( 5 0.0034 B ¯ 13 ) + w 4 ( 3.5 0.0086 B ¯ 14 ) + w 5 ( 2 0.034 B ¯ 15 ) + w 6 ( 3.5 0.086 B ¯ 16 ) + s 3 x 7 ( 1.5 0.0086 B ¯ 11 ) , ( 18 - 3 ) , B 11 = l n ( ( a 11 ) / ( 1 a 11 ) ) , B 12 = l n ( ( a 12 ) / ( 1 a 12 ) ) , ( 18 - 4 ) , B 13 = l n ( ( a 13 ) / ( 1 a 13 ) ) , B 14 = l n ( ( a 14 ) / ( 1 a 14 ) ) , ( 18 - 5 ) , B 15 = l n ( ( a 15 ) / ( 1 a 15 ) ) , B 16 = l n ( ( a 16 ) / ( 1 a 16 ) ) , ( 18 - 6 ) , B 21 = l n ( ( a 21 ) / ( 1 a 21 ) ) , B 22 = l n ( ( a 22 ) / ( 1 a 22 ) ) , ( 18 - 7 ) , B 23 = l n ( ( a 23 ) / ( 1 a 23 ) ) , B 24 = l n ( ( a 24 ) / ( 1 a 24 ) ) , ( 18 - 8 ) , B 25 = l n ( ( a 25 ) / ( 1 a 25 ) ) , B 26 = l n ( ( a 26 ) / ( 1 a 26 ) ) , ( 18 - 9 ) , B ¯ 11 = l n ( ( a ¯ 11 ) / ( 1 a ¯ 11 ) ) , B ¯ 11 = l n ( ( a ¯ 12 ) / ( 1 a ¯ 12 ) ) , ( 18 - 10 ) , B ¯ 13 = l n ( ( a ¯ 13 ) / ( 1 a ¯ 13 ) ) , B ¯ 14 = l n ( ( a ¯ 14 ) / ( 1 a ¯ 14 ) ) , ( 18 - 11 ) , B ¯ 15 = l n ( ( a ¯ 15 ) / ( 1 a ¯ 15 ) ) , B ¯ 16 = l n ( ( a ¯ 16 ) / ( 1 a ¯ 16 ) ) , ( 18 - 12 ) .
D M U 1 is evaluated DMU in this model.

2.5.4. A Numerical Example of Determined DEA Model for Normal Distribution

Let C r j = μ r B r j , for α r j [ 0 , 1 ] , 1 α r j [ 0 , 1 ] , so B r j = l n ( ( α r j ) / ( 1 α r j ) ) , i.e., B i j have no upper and lower constraints. Similarly, C ¯ i j = ω i B ¯ i j , for α ¯ i j [ 0 , 1 ] , 1 α ¯ i j [ 0 , 1 ] , so B ¯ i j = l n ( ( α ¯ r j ) / ( 1 ( α ¯ r j ) ) ) , i.e., B ¯ i j have no upper and lower constraints. The equivalent certain DEA model that obeys normal distribution is transformed into a linear Model (19) as follows:
Z j 0 19 = M i n . { x 7 0.000001 ( s 1 + s 2 + s 3 ) } , s . t . x 1 4 0.034 C 21 x 2 6 0.034 C 22 x 3 32 0.138 C 23 x 4 11.5 0.078 C 24 x 5 28 0.034 C 25 x 6 25.5 0.086 C 26 s 1 0.034 C 21 4 , ( 19 - 1 ) , x 1 7.5 0.0086 C 31 x 2 7 0.034 C 32 x 3 33.5 0.078 C 33 x 4 4.5 0.0086 C 34 x 5 20.5 0.0086 C 35 x 6 87.5 0.2154 C 36 + s 2 0.0086 C 31 7.5 , ( 19 - 2 ) , x 1 1.5 0.0086 ( C ¯ 11 ) + x 2 2.5 0.0086 ( C ¯ 12 ) + x 3 5 0.0034 ( C ¯ 13 ) + x 4 3.5 0.0086 ( C ¯ 14 ) + x 5 2 0.034 ( C ¯ 15 ) + x 6 3.5 0.086 ( C ¯ 16 ) + s 3 x 7 1.5 0.0086 ( C ¯ 11 ) , ( 19 - 3 ) .
D M U 1 is evaluated DMU in this model. Model (19) is solved by Lingo 11 software and the optimal solution is shown in Table 8.
As can be seen from Table 8, D M U 1 D M U 6 are non-efficient and the efficiency is lower than that of the original interval model but higher than that of the linear model. The efficiency from D M U 1 to D M U 6 presents an increasing trend.

2.5.5. Results Analysis of Three Uncertain Intervals Distribution DEA Model

The relative efficiency figures of the original Models (7)–(9), the linear Model (15) and the normal distribution Model (19) are as follows.
The comparison analysis of the efficiency value in Figure 3 shows that efficiency sorting of DMUs is the same, but the value of the same DMU in different models is different. At the same time, the results of the normal distribution model are much better than those of the uncertain model. Therefore, we use the uncertain normal distribution model to conduct empirical analyses in the following section.

2.6. Research on the Manufacturing Output Performance of Yangtze River Delta Region Based on Environmental Pollution Constraints

For this section, we select input–output indicators of the manufacturing industry in 16 cities in the Yangtze River Delta region, China, from 2005 to 2014 to conduct relative efficiency empirical analyzes and construct uncertain DEA models based on normal distribution. We assess and rank these cities according to the relative efficiency value and conduct a corresponding analysis. Furthermore, we validate the practicability and rationality of the uncertain normal distribution DEA model.

2.6.1. Concept Definition

The geographical range of the Yangtze River Delta can be measured in narrow and generalized terms. In the narrow sense, the Yangtze River Delta encompasses 16 member cities of the Yangtze River Delta Urban Economic Coordination Committee, including the cities of Shanghai, Nanjing, Wuxi, Changzhou, Suzhou, Nantong, Yangzhou, Zhenjiang and Taizhou in Jiangsu Province and Hangzhou, Ningbo, Jiaxing, Huzhou, Shaoxing, Zhoushan and Taizhou in Zhejiang Province. In the general sense, the Yangtze River Delta refers to all regions in Shanghai, Jiangsu Province and Zhejiang Province. This study takes the narrow sense as research objectives.
Environmental quality is an intrinsic property that exists objectively in an environmental system. The deterioration of environmental quality includes ecological damage and environmental pollution. The destruction of the ecological environment is usually due to human activities, which directly or indirectly damage biological communities and non-biological elements ecosystems, and ultimately leads to the deterioration of the human living environment. Environmental pollution is the discharge of pollutants from human activities into the environment, exceeding the capacity and the self-purification capacity of the environment, so that the composition of the environment is altered, leading to the deterioration of environmental quality, affecting and undermining people’s normal production and living conditions. Pollutants mainly refer to substances that are discharged into the atmosphere, water and soil. In this paper, changes in environmental quality mainly refer to environmental pollution with less consideration of eco-environment content, so that the content of the study is targeted.
Manufacturing refers to the industry that transforms manufacturing resources (materials, energy, equipment, tools, funds, technology, information and manpower), in accordance with market demand through the manufacturing process, into large tools and useful industrial and consumer products. Manufacturing embodies a country’s productivity level, which is an important factor in distinguishing between developing and developed countries. It occupies an important share in the national economy of the developed world. The manufacturing industry referred to in this paper starts with the agro-food processing industry and ends with metal products, machinery and equipment repairing (including 31 industries: agricultural and foodstuff processing industry; food manufacturing; liquor, beverage and refined tea manufacturing; the tobacco industry; textiles; leather, fur, feather and garment industry; wood processing and footwear; bamboo, rattan, palm, straw products; furniture manufacturing; printing and recording media reproduction industry; education, engineering, sports and entertainment products manufacturing; petroleum processing, coking and nuclear fuel processing; chemical raw materials and chemical manufacturing; pharmaceutical manufacturing; rubber and plastics products industry; nonmetal mineral products industry; ferrous metals smelting and calendering processing; nonferrous metals smelting and calendering processing; metal products industry; general equipment manufacturing; special equipment manufacturing; automotive manufacturing; railways, ships, aerospace and other transport equipment manufacturing; electrical machinery and equipment manufacturing; computer, communications and other electronic equipment manufacturing; instrumentation manufacturing; other manufacturing industries; waste resources comprehensive utilization industry; metal products, machinery and equipment repair industry).

2.6.2. Indicators Selection and Data Resources

Referring to the literature [69,70,71,72,73], this paper constructs environmental indicators shown in Table 9. With regards to these indicators, water pollution is characterized by the emission of waste water and air pollution is characterized by soot and sulphur dioxide emission. The treatment of solid wastes is characterized by comprehensive utilization of waste and, lastly, energy consumption is characterized by comprehensive energy consumption.
According to the actual development situation and data feasibility of the manufacturing industry in the Yangtze River Delta, this paper studies the output–input performance of this manufacturing industry based on the constraint of environmental pollution. The 16 cities in the Yangtze River Delta include Shanghai, eight cities in Jiangsu Province, namely Nanjing, Wuxi, Changzhou, Yangzhou, Suzhou, Zhenjiang, Nantong and Xuzhou, and seven cities in Zhejiang Province, namely Hangzhou, Ningbo, Shaoxing, Jiaxing, Zhoushan, Quzhou and Taizhou.
The data in this paper are gathered from the annual statistic yearbook of each city. The waste utilization rate, waste water emission, soot emission and sulphur dioxide emission are collected with industry data calibre. The remaining indicators are collected with manufacturing calibre (Shanghai comprehensive energy consumption uses industry-calibre), minimizing the regional differences due to factors such as population and administrative area. Except for the comprehensive utilization of waste, the remaining indicators are based on the total population of each city as the denominator for processing. All monetary data is deflated in real prices based on constant price value index in 2013.

2.6.3. Indicator Elaboration

This study selects average annual number of employees (AANE), completed investment in fixed assets (CIFA), and comprehensive energy consumption (CEC) of each city, industrial waste water discharge (IWWD), industrial soot emission (ISE), industrial sulphur dioxide emission (ISD) as the input indicators and comprehensive utilization of waste(CUW), and manufacturing gross output value (GOV) as the output indicators [74,75,76,77].

2.7. Empirical Data

We process the data collected from 2005 to 2014 and the corresponding means and variances are transformed into normal distribution data, as shown in Table 10, Table 11, Table 12 and Table 13.

2.7.1. Empirical Model

We take the interval data (Table 10, Table 11, Table 12 and Table 13) of the output indicators of the DMUs that obey the normal distribution as the empirical data, which resemble the deterministic DEA Model (17) of normal distribution. The equivalent nonlinear equilibrium output Model (20) is constructed as follows:
Z j 0 20 = M i n . { x 7 0.000001 ( s 1 + s 2 + s 3 + s 4 ) } , s . t . x 1 ( 96 + A 71 ) x 2 ( 87 + 4 A 72 ) x 3 ( 95 + 2 A 73 ) x 4 ( 98 + A 74 ) x 5 ( 94 + 4 A 75 ) x 6 ( 98 + 0.2 A 76 ) x 7 ( 92 + 5.3 A 77 ) x 8 ( 97 + 2 A 78 ) x 9 ( 98 + 2 A 79 ) x 10 ( 94 + A 7.10 ) x 11 ( 88 + 2 A 7.11 ) x 12 ( 88 + 5 A 7.12 ) x 13 ( 96 + 2 A 7.13 ) x 14 ( 98 + A 7.14 ) x 15 ( 96 + A 7.15 ) x 16 ( 95 + A 7.16 ) + s 1 96 + 1 A 71 x 1 ( 25190 + 3236 A 81 ) x 2 ( 8199 + 1779 A 82 ) x 3 ( 11238 + 1818 A 83 ) x 4 ( 6603 + 1589 A 84 ) x 5 ( 5038 + 1543 A 85 ) x 6 ( 21519 + 4073 A 86 ) x 7 ( 4040 + 1277 A 87 ) x 8 ( 6888 + 1920 A 88 ) x 9 ( 5009 + 2021 A 89 ) x 10 ( 9478 + 1333 A 8.10 ) x 11 ( 9098 + 1557 A 8.11 ) x 12 ( 6188 + 1213 A 8.12 ) x 13 ( 4360 + 915 A 8.13 ) x 14 ( 844 + 227 A 8.14 ) x 15 ( 950 + 241 A 8.15 ) x 16 ( 2894 + 374 A 8.16 ) + s 2 25190 + 3236 A 81 .
D M U 1 is evaluated DMU in this model. We subsequently take the interval data (Table 10, Table 11, Table 12 and Table 13) of the input indicators of the DMU that obeyed the normal distribution as the empirical data, which resemble the deterministic DEA Model (17) with normal distribution. The equivalent nonlinear equilibrium input Model (21) is constructed as follows:
s . t . x 1 ( 105586 5047 A 11 ) + x 2 ( 153931 46165 A 12 ) + x 3 ( 110321 18084 A 13 ) + x 4 ( 102264 24189 A 14 ) + x 5 ( 69503 19505 A 15 ) + x 6 ( 152959 28544 A 16 ) + x 7 ( 56470 17558 A 17 ) + x 8 ( 118301 28312 A 18 ) + x 9 ( 88585 30323 A 19 ) + x 10 ( 58528 8477 A 1.10 ) + x 11 ( 67685 11423 A 1.11 ) + x 12 ( 64686 12283 A 1.12 ) + x 13 ( 58441 9775 A 1.13 ) + x 14 ( 11634 2686 A 1.14 ) + x 15 ( 20687 4221 A 1.15 ) + x 16 ( 35372 6469 A 1.16 ) + s 3 x 7 ( 105586 5047 A 11 ) x 1 ( 26847 993 A 21 ) + x 2 ( 6923 584.6 A 22 ) + x 3 ( 12475 733 A 23 ) + x 4 ( 9644 + 1936 A 24 ) + x 5 ( 6258 795 A 25 ) + x 6 ( 29314 3079 A 26 ) + x 7 ( 4535 434 A 27 ) + x 8 ( 7884 942 A 28 ) + x 9 ( 5075 879 A 29 ) + x 10 ( 11980 516 A 2.10 ) + x 11 ( 15843 801 A 2 . 11 ) + x 12 ( 7994 280 A 2.12 ) + x 13 ( 8637 338 A 2.13 ) + x 14 ( 1752 104 A 2.14 ) + x 15 ( 1336 121 A 2.15 ) + x 16 ( 8061 1873 A 2.16 ) + s 4 x 7 ( 26847 993 A 21 ) x 1 ( 10094 629 A 31 ) + x 2 ( 7070 570 A 32 ) + x 3 ( 2020 128 A 33 ) + x 4 ( 15854 + 196 A 34 ) + x 5 ( 606 69 A 35 ) + x 6 ( 3077 257 A 36 ) + x 7 ( 578 57 A 37 ) + x 8 ( 694 55 A 38 ) + x 9 ( 1125 109 A 39 ) + x 10 ( 1970 88 A 3.10 ) + x 11 ( 1308 175 A 3.11 ) + x 12 ( 1195 76 A 3.12 ) + x 13 ( 912 128 A 3.13 ) + x 14 ( 62 24 A 3.14 ) + x 15 ( 666 40 A 3.15 ) + x 16 ( 168 7 A 3.16 ) + s 5 x 7 ( 10094 629 A 31 ) x 1 ( 221800 7898 A 41 ) + x 2 ( 33390 4964 A 42 ) + x 3 ( 36454 6652 A 43 ) + x 4 ( 28174 + 5445 A 44 ) + x 5 ( 9652 908 A 45 ) + x 6 ( 66976 3596 A 46 ) + x 7 ( 10167 958 A 47 ) + x 8 ( 16250 1204 A 48 ) + x 9 ( 9746 1475 A 49 ) + x 10 ( 63888 10972 A 4.10 ) + x 11 ( 17558 1354 A 4.11 ) + x 12 ( 28311 1839 A 4.12 ) + x 13 ( 18886 1596 A 4.13 ) + x 14 ( 1876 115 A 4.14 ) + x 15 ( 14209 793 A 4.15 ) + x 16 ( 5556 356 A 4.16 ) + s 6 x 7 ( 221800 7898 A 41 ) x 1 ( 58470 15487 A 51 ) + x 2 ( 48790 10931 A 52 ) + x 3 ( 48516 4426 A 53 ) + x 4 ( 38306 + 15207 A 54 ) + x 5 ( 12879 1704 A 55 ) + x 6 ( 50021 9908 A 56 ) + x 7 ( 27546 7248 A 57 ) + x 8 ( 40528 3253 A 58 ) + x 9 ( 32752 6556 A 59 ) + x 10 ( 47515 15058 A 5.10 ) + x 11 ( 26425 2031 A 5.11 ) + x 12 ( 29885 12655 A 5.12 ) + x 13 ( 24113 2222 A 5.13 ) + x 14 ( 12252 3256 A 5.14 ) + x 15 ( 666 40 A 5.15 ) + x 16 ( 10486 1802 A 5.16 ) + s 7 x 7 ( 58470 15487 A 51 ) x 1 ( 264620 46622 A 61 ) + x 2 ( 128380 8251 A 62 ) + x 3 ( 104328 13394 A 63 ) + x 4 ( 56916 + 10867 A 64 ) + x 5 ( 67284 11246 A 65 ) + x 6 ( 197259 20239 A 66 ) + x 7 ( 74044 9279 A 67 ) + x 8 ( 75753 7911 A 68 ) + x 9 ( 10918 11980 A 69 ) + x 10 ( 89183 14146 A 6.10 ) + x 11 ( 150393 19314 A 6.11 ) + x 12 ( 63048 3237 A 6.12 ) + x 13 ( 89829 10235 A 6.13 ) + x 14 ( 20768 3589 A 6.14 ) + x 15 ( 35248 6020 A 6.15 ) + x 16 ( 48820 10539 A 6.16 ) + s 8 x 7 ( 264620 46622 A 61 ) s 1 + s 2 + s 3 + s 4 + s 5 + s 6 + s 7 + s 8 = 1 .
The constraints of the empirical model are as follows:
s . t . A i 1 = l n ( ( a i 1 ) / ( 1 a i 1 ) ) , A i 2 = l n ( ( a i 2 ) / ( 1 a i 2 ) ) A i 3 = l n ( ( a i 3 ) / ( 1 a i 3 ) ) , A i 4 = l n ( ( a i 4 ) / ( 1 a i 4 ) ) A i 5 = l n ( ( a i 5 ) / ( 1 a i 5 ) ) , A i 6 = l n ( ( a i 6 ) / ( 1 a i 6 ) ) A i 7 = l n ( ( a i 7 ) / ( 1 a i 7 ) ) , A i 8 = l n ( ( a i 8 ) / ( 1 a i 8 ) ) A i 9 = l n ( ( a i 9 ) / ( 1 a i 9 ) ) , A ( i 10 ) = l n ( ( a ( i . 10 ) / ( 1 a ( i . 10 ) ) ) ) A ( i . 11 ) = l n ( ( a ( i . 11 ) / ( 1 a ( i . 11 ) ) ) ) A ( i . 12 ) = l n ( ( a ( i . 12 ) / ( 1 a ( i . 12 ) ) ) ) A ( i . 13 ) = l n ( ( a ( i . 13 ) / ( 1 a ( i . 13 ) ) ) ) A ( i . 14 ) = l n ( ( a ( i . 14 ) / ( 1 a ( i . 14 ) ) ) ) A ( i . 15 ) = l n ( ( a ( i . 15 ) / ( 1 a ( i . 15 ) ) ) ) A ( i . 16 ) = l n ( ( a ( i . 16 ) / ( 1 a ( i . 16 ) ) ) ) , i = 1 , 2 , , 9 .
By combining Models (20)–(22), we obtain the equivalent nonlinear equilibrium DEA model of normal distribution.

2.7.2. The Results of the Empirical Model

From Figure 4 and Table 14, the following can be seen:
  • Under the constraint of environmental pollution, the top four manufacturing input–output performance cities are found to be Yangzhou, Taizhou, Nantong and Shanghai, of which Yangzhou and Taizhou’s energy consumption and negative output are lower, but the per capital positive outputs are insufficient, belonging to the low-energy-consumption, low-pollution and low-output industrial structure. The energy consumption and negative output of Nantong and Shanghai are higher, but the per capital positive outputs have advantages in that they belong to the industrial structure with high energy, high pollution and high yield. However, Shanghai is at a stage of industrial transformation and the focus gradually shift from traditional to high-technology manufacturing. Nonetheless, it is still in the leading position in terms of scale and technology, so that the manufacturing output performance advantage is maintained under environmental pollution constraints.
  • Under the constraint of environmental pollution, Jiaxing, Quzhou and Shaoxing rank as the last three regions in manufacturing input–output performance. Their inputs and negative outputs are high, while the positive outputs are at a disadvantage, belonging to the high-input, high-pollution and low-output industrial structure.
  • Under the restriction of environmental pollution, although their efficiency is higher, Ningbo and Zhoushan have a slight advantage in terms of positive outputs, but the negative outputs and environmental costs of manufacturing development are higher and these regions consequently belonged to the high-pollution development model.
  • The Changzhou manufacturing development model is different from that of other cities. Regarding inputs, energy consumption and positive or negative outputs, Changzhou’s is significantly higher than those of other cities. Although its performance is top-ranking, the development model’s high yield is accompanied by high input, high energy consumption and high pollution. Therefore, such a model is unhealthy.
  • In general, except for Taizhou and Zhoushan, from the aspect of environmental pollution control of manufacturing input–output performance, cities in Zhejiang Province rank lower than other cities; there is a gap in the development level of the manufacturing industry between Shanghai and Jiangsu Province.
  • As for the Yangtze River coast, except for Shanghai, Yangzhou and Nantong, the manufacturing output performance under the constraint of environmental pollution of regions in the south is lower than that of regions in the north. Although the manufacturing levels in the south are high, the environmental costs are also higher.

3. Conclusions

In practice, input and output variables are usually interval data in the complex efficiency evaluation. The traditional approach is to consider the upper and lower limits of the interval data and many researchers consequently constructed models with upper and lower limits for evaluation, which inevitably leads to the lack of some important information. This study attempted to overcome these defects by studying the problems of input–output efficiency evaluation with interval data through two research angles.
The first perspective did not consider the distribution of interval data. Introducing the tuning parameters, α, the interval data were transformed into determination data. By using optimization methods, the input–output performance evaluation value was acquired. The advantage of this method was that the overall interval information could be considered more comprehensive and was conducive to the improvement of the non-efficient decision units.
The second perspective did consider the distribution characteristics of interval data. Based on the uncertainty theory proposed by Liu [68], the DEA evaluation model with interval data was constructed from dimensions of linear uniform distribution and normal distribution with uncertainty:
i
Linear uniform distribution model: According to the uncertain optimization theory, the DEA uncertainty evaluation model for linear uniform distribution was transformed into a deterministic nonlinear DEA evaluation model, which was further transformed into a deterministic linear DEA evaluation model.
ii
Uncertainty model of normal distribution: Based on uncertain optimization theory, we adopted the ‘ 4 σ ’ rule. First, we transformed the input–output interval indicators into an approximate normal distribution of the input–output DEA evaluation efficiency uncertainty optimization model. Then, we transformed it into an equivalent normal distribution of input–output DEA deterministic evaluation model.
iii
The advantage of these two kinds of transformation is that it can further reduce the computational cost of the model optimization solution and obtain the optimal solution of the uncertain model.
iv
The results of empirical analysis showed that the optimal solution of the normal distribution model was better than that of the linear distribution. Therefore, based on the normal distribution DEA evaluation model, this paper studied the input–output performance evaluation of the manufacturing industry in the Yangtze River Delta under the constraint of environmental pollution.
The paper focused on the regions of 16 member cities of Yangtze River Delta Regional Urban Economic Coordination Committee. We collected the manufacturing data of these cities from 2005 to 2014. According to the means and variances of these data, they were transformed into approximate normal distribution data indicators. Using the previously constructed normal distribution input–output DEA evaluation efficiency deterministic optimization model, we analyzed each city’s sorted input–output performance results. We consequently offered corresponding countermeasures and suggestions according to the differences in industrial structure in different cities.

Supplementary Materials

The following are available online at www.mdpi.com/2071-1050/9/12/2080/s1.

Acknowledgments

The authors would like to thank the editor and all anonymous referees for their valuable comments and suggestions. This research was supported by the National Natural Science Foundation of China ( 71171115 , 71571104 ) , the Six Talent Peaks Project in Jiangsu Province (2014-JY-014), and the Natural Science Foundation of Jiangsu ( B K 20141481 ) , ’333 Project’ research projects funded project (BRA2017456).

Author Contributions

Xiaoqing Chen was responsible for the model development and first draft, and finalized the manuscript. Zaiwu Gong conceived the main idea of the experiments and supervised the whole process of model development and manuscript drafting. All coauthors made significant contributions to the research contained in this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DEAData Envelopment Analysis
DMUdecision-making units
YRDYangtze River Delta
TPTotal Population
CUWcomprehensive utilization of waste
AANEaverage annual number of employees
CIFAcompleted investment in fixed assets
CECcomprehensive energy consumption
IWWDindustrial waste water discharge
ISEindustrial soot emission
ISDindustrial sulphur dioxide emission
GOVmanufacturing gross output value

References

  1. Fujii, H.; Managi, S. Wastewater management efficiency and determinant factors in the Chinese industrial sector from 2004 to 2014. Water 2017, 9, 586. [Google Scholar] [CrossRef]
  2. Fujii, H.; Cao, J.; Managi, S. Decomposition of Productivity Considering Multi-environmental Pollutants in Chinese Industrial Sector. Rev. Dev. Econ. 2015, 19, 75–84. [Google Scholar] [CrossRef]
  3. Fujii, H.; Cao, J.; Managi, S. Firm-level environmentally sensitive productivity and innovation in China. Appl. Energy 2016, 184, 915–925. [Google Scholar] [CrossRef]
  4. Yang, Q.; Kaneko, S.; Fujii, H.; Yoshida, Y. Do exogenous shocks better leverage the benefits of technological change in the staged elimination of differential environmental regulations? Evidence from China’s cement industry before and after the 2008 Great Sichuan Earthquake. J. Clean. Prod. 2017, 164, 1167–1179. [Google Scholar] [CrossRef]
  5. 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]
  6. Lu, K.; Nie, C.L. Application of Data Envelopment Analysis in Analysis of Maintenance Support Efficiency. Mod. Def. Technol. 2017, 1, 167–172. [Google Scholar]
  7. Li, J.J. Research on innovation performance of Listed Companies in Anhui Province—Based on DEA analysis. J. Chifeng Univ. 2016, 24, 78–80. [Google Scholar]
  8. Huang, H.S.; Yang, H.G. An analysis of the effectiveness of investment in Chengdu based on DEA model. China Mark. China Mark. 2016, 47, 22–25. [Google Scholar]
  9. Tan, J.; Luo, Z.Y.; Xu, G.W. Evaluation on Input Performance of Technological Innovation Based on CCR, BCC and SE-DEA Models. Sci. Technol. Econ. 2017, 1, 36–40. [Google Scholar]
  10. Kwon, D.S.; Cho, J.H.; Sohn, S.Y. Comparison of technology efficiency for CO2 emissions reduction among European countries based on DEA with decomposed factors. J. Clean. Prod. 2017, 151, 109–120. [Google Scholar] [CrossRef]
  11. Toloo, M. On finding the most BCC-efficient DMU: A new integrated MIP-DEA model. Appl. Math. Model. 2012, 36, 5515–5520. [Google Scholar] [CrossRef]
  12. Guo, I.L.; Lee, H.S.; Lee, D. An integrated model for slack-based measure of super-efficiency in additive DEA. Omega 2017, 67, 160–167. [Google Scholar] [CrossRef]
  13. Gutierrez, E.; Aguilera, E.; Lozano, S.; Guzmán, G.I. A two-stage DEA approach for quantifying and analyzing the inefficiency of conventional and organic rain-fed cereals in Spain. J. Clean. Prod. 2017, 149, 335–348. [Google Scholar] [CrossRef]
  14. Li, L.B.; Liu, B.L.; Liu, W.L.; Chiu, Y.H. Efficiency evaluation of the regional high-tech industry in China: A new framework based on meta-frontier dynamic DEA analysis. Socio-Econ. Plan. Sci. 2017. [Google Scholar] [CrossRef]
  15. Aryana, B. New version of DEA compressor for a novel hybrid gas turbine cycle: TurboDEA. Energy 2016, 111, 676–690. [Google Scholar] [CrossRef]
  16. Zhu, J.X.; Sun, Y.H.; Yang, L.M. Research on Stochastic DEA Model for Undesirable Outputs Evaluation. J. Qual. Econ. 2015, 3, 73–77. [Google Scholar]
  17. Shwartz, M.; Burgess, J.F.; Zhu, J. A DEA based composite measure of quality and its associated data uncertainty interval for health care provider profiling and pay-for-performance. Eur. J. Oper. Res. 2016, 253, 489–502. [Google Scholar] [CrossRef]
  18. Xiang, X.D.; Zhang, Z.Y. Interval DEA supply chain partner evaluation model. Stat. Decis. 2016, 21, 55–58. [Google Scholar]
  19. Chi, G.T.; Du, Y.Q.; He, X.G.; Liu, J.B. Research on Loan Pricing Model Based on the Interval DEA. Oper. Res. Manag. Sci. 2015, 5, 189–196. [Google Scholar]
  20. Fan, J.P.; Yue, W.Z.; Wu, M.Q. Dealing with interval DEA based on error propagation and entropy. Syst. Nato Pract. 2016, 5, 1293–1303. [Google Scholar]
  21. Chen, X.Q.; Gong, Z.W. DEA Efficiency of Energy Consumption in China’s Manufacturing Sectors with Environmental Regulation Policy Constraints. Sustainability 2017, 9, 210. [Google Scholar] [CrossRef]
  22. Olfat, L.; Amiri, M.; Soufi, J.B.; Pishdar, M. A dynamic network efficiency measurement of airports performance considering sustainable development concept: A fuzzy dynamic network-DEA approach. J. Air Transp. Manag. 2016, 57, 272–290. [Google Scholar] [CrossRef]
  23. Aggelopoulos, E.; Georgopoulos, A. Bank branch efficiency under environmental change: A bootstrap DEA on monthly Profit and Loss accounting statements of Greek retail branches. Eur. J. Oper. Res. 2017, 261, 1170–1188. [Google Scholar] [CrossRef]
  24. Vlontzos, G.; Pardalos, P.M. Assess and prognosticate green house gas emissions from agricultural production of EU countries, by implementing, DEA Window analysis and artificial neural networks. Renew. Sustain. Energy Rev. 2017, 76, 155–162. [Google Scholar] [CrossRef]
  25. Cooper, W.W.; Huang, Z.; Lelas, V.; Li, S.X.; Olesen, O.B. Chance constrained programming formulations for stochastic characterizations of efficiency and dominance in DEA. J. Prod. Anal. 1998, 9, 53–79. [Google Scholar] [CrossRef]
  26. Ding, Q.Y.; Wang, L. Empirical study on scale economy of securities companies: Based on methods of DEA and translog cost function. J. Shenyang Univ. Technol. 2014, 6, 529–534. [Google Scholar]
  27. Zhang, Z.H. Influencing Factors and Performance Analysis of Land Transfer in Shanxi; Northwest University: Xian, China, 2012. [Google Scholar]
  28. Duan, Z.P. Analysis on Energy Utilization-Efficiency of the Main Industries in Jiangxi Province; Nanchang University: Nanchang, China, 2010. [Google Scholar]
  29. Alzamora, R.M.; Apiolaza, L.A. A DEA approach to assess the efficiency of radiata pine logs to produce New Zealand structural grades. J. For. Econ. 2013, 19, 221–233. [Google Scholar] [CrossRef]
  30. Lin, B.; Xie, C. Energy substitution effect on transport industry of China-based on trans-log production function. Energy 2014, 67, 213–222. [Google Scholar] [CrossRef]
  31. Bounneche, M.D.; Boubchir, L.; Bouridane, A.; Nekhoul, B.; Ali-Chérif, A. Multi-spectral palmprint recognition based on oriented multiscale log-Gabor filters. Neurocomputing 2016, 205, 274–286. [Google Scholar] [CrossRef]
  32. Zhang, H.X. Study on chance constrained stochastic DEA model based on exponential distribution. Univ. Electron. Sci. Technol. China 2016, 8, 356–369. [Google Scholar]
  33. Peng, Y.; Chen, S.Y.; Sheng, W.W.; Pu, Y. Efficiency Based on the Two-stage Correlative DEA Malmquist Index on Western Regional Technical Innovation. Oper. Res. Manag. Sci. 2013, 3, 162–168. [Google Scholar]
  34. Li, G.C.; Fan, L.X.; Cheng, G.; Feng, Z.C. Total factor productivity growth in agriculture: A new window based DEA productivity index. J. Agrotechnol. 2013, 5, 4–17. [Google Scholar]
  35. Wang, W.G.; Ma, Y.Y. The Efficiency of Regional Logistics Industry in China-Based on the Three-stage DEA Model Using Malmquist-luenberger Index. Syst. Eng. 2012, 3, 66–75. [Google Scholar]
  36. Dai, X.; Kuosmanen, T. Best-practice benchmarking using clustering methods: Application to energy regulation. Omega 2014, 42, 179–188. [Google Scholar] [CrossRef]
  37. Jeang, A. Robust DEA methodology via computer model for conceptual design under uncertainty. J. Intell. Manuf. 2017. [Google Scholar] [CrossRef]
  38. Yuan, J.; Wu, Z.W. Research on Transfer Efficiency of High-Speed Rail Station Based on DEA Model. Luoyang Inst. Sci. Technol. 2017, 1, 48–51. [Google Scholar]
  39. Li, S.H.; Xiong, C. Analysis on the Efficiency and Public Satisfaction of Local Government Fiscal Expenditure on Social Security. Chin. Public Admin. 2016, 2, 104–111. [Google Scholar]
  40. Hu, X.W.; Wei, Y.B. Evaluation of Urban Public Transport Service Satisfaction Degree in Cold Regions’ Winter Based on DEA Method. J. Beijing Univ. Technol. 2015, 10, 1566–1573. [Google Scholar]
  41. Dai, Q.Z.; Lei, X.Y.; Li, Y.J.; Xie, Q.W. Allocating emissions based on DEA and maximizing the whole satisfaction degree. Syst. Nato Pract. 2014, 4, 917–924. [Google Scholar]
  42. Shabanpour, H.; Yousefi, S.; Saen, R.F. Future planning for benchmarking and ranking sustainable suppliers using goal programming and robust double frontiers DEA. Transp. Res. Part D Transp. Environ. 2017, 50, 129–143. [Google Scholar] [CrossRef]
  43. Nazari-Shirkouhi, S.; Keramati, A. Modeling customer satisfaction with new product design using a flexible fuzzy regression-DEA algorithm. Appl. Math. Model. 2017, 10, 1566–1573. [Google Scholar]
  44. Xi, W.; Li, A. The measurement of industrial capital service and the difference of production efficiency. Price Theory Pract. 2016, 12, 155–158. [Google Scholar]
  45. Liang, S.; Zeng, X.T.; Li, Y.P.; Qiao, X.L. The Effect Evaluation of Operation and Management on Soil and Water Conservation Project Based on Entropy Fuzzy Credibility-constrained DEA Method. China Rural Water Hydropower 2015, 8, 94–97. [Google Scholar]
  46. Wang, M.Q.; Li, Y.J. Non-radial Fuzzy DEA Model Based on Credibility Measure. Fuzzy Syst. Math. 2010, 6, 110–116. [Google Scholar]
  47. Fasanghari, M.; Amalnick, M.S.; Anvari, R.T.; Razmi, J. A novel credibility-based group decision making method for Enterprise Architecture scenario analysis using Data Envelopment Analysis. Appl. Soft Comput. 2015, 32, 347–368. [Google Scholar] [CrossRef]
  48. Azadeh, A.; Kokabi, R. Z-number DEA: A new possibilistic DEA in the context of Z-numbers. Adv. Eng. Inform. 2016, 30, 604–617. [Google Scholar] [CrossRef]
  49. Salvador-Carulla, L.; Torres-Jiménez, M.; GarcÍa-Alonso, C.R.; Fernández-Rodríguez, V. Evaluation of system efficiency using the Monte Carlo DEA: The case of small health areas. Eur. J. Oper. Res. 2015, 242, 525–535. [Google Scholar]
  50. Perelman, S.; SantÍn, D. How to generate regularly behaved production data? A Monte Carlo experimentation on DEA scale efficiency measurement. Eur. J. Oper. Res. 2009, 199, 303–310. [Google Scholar] [CrossRef]
  51. Cordero, J.M.; Santín, D.; Sicilia, G. Testing the accuracy of DEA estimates under endogeneity through a Monte Carlo simulation. Eur. J. Oper. Res. 2015, 244, 511–518. [Google Scholar] [CrossRef]
  52. Oh, S.C.; Shin, J. The impact of mismeasurement in performance benchmarking: A Monte Carlo comparison of SFA and DEA with different multi-period budgeting strategies. Eur. J. Oper. Res. 2015, 240, 518–527. [Google Scholar] [CrossRef]
  53. Andor, M.; Hesse, F. The StoNED age: The departure into a new era of efficiency analysis ? A monte carlo comparison of StoNED and the ”oldies”(SFA and DEA). J. Prod. Anal. 2014, 41, 85–109. [Google Scholar] [CrossRef]
  54. Azadeh, A.; Ahvazi, M.P.; Haghighii, S.M.; Keramati, A. Simulation optimization of an emergency department by modeling human errors. Simul. Model. Pract. Theory 2015, 67, 117–136. [Google Scholar] [CrossRef]
  55. Wu, Q. An Applied Study of Financial Distress Warning of Listed Companies Base on Logostic Model and Bayesian Network; Jilin University: Changchun, China, 2011. [Google Scholar]
  56. Jia, J.X. Dynamic selection method of manufacturers to customers—Application of DEA method and Bayesian correction method. Mod. Manag. Sci. 2006, 4, 39–40. [Google Scholar]
  57. Tsionas, E.G.; Papadakis, E.N. A Bayesian approach to statistical inference in stochastic DEA. Omega 2010, 38, 309–314. [Google Scholar] [CrossRef]
  58. Chiang, T.A.; Che, Z.H. A fuzzy robust evaluation model for selecting and ranking NPD projects using Bayesian belief network and weight-restricted DEA. Expert Syst. Appl. 2010, 37, 7408–7418. [Google Scholar] [CrossRef]
  59. Tsionas, E.G.; Papadakis, E.N. Combining DEA and stochastic frontier models: An empirical Bayes approach. Eur. J. Oper. Res. 2003, 147, 499–510. [Google Scholar] [CrossRef]
  60. Mitropoulos, P.; Talias, M.A.; Mitropoulos, I. Combining stochastic DEA with Bayesian analysis to obtain statistical properties of the efficiency scores: An application to Greek public hospitals. Eur. J. Oper. Res. 2015, 243, 302–311. [Google Scholar] [CrossRef]
  61. Khodabakhshi, M.; Asgharian, M. An input relaxation measure of efficiency in stochastic data envelopment analysis. Appl. Math. Model. 2009, 33, 2010–2023. [Google Scholar] [CrossRef]
  62. Eslami, R.; Khodabakhshi, M.; Jahanshahloo, G.R.; Lotfi, F.H.; Khoveyni, M. Estimating most productive scale size with imprecise-chance constrained input-output orientation model in data envelopment analysis. Comput. Ind. Eng. 2012, 63, 254–261. [Google Scholar] [CrossRef]
  63. Zhou, X.; Pedrycz, W.; Kuang, Y.; Zhang, Z. Type-2 fuzzy multi-objective DEA model: An application to sustainable supplier evaluation. Appl. Soft Comput. 2016, 46, 424–440. [Google Scholar] [CrossRef]
  64. Yu, X.W.; Lei, M.; Wang, Q.W.; Deng, J. Chinese Commercial Bank’s Efficiency (2005–2011): Based on Time-Series Regression and Chance-constrained DEA model. China J. Manag. Sci. 2012, S1, 356–362. [Google Scholar]
  65. Fujii, H.; Managi, S. An evaluation of inclusive capital stock for urban planning. Ecosyst. Health Sustain. 2015, 2. [Google Scholar] [CrossRef]
  66. Zarepisheh, M.; Soleimani-Damaneh, M. A dual simplex-based method for determination of the right and left returns to scale in DEA. Eur. J. Oper. Res. 2009, 194, 585–591. [Google Scholar] [CrossRef]
  67. Yagi, M.; Fujii, H.; Hoang, V.; Managi, S. Environmental efficiency of energy, materials, and emissions. J. Environ. Manag. 2015, 161, 206–218. [Google Scholar] [CrossRef] [PubMed]
  68. Liu, B. Uncertainty theory. In Uncertainty Theory; Springer: Berlin/Heidelberg, Germany, 2007; pp. 205–234. [Google Scholar]
  69. Zhao, H.X.; Jiang, X.W.; Cui, J.X. Shifting Path of Industrial Pollution Gravity Centers and Its Driving Mechanism in Pan-Yangtze River Delta. Environ. Sci. 2014, 11, 504–512. [Google Scholar]
  70. Song, J.; Wang, H.X.; Liu, S.Y. Quantitative assessment of stress of economic development to environment using ecological stress index. China J. WIT 2014, 22, 368–374. [Google Scholar] [CrossRef]
  71. Wang, Q.; Liu, Y.L.; Wu, Y.Y.; Wen, Y.L. Analysis of Spatial Characteristics of Emission Intensity of the Main Pollutants in China. Int. Environ. Sci. 2015, 3, 57–61. [Google Scholar]
  72. Chen, Z.H.; Lei, Z.J. The spatial-temporal characteristics and economic drivers of environmental pollution changes in China. Geogr. Res. 2015, 11, 2165–2178. [Google Scholar]
  73. Cao, Z.L.; Yang, J. Manufacturing Level Estimation and Trend Analysis of Environmental Pollution in China. Econ. Geol. 2013, 4, 107–113. [Google Scholar]
  74. Wang, L.; Gong, Z.; Gao, G.; Wang, C. Can energy policies affect the cycle of carbon emissions? Case study on the energy consumption of industrial terminals in Shanghai, Jiangsu and Zhejiang. Ecol. Indic. 2017, 83, 1–12. [Google Scholar] [CrossRef]
  75. Fujii, H.; Managi, S. Optimal production resource reallocation for CO2 emissions reduction in manufacturing sectors. Glob. Environ. Chang. 2015, 35, 505–513. [Google Scholar] [CrossRef]
  76. Lee, B.L.; Wilson, C.; Pasurka, C.A.; Fujii, H.; Managi, S. Sources of airline productivity from carbon emissions: An analysis of operational performance under good and bad outputs. J. Prod. Anal. 2017, 47, 223–246. [Google Scholar] [CrossRef]
  77. Kumar, S.; Fujii, H.; Managi, S. Substitute or complement? Assessing renewable and nonrenewable energy in OECD countries. Appl. Econ. 2015, 47, 1438–1459. [Google Scholar] [CrossRef]
Figure 1. Linear uncertainty distribution.
Figure 1. Linear uncertainty distribution.
Sustainability 09 02080 g001
Figure 2. Normal uncertainty distribution.
Figure 2. Normal uncertainty distribution.
Sustainability 09 02080 g002
Figure 3. DEA efficiency comparisons of three kind of uncertainty distribution.
Figure 3. DEA efficiency comparisons of three kind of uncertainty distribution.
Sustainability 09 02080 g003
Figure 4. DEA efficiency comparison of different DMUs.
Figure 4. DEA efficiency comparison of different DMUs.
Sustainability 09 02080 g004
Table 1. Input and output data of DMUs.
Table 1. Input and output data of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
x 1 [1,2][2,3][4,6][3,4][1,3][3,4]
y 1 [3,5][5,7][30,34][10,13][27,29][25,26]
y 2 [7,8][6,8][32,35][4,5][20,21][85,90]
Table 2. Relative efficiency of DMUs.
Table 2. Relative efficiency of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
θ 0 0.63450.41480.98960.481511
s 1 0 000000
s 1 + 0 000000
s 2 + 0 0004.629600
Table 3. Input and output data of DMUs.
Table 3. Input and output data of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6DMU7
x 1 [1,2][2,3][4,6][3,4][1,3][3,4][0.32,0.95]
y 1 [3,5][5,7][30,34][10,13][27,29][25,26][5,7]
y 2 [7,8][6,8][32,35][4,5][20,21][85,90][8,9]
Table 4. Transformation data of efficient DMUs.
Table 4. Transformation data of efficient DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
x 1 [0.32,0.95][0.3472,0.762][1.5657,3.5449][0.4815,0.963][1,3][3,4]
y 1 [5,7][7,9][34,38][13,16][27,29][25,26]
y 2 [8,9][8,10][35,38][5.2346,6.2346][20,21][85,90]
θ 0 111111
s 1 0 002.4763000
s 1 + 0 03.52650000
s 2 + 0 0004.629600
Table 5. Input and output data of DMUs.
Table 5. Input and output data of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
x 1 [ 1 , 2 ] [ 2 , 3 ] [ 4 , 6 ] [ 3 , 4 ] [ 1 , 3 ] [ 3 , 4 ]
y 1 [ 3 , 5 ] [ 5 , 7 ] [ 30 , 34 ] [ 10 , 13 ] [ 27 , 29 ] [ 25 , 26 ]
y 2 [ 7 , 8 ] [ 6 , 8 ] [ 32 , 35 ] [ 4 , 5 ] [ 20 , 21 ] [ 85 , 90 ]
Table 6. Relative efficiency of DMUs.
Table 6. Relative efficiency of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
θ 0 0.12310.0810.2230.0860.31480.7908
s 1 0 000000
s 1 + 0 0003.241400
s 2 + 0 000000
Table 7. Input and output data of DMUs.
Table 7. Input and output data of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
x 1 N [ 1.5 , 0.016 ] N [ 2.5 , 0.016 ] N [ 5 , 0.063 ] N [ 3.5 , 0.016 ] N [ 2 , 0.063 ] N [ 3.5 , 0.016 ]
y 1 N [ 4 , 0.063 ] N [ 6 , 0.063 ] N [ 32 , 0.25 ] N [ 11.5 , 0.141 ] N [ 28 , 0.063 ] N [ 25.5 , 0.016 ]
y 2 N [ 7.5 , 0.016 ] N [ 7 , 0.063 ] N [ 33.5 , 0.141 ] N [ 4.5 , 0.016 ] N [ 20.5 , 0.016 ] N [ 87.5 , 0.391 ]
Table 8. Relative efficiency of DMUs.
Table 8. Relative efficiency of DMUs.
DMUsDMU1DMU2DMU3DMU4DMU5DMU6
θ 0 0.26120.19320.49930.22940.56180.9967
s 1 0 000000
s 1 + 0 0003.864700
s 2 + 0 000000
Table 9. Explanation of empirical indicators.
Table 9. Explanation of empirical indicators.
General IndicatorsSub-IndicatorsUnit
Input IndicatorsCapitalCIFA(Million Dollar)
Labor ForceAANE(Hundreds of People)
EnergyCEC( Tons of Standard Coal)
Waste Water EmissionsIWWD(Million Tons)
Industrial Soot EmissionsISE(Tons)
Sulfur Dioxide EmissionsISD(Tons)
Output IndicatorsEnvironmental RegulationCUW(percentage)
Gross Output ValueGOV(Billion Yuan)
Table 10. Empirical indicators.
Table 10. Empirical indicators.
DMUsShanghaiNanjingWuxiChangzhou
Input IndicatorsCIFAN(105586,5047)N(153931,46165)N(110321,18084)N(102264,24189)
AANEN(26847,993)N(6923,585)N(12475,733)N(9644,1936)
CECN(10094,6295)N(7070,570)N(2020,128)N(1585,196)
IWWDN(221800,7898)N(33390,4964)N(36454,6652)N(28174,5445)
ISEN(58470,15487)N(48790,10931)N(48516,4426)N(38306,15207)
ISDEN(264620,46622)N(128380,8251)N(104328,13394)N(56916,10867)
Output IndicatorsCUWN(96,1)N(87,4)N(95,2)N(98,1)
GOVN(25190,3236)N(8199,1779)N(11238,1818)N(6603,1589)
Table 11. Empirical indicators.
Table 11. Empirical indicators.
DMUsYangzhouSuzhouZhenjiangNantong
Input IndicatorsCIFAN(69530,19505)N(152959,28544)N(56470,17558)N(118301,28312)
AANEN(6258,795)N(29314,3079)N(4535,434)N(7884,942)
CECN(606,69)N(3077,257)N(578,57)N(694,55)
IWWDN(9652,908)N(66976,3596)N(10167,958)N(16250,1204)
ISEN(12879,1704)N(50021,9908)N(27546,7248)N(40528,3253)
ISDEN(67284,11246)N(197259,20239)N(74044,9279)N(75753,7911)
Output IndicatorsCUWN(94,4)N(98,0.2)N(92.5,3)N(97,2)
GOVN(5038,1543)N(21519,4073)N(4040,1277)N(6888.44,1920)
Table 12. Empirical indicators.
Table 12. Empirical indicators.
DMUsXuzhouHangzhouNingboShaoxing
Input IndicatorsCIFAN(88585,30323)N(58528,8477)N(67685,11423)N(64686,12283)
AANEN(5075,879)N(11980,516)N(15843,801)N(7994,280)
CECN(1125,109)N(1970,88)N(1308,175)N(1195,76)
IWWDN(9746,1475)N(63888,10972)N(17558,1354)N(28311,1839)
ISEN(32752,6556)N(47515,15058)N(26425,2031)N(22383,3951)
ISDEN(109118,11980)N(89183,14146)N(150393,19314)N(63048,3237)
Output IndicatorsCUWN(98,2)N(94,1)N(88,2)N(88,5)
GOVN(5009,2021)N(9478,1333)N(9098,1557)N(6188,1213)
Table 13. Empirical indicators.
Table 13. Empirical indicators.
DMUsJiaxingZhoushanQuzhouTaizhou
Input IndicatorsCIFAN(59441,9775)N(11634,2686)N(20687,4221)N(35372,6469)
AAEN(8637,338)N(1752,104)N(1336,121)N(8061,1873)
CECN(912,128)N(62,24)N(666,40)N(168,7)
IWWDN(18886,1596)N(1876,115)N(14209,793)N(5556,356)
ISEN(24113,2222)N(12252,3256)N(29885,12655)N(10486,1802)
ISDEN(89829,10235)N(20768,3589)N(35248,6020)N(48820,10539)
Output IndicatorsCUWN(96,2)N(98,1)N(96,1)N(95,1)
GOVN(4360,915)N(844,227)N(920,241)N(2894,374)
Table 14. Empirical results.
Table 14. Empirical results.
DMUsShanghaiNanjingWuxiChangzhouYangzhouSuzhouZhenjiangNantong
θ 0 0.94990.57100.52970.80710.99870.47610.52650.9851
DMUsXuzhouHangzhouNingboShaoxingJiaxingZhoushanQuzhouTaizhou
θ 0 0.68380.70870.78750.48270.39870.94290.34350.9570

Share and Cite

MDPI and ACS Style

Gong, Z.; Chen, X. Analysis of Interval Data Envelopment Efficiency Model Considering Different Distribution Characteristics—Based on Environmental Performance Evaluation of the Manufacturing Industry. Sustainability 2017, 9, 2080. https://doi.org/10.3390/su9122080

AMA Style

Gong Z, Chen X. Analysis of Interval Data Envelopment Efficiency Model Considering Different Distribution Characteristics—Based on Environmental Performance Evaluation of the Manufacturing Industry. Sustainability. 2017; 9(12):2080. https://doi.org/10.3390/su9122080

Chicago/Turabian Style

Gong, Zaiwu, and Xiaoqing Chen. 2017. "Analysis of Interval Data Envelopment Efficiency Model Considering Different Distribution Characteristics—Based on Environmental Performance Evaluation of the Manufacturing Industry" Sustainability 9, no. 12: 2080. https://doi.org/10.3390/su9122080

APA Style

Gong, Z., & Chen, X. (2017). Analysis of Interval Data Envelopment Efficiency Model Considering Different Distribution Characteristics—Based on Environmental Performance Evaluation of the Manufacturing Industry. Sustainability, 9(12), 2080. https://doi.org/10.3390/su9122080

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