Data envelopment analysis (DEA), as a useful management and decision tool, has been widely used since it was first invented by Charnes et al. in 1978. On the one hand, the DEA models need accurate inputs and outputs data. On the other hand, in many situations, inputs and outputs are volatile and complex so that they are difficult to measure in an accurate way. The conflict leads to the researches of uncertain DEA models. This paper will consider DEA in uncertain environment, thus producing a new model based on uncertain measure. Due to the complexity of the new uncertain DEA model, an equivalent deterministic model is presented. Finally, a numerical example is presented to illustrate the effectiveness of the uncertain DEA model.

1. Introduction

Data envelopment analysis is a mathematical programming technique that measures the relative efficiency of decision making units with multiple inputs and outputs, which was initialized by Charnes et al. [1]. This was followed by variety of theory research work, including Banker et al. [2], Charnes et al. [3], Petersen [4], and Tone [5]. More DEA papers can refer to Seiford [6] in which 500 references are documented.

The original DEA models assume that inputs and outputs are measured by exact values. However, in many situations, such as in a manufacturing system, a production process, or a service system, inputs and outputs are volatile and complex so that they are difficult to measure in an accurate way. Thus many researchers tried to model DEA with various uncertain theories. Probability theory is the earliest theory which was used to establish the stochastic DEA models. Sengupta [7] generalized the stochastic DEA model using the expected value. Banker [8] incorporated statistical elements into DEA, thus developing a statistical method. Many papers [913] have employed the chance-constrained programming to DEA in order to accommodate stochastic variations in data. Fuzzy theory is another theory which was used to deal with the uncertainty in DEA. As one of the DEA initiators, Cooper et al. [1416] introduced how to deal with imprecise data such as bounded data, ordinal data, and ratio bounded data in DEA. Kao and Liu [17] developed a method to find the membership functions of the fuzzy efficiency scores when some inputs or inputs are fuzzy numbers. Entani et al. [18] proposed a DEA model with an interval efficiency by the pessimistic and the optimistic values. Many researchers have introduced possibility measure [19] into DEA [20, 21].

A lot of surveys showed that human uncertainty does not behave like fuzziness. For example, we say “the input is about 10.” Generally, we employ fuzzy variable to describe the concept of “about 10;" then there exists a membership function, such as a triangular one (9, 10, 11). Based on this membership function, we can obtain that “the input is exactly 10" with possibility measure 1. On the other hand, the opposite event of “not exactly 10" has the same possibility measure. The conclusion that “not 10" and “exactly 10" have the same possibility measure is not appropriate. This inspired Liu [22] to found an uncertainty theory which has become a branch of axiomatic mathematics for modeling human uncertainty. This paper will apply the uncertainty theory to DEA to deal with human uncertainty, thus producing some uncertain DEA models.

In this paper, we will assume the inputs and outputs are uncertain variables and propose some uncertain DEA models. The rest of this paper is organized as follows. Section 2 will introduce some basic concepts and properties about uncertain variables. Then an uncertain DEA model as well as its equivalent crisp model will be presented in Section 3. Finally, a numerical example will be given to illustrate the uncertain DEA model in Section 4.

2. Preliminaries

Uncertainty theory was founded by Liu [22] in 2007 and refined by Liu [23] in 2010. As extensions of uncertainty theory, uncertain process, and uncertain differential equations [24], uncertain calculus [25] were proposed. Besides, uncertain programming was first proposed by Liu [26] in 2009, which wants to deal with the optimal problems involving uncertain variable. This work was followed by an uncertain multiobjective programming, an uncertain goal programming [27], and an uncertain multilevel programming [28]. Since that, uncertainty theory was used to solve variety of real optimal problems, including finance [2931], reliability analysis [32, 33], graph [34, 35], and train scheduling [36, 37]. In this section, we will state some basic concepts and results on uncertain variables. These results are crucial for the remainder of this paper.

Let be a nonempty set, and let be a -algebra over . Each element  Ł is assigned a number . In order to ensure that the number has certain mathematical properties, Liu [22] presented the four axioms.

Axiom 1. for the universal set .

Axiom 2. for any event .

Axiom 3. For every countable sequence of events , we have

Axiom 4. Let be uncertainty spaces for . Then the product uncertain measure is an uncertain measure satisfying where are arbitrarily chosen events from for , respectively.

If the set function satisfies the first three axioms, it is called an uncertain measure.

Definition 1 (see Liu [22]). Let be a nonempty set, let be a -algebra over , and let be an uncertain measure. Then the triplet is called an uncertainty space.

Definition 2 (see Liu [22]). An uncertain variable is a measurable function from an uncertainty space to the set of real numbers; that is, for any Borel set of real numbers, the set is an event.

Definition 3 (see Liu [22]). The uncertainty distribution of an uncertain variable is defined by for any real number .

Example 4. The linear uncertain variable has an uncertainty distribution

Example 5. An uncertain variable is called zigzag if it has a zigzag uncertainty distribution denoted by , where are real numbers with .

Definition 6 (see Liu [25]). The uncertain variables are said to be independent if for any Borel sets .

Definition 7 (see Liu [23]). An uncertainty distribution of an uncertain variable is said to be regular if its inverse function exists and is unique for each . In this case, the inverse function is called the inverse uncertainty distribution of .

Example 8. The inverse uncertainty distribution of a zigzag uncertain variable is

Theorem 9 (see Liu [23]). Let be independent uncertain variables with regular uncertainty distributions , respectively. If is a strictly increasing function, then is an uncertain variable with inverse uncertainty distribution

Example 10. Let be an uncertain variable with regular uncertainty distribution . Since is a strictly increasing function for any constants and , the inverse uncertainty distribution of is

