Next Article in Journal
A Brain-Inspired Hyperdimensional Computing Approach for Classifying Massive DNA Methylation Data of Cancer
Next Article in Special Issue
Efficient Rule Generation for Associative Classification
Previous Article in Journal
A Class of Spline Functions for Solving 2-Order Linear Differential Equations with Boundary Conditions
Previous Article in Special Issue
Diagnosis in Tennis Serving Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Ensemble and Dimensionality Reduction DEA Models Based on Entropy Criterion

by
Parag C. Pendharkar
Information Systems School of Business Administration, Pennsylvania State University, Harrisburg 777 West Harrisburg Pike, Middletown, PA 17057, USA
Submission received: 28 July 2020 / Revised: 14 September 2020 / Accepted: 14 September 2020 / Published: 16 September 2020
(This article belongs to the Special Issue Algorithms in Decision Support Systems)

Abstract

:
Dimensionality reduction research in data envelopment analysis (DEA) has focused on subjective approaches to reduce dimensionality. Such approaches are less useful or attractive in practice because a subjective selection of variables introduces bias. A competing unbiased approach would be to use ensemble DEA scores. This paper illustrates that in addition to unbiased evaluations, the ensemble DEA scores result in unique rankings that have high entropy. Under restrictive assumptions, it is also shown that the ensemble DEA scores are normally distributed. Ensemble models do not require any new modifications to existing DEA objective functions or constraints, and when ensemble scores are normally distributed, returns-to-scale hypothesis testing can be carried out using traditional parametric statistical techniques.

1. Introduction

Data envelopment analysis (DEA) is a prominent technique for the non-parametric relative efficiency analysis of a set of decision-making units (DMUs) drawn from a similar production process [1]. DEA models are used in both operation research and data mining literature [2]. Some of the traditional properties of production functions, such as the monotonicity and convexity of the inputs and outputs, that are fundamental in DEA models are often found to be attractive in some data mining models where datasets are noisy and model resistance to learning noise is necessary [3]. An important aspect of DEA models is the reliability of DMU efficiency scores. It is generally accepted that the DEA efficiency estimates are reliable when the sample size is large [4]. Since the reliability of the DEA scores is dependent on the sample size, Cooper et al. [5] have suggested the following rule for the minimum number (n) of DMUs for reliable DEA analysis (each DMU has m inputs and s outputs):
n   m a x { 3 ( m + s ) ,   m × s }
For small-size datasets, where violations of the minimum number of DMUs specified by Equation (1a) frequently occur, dimensionality reduction (also known as variable reduction or variable selection) approaches are frequently used to select a subset of variables to satisfy Equation (1a). A variety of variable selection approaches are available in the literature. Among these variable selection approaches are statistical [6], regression [7], efficiency contribution measure [8], bootstrapping [9], hypothesis testing [10], variable aggregation [11] and statistical experiment designs [12]. Variable selection approaches are criticized extensively for applying parametric procedures and linear relationship assumptions for selecting variables to determine an unknown non-linear and non-parametric efficiency frontier. Nataraja and Johnson [13] provide a good description of some of these procedures and their pros and cons.
Pendharkar [14] proposed a competing approach to the dimensionality reduction/variable selection problem called the ensemble DEA. In his approach, traditional DEA analysis is conducted for all possible input and output combinations, and the efficiency scores of each DEA model for each DMU are averaged as an ensemble efficiency score for a DMU. Drawing from machine learning literature, Pendharkar [14] showed that the ensemble efficiency score is a reliable estimate of the “true” efficiency of a DMU. Even for small datasets, certain combinations of inputs will satisfy the criterion set by Equation (1a), while others will violate it, but the average ensemble score will be closer to the true efficiency of the DMU and will be reliable. Pendharkar [14] also proposed an exhaustive search procedure to generate all possible input and output combinations, and proposed a formula to compute the number of unique DEA models that need to be run to compute an average ensemble score. This number N of unique DEA models may be computed using the following formula:
N = ( i = 1 m ( m i ) ) × ( i = 1 s ( s i ) ) = ( 2 m 1 ) × ( 2 s 1 ) .
Using Banker et al.’s [15] variable-returns-to-scale (VRS) DEA BCC model, and data and models obtained from a few studies in the literature, Pendharkar [14] showed that the ensemble DEA model provides a better ranking of DMUs than the models proposed in a few studies from the literature.
This research investigates the additional properties and statistical distribution of the ensemble DEA model scores. It is shown that there are added benefits of ensemble efficiency scores. In particular, the ensemble efficiency scores maximize entropy, meaning that the DMU ranking distribution generated by the ensemble efficiency scores has a lower bias when compared to some competing radial and non-radial variable selection models recently reported in the literature, and second, the ensemble efficiency scores may be normally distributed under certain restrictive assumptions. The normal distribution of the efficiency score feature is particularly attractive because returns-to-scale hypothesis testing may be conducted by using traditional difference-in-means parametric statistical procedures. Both of these features are tested using data and models reported in a published study [16]. The rest of the paper is organized as follows: In Section 2, the basic DEA radial and non-radial models, ensemble DEA model and Entropy criterion for comparing different DEA models are described. In Section 3, using Iranian gas company data, the results of ensemble DEA models are compared with the results of variable selection models used in Toloo and Babaee’s [16] study. Additionally, in Section 3, the properties of the ensemble DEA scores are investigated in terms of the entropy criterion and their statistical distributions. In Section 4, the paper concludes with a summary and directions for future research.

