# An Integrated Slacks-Based Measure of Super-Efficiency with Input Saving and Output Surplus Scaling Factors and its Application in Paper Chemical Mills

Authors: Dong Guo; Zheng-Qun Cai
Journal: Journal of Chemistry
 (2020)

---

## Abstract

Data envelopment analysis (DEA) as a nonparametric programming approach has been widely extended and applied in many areas. Conventional DEA models can well measure the efficiency of inefficient decision-making units (DMUs) but cannot further discriminate the efficient DMUs. A lot of methods are proposed to address this problem. One of the most important methods is the slacks-based measure of super-efficiency model (S-SBM model) developed by Tone in 2002. However, the projection for a DMU on the efficient frontier identified by S-SBM model may not be strongly Pareto-efficient that makes the super-efficiency score misestimated. This paper revises the usual slacks-based measure of super-efficiency by incorporating input saving and output surplus scaling factors into the objection function for measuring DMUs. We integrate SBM model and S-SBM model effectively and yield input saving and output surplus scaling factors as well as input and output slacks under only one integrated model. According to the study, the projection reference point identified by our method is strongly Pareto-efficient. Meanwhile, how each decision variable influences the efficiency score for a specific DMU is revealed and illustrated through two numerical examples and an empirical study in paper chemical mills.


---

## Body

## 1. Introduction

Over the past several decades, the data envelopment analysis (DEA) initially proposed by Charnes et al. has proved to be an effective data-oriented programming method for measuring the efficiencies of a group of homogenous decision making units (DMUs) [1]. The technique of DEA depicts a best-practice production efficient frontier formed by observed DMUs and provides a benchmark or reference point on this frontier for each DMU to compute its efficiency score. Using DEA does not assume any production function or presuppose any specific weight restriction to inputs and outputs. So in the recent years, DEA has been widely extended and applied in many fields [2]. The efficiency value obtained by the classic CCR model indicates how efficiently a DMU has performed when compared with other DMUs so as to determine its efficient level within the group of all DMUs. The CCR model works as a radial model and has both input- and output-orientation styles which permits the DMU under assessment to proportionally reduce all its inputs for producing its given outputs, or to proportionally expand all its outputs by using its given inputs.However, when the evaluation target group includes a considerable number of DMUs, more than one unit will always get the same efficient score of unity. In this case, the CCR model cannot further differentiate the efficiency performances of these DMUs and cannot provide more recognizing information on them. For example, it does not detect whether the evaluated DMU is weakly efficient. As Chen pointed out, the “radial” efficiency model may make some DMU measured against a weakly efficient point on the efficient frontier in the production possibility set [3]. Specifically, the weakly efficient reference point for the current evaluated DMU may still have a positive amount of input excesses or output shortfalls, for it is not the strongly Pareto-efficient reference point. From the view of DEA, the efficiency score obtained by the weakly efficient reference point makes the evaluated DMU misestimated with respect to its strongly Pareto-efficient reference point on the efficient frontier.So far, many methods have been developed and studied in order to enhance the cognition levels and discrimination abilities to distinguish DMUs, such as the cross-efficiency technique [4, 5], the benchmark ranking method, and others [6]. These newly developed methods are mainly built to solve problems resulting from the original CCR model in a certain aspect. Therein, Andersen and Petersen creatively developed the first radial super-efficiency model under the assumption of constant returns to scale (CRS) to reassess the efficient DMUs under CCR model [7]. They exclude the DMUs being evaluated from the reference set by the envelopment linear program so as to retrieve the called super-efficiency scores for those efficient DMUs, while, for the variable returns to scale (VRS) super-efficiency model, Seiford and Zhu found that the problem of infeasibility may occur [8]. Chen further ascribed the sources of super-efficiency for an efficient DMU to its achieved positive amounts of input saving and output surplus regarding its efficient reference point on the super-efficiency frontier [9]. Cook et al. derived a revised VRS super-efficiency model which could generate optimal solutions for some efficient DMUs that is infeasible in the original VRS super-efficiency model [10]. The resulting super-efficiency scores gained by their model could be described from both inputs and outputs aspects to some degree. Later, Lee et al. introduced a two-stage process to address the infeasibility issue under VRS [11]. Next, this two-stage process was merged into a single linear program by the work of Chen and Liang [12]. More recently, Lee and Zhu settled the infeasibility caused by zero input data and decompose the acquired super-efficiency score into three indices: radial efficiency index, input saving index, and output surplus index [13].The above super-efficiency models are all of a radial type and they all have both input- and output-oriented forms. Tone built a popular nonradial model named slacks-based measure (SBM), which uses a fractional objection function depending on input and output slacks instead of a simple radial efficiency variable [14]. The SBM model only has one style, for there is no distinction between input-orientation and output-orientation under SBM. The efficiency score computed by the SBM model is also between 0 and 1. Compared to the traditional radial super-efficiency DEA models under CRS, the SBM model avoids the problem of infeasibility. Moreover, input slacks and output slacks in the SBM model can be utilized to detect the input excesses and output shortfalls of a given DMU, respectively. However, for SBM-efficient DMUs, Tone designed a super-efficiency model (S-SBM model) to examine their super-efficiency scores in order to allow SBM-efficient DMUs to be also ranked and compared [15]. Liu and Chen developed a HypoSBM model to distinguish the worst-performance DMUs from the bad ones [16]. Du et al. extended the SBM super-efficiency model to the additive slacks-based DEA model [17]. Fang et al. established a two-stage process, which was an alternative disposal treatment for the SBM method proposed by Tone [18]. They demonstrated that their two-stage approach generated the same results as Tone’s models. Chen indicated that there exists a discontinuous gap between the SBM score and S-SBM score for a DMU who has a weakly reference point under S-SBM model [3]. In his study, an ambidextrous joint computation model (J-SBM model) for slacks-based measure was provided.However, J-SBM model has two noted shortcomings. First, it fails to represent all input saving and output surplus scaling factors explicitly for each DMU. Second, it has a more complicated operational procedure and needs a three-stage computational process. The current paper further investigates this topic and presents a modified slacks-based measure of super-efficiency with input saving and output surplus scaling factors to reevaluate these DMUs. Our approach is devoted to detecting the specific scaling factors for each input and output as well as input and output slacks explicitly by one integrated model. Meanwhile, the phenomenon of discontinuousness between SBM score and S-SBM score for the same DMU can be eliminated. The results indicate that the projection obtained through our proposed model is strongly Pareto-efficient. And in this approach, we can see clearly how each decision variable influences the final efficiency score for a specific DMU.The structure of this paper is organized as follows. Section2 briefly reviews several kinds of slacks-based models. Section 3 presents a modified slacks-based measure of super-efficiency with input saving and output surplus scaling factors to determine DMUs’ efficiency scores. In Section 4, two numerical examples are applied to compare our approach with the previous models. Section 5 applies our approach to an empirical example where the performance of 32 paper chemical mills in China is evaluated. The main conclusions and remarks are given in Section 6.


