Estimation of Congestion in Free Disposal Hull Models Using Data Envelopment Analysis

Abstract

This paper deals with evaluating congestion in free disposal hull (FDH) models. There are several approaches in data envelopment analysis (DEA) literatures which discuss the theory and application of congestion. However, almost all of these approaches considered convex DEA technologies. So, in the case of nonconvex technologies, including FDH technology, this field is almost nil. This paper makes an attempt to fill in this void. To do so, this study provides a pairwise comparisons-based algorithm to evaluate congestion in FDH model. This algorithm identifies the sources of congestion and estimates its amounts. It is also capable of detecting the losses amounts of output due to congestion. The validity of the proposed model is demonstrated using some numerical and empirical examples.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Academic Editor:Mohammad Khodabakhshi

Department of Mathematics, Islamic Azad University, Science and Research Branch, Tehran 1477893855, Iran

Received 21 May 2014; Revised 14 July 2014; Accepted 20 July 2014; 15 October 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Evaluation of decision making units (DMUs) is an important task especially from a managerial point of view. DEA is a nonparametric and mathematical programming based approach for evaluating the performance of a set of homogeneous DMUs using multiple inputs to produce multiple outputs. In performance analysis, in particular in DEA, the concept of congestion plays a seminal role in theory and application. Congestion is a special phenomenon in the production process which is defined in economics where outputs are reduced due to excessive amount of inputs or an increase in one or more outputs results in a reduction in one or more inputs. For an actual example of congestion in a coal mine where a large crowd of the miners are working in a narrow underground, the amount of minerals excavated will be reduced [1].

Heretofore, various approaches have been presented in DEA for the treatment of congestion. The concept of congestion was first introduced in the literatures by Färe and Grosskopf [2] in the context of DEA. Subsequently an operationally implementable form was given by Färe et al. [3] and Cooper et al. [4-6]. Afterwards, Tone and Sahoo [7] developed a new slack-based approach to evaluate the scale elasticity in the presence of congestion with a unified framework. Wei and Yan [8] used DEA output additive models and proposed a necessary and sufficient condition for existence of congestion. Jahanshahloo and Khodabakhshi [9, 10] provided an approach of input congestion based on the relaxed combinations of inputs. Later on, Khodabakhshi [11] provided a one-model approach of input congestion based on input relaxation model. Also Khodabakhshi [12] proposed a method to detect the input congestion in the stochastic DEA. To see more references about this approach, the readers are referred to [13, 14]. Jahanshahloo et al. [15] and Khodabakhshi et al. [16] proposed some methods for computing the congestion in DEA models with production trade-offs and weight restrictions. Sueyoshi and Sekitani [17] proposed a modified approach which is able to measure congestion under the occurrence of multiple solution. There exist some papers which reviewed congestion papers, as that of Khodabakhshi et al. [18].

All of the above-mentioned investigations deal with congestion in convex technologies. In convex models, the targets resulting from efficiency assessment correspond to the points on the continuous efficiency frontiers. This means that DMUs might be compared with unreal DMUs which sometimes is meaningless in real life, for example, when we want to evaluate the efficiency of various car engines.

FDH models were first formulated by Deprins et al. [19]. The PPS of FDH model is made by deterministic (or observed activities) and free disposability postulates. So the PPS of FDH model is nonconvex. One appealing characteristic of FDH model due to nonconvexity nature of FDH efficiency frontier is that, in FDH model, targets correspond to observed units which is more compatible with real life because, in some circumstances, the observed unit is more comfortable when compared with a real unit rather than with a virtual one. As can be seen from the foregoing, there are several methods for evaluating congestion in convex DEA models, but for FDH models, although there are a few papers which are concerned with the field of estimation returns to scale (RTS), see, for example, [20-24], methods to estimate congestion can be hardly found. Therefore a new scheme is required to deal with congestion in FDH models.

In this paper, we first present definitions of output efficiency for DMUs under a series of DEA output additive models. Then, using these definitions, we develop a necessary and sufficient condition for existence of congestion in FDH model. Afterwards, we provide a polynomial time algorithm based on pairwise comparisons which evaluates congestion for DMUs using certain differences of inputs and outputs. This algorithm simply identifies the sources of congestion and estimates its amounts for congested DMUs.