2. DEA Preliminaries, Ensemble DEA Model, Entropy Criterion for DEA Model Comparisons and Statistical Distribution of Ensemble Scores

The basic DEA model assumes n DMUs, with each DMU consisting of m different inputs that produce s different outputs. The input and output vectors are semi-positive, and for DMUj (j = 1, …, n), the space for the input and output vectors (xj,yj) ϵ + m + s . For a DMUo, its relative efficiency may be computed by using the linear programming model under the constant returns-to-scale assumption. This efficiency is computed by solving the following model:
max r = 1 s u r y r o ,
subject to:
i = 1 m v i x i o = 1
r = 1 s u r y r j i = 1 m v i x i j 0       f o r   a l l   j = 1 , , n
v i , u r ε   f o r   a l l   i = 1 , , m   a n d   r = 1 , , s
where vi and ur are the weights associated with the ith input and jth output, respectively. The constant ε > 0 is infinitesimally non-Archimedean. The model (2a)–(2d) is often called the primary CCR model [1], and its dual is written as follows:
m i n i m i z e   θ ε ( i = 1 m s i + r = 1 s s r + ) ,  
subject to:
j = 1 n λ j x i j + s i = θ x i o ,       i = 1 , , m
j = 1 n λ j y r j s r + = y r o ,       r = 1 , , s ,   and
λ j , s i , s r + 0   f o r   a l l   i = 1 , , m ; j = 1 , , n ; r = 1 , , s
The VRS BCC model augments the system (2e)–(2h) by adding the following constraint:
j = 1 n λ j = 1
The aforementioned models are radial DEA models that are criticized for not providing input or output projections (for inefficient DMUs) that satisfy Pareto optimality conditions [17]. Fare and Lovell [18] independently proposed radial DEA models that allow for input or output reductions at variable rates. The radial version of the CCR model is mathematically represented in the following dual form:
m i n i m i z e   1 m i = 1 m θ i
subject to:
j = 1 n λ j x i j θ i x i o ,       i = 1 , , m
j = 1 n λ j y r j y r o ,         r = 1 , , s
θ i 1 ,       i = 1 , , m    
λ i 0 ,       j = 1 , , n
Pendharkar [14] proposed an ensemble DEA model based on the popularity of ensemble models in machine learning literature. The ensemble DEA model requires an exhaustive search procedure using a binary vector z whose components indicate whether an input or output is considered in performing DEA analysis. The dimension of this binary vector is (m + s). Figure 1 illustrates the z vector and exhaustive search tree for two-input-and-one-output datasets. The exhaustive tree is pruned (dotted edges) for models that have either no inputs or no outputs. DEA analysis is then conducted on the remaining models, and the efficiency results of each model for each DMU are averaged and used as ensemble DEA scores. To illustrate the ensemble DEA approach on a two-input-and-one-output dataset, a CCR DEA analysis using partial Cobb–Douglas production function data on US economic growth between 1899 and 1910 [19] is conducted. Table 1 illustrates the results of our DEA analysis and resulting ensemble scores. The two inputs were labor in person-hours worked per year and the amount of capital invested. The output was the total annual production. The results of the analysis show that the traditional DEA with z = [111] does not provide unique rankings (for the years 1901 and 1902 receive the same efficiency score), but the ensemble DEA model provides unique DMU rankings. Pendharkar’s [14] study provides a theoretical basis for the reliability of ensemble DEA scores.
The maximum entropy (ME) principle has been applied to DEA DMU ranking distribution [20] and model comparisons [21]. The ME principle measures the DMU ranking bias by using a more general family of distributions [22]. Several statistical distributions can be characterized as ME densities [23]. The ME distributions are the least biased distributions obtained by imposing moment constraints that are inherent in the data [21]. To obtain the ME for a given set of DMUs and their efficiencies, normalized ranks are first obtained by computing θ i * i = 1 n θ i * , for each DMU, and then computing the ME for a certain model z as follows:
ME z = i = 1 n ( θ i * i = 1 n θ i * ) l n ( θ i * i = 1 n θ i * )
The ME for the DEA models in Table 1 are ME111 = 2.4768, ME101 = 2.4775 and ME011 = 2.4757. The model with labor as an input and production as an output (z = [101]) has the highest entropy and has the least bias, with a maximum difference between DMU efficiencies for closely ranked DMUs for the years 1901 and 1902. The ensemble entropy is 2.4769, and since it is an average of all z-vector combinations, the comparison benchmark for ensemble entropy is the model with z = [111]. The ensemble entropy is higher than the benchmark. The highest possible entropy value or upper bound (UB) for a model is given by the following expression:
ME UB = n × ( ( 1 n ) l n ( 1 n ) )
The MEUB for the data in Table 1 is 2.485, and the ensemble entropy is very close to the maximum value. It is important to note that obtaining the maximum value is not always desirable, but it provides a theoretical benchmark estimate for a completely unbiased normalized DMU score distribution.
To compute ensemble efficiency scores, an n × m matrix E of DEA efficiency scores is necessary. The rows of such a matrix are the numbers of DMUs, and the columns are the numbers of models given by the numbers of eligible models considered in computing ensemble efficiency scores. This number of eligible models will have an upper bound given by N, computed using Equation (1b). The elements of this matrix will be efficiency scores for each DMU computed for a given model identified by column number. Figure 2 illustrates a five-DMU-and-five-model matrix. The ensemble efficiency score ( θ i E ) for each DMU is computed using the following formula:
θ i E = j = 1 m θ i j * m
A few observations can be made about any row i ϵ {1, …, n} of the ensemble efficiency score matrix. First, all the elements of a given row are an independent computation of efficiency scores by the same DMU under a different model number with its unique set of input(s) and output(s). Second, in all the elements of a given row, the DMU is maximizing its efficiency given its model constraints. Thus, each row represents independent evaluations by a DMU under the maximum decisional efficiency (MDE) principle [24]. The MDE principle was introduced by Troutt [25] to develop a function to the aggregate the performance of multiple decision-makers. The underlying assumption of the MDE principle is that all decision-makers seek to maximize their decisional efficiencies. Troutt [26] later used the MDE approach to rank DMUs and showed that DMUs deemed efficient under MDE are also efficient when ranked using the DEA. For a linear aggregator function, such as the one used in Equation (2j), Troutt [26] illustrated that the decisional efficiencies θ can be described by the following probability density function (pdf):
g ( θ ) = c α e α θ ,   α   >   0   and   θ     [ 0 , 1 ]
The pdf in (2k) is monotone, increasing on its interval with a mode at θ = 1 (see Figure 5 for illustration). Using the laws of probability, the value of cα = α (eα − 1)−1. Since each element in a given row of the ensemble efficiency score matrix is an independent evaluation by a decision-maker (i.e., a DMU in an ensemble model) trying to maximize its decisional efficiency θ i j * for j = {1, …, m}, the probability density function for each row (DMU) can be written as:
g ( θ i ) = c α i e α i θ i ,   α i   >   0   and   θ i   ϵ   [ 0 , 1 ]
The central limit theorem mentions that the cumulative distribution functions (cdfs) of the sums of independently identically distributed random variables asymptotically converge to a Gaussian cdf. The ensemble efficiency scores are normalized sums of independent efficiency assessments that will be distributed with a pdf given by (2l). These sums can be considered independent and identically distributed if α1 = α2 = … = αn. Under the restrictive assumption that α1 = α2 = … = αn, the ensemble efficiency scores are guaranteed to asymptotically converge to a normal distribution by the central limit theorem. In practice, however, the ensemble efficiency scores are not entirely random or perfectly identically distributed (due to the slight likely variation of Equation (2l)’s αi parameters for each row), and each ensemble model does introduce a degree of mild randomization. For mild differences in the row pdf parameters αi, where α1α2 ≈ … ≈ αn, the ensemble efficiency scores are likely to be normally distributed. A reader may note that under ideal conditions, where α1 = α2 = … = αn and individual DMU scores follow Equation (2l)’s distribution, the entropy of the ensemble scores will be highest and close to the highest upper bound mentioned in Equation (2i) because the distribution in Equation (2i) has a mode of 1 (see Figure 5). Thus, it may be argued that the likelihood of normality of the ensemble scores increases when the entropy of the ensemble scores is closer to its upper bound given by Equation (2i). It is important to note that an entropy equal to the exact value of the upper bound given by Equation (2i) is undesirable because at that value, the distribution is a uniform distribution where all the DMUs are fully efficient for all the models. The entropy of the pdf in Equation (2k) is maximized on the interval [0, 1] when the mean of the distribution is greater than 0.5 [27]. Additionally, another important aspect of the distribution of the ensemble efficiency scores is that both the rows and columns of ensemble efficiency scores (Figure 2) play a role in the pdf of the ensemble efficiency scores because the rows represent sampling from the MDE distributions and the columns represent sampling from the distribution of the sums of independent variables. Larger sample sizes increase the statistical reliability and robustness of the results.