## 2. Preliminaries

Suppose there aren DMUs, {DMUkk=1,2,…,n}. Let xk=x1k,…,xmk and yk=y1k,…,ysk denote the input and output vectors of the kth DMU. The ith input of the kth DMU is denoted as xik and the rth output of the kth DMU is denoted as yrk, respectively. λj is the intensity coefficient for the kth DMU which means its contribution to forming the efficient frontier. Assume that all input and output data are positive.


### 2.1. Tone’s SBM Measure

Tone developed the following SBM model to evaluate the relative efficiency ofDMUk [14], where input and output slacks are denoted, respectively, as zi−i=1,…,m and zr+r=1,…,s. Unlike the CCR model, the objective function in SBM model has a fractional form which directly includes input and output slacks.(1)ρk∗=min1−1/m∑i=1mzi−/xik1+1/s∑r=1szr+/yrk,s.t.∑j=1nλjxij=xik−zi−,i=1,…,m,∑j=1nλjyrj=yrk+zr+,r=1,…,s,λj≥0,j=1,…,n,zi−≥0,i=1,…,m,zr+≥0,r=1,…,s,where DMUk is called SBM-efficient if and only if zi−=zr+=0 for all i and r, that is, ρk∗=1; otherwise, it is called SBM-inefficient. In order to further discriminate SBM-efficient DMUs with the same SBM efficiency score of 1, Tone introduced a SBM super-efficiency model which was referred to as S-SBM model in the following formula [15].(2)δk∗=min1/m∑i=1mx¯ik/xik1/s∑r=1sy¯rk/yrk,s.t.∑j=1j≠knλjxij≤x¯ik,i=1,…,m,∑j=1j≠knλjyrj≥y¯rk,r=1,…,s,x¯ik≥xik,i=1,…,m,y¯rk≤yrk,r=1,…,s,λj,y¯rk≥0,j≠k,j=1,…,n,i=1,…,s.It should be noticed that the S-SBM model can only be used for the DMU whose SBM efficiency scoreρ∗=1. Through model (2), these SBM-efficient DMUs get its super-efficiency scores. However, model (2) cannot discriminate SBM-inefficient DMUs for they will get the same efficiency score of 1.


### 2.2. Fang’s Models

Fang et al. provided a two-stage process which brings in the same efficiency scores as those obtained by Tone’s two models [18]. In the first stage, they replace x¯ik, y¯rk with xik+wi−, yrk+wr+ to form model (3) for detecting both input savings and output surpluses for all DMUs first.(3)δk∗=min1+1/m∑i=1mwi−/xik1−1/s∑r=1swr+/yrk,s.t.∑j=1j≠knλjxij≤xik+wi−,i=1,…,m,∑j=1j≠knλjyrj≥yrk−wr+,r=1,…,s,λj≥0,j≠k,j=1,…,n,wi−≥0,i=1,…,m,wr+≥0,wr+≤yrk,r=1,…,s.Through Model (3), the optimal wi−∗ and wr+∗ for each input and output for all DMUs can be obtained. For SBM-efficient DMUs in model (1), there will exist at least one i or r, so that wi−∗>0 or wr+∗>0 in model (3), while, for SBM-inefficient DMUs in model (1), they have no positive input saving and output surplus for each corresponding input and output at all. So these DMUs will get wi−=wr+=0 for all i and r. Then, plug wi−∗ and wr+∗ into the following model, and a slacks-based measure is reconstructed as model (4) exhibits.(4)ρk∗=min1−1/m∑i=1msi−/xik1+1/s∑r=1ssr+/yrk,s.t.∑j=1j≠knλjxij=xik+wi−∗−si−,i=1,…,m,∑j=1j≠knλjyrj=yrk−wr+∗+sr+,r=1,…,s,λj≥0,j≠k,j=1,…,n,si−≥0,i=1,…,m,sr+≥0,r=1,…,s.Obviously, model (3) is equivalent to model (2). So model (3) also determines the super-efficiency scores for SBM-efficient DMUs. The optimal SBM efficiency scores and the optimal slacks si−∗,sr+∗ for SBM-inefficient DMUs can be computed based on the optimal wi−∗,wi+∗ by model (4).At last, the final efficiency score forDMUk based on SBM is defined as(5)φk∗=1+1/m∑i=1mwi−∗/xik1−1/s∑r=1swr+∗/yrk,if1+1/m∑i=1mwi−∗/xik1−1/s∑r=1swr+∗/yrk>1,1−1/m∑i=1msi−∗/xik1+1/s∑r=1ssr+∗/yrk,otherwise.