The rest of the paper unfolds as follows. In the next section, FDH model and some of its properties and definitions will be presented to facilitate later discussions. In Section 3, we present a method with many computational advantages for evaluating congestion in FDH model. The validity of the proposed model is demonstrated using three numerical examples in Section 4. Finally, Section 5 gives the conclusion of this paper.

2. Preliminaries

In this section we first briefly describe some characteristic property of FDH model. Consider n DMUs where each _{DMU j} ( j = 1 , ... , n ) utilizes m inputs _{x i j} ( i = 1 , ... , m ) to produce s outputs _{y r j} ( r = 1 , ... , s ) . Let _{x j} = ^{( _{x 1 j} , ... , _{x m j} ) T} and _{y j} = ^{( _{y 1 j} , ... , _{y s j} ) T} . We will also assume that _{x j} ...5; 0 , _{x j} ...0; 0 and _{y j} ...5; 0 , _{y j} ...0; 0 . The production possibility set T is represented as [figure omitted; refer to PDF]

Deprins et al. [19] have deduced the following production possibility set. This set is denoted by _{T FDH} , regarding the assumptions of deterministic and free disposability of the production technology: [figure omitted; refer to PDF] The additive FDH model to evaluate the efficiency of a special _{DMU p} ( p ∈ { 1 , ... , n } ) under the _{T FDH} is as follows: [figure omitted; refer to PDF]

Definition 1 (FDH efficiency).

Consider model (3). If the optimal objective value is zero, then _{DMU p} is said to be FDH efficient.

It is worth noting that different to CCR and BCC models, the FDH model does not operate with the convexity assumption. Therefore, this model has a discrete nature which causes the efficient target point for an inefficient DMU simply to be assigned as a point among only actually observed DMUs. Thus, the efficiency analysis is done relative to the other given DMUs instead of a hypothetical efficiency frontier. This has the advantage that the achievement goal for an inefficient DMU given by its efficient target point will be more credible than in cases of CCR and BCC models.

Definition 2 (FDH output efficiency).

Consider the following model. If _{Z FDH} = 0 , then _{DMU p} is said to be FDH output efficient: [figure omitted; refer to PDF]

Definition 3 (congestion).

Evidence of congestion is present in the performance of any DMU, when a decrease in one or more inputs is associated with increases that are maximally possible in one or more outputs without worsening other inputs or outputs. Conversely, congestion is said to occur when some of the outputs that are maximally possible are reduced by increasing one or more inputs without improving any other inputs or outputs [25].

A very restrictive form of the above definition yields the definition of strong congestion as follows.

Definition 4 (strong congestion).

If a proportionate reduction in all inputs of a DMU warrants an increase in all maximally possible outputs, then strong congestion occurs [7].

Definition 5 (technical efficiency).

Efficiency is achieved by _{DMU 0} if and only if it is not possible to improve some of its inputs or outputs without worsening some of its other inputs or outputs [25].

Definition 6 (technical inefficiency).

Technical inefficiency is said to be present in the performance of _{DMU 0} when the evidence shows that it is possible to improve some input or output without worsening some other inputs or outputs [25].

3. Congestion in FDH Model

In _{T FDH} , the efficiency surface is a staircase based on those given DMUs that are not dominated by other given DMUs. Figure 1 describes an illustrative example of _{T FDH} which is made by eight DMUs denoted by A, B, ..., H with one input and one output.

Figure 1: _{T FDH} .

[figure omitted; refer to PDF]

It should be noted that evaluating congestion in customary models for convex PPS has been studied on _{T NEW} , which is a PPS without input disposability postulate. Let us denote _{T NEW} corresponding to _{T FDH} as _{T NFDH} , which can be defined as follows: [figure omitted; refer to PDF] Figure 2 exhibits _{T NFDH} for the example in Figure 1. As seen from Figure 2, _{T NFDH} has a discrete nature and so the study of congestion on it is complicated. So to overcome this difficulty we introduce a new set as follows: [figure omitted; refer to PDF] Apparently, the set of ^{FDH - 1} is gained by reversing the sign of input inequalities in _{T FDH} . ^{FDH - 1} set corresponding to the example in Figure 1 is illustrated in Figure 3.

Figure 2: _{T NFDH} .

[figure omitted; refer to PDF]

Figure 3: FD ^{H - 1} .

[figure omitted; refer to PDF]

We use the following model to deal with the congestion phenomenon in FDH model: [figure omitted; refer to PDF] We call the above model " ^{FDH - 1} output additive model."