3. Comparing Variable Selection Models and Ensemble Model Using Gas Company Data and Entropy Criterion

For small datasets, many input or output variables are aggregated so that the selected variables satisfy the heuristic given in Equation (1a). There are two problems with all the variable selection approaches. First, they use an artificial criterion to select variables for a non-linear and non-parametric approach. Any artificial/subjective criterion will make some assumptions that are harder to justify. Second, these techniques have several selection parameters and thresholds that often lead to inconsistencies in applying these techniques. For example, Toloo and Babaee [16] illustrate three problems with a variable selection approach and suggested an improved approach. By contrast, the ensemble DEA approach does not make any assumptions, and for small datasets, trying out different input and output combinations and aggregating efficiency scores provide more reliable efficiency estimates than variable selection models. Part of the reason for the stability of ensemble DEA efficiency scores is that, even for small datasets, some DEA models in an ensemble will always satisfy the heuristic given in Equation (1a), which will increase the reliability of the ensemble efficiency scores due to model averaging. This stability of ensemble efficiency scores is illustrated by comparing ensemble scores with the results of models from Toloo and Babaee’s [16] study and using the entropy criterion.
To compare the results, the dataset from Toloo and Babaee’s [16] study is used. The dataset consists of three inputs and four outputs from an Iranian gas company. The inputs are budget (x1), staff (x2) and cost (x3). The outputs are customers (y1), the length of the gas network (y2), the volume delivered (y3) and gas sales (y4). Table 2 lists these data. Table 3 lists the efficiency scores of the ensemble DEA with the CCR and BCC models and models used by Toloo and Babaee [16]. Using formula (1b), a total of 105 unique DEA models were used to compute the DEA ensemble efficiency score.
The entropies of the Ensemble CCR, Ensemble BCC, Non-Radial and Radial models were 2.616, 2.621, 2.615 and 2.599, respectively. The MEUB from Equation (2i) is 2.639. Comparing the Ensemble CCR with the Non-Radial and Radial CCR models shows that the Ensemble CCR model has a higher entropy. Only the VRS Ensemble BCC model has a higher entropy than the Ensemble CCR model. The standard deviations of the Ensemble BCC model are mostly lower than the CCR model’s as well. More importantly, the Ensemble CCR model generates unique rankings for the DMUs, whereas the Non-Radial and Radial models generate a tie for three DMUs. The Ensemble BCC model also generates unique rankings, but the differences occur at the third decimal place. The Ensemble BCC efficiency scores for DMU 10, 12 and 13 were 0.960, 0.959 and 0.962, respectively.
Figure 3 and Figure 4 illustrate the numbers of models (out of 105 total models) where a DMU was fully efficient. These figures are useful for understanding to what extent the assumption α1α2 ≈ … ≈ αn was satisfied for the theoretical normal distribution of the ensemble efficiency scores. For these parameters to be similar, the expectation is that a similar number of fully efficient DMUs should exist across all models. Clearly, some DMUs are never fully efficient under any of 105 models and the assumption of identical distributions is violated. While the assumption is violated, Figure 4 illustrates that some DMUs, e.g., 1, 8, 10, 12 and 13, have a somewhat similar number of fully efficient DMUs to others. These ensemble scores of these DMUs may be considered as normalized random sums generated from identical distributions (such as Distribution 1). All of these DMUs have ensemble efficiency scores greater than 0.95. Similarly, DMUs 5, 6 and 11, in Figure 4, have no fully efficient scores, and these may also be considered as random normalized sums generated from identically distributed pdfs (such as Distribution 2).
The ensemble scores for this dataset appear to be random normalized sums from two or more pdfs of the forms given in Equation (2k). Given that these are independent random normalized sums, it can be easily shown that the product of two or more independent MDE pdfs is also an MDE pdf. Figure 5 illustrates two sample MDE pdfs for two different values of alpha. The entropy of an MDE pdf is maximized when the mean of a distribution is greater than 0.5 [27]. For the ensemble BCC model, from Table 3, this criterion is satisfied. The lowest value of the ensemble BCC score is 0.57, which is greater than the mean of 0.5 required to maximize entropy and higher than the lowest values for the efficiency scores for the radial, non-radial and ensemble CCR models. As a result, the ensemble BCC model appears to maximize its entropy slightly better than the ensemble CCR model.
While ensemble scores have a minor violation of an identical distribution for some DMUs, a formal test of the normality of the distribution of the ensemble efficiency scores was conducted. Table 4 illustrates the results of these tests. The Shapiro–Wilk statistic for the Ensemble CCR model is 0.944, and that for the Ensemble BCC model is 0.876, which, at 14 degrees of freedom, are non-significant, consistent the null hypothesis that the efficiency score distribution is normally distributed at the 95% level of statistical significance.
A paired sample t-test for the difference in mean efficiency scores between the Ensemble CCR and the Ensemble BCC models gives a |t|-value of 3.524, which is significant at the 99% level of statistical significance (df = 13), indicating that a variable returns-to-scale relationship exists between inputs and outputs. The normality of the ensemble efficiency score distributions increases the power of parametric statistical tests.