### 2.3. Chen’s Model

Chen investigated a data set in Tone's which can be seen in Table1 and found the reference points generated by model (2) or (3) might not be Pareto-efficient [3]. This resulted in a discontinuous gap between the SBM score and S-SBM score. The author indicated that the issue of discontinuity makes troubles to rationalize the efficiency score. So he established an ambidextrous joint computation model (6) (J-SBM model) to find the Pareto-efficient point for each DMU.(6)minϕk=JSBMkxJSBMkx−Mb1+1−b1b2,s.t.JSBMkx=1−1mb1∑i=1msi−xik−1−b1b2∑i=1msi−xik+1−b11−b2∑i=1ms˜i−xik,JSBMky=1+1sb1∑r=1ssr+yrk−1−b1b2∑r=1ssr+yrk+1−b11−b2∑r=1ss˜r+yrk,6.1b1∑j=1,j≠knλjxij=b1xik−si−,i=1,…,m,b1∑j=1,j≠knλjyrj=b1yrk+sr+,r=1,…,s,6.21−b1b2∑j=1,j≠knλjxij=1−b1b2xik+si−,i=1,…,m,1−b1b2∑j=1,j≠knλjyrj=1−b1b2yrk−sr+,r=1,…,s,6.31−b11−b2∑j=1,j≠knλjxij=1−b11−b2xik−s˜i−,i=1,…,m,1−b11−b2∑j=1,j≠knλjyrj=1−b11−b2yrk+s˜r+,r=1,…,s,λj≥0,j=1,…,n,si−≥0,s˜i−free fori=1,…,m,sr+≥0,s˜r+free forr=1,…,s,b1,b2∈0,1,Mis a large enough positive number.Table 1
Data set 1 from Tone [15].
DMU
x
1
x
2
y
1
Tone’s or Fang’s approach
Chen’s approach
Our approach
S-SBM
SBM
J-SBM
Integrated model (10)
δ
k
∗
ρ
k
∗
φ
k
∗
ϕ
k
∗
ψ
k
∗
A
4
3
1
1
0.833
0.833
0.833
0.833
B
7
3
1
1
0.619
0.619
0.619
0.619
C
8
1
1
1.125
1
1.125
1.125
1.125
D
4
2
1
1.25
1
1.25
1.25
1.25
E
2
4
1
1.5
1
1.5
1.25
1.25
F
10
1
1
1
0.9
0.9
0.9
0.9
G
12
1
1
1
0.833
0.833
0.833
0.833In model (6), b1 and b2 are two binary variables that are used to control which one of the three kinds of constraint conditions (6.1), (6.2), and (6.3) is chosen. For the SBM-inefficient DMUs, b1=b2=1, constraint condition 6.1 is active; now model (6) works as the SBM model. For SBM-efficient DMUs, model (6) first acts as the S-SBM model under active constraint condition 6.2 when b1=0,b2=1. In the meantime, if the super-efficiency reference point is not Pareto-efficient for the evaluated DMU, model (6) will activate constraint condition 6.3, and the corresponding super-efficiency score will be corrected for the DMU. Additionally, Chen further proved that the reference points for all DMUs under model (6) are Pareto-efficient [3].As can be seen from the above procedures in operating model (6), the constraint condition 6.3 is set intentionally to correct the misestimated efficiency score due to the less Pareto-efficient reference point caused by the S-SBM model (2) or (3). For these DMUs who have a weakly efficient reference point under model (2) or (3), however, the achievement of its super-efficiency needs a three-stage process. For instance, when DMUE in Table 1 gets its super-efficiency of 1.25 under the constraint condition 6.3, the constraint conditions (6.1), (6.2) have been inspected before.


## 2.1. Tone’s SBM Measure

Tone developed the following SBM model to evaluate the relative efficiency ofDMUk [14], where input and output slacks are denoted, respectively, as zi−i=1,…,m and zr+r=1,…,s. Unlike the CCR model, the objective function in SBM model has a fractional form which directly includes input and output slacks.(1)ρk∗=min1−1/m∑i=1mzi−/xik1+1/s∑r=1szr+/yrk,s.t.∑j=1nλjxij=xik−zi−,i=1,…,m,∑j=1nλjyrj=yrk+zr+,r=1,…,s,λj≥0,j=1,…,n,zi−≥0,i=1,…,m,zr+≥0,r=1,…,s,where DMUk is called SBM-efficient if and only if zi−=zr+=0 for all i and r, that is, ρk∗=1; otherwise, it is called SBM-inefficient. In order to further discriminate SBM-efficient DMUs with the same SBM efficiency score of 1, Tone introduced a SBM super-efficiency model which was referred to as S-SBM model in the following formula [15].(2)δk∗=min1/m∑i=1mx¯ik/xik1/s∑r=1sy¯rk/yrk,s.t.∑j=1j≠knλjxij≤x¯ik,i=1,…,m,∑j=1j≠knλjyrj≥y¯rk,r=1,…,s,x¯ik≥xik,i=1,…,m,y¯rk≤yrk,r=1,…,s,λj,y¯rk≥0,j≠k,j=1,…,n,i=1,…,s.It should be noticed that the S-SBM model can only be used for the DMU whose SBM efficiency scoreρ∗=1. Through model (2), these SBM-efficient DMUs get its super-efficiency scores. However, model (2) cannot discriminate SBM-inefficient DMUs for they will get the same efficiency score of 1.