Example 11. Let be independent uncertain variables with regular uncertainty distributions , respectively. Since is a strictly increasing function, the sum is an uncertain variable with inverse uncertainty distribution

Theorem 12 (see Liu [23]). Assume the constraint function is strictly increasing with respect to and strictly decreasing with respect to . If are independent uncertain variables with uncertainty distributions , respectively, then the chance constraint holds if and only if

3. DEA Model

In many situations, inputs and outputs are volatile and complex so that they are difficult to measure in an accurate way. This inspired many researchers to apply probability to DEA. As we know, probability or statistics needs a large amount of historical data. In the vast majority of real cases, the sample size is too small (even no sample) to estimate a probability distribution. Then we have to invite some domain experts to evaluate their degree of belief that each event will occur. This section will give some researches to empirical uncertain DEA using the theory introduced in Section 2. The new symbols and notations are given as follows: : the th DMU,   ;DMU0: the target DMU;   : the uncertain inputs vector of ,   ; : the uncertainty distribution of ,   ,   ; : the inputs vector of the target DMU0; : the uncertainty distribution of ,   ; : the uncertain outputs vector of ,   ; : the uncertainty distribution of ,   ,   ; : the outputs vector of the target DMU0; : the uncertainty distribution of ,   .

3.1. Uncertainty Distributions of Inputs and Outputs

Liu and Ha [38] proposed a questionnaire survey for collecting expert’s experimental data. It is based on expert’s experimental data rather than historical data. The starting point is to invite one expert who is asked to complete a questionnaire about the meaning of an uncertain input (output) like “How many is the input (output)."

We first ask the domain expert to choose a possible value that the uncertain input may take and then quiz him,“How likely is less than or equal to ?"

Denote the expert’s belief degree by . An expert’s experimental data is thus acquired from the domain expert.

Repeating the above process, we can obtain the following expert’s experimental data: that meet the following consistence condition (perhaps after a rearrangement):

Based on those expert’s experimental data, Liu and Ha [38] suggested an empirical uncertainty distribution,

Assume there are domain experts and each produces an uncertainty distribution. Then we may get uncertainty distributions . The Delphi method was originally developed in the 1950s by the RAND Corporation based on the assumption that group experience is more valid than individual experience. Wang et al. [39] recast the Delphi method as a process to determine the uncertainty distribution. The main steps are listed as follows.

Step 1. The domain experts provide their expert’s experimental data,

Step 2. Use the th expert’s experimental data to generate the th expert’s uncertainty distribution .

Step 3. Compute , where are convex combination coefficients.

Step 4. If are less than a given level , then go to Step 5. Otherwise, the th expert receives the summary ( and reasons) and then provides a set of revised expert’s experimental data. Go to Step 2.

Step 5. The last is the uncertainty distribution of the input (output).

3.2. Uncertain DEA Model

Similar to traditional DEA model [3], the objective of the uncertain DEA model is to maximize the total slacks in inputs and outputs subject to the constraints. Then the uncertain DEA model can be given as follows:

Definition 13 ( -efficiency). DMU0 is -efficient if and are zero for and , where and are optimal solutions of (21).

Since the uncertain measure is involved, this definition is different from traditional efficiency definition. For instance, as determined by the choice of , there is a risk that DMU0 will not be efficient even when the condition of Definition 13 is satisfied.

Since is one of the , we can always get a solution with ,    , and all slacks zero. Thus this uncertain DEA model has feasible solution and the optimal value for all .

3.3. Deterministic Model

Model (21) is an uncertain programming model, which is too complex to compute directly. This section will give its equivalent crisp model to simplify the computation process.

Theorem 14. Assume that are independent uncertain inputs with uncertainty distribution for each , , and are independent uncertain outputs with uncertainty distribution for each , . Then holds if and only if

Proof. Without loss of generality, let and ; then we will consider the equation
Rewrite (24) as
Since is an uncertain variable which is decreasing with respect to , its inverse uncertainty distribution is
For each ,   is an uncertain variable whose inverse uncertainty distribution is
It follows from the operational law that the inverse uncertainty distribution of the sum is From which we may derive the result immediately for and . Similarly, we can get other results.
Following Theorem 14, the uncertain DEA model can be converted to the crisp model as follows: which is a linear programming model. Thus it can be easily solved by many traditional methods.

4. A Numerical Example

This example wants to illustrate the uncertain DEA model. For simplicity, we will only consider five DMUs with two inputs and two outputs which are all zigzag uncertain variables denoted by . Table 1 gives the information of the DMUs.

For illustration, let DMU1 be the target DMU; then the uncertain DEA model (29) can be written as

Table 2 shows the results of evaluating DMUs with confidence level . The results can be interpreted in the following way: DMU1 and DMU3 are inefficient, whereas DMU2, DMU4, and DMU5 are efficient. Moreover, DMU3 is more efficient than DMU1 from the total slacks , since they are both inefficient.

Uncertain efficiencies obtained from model (30) for different confidence levels are shown in Table 3. DMU1 is inefficient at all confidence levels, whereas DMU2, DMU4, and DMU5 are always efficient at all levels. It can be seen that the number of the efficient DMUs is affected by the confidence level . The higher the confidence level is, the bigger the number of efficient DMUs is. This phenomena indicate that uncertain DEA is more complex than the traditional DEA because of the inherent uncertainty contained in inputs and outputs.

5. Conclusion

Due to its widely practical used background, data envelopment analysis (DEA) has become a pop area of research. Since the data cannot be precisely measured in some practical cases, many papers have been published when the inputs and outputs are uncertain. This paper has given some researches to uncertain DEA model. A new DEA model as well as its equivalent deterministic model was presented. For illustration, a numerical example was designed.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.


This work was supported by National Natural Science Foundation of China (nos. 71201005 and 61104132).