To see what is involved, we note that the input (like the output) constraints take the form ^{∑ j = 1 n} ... _{λ j x i j} ...5; _{x i p} . Hence, in this adaptation of additive models, the objective is to maximize the outputs without reducing any of the inputs.

Definition 7 ( ^{F D H - 1} output efficiency).

Consider the model (7). If _{Z FD ^{H - 1}} = 0 , then _{DMU p} is said to be FD ^{H - 1} output efficient.

Lemma 8.

_{D M U p} is ^{F D H - 1} output efficient if and only if the following system has no solution: [figure omitted; refer to PDF]

Proof.

It is clear using Definition 7.

Definition 9 (congestion in FDH model).

Let _{DMU p} = ( _{x p} , _{y p} ) be FD ^{H - 1} output efficient; if there exists _{DMU k} = ( _{x k} , _{y k} ) , such that _{x k} ...4; _{x p} , _{x k} ...0; _{x p} and _{y k} ...5; _{y p} , _{y k} ...0; _{y p} , then _{DMU p} has evidence of congestion.

Based upon Definition 9, units F, G, and H in Figure 3 have evidence of congestion and unit C is technically inefficient.

Definition 10 (strong congestion in FDH model).

Let _{DMU p} = ( _{x p} , _{y p} ) be congested in FDH model; if there exists _{DMU k} = ( _{x k} , _{y k} ) , such that _{x k} < _{x p} and _{y k} > _{y p} , then _{DMU p} has evidence of strong congestion.

Lemma 11.

Let _{D M U p} be ^{F D H - 1} output efficient; then _{D M U p} has evidence of congestion if and only if the following system has a solution: [figure omitted; refer to PDF]

Proof.

Let _{DMU p} has evidence of congestion, so from Definition 9, there exists _{DMU k} = ( _{x k} , _{y k} ) , such that _{x k} ...4; _{x p} , _{x k} ...0; _{x p} and _{y k} ...5; _{y p} , _{y k} ...0; _{y p} . Thus _{λ k} = 1 and _{λ j} = 0 ( j = 1 , ... , n , j ...0; k ) is a solution of (9).

Conversely, suppose that λ ~ = ( 0 , ... , 0,1 , 0 , ... , 0 ) , whose qth component is one, is a solution of (9). So, we have _{x q} ...4; _{x p} and _{y q} ...5; _{y p} , _{y q} ...0; _{y p} . Also, certainly _{x q} ...0; _{x p} , since, according to the assumption of lemma, _{DMU p} is FD ^{H - 1} output efficient, so _{x q} = _{x p} contradicts Lemma 8. Hence, by Definition 9, _{DMU p} has evidence of congestion.

Lemma 12.

_{D M U p} is not FDH output efficient if and only if the following linear system has a solution: [figure omitted; refer to PDF]

Proof.

Using definition of FDH output efficiency, the proof is completed.

We now present the main result of the proposed method.

Theorem 13.

Let _{D M U p} be ^{F D H - 1} output efficient; then _{D M U p} has evidence of congestion if and only if _{D M U p} is not FDH output efficient.

Proof.

Using Lemmas 11 and 12, the proof is completed.

Now, using Theorem 13, we can provide the following procedure to evaluate congestion in FDH model.

(1) Solve model (7) corresponding to ( _{x p} , _{y p} ) ; let ( ^{λ *} , ^{s + *} ) be the optimal solution of it. Let _{y ^ p} = _{y p} + ^{s + *} . It is evident that ( _{x p} , _{y ^ p} ) is FD ^{H - 1} output efficient.

(2) Solve model (4) for ( _{x p} , _{y ^ p} ) .

(3) If _{Z FDH} > 0 , then _{DMU p} is congested.

Remark 14.

Models (4) and (7) are mixed-integer programming, but we can simply show that it does not need any mathematical programming problem to solve. Indeed, an enumeration algorithm based on pairwise comparisons, similar to Tulken's enumeration algorithm for the case of radial FDH model [26], can be used.

Now, based upon foregoing procedure and Remark 14, we propose the following algorithm. The proposed algorithm includes two parts. In Part (a), we recognize the existence of congesting in performance of _{DMU p} and in Part (b), if _{DMU p} is recognized to be congested in Part (a), the amount of congestion for each input as well as the reduction amount of each output due to congestion will be estimated.