## 2.2. Fang’s Models

Fang et al. provided a two-stage process which brings in the same efficiency scores as those obtained by Tone’s two models [18]. In the first stage, they replace x¯ik, y¯rk with xik+wi−, yrk+wr+ to form model (3) for detecting both input savings and output surpluses for all DMUs first.(3)δk∗=min1+1/m∑i=1mwi−/xik1−1/s∑r=1swr+/yrk,s.t.∑j=1j≠knλjxij≤xik+wi−,i=1,…,m,∑j=1j≠knλjyrj≥yrk−wr+,r=1,…,s,λj≥0,j≠k,j=1,…,n,wi−≥0,i=1,…,m,wr+≥0,wr+≤yrk,r=1,…,s.Through Model (3), the optimal wi−∗ and wr+∗ for each input and output for all DMUs can be obtained. For SBM-efficient DMUs in model (1), there will exist at least one i or r, so that wi−∗>0 or wr+∗>0 in model (3), while, for SBM-inefficient DMUs in model (1), they have no positive input saving and output surplus for each corresponding input and output at all. So these DMUs will get wi−=wr+=0 for all i and r. Then, plug wi−∗ and wr+∗ into the following model, and a slacks-based measure is reconstructed as model (4) exhibits.(4)ρk∗=min1−1/m∑i=1msi−/xik1+1/s∑r=1ssr+/yrk,s.t.∑j=1j≠knλjxij=xik+wi−∗−si−,i=1,…,m,∑j=1j≠knλjyrj=yrk−wr+∗+sr+,r=1,…,s,λj≥0,j≠k,j=1,…,n,si−≥0,i=1,…,m,sr+≥0,r=1,…,s.Obviously, model (3) is equivalent to model (2). So model (3) also determines the super-efficiency scores for SBM-efficient DMUs. The optimal SBM efficiency scores and the optimal slacks si−∗,sr+∗ for SBM-inefficient DMUs can be computed based on the optimal wi−∗,wi+∗ by model (4).At last, the final efficiency score forDMUk based on SBM is defined as(5)φk∗=1+1/m∑i=1mwi−∗/xik1−1/s∑r=1swr+∗/yrk,if1+1/m∑i=1mwi−∗/xik1−1/s∑r=1swr+∗/yrk>1,1−1/m∑i=1msi−∗/xik1+1/s∑r=1ssr+∗/yrk,otherwise.


## 2.3. Chen’s Model

Chen investigated a data set in Tone's which can be seen in Table1 and found the reference points generated by model (2) or (3) might not be Pareto-efficient [3]. This resulted in a discontinuous gap between the SBM score and S-SBM score. The author indicated that the issue of discontinuity makes troubles to rationalize the efficiency score. So he established an ambidextrous joint computation model (6) (J-SBM model) to find the Pareto-efficient point for each DMU.(6)minϕk=JSBMkxJSBMkx−Mb1+1−b1b2,s.t.JSBMkx=1−1mb1∑i=1msi−xik−1−b1b2∑i=1msi−xik+1−b11−b2∑i=1ms˜i−xik,JSBMky=1+1sb1∑r=1ssr+yrk−1−b1b2∑r=1ssr+yrk+1−b11−b2∑r=1ss˜r+yrk,6.1b1∑j=1,j≠knλjxij=b1xik−si−,i=1,…,m,b1∑j=1,j≠knλjyrj=b1yrk+sr+,r=1,…,s,6.21−b1b2∑j=1,j≠knλjxij=1−b1b2xik+si−,i=1,…,m,1−b1b2∑j=1,j≠knλjyrj=1−b1b2yrk−sr+,r=1,…,s,6.31−b11−b2∑j=1,j≠knλjxij=1−b11−b2xik−s˜i−,i=1,…,m,1−b11−b2∑j=1,j≠knλjyrj=1−b11−b2yrk+s˜r+,r=1,…,s,λj≥0,j=1,…,n,si−≥0,s˜i−free fori=1,…,m,sr+≥0,s˜r+free forr=1,…,s,b1,b2∈0,1,Mis a large enough positive number.Table 1
Data set 1 from Tone [15].
DMU
x
1
x
2
y
1
Tone’s or Fang’s approach
Chen’s approach
Our approach
S-SBM
SBM
J-SBM
Integrated model (10)
δ
k
∗
ρ
k
∗
φ
k
∗
ϕ
k
∗
ψ
k
∗
A
4
3
1
1
0.833
0.833
0.833
0.833
B
7
3
1
1
0.619
0.619
0.619
0.619
C
8
1
1
1.125
1
1.125
1.125
1.125
D
4
2
1
1.25
1
1.25
1.25
1.25
E
2
4
1
1.5
1
1.5
1.25
1.25
F
10
1
1
1
0.9
0.9
0.9
0.9
G
12
1
1
1
0.833
0.833
0.833
0.833In model (6), b1 and b2 are two binary variables that are used to control which one of the three kinds of constraint conditions (6.1), (6.2), and (6.3) is chosen. For the SBM-inefficient DMUs, b1=b2=1, constraint condition 6.1 is active; now model (6) works as the SBM model. For SBM-efficient DMUs, model (6) first acts as the S-SBM model under active constraint condition 6.2 when b1=0,b2=1. In the meantime, if the super-efficiency reference point is not Pareto-efficient for the evaluated DMU, model (6) will activate constraint condition 6.3, and the corresponding super-efficiency score will be corrected for the DMU. Additionally, Chen further proved that the reference points for all DMUs under model (6) are Pareto-efficient [3].As can be seen from the above procedures in operating model (6), the constraint condition 6.3 is set intentionally to correct the misestimated efficiency score due to the less Pareto-efficient reference point caused by the S-SBM model (2) or (3). For these DMUs who have a weakly efficient reference point under model (2) or (3), however, the achievement of its super-efficiency needs a three-stage process. For instance, when DMUE in Table 1 gets its super-efficiency of 1.25 under the constraint condition 6.3, the constraint conditions (6.1), (6.2) have been inspected before.