4. Summary, Conclusions and Directions for Future Work

A significant amount of research in the DEA literature has focused on dimensionality reduction/variable selection techniques for small datasets. These techniques are often criticized and have their limitations, with no clear way of selecting which technique is the best. A better approach would be to use an ensemble DEA score that does not make any additional assumptions and provides models that have high entropy values and normally distributed scores under restrictive assumptions. Pendharkar [14], in his study, has already provided a theoretical foundation for the reliability of ensemble DEA scores. The added benefit of ensemble DEA scores is that they provide unique DMU rankings.
The normality of ensemble DEA scores is not guaranteed unless the ensemble DEA scores are normalized sums generated from independent identically distributed MDE pdfs. This assumption may not be strictly satisfied in most real-world datasets, but the current study shows that minor deviation from this assumption may be tolerated because the entropy of all MDE pdfs is maximized when normalized sums have a value greater than 0.5. This means that, typically, the differences in means between the underlying pdfs (Equation (2l)) for ensemble entropy scores will be less than 0.5, and, while these pdfs may not be identically distributed, the means of these distributions will be close, resulting in the likely normal distribution of ensemble scores in most real-world cases. The normality of ensemble DEA scores allows for the application of traditional statistical tests for return-of-scales hypothesis tests. Traditional DEA hypothesis-testing methods are not perfect and are known to be slightly biased [28]. Future research may focus on comparing ensemble DEA-based hypothesis testing with traditional DEA hypothesis testing to identify which method provides reliable results.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. Pendharkar, P.C. Data envelopment analysis models for probabilistic classification. Comput. Ind. Eng. 2018, 119, 181–192. [Google Scholar] [CrossRef]
  3. Pendharkar, P. A data envelopment analysis-based approach for data preprocessing. IEEE Trans. Knowl. Data Eng. 2005, 17, 1379–1388. [Google Scholar] [CrossRef]
  4. Ruggiero, J. A new approach for technical efficiency estimation in multiple output production. Eur. J. Oper. Res. 1998, 111, 369–380. [Google Scholar] [CrossRef]
  5. Cooper, W.W.; Seiford, L.M.; Tone, K. Data Envelopment Analysis: A Comprehensive Text with Models, Applications, References and DEA-Solver Software, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  6. Ueda, T.; Hoshiai, Y. Application of principal component analysis for parsimonious summarization of DEA inputs and/or outputs. J. Oper. Res. Soc. Jpn. 1997, 40, 466–478. [Google Scholar] [CrossRef] [Green Version]
  7. Ruggiero, J. Impact assessment of input omisson on DEA. J. Inf. Technol. Decis. Mak. 2005, 4, 359–368. [Google Scholar] [CrossRef]
  8. Wagner, J.M.; Shimshak, D.G. Stepwise selection of variables in data envelopment analysis: Procedures and managerial perspectives. Eur. J. Oper. Res. 2007, 180, 57–67. [Google Scholar] [CrossRef]
  9. Simar, L.; Wilson, P.W. Testing restrictions in nonparameteric efficiency models. Commun. Stat. 2001, 30, 159–184. [Google Scholar] [CrossRef]
  10. Banker, R.D. Hypothesis tests using data envelopment analysis. J. Prod. Anal. 1996, 7, 139–159. [Google Scholar] [CrossRef]
  11. Amirteimoori, A.R.; Despotis, D.K.; Kordrostami, S. Variable reduction in data envelopment analysis. Optimization 2012, 63, 735–745. [Google Scholar] [CrossRef]
  12. Morita, H.; Avkiran, N.K. Selecting inputs and outputs in data envelopment analysis by designing statistical experiments. J. Oper. Res. Soc. Jpn. 2009, 52, 163–173. [Google Scholar] [CrossRef] [Green Version]
  13. Nataraja, N.R.; Johnson, A.L. Guidelines for using variable selection techniques in data envelopment analysis. Eur. J. Oper. Res. 2011, 215, 662–669. [Google Scholar] [CrossRef]
  14. Pendharkar, P.C. Ensemble based ranking of decision making units. Inf. Syst. Oper. Res. 2013, 51, 151–159. [Google Scholar] [CrossRef]
  15. Banker, R.D.; Charnes, A.; Cooper, W.W. Some models for estimating technical and scale inefficiencies in data envelopment analysis. Manag. Sci. 1984, 30, 1078–1092. [Google Scholar] [CrossRef] [Green Version]
  16. Toloo, M.; Babaee, S. On variable reductions in data envelopment analysis with an illustrative application to a gas company. Appl. Math. Comput. 2015, 270, 527–533. [Google Scholar] [CrossRef]
  17. Ray, S.C. Data Envelopment Analysis: Theory and Techniques for Economics and Operations Research; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  18. Fare, R.; Lovell, C.A.K. Measuring the technical efficiency. J. Econ. Theory 1978, 19, 150–162. [Google Scholar] [CrossRef]
  19. Stewart, J. Multivariable Calculus: Concepts and Contexts, 3rd ed.; Thomson Learning: London, UK, 2005. [Google Scholar]
  20. Soleimani-Damaneh, M.; Zarepisheh, M. Shannon’s entropy for combining the efficiency results of different DEA models: Method and Application. Expert Syst. Appl. 2009, 36, 5146–5150. [Google Scholar] [CrossRef]
  21. Xie, Q.; Dai, Q.; Li, Y.; Jiang, A. Increasing the discriminatory power of DEA Using Shannon’s Entropy. Entropy 2014, 16, 1571–1585. [Google Scholar] [CrossRef] [Green Version]
  22. Park, S.Y.; Bera, A.K. Maximum entropy autoregressive conditional heteroskedasticity model. J. Econom. 2009, 150, 219–230. [Google Scholar] [CrossRef]
  23. Kagan, A.M.; Linik, Y.V.; Rao, C.R. Characterization Problems in Mathematical Statistics; Wiley: New York, NY, USA, 1973. [Google Scholar]
  24. Pendharkar, P.C. Cross efficiency evaluation of decision-making units using the maximum decisional efficiency principle. Comput. Ind. Eng. 2020, 145, 106550. [Google Scholar] [CrossRef]
  25. Troutt, M.D. Maximum decisional efficiency estimation principle. Manag. Sci. 1995, 41, 76–82. [Google Scholar] [CrossRef]
  26. Troutt, M.D. Derivation of the maximin efficiency ratio model from the maximum decisional efficiency principle. Ann. Oper. Res. 1997, 73, 323–338. [Google Scholar] [CrossRef]
  27. Troutt, M.D.; Zhang, A.; Tadisina, S.K.; Rai, A. Total factor efficiency/productivity ratio fitting as an alternative to regression and canonical correlation models for performance data. Ann. Oper. Res. 1997, 74, 289–304. [Google Scholar] [CrossRef]
  28. Banker, R.D. Maximum likelihood, consistency and data envelopment analysis: A statistical foundation. Manag. Sci. 1993, 39, 1265–1273. [Google Scholar] [CrossRef]
