ABSTRACT
We present a novel approach to the inverse data envelopment analysis (DEA) process, based on the slacks-based measure (SBM) model in the presence of undesirable outputs. In the first step, we determine the minimum undesirable outputs levels that a decision-making unit (DMU) can produce, according to the level of inputs and desired outputs, based on weak disposability. We then introduce counterpart (hypothetical) units corresponding to each DMU. In the second step, we determine the optimal level of inputs and outputs for the newly created units, based on the SBM model in the presence of undesirable outputs. We present two new criteria models based on the SBM model, to compare the efficiency of the newly created unit in the inverse DEA process with the original unit. We demonstrate that the efficiency scores obtained from these two models are equal. Using the SBM-based inverse non-radial DEA process, we can use all inefficiency slacks related to all input and undesirable output components in the inverse non-radial DEA process to estimate inputs and outputs simultaneously.
Keywords:
Data Envelopment Analysis; inverse DEA; undesirable output; efficiency; non-radial DEA
1 INTRODUCTION
According to the World Environment Organization, adverse waste products and social activities, such as air pollutants and hazardous waste, are increasingly being produced in the environment and pose a serious threat to it. For example, the number of defective goods, amount of pollution and waste, or release of CO2 in the production process are all undesirable and should be reduced.
Therefore, evaluating the performance of organizations that produce undesirable outputs without considering these outputs is not a logical or correct evaluation, and these outputs should be included in the performance evaluation (Scheel, 2001; Sueyoshi and Goto, 2010SUEYOSHI T & GOTO M. 2010. Should the US clean air act include CO2 emission control? Examination by data envelopment analysis. Energy Policy, 38(10): 5902-5911., 2012SUEYOSHI T & GOTO M. 2012. DEA radial and non-radial models for unified efficiency under natural and managerial disposability: Theoretical extension by strong complementary slackness conditions. Energy Economics, 34: 700-713.). In this regard, the DEA technique, which was first introduced by Charnes et al. (1978CHARNES A, COOPER W & RHODES E. 1978. Measuring the efficiency of decision making units. European Journal of Operational Research, 2(6): 429-444.), is a suitable technique. When evaluating the performance of a DMU based on DEA models, we can use input or output-oriented models. In the input-oriented models, by decreasing the inputs (or increasing the outputs), we depict the inefficient unit on the efficiency frontier. In these models, the amount of outputs (or inputs) does not change. However, if we have an undesirable output among the output components, increasing these outputs will reduce the efficiency of the evaluated DMU. The DEA models are important for dealing with undesirable outputs. The DEA technique uses the Pareto optimality property in measuring efficiency. The efficient unit is denoted as non-dominated. In this method, if the DMU can maximize the production of one desirable output, minimize one of the undesirable outputs, or minimize one of the inputs, the DMUs can be regarded as efficient without taking into account the number of other parameters linked with it. Only if a DMU does not produce any undesirable outputs can it be considered efficient. Any amount of undesirable output generated should decrease the efficiency measures. However, since the production of desirable outputs normally comes with the generation of undesirable outputs, a certain quantity of undesirable outputs should be permitted because they will be generated until the desirable outputs are achieved. In this method, if a DMU generates a moderate amount of undesirable output, it still has a chance to be judged as efficient. The DMU is defined as inefficient if an excessive level of undesirable outputs is produced (Liu et al., 2010LIU W, MENG W, LI X & ZHANG D. 2010. DEA models with undesirable inputs and outputs. Annals of Operations Research, 173: 177-194.; Kao and Hwang, 2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.).
DEA models with undesirable inputs and outputs have been frequently discussed in DEA literature. For instance, Shi et al. (2010SHI G, BI J & WANG J. 2010. Chinese regional industrial energy efficiency evaluation based on a DEA model of fixing non-energy inputs. Energy Policy, 38(10): 6172-6179.) considered undesirable outputs for assessing energy efficiency in the Chinese manufacturing industry, and Yeh et al. (2010YEH T, CHEN T & LAI P. 2010. A comparative study of energy utilization efficiency between Taiwan and China. Energy Policy, 38(5): 2386-2394.) compared the total factors of energy efficiency in China. Sueyoshi and Goto (2010SUEYOSHI T & GOTO M. 2010. Should the US clean air act include CO2 emission control? Examination by data envelopment analysis. Energy Policy, 38(10): 5902-5911.) proposed a new DEA model to calculate the efficiency of electric fossil fuels, taking into account the CO2 produced by each production unit. Chen et al. (2017CHEN Z, WANKE P, ANTUNES J & ZHANG N. 2017. Chinese airline efficiency under CO2 emissions and flight delays: A stochastic network DEA model. Energy Economics, 68: 89-108.) proposed a stochastic network DEA model for evaluating Chinese airline efficiency under CO2 emissions and flight delays. Izadikhah and Saen (2018IZADIKHAH M & SAEN R. 2018. Assessing sustainability of supply chains by chance-constrained two-stage DEA model in the presence of undesirable factors. Computers & Operations Research, 100: 343-367.) developed a chance-constrained two-stage DEA model to present undesirable factors for evaluating the sustainability of supply chains. Ren et al. (2020REN J, GAO B, ZHANG J & CHEN C. 2020. Measuring the energy and carbon emission efficiency of regional transportation systems in China: Chance-constrained DEA models. Mathematical Problems in Engineering, Available at: https://doi.org/10.1155/2020/9740704.
https://doi.org/10.1155/2020/9740704...
) proposed a new chance-constrained DEA model to measure the energy and carbon emission efficiency of regional transportation systems in China.
One of the convenient methods proposed for dealing with undesirable outputs is the approach introduced by Kao and Hwang (2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.). They provided a suitable model for measuring the effect of undesirable outputs on efficiency measurement by addressing the minimum level of such outputs that can be allowed for each DMU. They used the results to create a new efficiency frontier of production by introducing the concept of a counterpart (hypothetical) unit corresponding to each of the DMUs. The efficiency of a production unit was presented based on the minimum level of the undesirable output, and the new frontier was used to calculate the efficiency scores of a manufacturing unit and compare with other units with the minimum level of undesirable outputs. (For more detail see: Anderson, Daim, and Kim, 2008ANDERSON T, DAIM T & KIM J. 2008. Technology forecasting for wireless communication. Technovation, 28: 602-614.; Inman, Anderson, and Harmon, 2006INMAN O, ANDERSON T & HARMON R. 2006. Predicting US jet fighter aircraft introductions from 1944 to 1982: A dogfight between regression and TFDEA. Technological Forecasting and Social Change, 73: 1178-1187.; Lim, Anderson, and Shott, 2015LIM DJ, ANDERSON T & SHOTT T. 2015. Technological forecasting of supercomputer development: The march to exascale computing. Omega, 51: 128-135.). The difference between these two efficiencies was then used to calculate the impact of creating too many undesirable outputs on a specific manufacturing unit’s efficiency. Data transformation technologies can produce skewed efficiency measurements (Färe and Grosskopf, 2004FÄRE R, GROSSKOPF S, LOVELL C & PASURKA C. 1989. Multilateral productivity comparisons when some outputs are undesirable: A nonparametric approach. The Review of Economics and Statistics, 71: 90-98., & Kao, 2020KAO C. 2020. Measuring efficiency in a general production possibility set allowing for negative data. European Journal Operational Research, 282(3): 980-988.). The approaches that use the source data to determine efficiency are provided by Kao and Hwang (2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.) for this purpose. The SBM model, one of Kao’s (2017KAO C. 2017. Network data envelopment analysis: Foundations and extensions. Switzerland: Springer International Publishing.) three types of direct approaches, can distinguish between the observed amount and the optimal one, making it suitable for analyzing the effect of creating a large level of undesirable outputs in measuring efficiency. Therefore, the SBM idea is employed in this study to create models in inverse DEA in the presence of undesirable outputs.
DEA models can be divided into radial and non-radial categories for evaluating efficiency (see: Gerami et al. 2022GERAMI J, MOZAFFARI M, WANKE P & CORREA H. 2022. Improving information reliability of non-radial value efficiency analysis: an additive slacks based measure approach. European Journal of Operational Research, 298(3): 967-978.). Radial models do not consider slacks in inputs and outputs, which are important components for evaluating efficiency. Radial approaches may not produce robust results in efficiency evaluation based on radial reduction of inputs and outputs. All inverse DEA models are radial models, and as a result, they may help DMs estimate the level of outputs (inputs) by ignoring the effects of slack variables corresponding to input and output components. However, radial models suffer from the limitation of being incapable of addressing the slacks in the model. Limited studies have used non-radial inverse DEA models. For instance, Jahanshahloo et al. (2014JAHANSHAHLOO G, HOSSEINZADEH LOTFI F, ROSTAMY-MALKHALIFEH M & GHOBADI S. 2014. Using enhanced Russell model to solve inverse data envelopment analysis problems. The Scientific World Journal, p. 1-10.) proposed an inverse DEA model based on the Enhanced Russell Model, assuming constant efficiency scores among dimensions. Zhang and Cui (2020ZHANG G & CUI G. 2020. A general inverse DEA model for non-radial DEA. Computers & Industrial Engineering, 142: 106368.) proposed a general inverse DEA model for non-radial DEA, assuming unchanged overall efficiency scores but not for every dimension. They presented a basic form of non-radial inverse DEA models, which was nonlinear in a special case for the SBM model. They used nonlinear optimization methods to solve the proposed model and one-dimensional search as a suitable method to solve it.
We have developed a non-radial SBM model for estimating optimal inputs and outputs in the presence of undesirable outputs. The efficiency scores of each DMU are first obtained using the non-radial SBM model, and then a non-radial SBM model is presented under the inverse DEA process for simultaneous estimation of optimal inputs and outputs. The proposed approach can deal with both non-radial measurements and undesirable outputs. In the first step, the optimal level of undesirable outputs is obtained for each DMU using the approach proposed by Kao and Hwang (2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.). In the second step, the inverse DEA process is performed using the non-radial SBM model to obtain the optimal levels of inputs and outputs for the created new units that have the same efficiency scores as their corresponding original units. We present two new criteria models for evaluating the efficiency of the created new units in the presence of counterpart units and show that the efficiency of the created new units in the inverse DEA process is equal to that of the original units.
2 LITERATURE REVIEW
2.1 The undesirable outputs in DEA
When additional input and output parameters are included in the evaluation, the DEA approach can yield higher efficiency scores for the DMUs being evaluated. This can lead to unreasonable results. Moreover, when considering undesirable outputs in the efficiency measurement of a DMU, the scores will either be equal to or higher than those that do not consider the undesirable outputs. To arrive at a more robust efficiency metric, we should develop traditional DEA models and provide suitable models. In many manufacturing processes, undesirable outputs may be generated unexpectedly along with the planned outputs. For instance, in the pulping process, wastewater or sewage is created, and in the combustion of fossil fuels, CO2 and NOx are formed. Although these undesirable outputs may not have a direct impact on the DMU generating them, they have a negative impact on the environment. The production of such outputs should be minimized to protect the environment. However, due to the costs of preventing the production of undesirable outputs, they are usually kept at the minimum standard level allowed by law. If the costs of preventing their production outweigh the costs of breaking the rules, the production units may be forced to control the creation of undesirable outputs. Undesirable outputs are often ignored in efficiency measurements, but understanding their impact on efficiency measurements is crucial to increase environmental awareness. There are numerous ways to determine the efficiency of a production unit, with the DEA framework being the most commonly used. This technique was first developed by Charnes et al. (1978CHARNES A, COOPER W & RHODES E. 1978. Measuring the efficiency of decision making units. European Journal of Operational Research, 2(6): 429-444.), who created efficiency evaluation models based on the constant returns to scale (CRS) technology. Banker et al. (1984) further developed this method for the variable returns to scale (VRS) technology. DEA measures the relative efficiency of a set of DMUs that use multiple inputs to generate multiple outputs.
Scheel (2001) categorized methods of measuring the efficiency of DMUs with the presence of undesirable outputs into two categories: direct and indirect. Direct methods use original data, while indirect methods involve data transformations. The concept of weak disposability is directly linked to direct approaches (Färe et al., 1989FÄRE R, GROSSKOPF S, LOVELL C & PASURKA C. 1989. Multilateral productivity comparisons when some outputs are undesirable: A nonparametric approach. The Review of Economics and Statistics, 71: 90-98.; Liu et al., 2010LIU W, MENG W, LI X & ZHANG D. 2010. DEA models with undesirable inputs and outputs. Annals of Operations Research, 173: 177-194.). Inverse input methods, additive inverse, translated inverse, and multiplicative inverse methods are examples of indirect methods. The direct approaches are classified by Zhao et al. (2008ZHOU P, ANG B & POH K. 2008. A survey of data envelopment analysis in energy and environmental studies. European Journal of Operational Research, 189: 1-18.) into the directional distance function (DDF), hyperbolic models (Kao, 2017KAO C. 2017. Network data envelopment analysis: Foundations and extensions. Switzerland: Springer International Publishing.), and the slacks-based measure (SBM). Vazquez-Rowe et al. (2010VAZQUEZ-ROWE I, IRIBARREN D, MOREIRA M & FEIJOO G. 2010. Combined application of life cycle assessment and data envelopment analysis as a methodological approach for the assessment of fisheries. The International Journal of Life Cycle Assessment, 15: 272-283.) used the SBM-DEA model in the presence of undesirable outputs to assess fisheries in Spain.
Undesirable outputs are considered in the evaluation of environmental efficiency, and Song et al. (2012SONG M, AN Q, ZHANG W, WANG Z & WU J. 2012. Environmental efficiency evaluation based on data envelopment analysis: A review. Renewable and Sustainable Energy Reviews, 16(7): 4465-4469.) classified the methods into three categories including disposability-related methods, data transformation, and input reversal. In an extensive review of papers on performance benchmarking in the presence of undesirable outputs, five approaches are classified by Dakpo et al. (2016DAKPO K, JEANNEAUX P & LATRUFFE L. 2016. Modelling pollution-generating technologies in performance benchmarking: Recent developments, limits and future prospects in the nonparametric framework. European Journal of Operational Research, 250: 347-359.), including data transformation, materials balance principles, free disposability of the inputs, weak disposability of the undesirable outputs, as well as the two sub-technologies from the biproportional production model of Murty et al. (2012MURTY S, RUSSELL R & LEVKOFF S. 2012. On modeling pollution-generating technologies. Journal of Environmental Economics and Management, 64: 117-135.), and the natural mammalian performance benchmarking model. Wu et al. (2013WU C, LI Y, LIU Q & WANG K. 2013. A stochastic DEA model considering undesirable outputs with weak disposability. Mathematical and Computer Modelling, 58(5-6): 980-989.) considered wasting water, emission of toxic gases, and the production of useless solid material as undesirable stochastic outputs in evaluating several provinces in China.
Dakpo et al. (2016DAKPO K, JEANNEAUX P & LATRUFFE L. 2016. Modelling pollution-generating technologies in performance benchmarking: Recent developments, limits and future prospects in the nonparametric framework. European Journal of Operational Research, 250: 347-359.) proposed the drawbacks of each method. They evaluated the efficiency of 136 winter wheat farms in Poland under six DEA models with the life cycle assessment approach. The methods included the ones that disregard the undesirable outputs, treat the undesirable outputs as inputs, impact rate, data transformation, You and Yan’s (2011YOU S & YAN H. 2011. A new approach in modelling undesirable output in DEA model. Journal of the Operational Research Society, 62(12): 2146-2156.) ratio model, as well as the slacks-based model. The comparison among the models showed that the slacks-based model provided better results in differentiating the performance among the firms. However, they did not discuss weak disposability. Pishgar-Komleh et al. (2020PISHGAR-KOMLEH S, ZYLOWSKI T, ROZAKIS S & KOZYRA J. 2020. Efficiency under different methods for incorporating undesirable outputs in an LCA+DEA framework: A case study of winter wheat production. Journal of Environmental Management, 260, Article 110138.) classified DEA models in the face of undesirable outputs as follows: a) ignoring undesirable outputs, b) treating undesirables as inputs to the DEA model, c) data transformation, d) impact rate, e) ratio model. They discussed several papers that used each method.
Cecchini et al. (2018CECCHINI L, VENANZI S, PIERRI A & CHIORRI M. 2018. Environmental efficiency analysis and estimation of CO2 abatement costs in dairy cattle farms in Umbria (Italy): a SBM-DEA model with undesirable output. Journal of Cleaner Production, 197: 895-907.) proposed a slack-based measure (SBM)-DEA model for analyzing the environmental efficiency of dairy cattle farms with undesirable outputs in Italy. Adenuga et al. (2018ADENUGA A, DAVIS J, HUTCHINSON G, DONNELLAN T & PATTON M. 2018. Modelling regional environmental efficiency differentials of dairy farms on the island of Ireland. Ecological Indicators, 95: 851-861.) used the directional output distance function DEA-based model to optimize both desirable and undesirable outputs on dairy farms in Ireland. Dong et al. (2018DONG G, WANG Z & MAO X. 2018. Production efficiency and GHG emissions reduction potential evaluation in the crop production system based on emergy synthesis and nonseparable undesirable output DEA: A case study in Zhejiang Province, China. PLOS ONE, 13(11): 0206680.) considered the SBM-DEA model in the presence of undesirable outputs to evaluate the resource use in crop production in a province in China. Angulo-Meza et al. (2019ANGULO-MEZA L, GONZALEZ-ARAYA M, IRIARTE A, REBOLLEDO-LEIVA R & MELLO J. 2019. A multiobjective DEA model to assess the eco-efficiency of agricultural practices within the CF+DEA method. Computers and Electronics in Agriculture, 161: 151-161.) proposed a multi-objective DEA model to reduce the CF of organic blueberry orchards while maintaining high production levels.
Pishgar-Komleh et al. (2020PISHGAR-KOMLEH S, ZYLOWSKI T, ROZAKIS S & KOZYRA J. 2020. Efficiency under different methods for incorporating undesirable outputs in an LCA+DEA framework: A case study of winter wheat production. Journal of Environmental Management, 260, Article 110138.) defined efficiency under different methods for incorporating undesirable outputs in an LCA+DEA framework. They applied the proposed approach to winter wheat production in Poland. Kao and Hwang (2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.) proposed an approach for measuring the effects of undesirable outputs on the efficiency of production units. They developed a concept for determining the minimum amount of undesirable outputs that a DMU is allowed to generate based on the assertion of weak disposability and proposed a suitable approach for measuring the effect of undesirable outputs on the efficiency measurement.
2.2 Inverse DEA
The inverse DEA models, compared to the traditional DEA models, achieve the rate of output variations by increasing inputs while maintaining the efficiency level of the unit under evaluation. The inverse optimization problem is useful for the decision-maker (DM) to obtain information regarding resource allocation. The traditional DEA models focus only on efficiency evaluation. The idea of the inverse DEA was first introduced by Zhang and Cui (1999ZHANG XS & CUI JC. 1999. A project evaluation system in the state economic information system of China an operations research practice in public sectors. International Transactions in Operational Research, 6(5): 441-452.) and then formulated by Wei et al. (2000WEI Q, ZHANG J & ZHANG X. 2000. An inverse DEA model for inputs/outputs estimate. European Journal of Operational Research, 121(1): 151-163.) based on multi-objective programming techniques. The inverse DEA models were then studied from various theoretical and practical perspectives. Hadi-Vencheh et al. (2008HADI-VENCHEH A, FOROUGHI A & SOLEIMANI-DAMANEH M. 2008. A DEA model for resource allocation. Economic Modelling, 25(5): 983-993.) proposed an inverse DEA model for estimating inputs. They initially increased the outputs and assumed that the efficiency of the unit under evaluation would not change. Zhang and Cui (2016ZHANG M & CUI JC. 2016. The extension and integration of the inverse DEA method. Journal of the Operational Research Society, 67: 1212-1220.) integrated twelve different scenarios into the inverse DEA models. Other applications of the inverse DEA include input estimation for resource allocation (Xiaoya & Jinchuan, 2008XIAOYA L & JINCHUAN C. 2008. A comprehensive DEA approach for the resource allocation problem based on scale economies classification. Journal of Systems Science and Complexity, 21(4): 540-557.) and inter-temporal application (Jahanshahloo et al., 2015JAHANSHAHLOO G, SOLEIMANI-DAMANEH M & GHOBADI S. 2015. Inverse DEA under inter-temporal dependence using multiple-objective programming. European Journal of Operational Research, 240(2): 447-456.). Ghobadi (2018GHOBADI S. 2018. Inverse DEA using enhanced Russell measure in the presence of fuzzy data. International Journal of Industrial Mathematics, 10(2): 165-180.) presented inverse DEA models with fuzzy data. Ghiyasi (2017GHIYASI M. 2017. Inverse DEA based on cost and revenue efficiency. Computers & Industrial Engineering, 114: 258-263.) developed inverse DEA models based on the concepts of cost and revenue efficiency. Other applications of the inverse DEA models include recent studies on sustainability (Hassanzadeh et al., 2018HASSANZADEH A, YOUSEFI S, SAEN R & HOSSEININIA S. 2018. How to assess sustainability of countries via inverse data envelopment analysis? Clean Technologies and Environmental Policy, 20: 29-40.), enterprise merger (Amin et al., 2017aAMIN G, EMROUZNEJAD A & GATTOUF S. 2017a. Minor and major consolidations in inverse DEA: Definition and determination. Computers & Industrial Engineering, 103: 193-200., bAMIN G, EMROUZNEJAD A & GATTOUF S. 2017b. Modelling generalized firms’ restructuring using inverse DEA. Journal of Productivity Analysis, 48: 51-61.), revenue target management (Lin, 2010LIN H. 2010. An efficiency-driven approach for setting revenue target. Decision Support Systems, 49(3): 311-317.), and pricing strategy (Frija et al., 2011FRIJA A, WOSSINK A, BUYSSE J, SPEELMAN S & HUYLENBROECK G. 2011. Irrigation pricing policies and its impact on agricultural inputs demand in Tunisia: A DEA-based methodology. Journal of Environmental Management, 92(9): 2109-2118.). Lertworasirikul et al. (2011LERTWORASIRIKUL S, P C & FANG S. 2011. Inverse data envelopment analysis model to preserve relative efficiency values: The case of variable returns to scale. Computers and Industrial Engineering, 61(4): 1017-1023.) proposed inverse DEA models in the VRS technology. Amin et al. (2017bAMIN G, EMROUZNEJAD A & GATTOUF S. 2017b. Modelling generalized firms’ restructuring using inverse DEA. Journal of Productivity Analysis, 48: 51-61.) developed a generalized inverse DEA model for restructuring, considering the number of original DMUs as pre-restrictions for producing new units. Lim (2016LIM D. 2016. Inverse DEA with frontier changes for new product target setting. European Journal of Operational Research, 254(2): 510-516.) proposed an inverse DEA model based on the frontier changes for new product target setting. He discusses an inverse DEA problem that considered expected changes in the production frontier in the future by adding the inverse optimization problem with a time series application of DEA to provide a tool for new product target setting practices. Wegener and Amin (2019WEGENER M & AMIN G. 2019. Minimizing greenhouse gas emissions using inverse DEA with an application in oil and gas. Expert Systems with Applications, 122: 369-375.) applied an inverse DEA model in the oil and gas industry. Emrouznejad et al. (2019EMROUZNEJAD A, YANG G & AMIN G. 2019. A novel inverse DEA model with application to allocate the CO2 emissions quota to different regions in Chinese manufacturing industries. Journal of the Operational Research Society, 70(7): 1079-1090.) proposed an inverse DEA model to allocate undesirable outputs. Ghiyasi (2019GHIYASI M. 2019. Novel criterion models in the inverse DEA problem. International Journal of Operational Research, 35(1): 20-36.) proposed a new criterion model, reducing the computational complications in the inverse DEA. Ghiyasi and Zhu (2020GHIYASI M & ZHU N. 2020. An inverse semi-oriented radial data envelopment analysis measure for dealing with negative data. IMA Journal of Management Mathematics, 31(4): 505-516.) proposed an inverse semi-oriented radial DEA in the presence of negative data. Amin and Ibn Boamah (2021AMIN G & IBNBOAMAH M. 2021. A two-stage inverse data envelopment analysis approach for estimating potential merger gains in the US banking sector. Managerial and Decision on Economics, 42(6): 1454-1465.) developed a two-stage inverse DEA approach for estimating potential merger gains in the US banking sector. Zeinodin and Ghobadi (2020ZEINODIN E & GHOBADI S. 2020. Merging decision-making units under inter-temporal dependence. IMA Journal of Management Mathematics, 31(2): 139-166.) examined the merger issue under the inter-temporal dependence structure. Gerami et al. (2023GERAMI J, MOZAFFARI M, WANKE P & HL C. 2023. Generalized inverse DEA model for firm restructuring based on value efficiency. IMA Journal of Management Mathematics, 34(3): 541-580.) proposed a generalized inverse DEA model for firm restructuring based on value efficiency.
3 MEASURING THE EFFICIENCY IN THE PRESENCE OF UNDESIRABLE OUTPUTS BASED ON THE SBM NON-RADIAL MODEL
Suppose we have n decision units as DMU j =(x j , y j ), j=1, ..., n. The input and output vectors corresponding to DMU j , j=1, ..., n are denoted as x j =x 1j , ..., x mj and y j =y 1j , ..., y sj .
Assume that , i=1, ..., m and , r=1, ..., s are the slacks of the inputs and outputs, respectively. j is the intensity variable corresponding to DMU j .
The SBM model for measuring the efficiency of DMU k as a unit under evaluation in the variable returns to scale technology was presented by Tone (2001TONE K. 2001. A slacks-based measure of efficiency in data envelopment analysis. European Journal of Operational Research, 130(3): 498-509.) as follows:
Suppose, (λ*, s x *, s y *) where and is an optimal solution of model (1).
Definition 1.DMUkis efficient in evaluation withmodel (1)if and only if the objective function score ofmodel (1)is equal to one or equivalent to, i=1, ..., m, , r=1, ..., s.
When there are undesirable outputs, such as uj=u1j , ..., u fj , j=1, ..., n,model (1)cannot be used to assess the DMUs’ efficiency. In this regard,Kao and Hwang (2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.) proposed a new approach for efficiency evaluation in the presence of undesirable outputs. They first introduced the concept of a counterpart (hypothetic) unit and presented the efficiency based on the minimum level of undesirable output.
To obtain the amount of efficiency corresponding to each of the DMUs on the basis of the projection of the unit under evaluation on the new frontier created based on the counterpart units (for more detail see:Entani et al. 2002ENTANI T, MAEDA Y & TANAKA H. 2002. Dual Models of Interval DEA and its extensions to interval data. European Journal of Operational Research, 136(1): 32-45.), they used this new frontier to calculate the efficiency scores that produce the least undesirable output. The difference between these two efficiencies is then used to calculate the impact of creating too many undesirable outputs on a production unit’s efficiency. The counterpart units have the same input and output levels as the original units, but their undesirable output level is lower than that of the original units. To determine the level of undesirable outputs from these DMUs, assume that, f=1, ..., h are slacks of the undesirable outputs, we can solve the following model (Kao and Hwang, 2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.):
Suppose (λ*, s u *) where is an optimal solution of model (2) . In , the level of desirable outputs remains constant, and DMU k is allowed to generate the lowest undesirable outputs as below:
On the basis of the frontier constructed from the counterpart DMUs, namely, j=1, ..., n, the efficiencies can be calculated through the following model (Kao and Hwang, 2021KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.):
Inmodel (3),, i=1, ..., m,, r=1, ..., s, and, f=1, ..., h, are slacks of the inputs, desirable and undesirable outputs, respectively. λj is the intensity variable corresponding to DMU j . Inmodel (3), the undesirable outputs are reduced relative to the desirable output under weak disposability. Suppose, (λ*, s x *, s u *, s y *) where,andis an optimal solution ofmodel (3).
Definition 2.DMUkis called efficient in evaluation withmodel (3)if and only if the objective function score ofmodel (3)is equal to one or equivalent, i=1, ..., m, , r=1, ..., s,, f=1, ..., h.
Now, based on model (3) , we present the input-oriented non-radial SBM model as follows:
Inmodel (3),, i=1, ..., m, and, f=1, ..., h, represent the slacks of the inputs and undesirable outputs, respectively. λj is the intensity variable corresponding to DMU j .Model (4)obtains the efficiency scores of the DMU under evaluation. Suppose, (λ*, s x *, s u *) is an optimal solution ofmodel (4)whereand.
Definition 3.DMUkis efficient in evaluation withmodel (4)if and only if the objective function score ofmodel (4)is equal to one, i.e., ψk =1 or equivalently,, i=1, ..., m and, f=1, ..., h.
4 SBM-BASED INVERSE NON-RADIAL DEA
We present a new approach based on the non-radial SBM model for dealing with the presence of undesirable outputs. We use the input-oriented SBM model, which calculates the efficiency of DMUs based only on the slack values of the inputs and undesirable outputs. Therefore, we first solve model (4) to obtain the efficiency and slacks of the inputs and undesirable outputs. Suppose (λ*, s x *, s u *) is an optimal solution of model (4) where and , To estimate the input, desirable output, and output components in the inverse DEA process, we created a new corresponding DMU k as follows, taking into account the level of inputs and outputs:
where ∆x k , ∆y k , ∆u k , show the deviations of the inputs, desirable and undesirable outputs of DMU k , respectively. In the inverse DEA process based on the SBM model, the optimal inputs and outputs of the newly created DMU, i.e. , are determined in such a way that this unit’s efficiency is the same as that of the corresponding unit, i.e. DMU k . To obtain the optimal inputs and outputs of the newly created unit, we solve the following model in the inverse non-radial DEA process:
In model (5), , i=1, ..., m, and , f=1, ..., h, are the slacks of the inputs and undesirable outputs, respectively. λj is the intensity variable corresponding to DMU j , and is the intensity variable corresponding to . Additionally, ∆x ik , i=1, ..., m, ∆y rk , r=1, ..., s, ∆u fk , f=1, ..., h, represent the deviations of the inputs, desirable and undesirable outputs of DMU k , respectively. The constraint x ik +∆x ik ≥0, u fk +∆u fk ≥0, i=1, ..., m, f=1, ..., h, y rk +∆y rk ≥0, r=1, ..., s, ensures that the optimal level the newly created unit’s inputs and outputs, i.e. , is non-negative.
If we present model (5) with the slack values corresponding to the desirable output components, the linearization process of is not easy. Therefore, we only present the input-oriented model (5) with the slack values corresponding to the input and undesirable output components. Model (5) is a nonlinear model used to determine the efficiency of the created new unit, i.e., . However, the optimal level of the new unit’s inputs and outputs, i.e., , is unclear, making it difficult to solve model (5) and determine its efficiency. Additionally, the optimal solution of this model may be zero. To determine the efficiency score of the created new unit, i.e., , we need to obtain the optimal level of inputs and outputs corresponding to it. However, due to the non-linear shape of model (5), we need to present a new linear model. This model assumes that the created new unit’s efficiency, i.e., , is the same as that of the corresponding original unit, i.e., DMU k , and that we obtain the efficiency score corresponding to DMU k from model (4). To do this, we need to present the model in a way that includes the optimal values of the slack variables regarding the input and undesirable output components resulting from model (4) corresponding to the unit under evaluation, i.e., DMU k . This ensures that the optimal values of the slack variables regarding the input and undesirable output components corresponding to the unit under evaluation, i.e., DMU k , from model (4) are equal to the ones of the created new unit, i.e., , from model (5) in a special case:
In this case, we can guarantee that the created new unit’s efficiency, i.e. is the same as the one of the corresponding original unit, i.e. DMU k , and with this assumption, the created new unit’s optimal inputs and outputs are obtained. Equation (6) is equivalent to the following equation:
We can now use the slacks values of the input and undesirable output components corresponding to the newly created unit, i.e., , in the SBM-based inverse non-radial DEA process to determine the optimal level of this unit’s inputs and outputs. We propose model (8) as follows:
Model (8) obtains the minimum level of inputs and undesirable outputs from the created new unit, i.e. , in a manner that the created new unit’s efficiency, i.e. , is equal to that of the corresponding unit, i.e. DMU k . In this model, we use the optimal values of the slack variables regarding the input and undesirable output components corresponding to the created new unit., i.e. from model (5), which we obtain in equation (7), according to the optimal values of the slack variables regarding the input and undesirable output components corresponding to the unit under evaluation, i.e. DMU k from model (4). We apply model (8) to determine the minimum level of inputs and undesirable outputs from the created new unit, i.e., . In model (8), the optimal level of inputs, desirable outputs, and undesirable outputs of the created new unit, i.e. , is determined based on the counterpart DMUs, namely, , j=1, ..., n.
Suppose such that , and is an optimal solution of model (8). In this case, the optimal levels of the created new unit’s inputs and outputs corresponding to i.e. DMUk based on model (8) are considered as follows:
show the optimal deviations of the inputs, desired outputs, and undesirable outputs of DMU k , respectively. In the following, to show that the created new unit’s efficiency, i.e. , in the inverse non-radial DEA process is the same as the one of the original unit corresponding to it, first we provide the following criterion model:
Model (9) obtains the created new unit’s efficiency score, i.e. , based on the slack variables of the input and undesirable output components in the structure of the SBM model. In Theorem (1), we show that the score of the optimal objective function obtained from model (9) in the evaluation and the optimal score obtained from model (4) in the evaluation DMU k are equal. In other words, the efficiency scores of the unit under evaluation, i.e. DMU k , and its corresponding unit in the inverse non-radial DEA process, i.e. , are equal.
Theorem 1. The optimal objective function scores of models (9) is less than or equal to the optimal objective function scores of models (4).
Proof. As stated earlier, to solve model (9), we need to first solve model (8) and use the optimal solution obtained from model (8) to solve model (9). Similarly, to solve model (8), we need to solve model (4) and use the optimal solution obtained from model (4) to solve model (8). Let (ψk , λ*, s x *, s u *) be a desired optimal solution of model (4). Based on this solution, we solve model (8) and let be the desired optimal solution of model (8). Next, based on this solution, we solve model (9) and let be the desired optimal solution of model (9).
To prove the theorem, we need to show that . Given that is a desired optimal solution of model (8), it satisfies in constraints of model (8). Therefore, we have
(λ*, s x *, s u *) is a desired optimal solution of model (4) which was used to solve model (8). Now we get a feasible solution for model (9). For this purpose, we put
According to the relations (10) and the constraints of model (9), it is clear that satisfies the constraints of model (9), and therefore is a feasible solution for model (9). The objective function score of model (9) for this feasible solution is as follows:
According to relations (11), the optimal objective function scores of models (9) is less than or equal to the optimal objective function scores of models (4), this is attributed to the fact that the objective function score of model (9) for a feasible solution is equal to the optimal objective function score of model (4), and both models are of the minimization type. Then .
Conversely, suppose is a desired optimal solution of model (9). Then it satisfies the constraints of model (9).
To solve model (9), we must first solve model (8) and use the optimal solution obtained from model (8) to solve model (9). In model (8), we replace the variables , i=1, ..., m, , f=1, ..., h, with these scores, then we will have
Based on the relations above, we can derive a special case where we can establish a relation between the scores and , i=1, ..., m, as well as a relation between the scores , , f=1, ..., h, we put
Given that is an optimal solution of model (4) which we used to solve model (8). Then we have
Based on the relations above, we can conclude that the optimal objective function scores of models (4) is less than or equal to the optimal objective function scores of models (9) because the objective function score of the model (4) for a feasible solution is equal to the optimal objective function score of model (9), and both models are of the minimization type. Therefore, we have .
From these results, we can deduce that . The proof of the theorem is complete. □
We will now calculate the efficiency score of the newly created unit, denoted as , based on model (5). Assume that is an optimal solution of model (8). We put the following values in model (5):
Therefore, the new criteria model for evaluating the efficiency of based on model (5) can be presented as follows:
Theorem 2. Suppose that and are the optimal solutions of models (12) and (9), respectively. Then, we have .
Proof. Let be an optimal solution of model (12), and assume that is not efficient based on model (12). In this case, we will have .
We can present the optimal solution by considering the constrains of model (12) as follows:
By placing , i=1, ..., m, , f=1, ..., h, , j=1, ..., n,
By considering the constraints in model (9) and relations (13), it is possible to obtain a feasible solution for model (9). The objective function score of model (9) for this feasible solution will be equal to . From the optimization point of view, since model (9) is a minimization type, its optimal objective function score will be we will therefore have:
On the other hand, suppose that (γ*, s x-cr *, s u-cr *) is an optimal solution of model (9). Now consider the constraints of model (9) for this optimal solution as follows:
By placing
Given the constraints in model (12) and the relationships in (15), we can obtain a feasible solution for model (12). The objective function score of model (12) for this feasible solution will be equal to . From the optimization point of view, considering that model (12) is a model of minimization type, and its optimal objective function score is equal to , we will therefore have:
According to the relations (14) and (16), we have:
Therefore, the optimal objective function scores of models (9) and (12) are equal.
If this unit is efficient in evaluating based on model (12), the objective function score of model (12) will be equal to one. In this case, model (12) is equivalent to the input-oriented SBM model or model (4), which we use to evaluate n+1 DMUs including the counterpart DMUs, namely , j=1, ..., n, and . Therefore, the objective function score of model (4) in the evaluation of this n+1 DMU will be equal to one. According to Theorem (1), the optimal objective function score obtained from models (4) and (9) are equal. Therefore, the optimal objective function score of model (9) is also equal to one, meaning is also efficient in evaluation with model (9). Therefore, if is efficient in evaluation with model (12), it also will be an efficient unit in evaluation with model (9). The optimal objective function scores of models (9) and (12) are then equal, and the proof of the theorem is complete. □
If the efficiency frontier remains unchanged in the presence of newly created DMUs, or in other words, the new units can be expressed as a combination of the original observed DMUs, then in model (12). Model (12) calculates the efficiency score of the created new unit, namely , based on model (5). We can say that model (12) is a criteria model for evaluating the efficiency of based on model (5).
Now, according to Theorems (1) and (2), the optimal objective function scores of models (4) and (12) are equal. This is because, according to Theorem (1), the optimal objective function scores of models (4) and (9) are equal, and according to Theorem (2), the optimal objective function scores of models (9) and (12) are equal. Therefore, the optimal objective function scores of models (4) and (12) are also equal. This means we can use either of the criteria models (9) and (12) to show that the created new unit’s efficiency in the inverse non-radial DEA process is equal to that of the original unit. Given that model (9) has a smaller number of variables compared to model (12), model (9) is proposed as a suitable criteria model. We now present the following algorithm as a three-step approach to estimate the optimal level of the newly created unit’s inputs and outputs in the SBM-based inverse non-radial DEA process as follows (Table 1):
An algorithm three-step approach for SBM-based inverse non-radial DEA in presence of undesirable output.
5 NUMERICAL EXAMPLE
This example presented below is intended to help the reader better understand the proposed approach. We consider a case of 5 DMUs, each having one input, one desirable output and one undesirable output, as shown in Table 2. The input, desirable output, and undesirable output are indicated by the symbols X,Y, and U, respectively.
First, using model (1), we obtain efficiency score for each DMU, as reported in Table 2. It can be seen that all units are inefficient except for units A, C and E. The SBM model divides the DMUs into two categories: strong efficient and inefficient. Model (1) calculates the efficiency score of each DMU without considering the undesirable output. To calculate the efficiency of each DMU in the presence of undesirable outputs, we solve model (2) and obtain the lowest undesirable output based on the input and desirable output values. To do this, we solve model (2) and introduce the counterpart (hypothetic) units corresponding to each of the DMUs, and we obtain the true efficiency scores and slacks of the inputs and undesirable outputs by the non-radial projection of these units on the frontier of efficiency that is created on the basis of the counterpart (hypothetic) units. These units have the same levels of inputs and desirable outputs as the original units, but their undesirable output level is less than the original units. These units have a minimum level of undesirable outputs based on the levels of inputs and desirable outputs. The results further report the input, desirable output, and undesirable output values corresponding to each of the counterpart DMUs. Model (3) is solved to obtain the corresponding efficiency score for each of the DMUs. This model obtains the efficiency score based on the slack variable of the inputs and undesirable outputs based on the counterpart (hypothetic) units. According to this model, none of the DMUs are efficient. Unit C has the highest efficiency score among the DMUs.
To perform the SBM-based inverse non-radial DEA process, we solve model (4) to obtain the efficiency scores and slacks of the inputs and undesirable outputs. This model achieves the true efficiency scores and slacks of each DMU by representing the units in the non-radial form on the frontier of efficiency that is created based on the counterpart (hypothetic) units. The results are given in the second to fourth columns of Table 3.
The result of models (4), (9), (12) in the inverse non-radial DEA process for data set in numerical example.
We can see that the results of models (3) and (4) are the same, and none of the DMUs are efficient. It is shown in Tables (2) and (3) that the efficiency scores of each DMU based on model (1) are greater than the ones obtained from model (3); however, this is not true from an optimization perspective. This is because we expect the efficiency of each DMU to deteriorate and decrease in the presence of undesirable outputs. Therefore, it is recommended that instead of solving model (1), we use models (3) and (4) to calculate the efficiency of DMUs.
As stated in the third section, in the inverse non-radial DEA process, we must solve model (8), through which the lowest inputs and undesirable outputs can be determined. Before solving model (8), we must first solve model (4) and obtain the efficiency scores and optimal slacks of the inputs and undesirable outputs of the unit under evaluation. As previously stated, model (8) obtains the minimum level of inputs and undesirable outputs from the created new unit in such a way that the created new unit’s efficiency is the same as the original unit. The third to fifth columns of Table 4 show the deviations of the input, desirable outputs, and undesirable output components of the DMUs, respectively. Also, the sixth to eighth columns of Table 4 show the respective new levels.
To show that the efficiency scores corresponding to the original DMUs and the ones of the created new units are equal, we can solve two criterion models (9) and (12). We observe that the efficiency scores of the created new unit are the same as those of the original units. Model (9) with the SBM model structure obtains the efficiency scores of the created new unit, and we obtain their input and undesirable output levels from model (8). Additionally, model (9) non-radially depicts the created new units on the frontier of efficiency that is created on the basis of the counterpart units. Next, we use model (12) as a criterion model to obtain the efficiency scores of the created new unit. The level of inputs and outputs of this new unit was determined based on model (8). The results are presented in the last column of Table 4, showing that the efficiency scores of the created new units based on model (12) are the same as those of the corresponding original units.
6 CASE STUDY
In this section, we apply our method to the dataset of the energy sector in the Iranian oil industry. The energy sector in Iran is of particular importance, which is why it is necessary to have accurate, comprehensive, and reliable statistics about it. Energy statistics are important tools for making correct and appropriate decisions and policies in the energy sector. These statistics must have the necessary coordination with the decision-making process of energy policies in the country. In recent years, the planning office and macroeconomics of electricity and energy of the Ministry of Energy have tried to provide comprehensive information and publish a set of basic information chronologically from 2005 to meet the need of officials, managers, researchers, and experts of the country for reliable statistical sources and also to create appropriate, quick, and easy access to energy information. In this paper, we evaluate the performance of this industry over 14 years during the period 2005-2018 using the DEA technique.
To evaluate the performance of this industry, we consider each year as a Decision Making Unit (DMU). We must determine the inputs and outputs corresponding to each DMU for performance evaluation. The correct selection of inputs and outputs can determine the efficiency values accurately. After determining the efficiency of each of the DMUs, one of the managers’ questions is how to develop DEA models to obtain the amount of changes in the inputs and outputs from the DMUs that do not change the efficiency of these units. Therefore, inverse DEA models are suggested. It is important to use models that correctly calculate the efficiency score and consider all inefficiency factors in all input and output components. Therefore, we can develop inverse DEA models in terms of input and output slacks based on the SBM model. One of the important issues in the evaluation of energy sectors is the presence of unwanted data that has a negative effect on the performance of these sectors. According to the above explanations, the models presented in this paper are among the most appropriate models in dealing with unfavorable outcomes. These models use the inverse DEA technique to evaluate energy sectors during different years based on non-radial models.
We used the opinions of senior managers of this industry to select inputs and outputs. Inputs and outputs were selected according to the indicators that have had the greatest impact on the performance of this industry over many years. These inputs and outputs are as follows. Each DMU has three inputs, two desirable outputs, and two undesirable outputs. The first input is the amount of drilled oil wells during that year, including the fields of exploration, development, repair, and supplementary wells, measured in meters. The second input shows the amount of gasoline imports during different years, measured in one thousand cubic meters. Due to the inability of the country’s refineries to produce the gasoline needed by the country, gasoline is exported, and the required gasoline is imported from other producing countries. The third input includes the amount of different petroleum products transported by various means, including pipelines, railway tanks, road tankers, petroleum products trucks, refueling vessels, and refueling vessels, which are measured in millions of tons per kilometer.
Outputs include desirable and undesirable outputs. Desirable outputs include two outputs, which are factors whose larger value increases the rate of efficiency and performance. Total crude oil extracted from oil wells is used in three ways. The first part is the total amount of crude oil provided to the country’s refineries in the form of petroleum products, from which other petroleum products and derivatives are produced and used for domestic consumption and in various industries, and the rest of the refinery products are also exported. The second part is the amount of crude oil that is exported directly as furnace oil. The third part is exported as crude oil to other countries.
The first desirable output is the total production of the country’s petroleum products, which are produced in ten different refineries. These refineries produce petroleum products, including gasoline, liquefied petroleum gas, heavy furnace oil, sulfur, gas condensate, solvents, jet fuel, industrial oil, gas oil, crude oil, engine oils, kerosene, and other products. Most of the petrochemical products in the country are used in various industries and households, and a limited part is exported to other countries. The total amount of crude oil supplied to different refineries is measured in cubic meters per day. The second desirable output is the amount of crude oil exported to other countries during the 14 years studied. The amount of crude oil exported is measured in million barrels per year.
In this study, the undesirable outputs include the emission of polluting and greenhouse gases resulting from the country’s energy production and consumption. Emission coefficients of pollutants and greenhouse gases used in the power plant sector from 2005 onwards, based on the studies conducted under the title of “Compilation of Pollution Atlas of Power Plants in 2005,” have been revised. In this study, nitrogen oxides and sulfur, which are the most significant pollutants and greenhouse gases emitted from the country’s energy production and consumption, are considered as undesirable outputs. The emissions of these gases (NOx, SO2) resulting from the country’s energy production and consumption in each of the studied years are considered as undesirable outputs. Therefore, the first and second undesirable outputs are the emissions of NOx and SO2 gases caused by the country’s energy production and consumption in each of the studied years, respectively. The dataset is provided in Table 5.
We first solve model (1) to obtain the efficiency score of each DMU. Model (1) calculates the efficiency score of each DMU without considering the undesirable output. We can see that Units 1, 3, 7, 8, 9, 11, 12, and 14 are efficient, while the others are inefficient. The optimal scores of the inefficient slacks related to the two input components and the two desirable output components are provided in the other columns of Table 6.
Now, to calculate the efficiency of each DMU in the presence of undesirable outputs, we solve model (2) and obtain the minimum level of undesirable output based on the levels of input and desirable output. As stated in the third section, we first solve model (2) and introduce a counterpart (hypothetical) unit for each of the DMUs. Table 7 displays the slack of undesirable outputs and the corresponding values of undesirable outputs for each of the counterpart DMUs.
Model (3) is used to obtain the corresponding efficiency for each DMU. We can observe that Units 1, 3, 9, 11, and 14 are efficient, while the others are inefficient. Table 8 illustrates the efficiency scores, input, desirable output, and undesirable output slacks corresponding to all DMUs based on model (3).
In order to perform the SBM-based inverse non-radial DEA process, we solve model (4) and obtain the efficiency scores and slacks. Table 6 shows the results of model (4). We can see that units 1, 3, 9, 11, and 14 are efficient and others are inefficient. The model (4) introduces the oil industry as efficient in 2005, 2007, 2013, 2015, and 2018. Over these years, due to the low amount of costs incurred in this sector, including the cost of operating wells and oil rigs, the cost of transporting crude oil to the country’s refineries, and the cost of imported gasoline, we observe favorable performance for the oil industry in the energy sector. However, during the years 2006, 2008, 2009, 2010, 2011, 2012, 2014, and 2016, the oil industry did not perform well in the energy sector due to rising costs. The model (4) introduces the oil industry as inefficient in these years.
The difference between models (3) and (4) is that the former measures the efficiency scores based on the counterpart (hypothetic) units. Model (3) is based on the SBM model in terms of the slack variables of the input, desirable output, and undesirable output components. Model (4) is similar to model (3), but model (4) does not consider the slacks of the desirable output components. We can say that model (3) is the input-oriented SBM model in the presence of undesirable outputs. In models (3) and (4), we consider the undesirable outputs, but in model (1), we do not consider the undesirable outputs. As can be seen in Tables (6), (7), and (8), the efficiency scores of each DMU based on model (1) are greater than their corresponding ones from models (3) and (4), which is not correct because we expect the efficiency score of each DMU to deteriorate in the presence of undesirable outputs. Therefore, it is recommended that instead of solving model (1), we use models (3) and (4) to calculate the efficiency of the DMUs.
As stated in the fourth section, in the inverse non-radial DEA process, we must solve model (8), through which the lowest inputs and undesirable outputs can be determined. Before solving model (8), we must first solve model (4) and obtain the efficiency scores and optimal slacks of the inputs and undesirable outputs of the unit under evaluation. As previously stated, in the inverse non-radial DEA process, model (8) obtains the minimum level of inputs and undesirable outputs from the created new unit in a manner that the created new unit’s efficiency is the same as that of the original unit. Table 9 shows the optimal objective function score and the deviation of the input, desirable outputs, and undesirable output components of the 14 DMUs based on model (8).
Table 10 shows the new levels of the input, desirable outputs, and undesirable output components corresponding to each of the DMUs derived from model (8). As described in the fourth section, the vector of the inputs, desirable outputs, and undesirable output of the created new unit in the inverse non-radial DEA process is as follows:
It should be noted that Table 5 lists the vector of the inputs, desirable outputs, and undesirable output corresponding to each of the original DMUs as DMU k =(x k , y k , u k ).
Now, in order to analyze the results of the models, we will examine the results of DMU 4, which is an inefficient unit. Based on the proposed algorithm, model (4) is solved to obtain the efficiency and the optimal levels of slacks of the inputs and undesirable outputs for DMU 4. The efficiency score is equal to 0.8168, the optimal levels of slacks of the inputs, desirable outputs, and undesirable outputs are equal to 0.8168, 58889.57, 0, 4896.477, 357445.3, and 720117.7 respectively. These values indicate the rate of improvement in the inputs and undesirable output components. They are obtained by a non-radial projection of these units on the frontier of efficiency, which is built on by the counterpart units on the basis of the model (4).
In the following, we solve model (8), based on which we obtain the input and output variations. The input variations for the first, second, and third inputs (the amount of drilled oil wells during 2008 year, the amount of gasoline imports during 2008 year, the amount of different petroleum products transported by various means during 2008 year) are 14434.02, 1512.7 and 1699.61 respectively. This shows that the first and second inputs can be increased and decreased by 14434.02, 1512.7, respectively, and third input can be decreased by 1699.61. In this case, the newly created unit, , has an efficiency score of 0.8168, which is equal to the efficiency score of the original unit, i.e., DMU 4. Similarly, the desirable outputs (the total production of the country’s petroleum products during 2008 year, the amount of crude oil exported to other countries during 2008 year) variations for the first and second desirable outputs are -9960 and -41.8, respectively. This shows that the first and second desirable outputs can be decreased by 9960 and 3723, respectively. In this case, the newly created unit, , has an efficiency score of 0.8168, which is the same as the one of the original unit, i.e. DMU 4. The undesirable outputs (the emissions of NO x gases caused by the country’s energy production and consumption during 2008 year, the emissions of SO 2 gases caused by the country’s energy production and consumption during 2008 year) variation is -241753 and -198194. This shows that the first and second undesirable outputs can be decreased by 241753 and 198194, respectively. With this amount of changes, the efficiency score of the newly created unit, , is the same as the one of the original unit, i.e. DMU 4. Of course, all the changes related to the inputs, desirable output and undesirable output components are done together. Then, the optimal levels of the newly created unit’s inputs and outputs corresponding to DMU 18 based on model (8) are considered as follows:
The efficiency scores of the counterpart (hypothetic) units corresponding to DMU 4, and based on models (9) and (12) are equal, which is 0.7657.
To show that the efficiency scores of the original DMU and the newly created unit are equal, we can solve the criterion models (9) and (12). We can observe that the efficiency scores of the new unit based on model (9) are equivalent to those of their original units based on model (4). Model (9) with the SBM model structure obtains the efficiency scores of the new unit, and we acquire their input and undesirable output levels from model (8). Additionally, model (9) displays the new unit non-radially on the frontier of efficiency that is created based on the counterpart units. Subsequently, we use the criterion model (12) to determine the efficiency score of the new unit. The level of inputs and outputs of this new unit was determined based on model (8). We can see that the efficiency scores of the new unit based on model (12) are the same as those of the corresponding unit based on model (4). We find that the efficiency scores obtained from models (4), (9), and (12) are the same. Therefore, we have demonstrated the validity of the theorems (1) and (2). Table 11 displays the results based on models (9) and (12).
Now, we compare the results obtained from our approach with the results of Hadi-Vencheh and Foroughi (2006HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.). In this way, we ignore undesirable output in the inverse DEA process. We also assume that in the radial model it is possible to decrease or increase inputs. At first, we apply input orientation model in the radial form. The results are given in Table 12.
The results of Hadi-Vencheh and Foroughi (2006HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.) in the input oriented.
In the following, we compare the results obtained from our approach in this paper with approach provided by Hadi-Vencheh and Foroughi (2006HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.). We consider models in the output orientation and radial form. Also, we ignore undesirable output in the inverse DEA process. The results are given in Table 13.
The results obtained from Hadi-Vencheh and Foroughi (2006HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.) in the output oriented.
In the comparison made between the approach presented in this paper and the approach presented by Hadi-Vencheh and Foroughi (2006HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.) based on the inverse DEA process, we did not consider the undesirable outputs, then the results are not comparable in any way, and the results are provided only so that the reader can see the difference between approaches based on radial and non-radial models. As can be seen, the input vector based on the approach presented in this paper and the approach of Hadi-Vencheh and Foroughi (2006HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.) are larger in some components and smaller in others according to the results in tables (10) and (12). We can provide a similar interpretation for the output vector according to the results in tables (10) and (13).
In general, the approaches provided in the field of inverse DEA are divided into two categories. The first category includes approaches that consider two or more DMUs. By merging the levels of inputs and outputs of these DMUs, a new unit is proposed; this unit has a predetermined target efficiency. However, in this paper, the inverse DEA process is done only based on the DMU under evaluation, and the target efficiency is predetermined. This efficiency score is the efficiency of the DMU under evaluation based on the DEA model. Also, by solving only one model, we can determine the optimal level of inputs and outputs from the unit under evaluation in such a way that the newly created unit has the same efficiency score as the corresponding original unit.
The second category includes approaches that have an orientation. In this approach, the decision-maker predetermines the target efficiency score. In this situation, we are faced with two types of models in the inverse DEA process. These models are input- and output-oriented. In the input-oriented approach, we obtain the optimal level of inputs according to the target efficiency score. Also, the output level is constant during the inverse DEA process. Similarly, in the output orientation, we obtain the optimal level of outputs according to the target efficiency score. The input level is constant during the inverse DEA process. However, our approach in this paper is different from these approaches, and the target efficiency score is the same as the efficiency score of the unit under evaluation. The model determines the optimal level of inputs and outputs simultaneously. In other words, the model determines the best level of inputs and outputs (undesirable and desirable) for the unit under evaluation based on the conditions of this unit.
Therefore, this paper introduces a new approach in the field of inverse DEA that is different from previous approaches.
7 CONCLUSION
In this paper, we present the SBM-based inverse non-radial DEA process in the presence of undesirable outputs. Firstly, we presented the counterpart (hypothetical) units corresponding to each of the DMUs. These units have the same input and desirable output levels as the corresponding original units, but they have the lowest level of undesirable outputs. The observed units are used on the basis of weak disposability in DEA. We obtain the true efficiency scores and slacks regarding the input and undesirable output components of each DMU by depicting the units non-radially on the frontier of efficiency created on the basis of the counterpart (hypothetical) units. Next, based on the efficiency score and slacks regarding the input and undesirable output components of each DMU, the SBM-based inverse non-radial DEA process is presented. We demonstrated that this model obtains the optimal level of input, desirable outputs, and undesirable output components for the created new unit in the inverse non-radial DEA process. We presented two new criteria models to compare the efficiency scores between the original unit and the created new unit and showed that these efficiency scores are equal based on the two criteria models presented in this paper. By using the SBM-based inverse non-radial DEA process, we can use all the inefficiency slacks related to all the input and undesirable output components in the inverse non-radial DEA process to estimate the inputs and outputs simultaneously. This paper proposes a new idea of creating a set of hypothetical counterpart DMUs of the original DMUs. These DMUs used the same amount of inputs to produce the same amount of desirable outputs, however generated the smallest amount of undesirable outputs calculated from the observed DMUs based on weak disposability. We used hypothetical DMUs to construct the production frontier in order to do inverse DEA process. We do inverse DEA process using hypothetical DMUs. The new level of inputs, desirable output and undesirable output of the original DMUs obtained by considering the new frontier created by these DMUs.
We can say that the advantages of our approach in this paper over the previous approach are as follows:
-
In our proposed approach in this paper for the inverse DEA process, we presented two criteria models in order to check whether the efficiency score of the new unit created and the initial unit corresponding to it are equal or not. The new unit created in the inverse DEA process has an efficiency score equal to that of the corresponding initial unit; however, this unit has input and undesirable output levels less than or equal to the input and undesirable output levels of the initial unit corresponding to this unit. Also, the undesirable output level of the new unit created is greater than or equal to the undesirable output level of the initial unit corresponding to itself.
-
This paper presents a new approach to the inverse DEA. By solving a model, we can determine the optimal level of inputs and outputs from the unit under evaluation in such a way that this unit has the same efficiency score as the corresponding initial unit. It can be said that the inverse DEA process presented in this paper is done only based on the unit under evaluation. The new unit created can be considered the target (benchmark) unit corresponding to the initial unit corresponding to this unit.
-
Most of the approaches provided in the inverse DEA are based on radial models. However, the our proposed approach in this paper is based on non-radial models. In non-radial models such as SBM, the target corresponding to the unit under evaluation is determined in the best way. These models obtain the efficiency score of the unit under evaluation based on all the inefficiency slack values corresponding to the input and output components. The image of the inefficient unit on the efficiency frontier is determined non-radially. These models obtain the optimal level of input and output components for the target unit. Also, considering that the inverse DEA process presented in this paper is based on non-radial models, the target unit corresponding to the initial unit is determined in the best way possible.
-
The approach presented in this paper provides the ability to deal with undesirable outputs in the inverse DEA process. In this regard, we first determine the lowest level of undesirable output. In the following, new units are introduced with the title of counterpart units corresponding to each of the original units. The counterpart unit has the same level of inputs and desirable outputs as the corresponding original unit; however, it has a Less amount of undesirable outputs than the original unit. We perform the inverse DEA process for original units based on counterpart units. In other words, the unit under evaluation is projected on the new frontier created by the counterpart units. Next, in the inverse DEA process, we determine the optimal level of input and output (desirable and undesirable) components from the target unit corresponding to the original unit. It can be said that the target unit has the best level of inputs and outputs (desirable and undesirable) compared to the corresponding initial unit, but it has the same efficiency score as this unit.
-
The models presented in this paper are without orientation. These models get the optimal level of inputs and outputs simultaneously. However, the models in the previous approaches were often presented as input- and output-oriented. In input-output-oriented models, we first place the target efficiency score in the inverse DEA model, and the new level of inputs (outputs) is determined by keeping the level of outputs (inputs) constant.
-
The inverse DEA model presented in this paper can be converted into a linear model with a simple transformation, and it can be easily solved with available software for solving linear programming models.
-
The approach presented in this paper determines the target unit corresponding to the unit under evaluation in the presence of undesirable outputs. The target unit can be considered a benchmark for the original unit under evaluation. Although this unit may be an inefficient one, it has optimal input and output levels compared to the original unit. In other words, we showed that, considering the existing conditions, the unit under evaluation can have a better performance than its initial performance and should improve the level of its inputs and outputs.
-
In the models presented in this paper, the level of inputs and outputs may decrease or increase, but the amount of input and output from the newly created unit is non-negative. And it is possible for the units to adjust the level of their inputs and outputs in order to reach the optimal level of their inputs and outputs.
-
From a managerial point of view, the target unit achieved for the unit under evaluation may not be able to be created in practice. As we know, the target unit corresponding to each of the DMUs based on DEA models may not be able to be created, and the optimal level of inputs and outputs cannot be reached in some cases. However, the DEA offers targets. Only if the target unit is an observed unit does it have input and output levels that can be created and are real. However, if the target unit is a virtual unit, it may not be able to be created. The models presented in this paper are also based on the SBM model in DEA and are not exempt from this rule. In DEA, we can achieve attainable targets only if we use FDH models. However, by changing the conditions of production, it is possible to achieve targets that can be produced in practice.
In future work, we can apply the proposed approach in the industrial sector, as one of the primary issues in the industrial sector is how to deal with the outcome of undesirable outputs among the data and how to determine and control the level of these outputs to protect the environment. One of the limitations of the research is that we may obtain targets based on inverse DEA models for DMUs that cannot be created in the real world. In terms of future studies, we can solve the proposed approach in this paper by considering the target efficiency of DMUs and obtaining the amount of deviation of the input and output components needed to reach the target efficiency level. We can also solve the models presented in this paper for other technologies such as CRS, IRS, DRS, and semi-additivity technologies. Additionally, we can apply the proposed approach to other data structures in DEA, such as the two-stage network structures (Mills et al., 2021MILLS E, ZENG K, LIU F & LI F. 2021. Modelling innovation efficiency, its micro-level drivers, and its impact on stock returns. Chaos, Solitons & Fractals, 152: 111303.), or use them in the presence of inaccurate data, such as fuzzy and probabilistic data (Valami, 2009VALAMI H. 2009. Cost efficiency with triangular fuzzy numbers in input prices: an application of DEA. Chaos, Solitons & Fractals, 42(3): 1631-1637.).
References
- ADENUGA A, DAVIS J, HUTCHINSON G, DONNELLAN T & PATTON M. 2018. Modelling regional environmental efficiency differentials of dairy farms on the island of Ireland. Ecological Indicators, 95: 851-861.
- ALEMDAR T & OREN M. 2006. Determinants of technical efficiency of wheat farming in Southeastern Anatolia, Turkey: a nonparametric technical efficiency analysis. Journal of Applied Sciences, 6: 827-830.
- AMIN G, EMROUZNEJAD A & GATTOUF S. 2017a. Minor and major consolidations in inverse DEA: Definition and determination. Computers & Industrial Engineering, 103: 193-200.
- AMIN G, EMROUZNEJAD A & GATTOUF S. 2017b. Modelling generalized firms’ restructuring using inverse DEA. Journal of Productivity Analysis, 48: 51-61.
- AMIN G & IBNBOAMAH M. 2021. A two-stage inverse data envelopment analysis approach for estimating potential merger gains in the US banking sector. Managerial and Decision on Economics, 42(6): 1454-1465.
- ANDERSON T, DAIM T & KIM J. 2008. Technology forecasting for wireless communication. Technovation, 28: 602-614.
- ANGULO-MEZA L, GONZALEZ-ARAYA M, IRIARTE A, REBOLLEDO-LEIVA R & MELLO J. 2019. A multiobjective DEA model to assess the eco-efficiency of agricultural practices within the CF+DEA method. Computers and Electronics in Agriculture, 161: 151-161.
- CECCHINI L, VENANZI S, PIERRI A & CHIORRI M. 2018. Environmental efficiency analysis and estimation of CO2 abatement costs in dairy cattle farms in Umbria (Italy): a SBM-DEA model with undesirable output. Journal of Cleaner Production, 197: 895-907.
- CHARNES A, COOPER W & RHODES E. 1978. Measuring the efficiency of decision making units. European Journal of Operational Research, 2(6): 429-444.
- CHEBIL A, FRIJA A & THABET C. 2015. Economic efficiency measures and its determinants for irrigated wheat farms in Tunisia: a DEA approach. New medit, 14: 32-38.
- CHEN Z, WANKE P, ANTUNES J & ZHANG N. 2017. Chinese airline efficiency under CO2 emissions and flight delays: A stochastic network DEA model. Energy Economics, 68: 89-108.
- DAKPO K, JEANNEAUX P & LATRUFFE L. 2016. Modelling pollution-generating technologies in performance benchmarking: Recent developments, limits and future prospects in the nonparametric framework. European Journal of Operational Research, 250: 347-359.
- DONG G, WANG Z & MAO X. 2018. Production efficiency and GHG emissions reduction potential evaluation in the crop production system based on emergy synthesis and nonseparable undesirable output DEA: A case study in Zhejiang Province, China. PLOS ONE, 13(11): 0206680.
- EMROUZNEJAD A, YANG G & AMIN G. 2019. A novel inverse DEA model with application to allocate the CO2 emissions quota to different regions in Chinese manufacturing industries. Journal of the Operational Research Society, 70(7): 1079-1090.
- ENTANI T, MAEDA Y & TANAKA H. 2002. Dual Models of Interval DEA and its extensions to interval data. European Journal of Operational Research, 136(1): 32-45.
- FRIJA A, WOSSINK A, BUYSSE J, SPEELMAN S & HUYLENBROECK G. 2011. Irrigation pricing policies and its impact on agricultural inputs demand in Tunisia: A DEA-based methodology. Journal of Environmental Management, 92(9): 2109-2118.
- FÄRE R, GROSSKOPF S, LOVELL C & PASURKA C. 1989. Multilateral productivity comparisons when some outputs are undesirable: A nonparametric approach. The Review of Economics and Statistics, 71: 90-98.
- GERAMI J, MOZAFFARI M, WANKE P & CORREA H. 2022. Improving information reliability of non-radial value efficiency analysis: an additive slacks based measure approach. European Journal of Operational Research, 298(3): 967-978.
- GERAMI J, MOZAFFARI M, WANKE P & HL C. 2023. Generalized inverse DEA model for firm restructuring based on value efficiency. IMA Journal of Management Mathematics, 34(3): 541-580.
- GHIYASI M. 2017. Inverse DEA based on cost and revenue efficiency. Computers & Industrial Engineering, 114: 258-263.
- GHIYASI M. 2019. Novel criterion models in the inverse DEA problem. International Journal of Operational Research, 35(1): 20-36.
- GHIYASI M & ZHU N. 2020. An inverse semi-oriented radial data envelopment analysis measure for dealing with negative data. IMA Journal of Management Mathematics, 31(4): 505-516.
- GHOBADI S. 2018. Inverse DEA using enhanced Russell measure in the presence of fuzzy data. International Journal of Industrial Mathematics, 10(2): 165-180.
- GOLANY B & ROLL Y. 1989. An application procedure for DEA. Omega, 17: 237-250.
- HADI-VENCHEH A & FOROUGHI A. 2006. A generalized DEA model for inputs/outputs estimation. Mathematical and Computer Modelling, 43(5): 447-457.
- HADI-VENCHEH A, FOROUGHI A & SOLEIMANI-DAMANEH M. 2008. A DEA model for resource allocation. Economic Modelling, 25(5): 983-993.
- HASSANZADEH A, YOUSEFI S, SAEN R & HOSSEININIA S. 2018. How to assess sustainability of countries via inverse data envelopment analysis? Clean Technologies and Environmental Policy, 20: 29-40.
- INMAN O, ANDERSON T & HARMON R. 2006. Predicting US jet fighter aircraft introductions from 1944 to 1982: A dogfight between regression and TFDEA. Technological Forecasting and Social Change, 73: 1178-1187.
- IZADIKHAH M & SAEN R. 2018. Assessing sustainability of supply chains by chance-constrained two-stage DEA model in the presence of undesirable factors. Computers & Operations Research, 100: 343-367.
- JAHANSHAHLOO G, HOSSEINZADEH LOTFI F, ROSTAMY-MALKHALIFEH M & GHOBADI S. 2014. Using enhanced Russell model to solve inverse data envelopment analysis problems. The Scientific World Journal, p. 1-10.
- JAHANSHAHLOO G, SOLEIMANI-DAMANEH M & GHOBADI S. 2015. Inverse DEA under inter-temporal dependence using multiple-objective programming. European Journal of Operational Research, 240(2): 447-456.
- KAO C. 2017. Network data envelopment analysis: Foundations and extensions. Switzerland: Springer International Publishing.
- KAO C. 2020. Measuring efficiency in a general production possibility set allowing for negative data. European Journal Operational Research, 282(3): 980-988.
- KAO C & HWANG S. 2021. Measuring the effects of undesirable outputs on the efficiency of production units. European Journal of Operational Research, 292(3): 996-1003.
- KOEIJER T, WOSSINK G, STRUIK P & RENKEMA J. 2002. Measuring agricultural sustainability in terms of efficiency: the case of Dutch sugar beet growers. Journal of Environmental Management, 66: 9-17.
- KUOSMANEN T & KORTELAINEN M. 2005. Measuring eco-efficiency of production with data envelopment analysis. Journal of Industrial Ecology, 9: 59-72.
- LERTWORASIRIKUL S, P C & FANG S. 2011. Inverse data envelopment analysis model to preserve relative efficiency values: The case of variable returns to scale. Computers and Industrial Engineering, 61(4): 1017-1023.
- LIM D. 2016. Inverse DEA with frontier changes for new product target setting. European Journal of Operational Research, 254(2): 510-516.
- LIM DJ, ANDERSON T & SHOTT T. 2015. Technological forecasting of supercomputer development: The march to exascale computing. Omega, 51: 128-135.
- LIN H. 2010. An efficiency-driven approach for setting revenue target. Decision Support Systems, 49(3): 311-317.
- LIU W, MENG W, LI X & ZHANG D. 2010. DEA models with undesirable inputs and outputs. Annals of Operations Research, 173: 177-194.
- LOZANO S, IRIBARREN D, MOREIRA M & FEIJOO G. 2009. The link between operational efficiency and environmental impacts: a joint application of life cycle assessment and data envelopment analysis. Science of the Total Environment, 407: 1744-1754.
- MILLS E, ZENG K, LIU F & LI F. 2021. Modelling innovation efficiency, its micro-level drivers, and its impact on stock returns. Chaos, Solitons & Fractals, 152: 111303.
- MOHSENI P, BORGHEI A & KHANALI M. 2018. Coupled life cycle assessment and data envelopment analysis for mitigation of environmental impacts and enhancement of energy efficiency in grape production. Journal of Cleaner Production, 197: 937-947.
- MURTY S, RUSSELL R & LEVKOFF S. 2012. On modeling pollution-generating technologies. Journal of Environmental Economics and Management, 64: 117-135.
- NABAVI-PELESARAEI A, RAFIEE S, MOHTASEBI S, HOSSEINZADEH-BANDBAFHA H & CHAU KW. 2017. Energy consumption enhancement and environmental life cycle assessment in paddy production using optimization techniques. Journal of Cleaner Production, 162: 571-586.
- PICAZO-TADEO A, GOMEZ-LIMON J & REIG-MARTÍNEZ E. 2011. Assessing farming ecoefficiency: a data envelopment analysis approach. Journal of Environmental Management, 92: 1154-1164.
- PISHGAR-KOMLEH S, ZYLOWSKI T, ROZAKIS S & KOZYRA J. 2020. Efficiency under different methods for incorporating undesirable outputs in an LCA+DEA framework: A case study of winter wheat production. Journal of Environmental Management, 260, Article 110138.
- REBOLLEDO-LEIVA R, ANGULO-MEZA L, IRIARTE A & GONZALEZ-ARAYA M. 2017. Joint carbon footprint assessment and data envelopment analysis for the reduction of greenhouse gas emissions in agriculture production. Science of the Total Environment, 593: 36-46.
- REN J, GAO B, ZHANG J & CHEN C. 2020. Measuring the energy and carbon emission efficiency of regional transportation systems in China: Chance-constrained DEA models. Mathematical Problems in Engineering, Available at: https://doi.org/10.1155/2020/9740704
» https://doi.org/10.1155/2020/9740704 - SHI G, BI J & WANG J. 2010. Chinese regional industrial energy efficiency evaluation based on a DEA model of fixing non-energy inputs. Energy Policy, 38(10): 6172-6179.
- SONG M, AN Q, ZHANG W, WANG Z & WU J. 2012. Environmental efficiency evaluation based on data envelopment analysis: A review. Renewable and Sustainable Energy Reviews, 16(7): 4465-4469.
- SUEYOSHI T & GOTO M. 2010. Should the US clean air act include CO2 emission control? Examination by data envelopment analysis. Energy Policy, 38(10): 5902-5911.
- SUEYOSHI T & GOTO M. 2012. DEA radial and non-radial models for unified efficiency under natural and managerial disposability: Theoretical extension by strong complementary slackness conditions. Energy Economics, 34: 700-713.
- SYP A, FABER A, BORZECKA-WALKER M & OSUCH D. 2015. Assessment of greenhouse gas emissions in winter wheat farms using data envelopment analysis approach. Polish Journal of Environmental Studies, 24.
- TONE K. 2001. A slacks-based measure of efficiency in data envelopment analysis. European Journal of Operational Research, 130(3): 498-509.
- ULLAH A, PERRET S, GHEEWALA S & SONI P. 2016. Eco-efficiency of cotton-cropping systems in Pakistan: an integrated approach of life cycle assessment and data envelopment analysis. Journal of Cleaner Production, 134: 623-632.
- VALAMI H. 2009. Cost efficiency with triangular fuzzy numbers in input prices: an application of DEA. Chaos, Solitons & Fractals, 42(3): 1631-1637.
- VAZQUEZ-ROWE I, IRIBARREN D, MOREIRA M & FEIJOO G. 2010. Combined application of life cycle assessment and data envelopment analysis as a methodological approach for the assessment of fisheries. The International Journal of Life Cycle Assessment, 15: 272-283.
- WEGENER M & AMIN G. 2019. Minimizing greenhouse gas emissions using inverse DEA with an application in oil and gas. Expert Systems with Applications, 122: 369-375.
- WEI Q, ZHANG J & ZHANG X. 2000. An inverse DEA model for inputs/outputs estimate. European Journal of Operational Research, 121(1): 151-163.
- WU C, LI Y, LIU Q & WANG K. 2013. A stochastic DEA model considering undesirable outputs with weak disposability. Mathematical and Computer Modelling, 58(5-6): 980-989.
- XIAOYA L & JINCHUAN C. 2008. A comprehensive DEA approach for the resource allocation problem based on scale economies classification. Journal of Systems Science and Complexity, 21(4): 540-557.
- YEH T, CHEN T & LAI P. 2010. A comparative study of energy utilization efficiency between Taiwan and China. Energy Policy, 38(5): 2386-2394.
- YOU S & YAN H. 2011. A new approach in modelling undesirable output in DEA model. Journal of the Operational Research Society, 62(12): 2146-2156.
- ZEINODIN E & GHOBADI S. 2020. Merging decision-making units under inter-temporal dependence. IMA Journal of Management Mathematics, 31(2): 139-166.
- ZHANG G & CUI G. 2020. A general inverse DEA model for non-radial DEA. Computers & Industrial Engineering, 142: 106368.
- ZHANG M & CUI JC. 2016. The extension and integration of the inverse DEA method. Journal of the Operational Research Society, 67: 1212-1220.
- ZHANG XS & CUI JC. 1999. A project evaluation system in the state economic information system of China an operations research practice in public sectors. International Transactions in Operational Research, 6(5): 441-452.
- ZHOU P, ANG B & POH K. 2008. A survey of data envelopment analysis in energy and environmental studies. European Journal of Operational Research, 189: 1-18.
Publication Dates
-
Publication in this collection
02 Sept 2024 -
Date of issue
2024
History
-
Received
10 Sept 2023 -
Accepted
25 Nov 2023