## 3. A Modified Slacks-Based Measure of Super-efficiency

In this section, we first establish an equivalent form of S-SBM model, which explicitly contains input saving and output surplus scaling factors. Note thatx¯ik≥xik for all i and y¯rk≤yrk for all r in model (2). Let x¯ik=xik+tixik,ti≥0 and y¯rk=yrk−βryrk≥0,0≤βr≤1, or wi−=xik+tixik,ti≥0 and wr+=yrk−βryrk≥0,0≤βr≤1; then model (2) or model (3) will become the following model (7).(7)δk∗=min1+1/m∑i=1mti1−1/s∑r=1sβr,s.t.∑j=1j≠knλjxij≤1+tixik,i=1,…,m,∑j=1j≠knλjyrj≥1−βryrk,r=1,…,s,λj≥0,j≠k,j=1,…,n,ti≥0,i=1,…,m,0≤βr≤1,r=1,…,s.The difference between models (7) and (3) is that model (7) can not only measure the input saving tixik and the output surplus βryrk, but also present the specific scaling factors ti for xik and βr for yrk of DMUk.Letti∗,βr∗ be the optimal solution of model (7). Based on model (7), the standard SBM model (1) can be revised as follows:(8)ρk∗=min1−1/m∑i=1msi−/xik1+1/s∑r=1ssr+/yrk,s.t.∑j=1j≠knλjxij=1+ti∗xik−si−,i=1,…,m,∑j=1j≠knλjyrj=1−βr∗yrk+sr+,r=1,…,s,λj≥0,j≠k,j=1,…,n,si−≥0,i=1,…,m,sr+≥0,r=1,…,s.Similarly, we apply model (7) to all DMUs before the utilization of model (8) on inefficient DMUs. We define the final efficiency score for DMUk based on SBM with the following piecewise function:(9)φk∗=1+1/m∑i=1mti∗1−1/s∑r=1sβr∗,if1+1/m∑i=1mti∗1−1/s∑r=1sβr∗>1,1−1/m∑i=1msi−∗/xik1+1/s∑r=1ssr+∗/yrk,otherwise,where si−∗,si+∗ is the optimal solution in model (8).The above two-stage process is identical to that of Fang et al. Model (7) may also suffer from the problem that the projection reference point for a specific DMU may not be Pareto-efficient. As Chen mentioned, the SBM and S-SBM efficiency scores are calculated in two separate models with two different objection functions that have two different projecting styles for the evaluated DMU [3]. The SBM scores are achieved by their input and output slacks, while the S-SBM scores depend only on the reference point (may not be Pareto-efficient) on the super-efficiency frontier. Therefore, the whole evaluation system is not an integrated one and is divided into two styles, which results in the problem of a discontinuous gap.In order to integrate the whole evaluation system, we maintain the basic sense of SBM model and S-SBM model into one model and ensure that the desired model has the uniform projecting way for the evaluated DMU. So, we develop the following modified SBM measure based on model (7).(10)SE=min1−1/m∑i=1msi−/xik−ti1+1/s∑r=1ssr+/yrk−βr+M×1+1/m∑i=1mti1−1/s∑r=1sβr,s.t.∑j=1j≠knλjxij=1+tixik−si−,i=1,…,m,∑j=1j≠knλjyrj=1−βryrk+sr+,r=1,…,s,λj≥0,j≠k,j=1,…,n,si−≥0,i=1,…,m,sr+≥0,r=1,…,s,ti≥0,i=1,…,m,0≤βr<1,r=1,…,s.The logic behind model (10) is based on two major attributes. First, the measure we modified as1−1/m∑i=1msi−/xik−ti/1+1/s∑r=1ssr+/yrk−βr, which not only includes the input and output slacks in SBM measure, but also includes input saving and output surplus scaling factors in S-SBM model. It is able to comprehensively explain the reason for the resulting efficiency scores; that is to say, when input saving scaling factors or output surplus scaling factors appear together with input or output slacks, they all contribute to the efficiency score.Hence, the efficiency score can be defined as(11)ψk∗=1−1/m∑i=1msi−∗/xik−ti∗1+1/s∑r=1ssr+∗/yrk−βr∗=1−1/m∑i=1msi−∗/xik+1/m∑i=1mti∗1+1/s∑r=1ssr+∗/yrk−1/s∑r=1sβr∗,where ti∗,βr∗,si−∗,si+∗ is the optimal solution of model (10).Second, we integrate SBM model and S-SBM model under the rule that the evaluated DMU must have the minimum super-efficiency score if it has at least oneti∗>0 or βr∗>0. So we add the controllable term of M×1+1/m∑i=1mti/1−1/s∑r=1sβr to the objective function so as to make ti and βr as small as possible for all i and r.Here,M (we set it 107 so as to make it large enough in the empirical study) is a predetermined large number used mainly for magnifying the effect of input saving and output surplus scaling factors on the objective function under the condition of minimum S-SBM score. Simultaneously, it can control the optimal values of decision variables. Note the constraints in model (10):(12)∑j=1j≠knλjxij=1+tixik−si−,i=1,…,m,∑j=1j≠knλjyrj=1−βryrk+sr+,r=1,…,s.A severe problem which may occur is that whenti>0 and si−>0 for a certain i (or βr>0 and sr+>0 for a certain r) simultaneously exist, there is an offsetting contradiction between tixik and si− (or βryrk and sr+). To avoid this kind of situation, the controllable term of M×1+1/m∑i=1mti/1−1/s∑r=1sβr added to the objective function is able to ensure that either ti∗>0 or si−∗>0 for a certain i (either βr∗>0 or sr+∗>0 for a certain r) and both of them are not in existence.Model (10) has close connections with the SBM model and S-SBM model.Theorem 1.
For any SBM-inefficientDMUk, ψk∗=ρk∗.Proof.
For SBM-inefficientDMUk, there must be ti∗=0 and βr∗=0 for all i and r in model (10). In this case, model (10) degenerates into SBM model (1).Theorem 2.
For SBM-efficientDMUk, if the reference point under S-SBM model for DMUk is Pareto-efficient, ψk∗=δk∗.Proof.
For SBM-inefficientDMUk, suppose the reference point of DMUk in S-SBM model is Pareto-efficient, there must be si−∗ and sr+∗=0 for all i and r in model (10). At this moment, model (10) degenerates into model (7) which is equivalent to S-SBM model (2).
However, for other SBM-efficient DMUs, if the reference point under S-SBM model is not Pareto-efficient, model (10) will revise the super-efficiency score obtained by S-SBM model (7) but just make sure that model (10) can identify the Pareto-efficient reference for these DMUs.Theorem 3.
The reference point identified by model (10) is Pareto-efficient.Proof.
For a SBM-inefficient DMU, model (10) serves as the standard SBM model (1); the reference point identified the standard SBM model as Pareto-efficient, and so does model (10).
For a SBM-efficient DMU, there at least exists one inputi∈1,…,m or one output r∈1,…,s such that at least one scaling factor ti∗>0 or βr∗>0. Meanwhile, the existence of input and output slacks allows the evaluated DMU to decrease its inputs and increase its outputs. Thus, the reference point identified by model (10) is as follows:(13)∑j=1j≠knλj∗xij=1+ti∗xik−si−∗,i=1,…,m,∑j=1j≠knλj∗yrj=1−βr∗yrk+sr+∗,r=1,…,s.
If1+ti∗xik,1−βr∗yrk is not a Pareto-efficient reference point, then it moves si−∗ leftward and sr+∗ upward to reach the Pareto-efficient reference points.
However, it should be noted that model (10) is a quite different optimization design because the objective function and constraints are quite different to the standard SBM and S-SBM models. Although the projection identified by model (10) is Pareto-efficient, for SBM-efficient DMUs, the projected reference point by model (10) may not coincide with that by S-SBM model or J-SBM model (see examples in Section 4).
Clearly, the modified measurementψk has the unit-invariant property.Theorem 4.
ψ
k is unit-invariant.The fractional programming (10) can be transformed into a linear programming problem through the following Charnes-Cooper transformation by setting w1=1/1+1/s∑r=1ssr+/yrk−βr, w2=1/1−1/s∑r=1sβr (w1≤w2), λj′=w1⋅λj, si′−=w1⋅si−, sr′+=w1⋅sr+, ti′=w1⋅ti, βr′=w1βr, ti″=w2⋅ti, βr″=w2βr; then model (10) will become the following linear programming:(14)minw1−1m∑i=1msi′−xik−ti′+M×w2+1m∑i=1mti″,s.t.w1+1s∑r=1ssr′+yrk−βr′=1,w2−1s∑r=1sβr″=1,∑j=1j≠knλj′xij=w1xik+ti′xik−si′−,i=1,…,m,∑j=1j≠knλj′yrj=w1yrk−βr′yrk+sr′+,r=1,…,s,λj′≥0,j≠k,j=1,…,n,si′−≥0,i=1,…,m,sr′+≥0,r=1,…,s,0≤ti′≤ti″,i=1,…,m,0≤βr′≤w1,βr′≤βr″,r=1,…,s,0≤w1≤w.Supposew1∗,ti′∗,βr′∗,si′−∗,sr′+ is the optimal solution of model (11), then ti∗=ti′∗/w1∗,βr∗=βr′∗/w1∗,si−∗=si′−∗/w1∗,sr+∗=sr′+∗/w1∗ is the optimal solution of model (10).