Figure 1. Exhaustive Search Tree for possible unique combinations of two-input-one-output datasets.
Figure 1. Exhaustive Search Tree for possible unique combinations of two-input-one-output datasets.
Algorithms 13 00232 g001
Figure 2. An illustration of 5 × 5 ensemble efficiency score matrix.
Figure 2. An illustration of 5 × 5 ensemble efficiency score matrix.
Algorithms 13 00232 g002
Figure 3. Number of times a DMU is fully efficient in Ensemble CCR models.
Figure 3. Number of times a DMU is fully efficient in Ensemble CCR models.
Algorithms 13 00232 g003
Figure 4. Number of times a DMU is fully efficient in Ensemble BCC models.
Figure 4. Number of times a DMU is fully efficient in Ensemble BCC models.
Algorithms 13 00232 g004
Figure 5. The maximum decisional efficiency (MDE) probability density function (pdf) for α = 5 and α = 10, respectively.
Figure 5. The maximum decisional efficiency (MDE) probability density function (pdf) for α = 5 and α = 10, respectively.
Algorithms 13 00232 g005
Table 1. Ensemble data envelopment analysis (DEA) scores for 1899–1910 US economic growth data.
Table 1. Ensemble data envelopment analysis (DEA) scores for 1899–1910 US economic growth data.
YearProductionLaborCapitalDEA Model EfficienciesEnsemble
Score
z = [111]z = [101]z = [011]
18991001001000.6810.6810.6650.676
19001011051070.7220.7220.6780.707
19011121101140.6930.6930.6890.692
19021221171220.6930.6810.6930.689
19031241221310.7200.7200.7140.718
19041221211380.7700.7700.7580.766
19051431251490.7930.7100.7930.765
19061521341630.8090.7300.8090.783
19071511401760.8360.7940.8360.822
19081261231851.0001.0001.0001.000
19091551431980.9210.8700.9210.904
19101591472080.9410.8910.9410.924
Table 2. The Iranian gas company dataset.
Table 2. The Iranian gas company dataset.
DMUx1x2x3y1y2y3y4
1177,430401528,32580141,67577,564201,529
2221,33810941,186,90580334,96044,136840,446
3267,80610791,323,32525124,46127,690832,616
4160,912444648,68581623,74445,882251,770
5177,214801909,53965436,40972,676443,507
6146,325686545,11517718,00019,839341,585
7195,138687790,34869531,22140,154233,822
8108,146152236,72260623,88937,770118,943
9165,663494523,89965225,16328,402179,315
10195,728503428,56695943,44063,701195,303
1187,050343298,696221968917,334106,037
12124,313129198,59856521,03230,24261,836
1367,545117131,64915210,39814,13946,233
1447,208165228,730211939113,50542,094
Table 3. The results of experiments.
Table 3. The results of experiments.
DMUEnsemble CCREnsemble BCCNon-Radial &Radial &
10.87 (0.15)0.95 (0.11)0.980.75
20.75 (0.30)0.77 (0.28)11
30.61 (0.36)0.62 (0.36)0.90.82
40.71 (0.19)0.8 (0.19)0.790.63
50.77 (0.22)0.82 (0.21)0.950.83
60.58 (0.27)0.64 (0.27)0.760.64
70.54 (0.16)0.57 (0.14)0.570.47
80.98 (0.08)0.99 (0.04)11
90.57 (0.14)0.6 (0.14)0.610.46
100.86 (0.18)0.96 (0.11)0.850.77
110.47 (0.12)0.63 (0.14)0.550.46
120.93 (0.15)0.96 (0.11)11
130.63 (0.13)0.96 (0.09)0.680.51
140.6 (0.24)0.86 (0.17)0.560.51
& Results taken from Toloo and Babaee’s [16] study.
Table 4. The results of normality tests.
Table 4. The results of normality tests.
Kolmogorov–SmirnovShapiro–Wilk
StatisticdfSig.StatisticdfSig.
Ensemble BCC0.196140.1490.876140.051
Ensemble CCR0.182140.2000.944140.477

Share and Cite

MDPI and ACS Style

Pendharkar, P.C. A Comparison of Ensemble and Dimensionality Reduction DEA Models Based on Entropy Criterion. Algorithms 2020, 13, 232. https://0-doi-org.brum.beds.ac.uk/10.3390/a13090232

AMA Style

Pendharkar PC. A Comparison of Ensemble and Dimensionality Reduction DEA Models Based on Entropy Criterion. Algorithms. 2020; 13(9):232. https://0-doi-org.brum.beds.ac.uk/10.3390/a13090232

Chicago/Turabian Style

Pendharkar, Parag C. 2020. "A Comparison of Ensemble and Dimensionality Reduction DEA Models Based on Entropy Criterion" Algorithms 13, no. 9: 232. https://0-doi-org.brum.beds.ac.uk/10.3390/a13090232

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