Proposed Algorithm

Part (a)

Step 1.

Calculate the optimal value of model (7) by the following equation: [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Step 2.

Let _{y ^ p} = _{y p} + ^{s + *} , where ^{s + *} = _{y q} - _{y p} . Obtain the optimal value of model (4) by [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Step 3.

If _{Z FDH} > 0 , then _{DMU p} is congested, so go to Part (b); furthermore, if there exist j ∈ _{D ^ p} such that _{x j} < _{x p} and _{y j} > _{y ^ p} , then, based on Definition 4, _{DMU p} is strongly congested. If _{Z FDH} = 0 , then _{DMU p} is not congested and stop.

Part (b)

Step 4.

Define _{K p} as follows: [figure omitted; refer to PDF] Then calculate [figure omitted; refer to PDF]

Step 5.

Define _{T p} as follows: [figure omitted; refer to PDF] For j ∈ _{T p} define ^{s i c *} as the amount of congestion in i th input of _{DMU p} and ^{s ^ r + *} as reduction amount of r th output due to congestion as follows: [figure omitted; refer to PDF] ^{α *} = ^{∑ i = 1 m} ... ^{s i c *} is the total amount of congestion in all inputs of _{DMU p} .

Corollary 15.

If _{D ^ p} = ∅ , congestion has no appearance at _{D M U p} .

Proof.

In this case, it is obvious that system (10) has no solution. So, by Lemma 12, ( _{x p} , _{y ^ p} ) is FDH output efficient and regarding Theorem 13 congestion has no appearance at ( _{x p} , _{y ^ p} ) . Since congestion is a frontier concept, then _{DMU p} has no congestion

Remark 16.

If Card ( { j ∈ _{D p} |" _{Z FD ^{H - 1}} = ^{∑ r = 1 s} ... ( _{y r j} - _{y r p} ) } ) > 1 , then the projection of _{DMU p} , ( _{x p} , _{y ^ p} ) , is not determined uniquely.

Remark 17.

Card ( _{T p} ) ...5; 1 . If Card ( _{T p} ) = 1 , then the amount of congestion can be uniquely determined; otherwise there are alternatives for the amount of congestion for _{DMU p} .

Using the above algorithm, we can evaluate the congestion for each DMU in FDH model without any mathematical programming problems and with only some certain pairwise differences of inputs and outputs with regard to (11), (16), and (18). Therefore, the following theorem is obviously true.

Theorem 18.

The proposal algorithm is a polynomial time algorithm.

4. Numerical Examples

In this section, we apply our proposed procedure to measure the congestion effect on two numerical examples and an empirical example.

Example 1.

We consider the illustrative example provided in Section 3 which includes eight DMUs, A, B, C, D, E, F, G, and H, with one output and one input each, as shown in Figure 1. The data set of DMUs as well as the results of Part (a) of the proposed algorithm is given in Tables 1 and 2 displays results of Part (b). As shown in Table 1, there is no congestion in DMUs A, B, C, D, and E and congestion has appeared in DMUs F, G, and H. The input congestion amount and reduction amount of output due to congestion for congested units, resulting from the proposed algorithm, have been provided in the two last columns of Table 2, respectively. Obviously, in the case of one input and one output, each congested unit has evidence of strong congestion.

Table 1: Data set and results of Part (a) of the proposed algorithm in Example 1.

DMUs	I	O	D p	Z FDH - 1	s + *	y ^ p	D ^ p	Z FDH	Status
A	2	1	{A, B, C, D, E, F}	3	3	4	{A}	0	No congestion
B	3	3	{B, D, E, F}	1	1	4	{B}	0	No congestion
C	5	2	{C, E, F, H}	2	2	4	{ }		No congestion
D	6	4	{D, E}	0	0	4	{D}	0	No congestion
E	7	4	{E}	0	0	4	{D}	0	No congestion
F	8	3	{F}	0	0	3	{B, D, E}	1	Congestion
G	8	1	{F, G, H}	2	2	3	{B, D, E, F}	1	Congestion
H	9	2	{H}	0	0	2	{B, C, D, E, F}	2	Congestion

Table 2: Results of Part (b) of the proposed algorithm in Example 1.

DMUs	k p	α *	T p	s c *	s ^ + *
F	{E, D}	1	{E}	1	1
G	{E, D}	1	{E}	1	1
H	{E, D}	2	{E}	2	2

Example 2.

We consider an example adopted from Tone and Sahoo ([7], page 756) which has been listed in Table 3 of our study. This example consists of four DMUs, A, B, C, and D, using two inputs and producing two outputs. The results of Part (a) of the proposed algorithm are given in Table 3. From Table 3 we can see that there is no congestion in DMUs A and B. As shown in Table 4, for unit C, the congestion amount of inputs is ( ^{s 1 c *} , ^{s 2 c *} ) = ( 0,1 ) and ( ^{s ^ 1 + *} , ^{s ^ 2 + *} ) = ( 0,1 ) is output losses of C due to congestion. For unit D, the congestion amount of inputs is ( ^{s 1 c *} , ^{s 2 c *} ) = ( 1,1 ) and ( ^{s ^ 1 + *} , ^{s ^ 2 + *} ) = ( 1,1 ) is output losses of D due to congestion. Therefore considering Definition 10, unit D is strongly congested.

Table 3: Data set and results of Part (a) of the proposed algorithm in Example 2.

DMUs	( I 1 , I 2 )	( O 1 , O 2 )	D p	Z FDH - 1	( s 1 + * , s 2 + * )	y ^ p	D ^ p	Z FDH	Status
A	(1, 1)	(1, 1)	{A, B, C, D}	2	(1, 1)	(2, 2)	{ }		No congestion
B	(2, 2)	(2, 2)	{B}	0	(0, 0)	(2, 2)	{B}	0	No congestion
C	(2, 3)	(2, 1)	{C}	0	(0, 0)	(2, 1)	{B, C}	1	Congestion
D	(3, 3)	(1, 1)	{D}	0	(0, 0)	(1, 1)	{A, B, C, D}	2	Congestion

Table 4: Result of Part (b) of the proposed algorithm in Example 2.

DMUs	K p	α *	T p	( s 1 c * , s 2 c * )	( s ^ 1 + * , s ^ 2 + * )
C	{B}	1	{B}	(0, 1)	(0, 1)
D	{B}	2	{B}	(1, 1)	(1, 1)

Example 3.

Table 5 presents the input and output data of the Chinese textile industry during 1981-1997 assembled by Cooper et al. [5]. Each year has been treated as a DMU with two inputs and one output: labor ( _{X 1} ) measured in units of 1000 persons, capital ( _{X 2} ) measured in units of one million Ren Min Be (Chinese monetary unit), and output ( Y ) measured in units of one million Ren Min Be too. Note that the capital and output values have been adjusted to a 1991 base period to eliminate the impact of price variations. The congestion amount regarding FDH technology using the proposed method as well as the congestion amounts using Cooper et al.'s method [6] in BCC technology which are shown by ( ^{s 1 c *} , ^{s 2 c *} ) is provided in Table 6. As can be seen from Table 6, congestion appeared in performance of DMUs 8, 9, 10, and 15 in _{T FDH} . Besides DMUs 2, 8, 9, 10, 12, 13, and 15 have been recognized to be congested in _{T BCC} . Comparison of these two computational results shows that the DMUs which are congested in _{T FDH} are among the ones that are congested in _{T BCC} . That is, the set of congested DMUs in _{T FDH} is a subset of set of congested DMUs in _{T BCC} . It seems reasonable because _{T FDH} ⊆ _{T BCC} .

Table 5: Data set in Example 3.

DMU = year	Labor ( X 1 )	Capital ( X 2 )	Output ( Y )
DMU 1 = 1981	389.00	19.86	856.02
DMU 2 = 1982	412.30	21.16	866.85
DMU 3 = 1983	423.50	17.08	956.04
DMU 4 = 1984	417.30	18.10	1082.94
DMU 5 = 1985	570	12.61	1273.20
DMU 6 = 1986	600.50	13.45	1230.72
DMU 7 = 1987	641.10	15.91	1410.66
DMU 8 = 1988	715.30	23.72	1728.16
DMU 9 = 1989	736.00	25.97	2109.57
DMU 10 = 1990	745.00	18.24	2291.08
DMU 11 = 1991	756.00	14.40	2533.27
DMU 12 = 1992	743.00	17.50	2899.16
DMU 13 = 1993	684.00	25.08	3520.74
DMU 14 = 1994	691.00	25.45	4949.93
DMU 15 = 1995	673.00	29.35	4604.00
DMU 16 = 1996	634.00	23.05	4722.29
DMU 17 = 1997	596.00	25.02	4760.28

Table 6: Comparison results of congestion treatment in BCC and FDH technologies in Example 3.

DMU = year	( s 1 c * , s 2 c * ) in FDH	( s 1 c * , s 2 c * ) in BCC	Results of the proposed approach in T FDH	Results of Cooper et al.'sapproach in T BCC [6]
DMU 1 = 1981	(0, 0)	(0, 0)	No congestion	No congestion
D M U 2 = 1982	(0, 0)	(0, 0.72)	Congestion	No congestion
DMU 3 = 1983	(0, 0)	(0, 0)	No congestion	No congestion
DMU 4 = 1984	(0, 0)	(0, 0)	No congestion	No congestion
DMU 5 = 1985	(0, 0)	(0, 0)	No congestion	No congestion
DMU 6 = 1986	(0, 0)	(0, 0)	No congestion	No congestion
DMU 7 = 1987	(0, 0)	(0, 0)	No congestion	No congestion
DMU 8 = 1988	(81.30, 0.67)	(65.39, 0)	Congestion	Congestion
DMU 9 = 1989	(45.00, 0.52)	(45.00, 0)	Congestion	Congestion
DMU 10 = 1990	(2.00, 0.74)	(43.16, 0)	Congestion	Congestion
DMU 11 = 1991	(0, 0)	(0, 0)	No congestion	No congestion
DMU 12 = 1992	(0, 0)	(30.72, 0)	No congestion	Congestion
DMU 13 = 1993	(0, 0)	(1.79, 0)	No congestion	Congestion
DMU 14 = 1994	(0, 0)	(0, 0)	No congestion	No congestion
DMU 15 = 1995	(7.70, 4.33)	(0, 3.98)	Congestion	Congestion
DMU 16 = 1996	(0, 0)	(0, 0)	No congestion	No congestion
DMU 17 = 1997	(0, 0)	(0, 0)	No congestion	No congestion

5. Conclusion

In this paper we proposed a method based on pairwise comparison to evaluate congestion in FDH model. The results of the study have been proved with some lemmas and theorems. Our proposed method is able to identify congestion in performance of DMUs and it can determine the amount of excessive inputs for congested DMUs based on the calculation of certain pairwise differences of inputs and outputs. It also is capable of detecting the losses amounts of output due to congestion. One of the advantages of this method is that in all stages there is no need to solve any mathematical programming problems and so a polynomial time algorithm to identify congestion in FDH model is provided. Hence, it is superior from a computational point of view. The numerical examples demonstrated the compatibility of the proposed approach and so can be developed in performance analysis and large practical projects.

Because of having low complexity of computation, this method can be developed for imprecise data for further research.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] W. W. Cooper, L. Seiford, K. Tone Introduction to Data Envelopment Analysis and its Uses: With DEA-Solver Software and References , Springer, New York, NY, USA, 2005.

[2] R. Färe, S. Grosskopf, "Measuring congestion in production," Zeitschrift für Nationalökonomie , vol. 43, no. 3, pp. 257-271, 1983.

[3] R. Färe, S. Grosskopf, C. A. K. Lovell The Measurement of Efficiency of Production , Kluwer-Nijhoff Publishing, Boston, Mass, USA, 1985.

[4] W. W. Cooper, R. G. Thompson, R. M. Thrall, "Introduction: extensions and new developments in DEA," Annals of Operations Research , vol. 66, pp. 3-45, 1996.

[5] W. W. Cooper, H. Deng, B. Gu, S. Li, R. M. Thrall, "Using DEA to improve the management of congestion in Chinese industries (1981-1997)," Socio-Economic Planning Sciences , vol. 35, no. 4, pp. 227-242, 2001.

[6] W. W. Cooper, L. M. Seiford, J. Zhu, "A unified additive model approach for evaluating inefficiency and congestion with associated measures in DEA," Socio-Economic Planning Sciences , vol. 34, no. 1, pp. 1-25, 2000.

[7] K. Tone, B. K. Sahoo, "Degree of scale economies and congestion: a unified DEA approach," European Journal of Operational Research , vol. 158, no. 3, pp. 755-772, 2004.

[8] Q. L. Wei, H. Yan, "Congestion and returns to scale in data envelopment analysis," European Journal of Operational Research , vol. 153, no. 3, pp. 641-660, 2004.

[9] G. R. Jahanshahloo, M. Khodabakhshi, "Suitable comb ination of inputs for improving outputs in DEA with determining input congestion :considering textile indus try of China," Applied Mathematics and Computation , vol. 151, no. 1, pp. 263-273, 2004.

[10] G. R. Jahanshahloo, M. Khodabakhshi, "Determining assurance interval for non-Archimedean element in the improving outputs model in DEA," Applied Mathematics and Computation , vol. 151, no. 2, pp. 501-506, 2004.

[11] M. Khodabakhshi, "A one-model approach based on relaxed combinations of inputs for evaluating input congestion in DEA," Journal of Computational and Applied Mathematics , vol. 230, no. 2, pp. 443-450, 2009.

[12] M. Khodabakhshi, "Chance constrained additive input relaxation model in stochastic data envelopment analysis," International Journal of Information and Systems Sciences , vol. 6, no. 1, pp. 99-112, 2010.

[13] M. Khodabakhshi, "An output oriented super-efficiency measure in stochastic data envelopment analysis: considering Iranian electricity distribution companies," Computers and Industrial Engineering , vol. 58, no. 4, pp. 663-671, 2010.

[14] M. Asgharian, M. Khodabakhshi, L. Neralic, "Congestion in stochastic data envelopment analysis: an input relaxation approach," International Journal of Statistics and Management System , vol. 5, no. 1, pp. 84-106, 2010.

[15] G. R. Jahanshahloo, M. Khodabakhshi, F. H. Lotfi, M. R. Goudarzi, "Computation of congestion in DEA models with productions trade-offs and weight restrictions," Applied Mathematical Sciences , vol. 5, no. 13-16, pp. 663-676, 2011.

[16] M. Khodabakhshi, R. M. Goudarzi, M. Y. Maryaki, M. Hajiani, "A one-model approach for computation of congestion with productions trade-offs and weight restrictions," International Journal of Applied Operational Research , vol. 3, no. 4, pp. 69-80, 2013.

[17] T. Sueyoshi, K. Sekitani, "DEA congestion and returns to scale under an occurrence of multiple optimal projections," European Journal of Operational Research , vol. 194, no. 2, pp. 592-607, 2009.

[18] M. Khodabakhshi, F. Hosseinzadeh Lotfi, K. Aryavash, "Review of input congestion estimating methods in DEA," Journal of Applied Mathematics , vol. 2014, 2014.

[19] D. Deprins, L. Simar, H. Tulkens, M. Marchand, P. Pestieau, H. Tulkens, "Measuring labor-efficiency in post offices," The Performance of Public Enterprises , pp. 243-267, Elsevier Science, Amsterdam, The Netherlands, 1984.

[20] K. Kerstens, P. V. Eeckaut, "Estimating returns to scale using non-parametric deterministic technologies: a new method based on goodness-of-fit," European Journal of Operational Research , vol. 113, no. 1, pp. 206-214, 1999.

[21] V. V. Podinovski, "On the linearisation of reference technologies for testing returns to scale in FDH models," European Journal of Operational Research , vol. 152, no. 3, pp. 800-802, 2004.

[22] M. Soleimani-damaneh, G. R. Jahanshahloo, M. Reshadi, "On the estimation of returns-to-scale in FDH models," European Journal of Operational Research , vol. 174, no. 2, pp. 1055-1059, 2006.

[23] M. Soleimani-damaneh, M. Reshadi, "A polynomial-time algorithm to estimate returns to scale in FDH models," Computers and Operations Research , vol. 34, no. 7, pp. 2168-2176, 2007.

[24] M. Soleimani-Damaneh, A. Mostafaee, "Stability of returns to scale in FDH models," European Journal of Operational Research , vol. 196, pp. 1223-1228, 2009.

[25] W. W. Cooper, H. Deng, Z. M. Huang, S. X. Li, "A one-model approach to congestion in data envelopment analysis," Socio-Economic Planning Sciences , vol. 36, no. 4, pp. 231-238, 2002.

[26] H. Tulkens, "On FDH efficiency analysis: some methodological issues and applications to retail banking, courts, and urban transit," Journal of Productivity Analysis , vol. 4, no. 1-2, pp. 183-210, 1993.