## 4. Illustration Examples

Two numerical examples from Tone are used to verify our approach through comparing with the SBM model, S-SBM model, and J-SBM model. In Table1, the efficiency scores derived by the SBM model, S-SBM model, and J-SBM model are presented and those identified through our approach are listed in the last column. Table 2 shows detailed optimal solutions under J-SBM model (6) and we rewrite them in the form of input and output surplus scaling factors so as to facilitate comparison. Table 3 presents detailed optimal solutions under our model (10). As shown in Tables 2 and 3, these DMUs get the same efficiency scores but may be a different optimal solution for each decision variable. For example, DMUD gets the super-efficiency score of 1.25 under J-SBM model (6), and the achievement of super-efficiency score seems reflected in its input saving scaling factors t1∗=0.234 and t2∗=0.266, while under our model (10) it not only owns input saving scaling factors t1∗=0.137 and t2∗=0.163, but also has the output surplus scaling factor β1∗=0.08. The reason is that due to the different projecting style of these two models, DMUD gets different reference points, which are both Pareto-efficient projected points. And for DMUE, J-SBM model (6) and our model (10) identified the same Pareto-efficient reference point. Our model also overcomes the problem of the discontinuous gap and finds its input saving scaling factor t1∗=1 and input slack s2−∗=2, so the rational super-efficiency score is(15)ψk∗=1−1/m∑i=1msi−∗/xik−ti∗1+1/s∑r=1ssr+∗/yrk−βr∗=1−1/2×1/2+1/2×11=1.25.Table 2
Results from J-SBM model (6).
DMU
x
1
x
2
y
1
φ
k
∗
t
1
∗
t
2
∗
β
1
∗
s
1
−
∗
s
2
−
∗
s
1
+
∗
A
4
3
1
0.833
0
0
0
0
1
0
B
7
3
1
0.619
0
0
0
2.895
0.948
0.026
C
8
1
1
1.125
0.25
0
0
0
0
0
D
4
2
1
1.25
0.234
0.266
0
0
0
0
E
2
4
1
1.25
1
0
0
0
2
0
F
10
1
1
0.9
0
0
0
2
0
0
G
12
1
1
0.833
0
0
0
4
0
0Table 3
Results from our model (10).
DMU
x
1
x
2
y
1
ψ
k
∗
t
1
∗
t
2
∗
β
1
∗
s
1
−
∗
s
2
−
∗
s
1
+
∗
A
4
3
1
0.833
0
0
0
0
1
0
B
7
3
1
0.619
0
0
0
1.706
0.353
0.324
C
8
1
1
1.125
0.25
0
0
0
0
0
D
4
2
1
1.25
0.137
0.163
0.08
0
0
0
E
2
4
1
1.25
1
0
0
0
2
0
F
10
1
1
0.9
0
0
0
2
0
0
G
12
1
1
0.833
0
0
0
4
0
0In Table4, we mainly compare the approach from Tone’s or Fang’s and ours using data set 2 from Tone. For the first four DMUs (1–4), two approaches gain the same efficiency scores. But for DMU5, our approach detects that it has both input saving scaling factor t1∗=0.4 and output surplus scaling factor β2∗=0.475. Besides, it has input slack s2−∗=1.9 and output slack s1+∗=0.4. This situation means that S-SBM model does not project DMU5 at a Pareto-efficient targeted point and its efficiency score should be revised. According to our approach, the efficiency score should be modified as(16)ψk∗=1−1/m∑i=1msi−∗/xik−ti∗1+1/s∑r=1ssr+∗/yrk−βr∗=1−1/2×1.711/4+1/2×0.5261+1/2×0.526−1/2×0.428=1.0001.Table 4
Data set 2 from Tone [14] and results from our approach.
DMU
x
1
x
2
y
1
y
2
Tone’s or Fang’s approach
Our approach
φ
k
∗
ψ
k
∗
t
1
∗
t
2
∗
β
1
∗
β
2
∗
s
1
−
∗
s
2
−
∗
s
1
+
∗
s
2
+
∗
1
4
3
2
3
0.7980
0.7980
0
0
0
0
0
0.357
0.714
0
2
6
3
2
3
0.5682
0.5682
0
0
0
0
0
0.643
2.286
0
3
8
1
6
2
1.3333
1.3333
0
0
0
0.5
0
0
0
0
4
8
1
6
1
0.6667
0.6667
0
0
0
0
0
0
0
1
5
2
4
1
4
1.4545
1.0001
0.526
0
0
0.428
0
1.711
0.526
0


## 5. An Empirical Study

After high speed development on industries and economics during the past several decades, China has accumulated a lot of air environment problems and the Chinese government is paying more and more attention to ecological and environmental assessment. Recent years have seen a lot of studies based on DEA to measure China’s environmental pollution problems. For example, Zhou et al. [19, 20] construct a set of DEA models with the integral and zero-sum gain constraints for calculating air quality and a new nonradial directional distance function to scale the performance of water use and wastewater emission, respectively.In this section, the method that we have developed is applied to examine the efficiencies and rankings of 32 paper chemical mills along the Huai River in China. In this empirical study, four input and output variables are considered to evaluate each mill’s performance. The inputs of each paper mill include labor and capital, and good output as paper products as well as bad output as biochemical oxygen demand (BOD). The detailed data set of these 32 paper chemical mills is shown in Table5.Table 5
The Data set and efficiency results of 32 paper mills.
Paper mills
Inputs
Bad output
Good output
ψ
k
∗
Rank
Labor (person)
Capital (¥10000)
BOD (ton)
Paper (ton)
1
1077
2959.9
21.4290
27582
0.9005
5
2
452
3589
19.8062
29514
0.8544
7
3
319
5901.9
12.3287
14700
0.3590
26
4
1075
4892.8
9.1559
22354
0.4337
16
5
813
4079.7
11.9146
20669
0.4328
17
6
850
5239.6
5.2037
8222
0.2019
32
7
1090
3022.8
3.6054
15066
0.5286
12
8
122
3173.1
3.7278
8066
0.4937
14
9
297
2277.4
8.0765
19125
1.0734
3
10
1047
1491.9
8.9060
7601
0.3141
30
11
1010
3940.1
4.8940
11579
0.3184
28
12
262
3236.5
4.0835
23216
1.4546
1
13
551
4448.6
4.8750
21698
0.6357
10
14
671
1789.7
4.5334
8127
0.3617
25
15
577
2310.9
6.1362
11549
0.4179
20
16
208
3398.2
7.0186
10295
0.4130
21
17
667
5331.9
28.4877
29881
0.4904
15
18
878
3450.4
13.1680
19076
0.4236
19
19
640
3109.8
6.1616
12176
0.3694
24
20
927
3345.2
1.4533
5187
0.3024
31
21
167
4329.7
22.5809
24005
0.8937
6
22
903
3855.2
26.3390
23085
0.4258
18
23
720
1908.3
3.0787
6545
0.3182
29
24
629
3468.2
21.7332
27599
0.6853
9
25
152
5571.7
5.3061
23748
1.2256
2
26
1010
4647.1
9.1360
17323
0.3489
27
27
578
2513.3
5.0049
10617
0.3898
23
28
384
2247.4
4.4373
9083
0.3968
22
29
166
3968.1
10.6127
15151
0.5997
11
30
894
1368.5
4.7758
9911
0.5089
13
31
143
5350.2
17.4677
25260
0.8538
8
32
879
2873.2
26.9134
27721
0.9464
4
Source: Anhui Environmental Protection Bureau, the Fuyang Environmental Protection Bureau, and the Huainan Environmental Protection Bureau.In papermaking industry, the good products are always produced with bad products and it is impossible to increase good outputs and decrease bad outputs meanwhile. First, the attribute of “the smaller the better” for bad outputs is consistent with inputs. Second, bad outputs cannot create any new profit and dealing with them (such as sewage treatment and air purification) always comes at a price which the mill should afford [21]. Therefore, here we simply treat the bad outputs as inputs from the cost point of view.We use model (10) to assess the performance of 32 paper chemical mills and the computed efficiency scores and ranking results are also displayed in Table 5. The last column gives their ranks according to efficiency scores obtained by our method. Of the 32 paper mills, only three mills, named mills 9, 12, and 25, are identified to obtain super-efficiency scores and they are ranked as the top three accordingly. For example, mill 9 and 12 get super-efficiency of 1.0734 and 1.4546, respectively, due to that the scaling factor of output surplus in one good output (paper) is 0.0734 and 0.4546 separately and there exist no slacks in all inputs and outputs. In contrast, 29 mills fail to get super-efficiency scores. For instance, mill 6 ranks dead last because it does not have any scaling factor of each input or output. Meanwhile, positive slacks are found to exist in two inputs and one bad output. Although mill 31 achieves the input saving scaling factor of 0.0229 in the first input (labor) and the output surplus scaling factor of 0.1077 in the good output (paper), its efficiency score is identified as 0.8538 and is ranked eighth finally. The reason for this is that there exists a slack amount of 13.7051 in the bad output (BOD), which brings more negative impact to efficiency scores than the positive effect brought by input and output scaling factors.As can be seen from Table5, the efficiency score or the super-efficiency score for each mill can be obtained by the integrated model (10). This course by our method prevents transformations among three submodels of the method by Chen. The final efficiency scores ψk∗k=1,…,32 are affected synthetically by all decision variables that include input saving and output surplus scaling factors, input and output slacks. So we can fully differentiate the performances of all 32 paper mills and can provide a complete ranking criterion for all DMUs to be compared with regard to these efficiency scores.


## 6. Conclusions and Remarks

The classic SBM model proposed by Tone projects the DMU under evaluation at a Pareto-efficient reference point on the production frontier, while the S-SBM model proposed by Tone may not apply and the reference point on the super-efficiency frontier by S-SBM model may not be Pareto-efficient. Chen found an approach through transformations among three submodels to handle this issue.In the present paper, we build a universal model to realize the integration of SBM model and S-SBM model. Our approach not only shows the specific scaling factors for each input and output of a specific DMU explicitly, but also provides input saving, output surplus, and slacks information simultaneously only in one model. And we can see clearly how the efficiency scores are obtained by the optimal value for each decision making variable. Thus, more recognizing information on DMUs is revealed via two numerical examples and an empirical study in paper chemical mills. The composite efficiency scores for 32 mills are represented and ranked by integrating the effect of all decision variables, which include input saving, output surplus scaling factors, and input and output slacks.The current paper also overcomes the problem of a discontinuous gap between SBM score and S-SBM score found by Chen. Model (10) is proved to make sure that for each DMU under evaluation it supplies a Pareto-efficient reference point on the efficient frontier. This is achieved by incorporating input and output slacks into S-SBM model and controlling optimal values for each variable through the use of big M. For one specific input (output), input saving scaling factor (output surplus scaling factor) and input slack (output slack) are not permitted to appear at the same time. Especially, for a certain DMU, we can judge whether the reference point in S-SBM model is Pareto-efficient through our model. That is, the reference point 1+ti∗xik,1−βr∗yrk for DMUk in S-SBM model (7) is not Pareto-efficient if and only if there exists si′−∗>0 (i′≠i∈1,…,m) or sr′+∗>0 (r′≠r∈1,…,m) under model (10).



---

*Source: 6161343-2020-02-08.xml*