# Performance Analysis in Production Systems with Uncertain Data: A Stochastic Data Envelopment Analysis Approach

Authors: Seyedeh Fatemeh Bagheri; Alireza Amirteimoori; Sohrab Kordrostami; Mansour Soufi
Journal: Complexity
 (2022)

---

## Abstract

The problem of determining an optimal benchmark to inefficient decision-making units (DMUs) is an important issue in the field of performance analysis. Previous methods for determining the projection points of inefficient DMUs have only focused on one objective and other features have been ignored. This paper attempts to determine the best projection point for each DMU when the inputs and outputs data are in stochastic form and presents an alternative definition for the best projection by considering three main aspects: technical efficient, minimal cost, and maximal revenue as much as possible. Considering the important role of the electricity industry in the economic growth of each country, a practical example has been implemented on 16 regional electricity companies in Iran in 9 consecutive periods. The efficiency score along with the projection points of the three technical models (BCC model of Banker et al. (1984)), cost, and stochastic revenue are compared with the projection point obtained from the model presented in this article, which simultaneously meets these three objectives, showing the improvement of companies’ performance.


---

## Body

## 1. Introduction

Data envelopment analysis (DEA) was first introduced by Charnes et al. [1] (CCR model) to calculate the relative efficiency of a set of homogeneous decision-making units (DMU) and extended by Banker et al. [2] (BCC model). Traditional DEA models (BCC and CCR models) rely on past information and they are developed in deterministic case, and hence the uncertainty and stochasticity of the data are ignored. In the real world, we are faced with many decisions, most of which are associated with uncertainty. Stochastic data envelopment analysis models to address these shortcomings were proposed by Charnes and Cooper [3]. In the last two decades, performance analysis in an uncertain environment has been studied frequently. See for instances Hossain et al. [4], Mo et al. [5], Rong et al. [6], Amirteimoori et al. [7], and Ghasemi et al. [8].Some of the differences between DEA and stochastic DEA (SDEA) models have been presented in Table1. SDEA method is based on stochastic programming methods, which is one of the branches of mathematical programming. Stochastic programming models are divided into E-model and P-model in terms of the type of objective function. E-model based stochastic programming models are used to obtain the expected value or mathematical expectation that the objective function of such models does not include a random variable. To obtain the highest probability of an event occurring, stochastic programming models based on P-model are used so that in the objective function of such models, a random variable is used possibly. Charnes and Cooper [9] proposed a stochastic programming model with potential constraints. They were then introduced by Land et al. [10] E-model-based stochastic programming model. Programming models based on E-model are used to obtain the expected value. The objective function of these models does not include random variables.Table 1
Differences between models.
DEA modelSDEA modelThe values of inputs and outputs are definite.The values of inputs and outputs are random.No accidental error is entered in the modelRandom error enters the model in the form of a random componentThe model offers only past performanceAllows you to predict performance in the futureThe efficiency frontier is very sensitive to small changes in inputs and outputs.The efficiency frontier is less sensitive to changes in inputs and outputs.The error level in such models is considered zero.Efficiency is defined according to the level of error in the model.In the final model, the correlation between inputs and outputs is not considered.The correlation between inputs and outputs is considered in such models.The efficiency measure in which, in addition to input and output data, the price of the input of the DMUs is also required and its goal is cost minimization, called cost efficiency, was introduced by Farrell [11]. Färe et al. [12] developed cost efficiency. By studying data envelopment analysis and stochastic cost functions, Cooper and Ton [13] identified some specific problems in stochastic frontier analysis (SFA). Stochastic frontier analysis method is a parametric method. Shiraz et al. [14] formulated the stochastic cost model in a nonlinear form.Färe et al. [12] have played a key role in the development of revenue efficiency. In addition to input and output data, they have also paid attention to the output of the decision-making unit aiming at maximizing revenue. Lin [15] proposed a method to set revenue targets at the efficient frontier. Shiraz et al. [14] also presented a random income model that is nonlinear.In data envelopment analysis, DMUs are divided into two categories: efficient and inefficient units. Inefficient units can be efficient by decreasing their inputs level and increasing their outputs level.Achieving optimal scale size (OSS) has always been of interest to researchers. Cesaroni and Giovannola [16] introduced a new definition of optimal scale size based on minimizing the cost of each unit so that average cost productivity combines scale productivity and allocation and generalizes economic scale size in productivity analysis. In the following, Haghighatpisheh et al. [17] proposed a new definition of optimal scale size that uses both of the cost of inputs and the revenue of outputs. This measure of average cost-revenue efficiency (ACRE) shows the ratio of profitability efficiency to average productivity.The process of industrialization of developed countries shows proper planning, optimal use of resources, and determining the appropriate pattern as the main goal in the development process of these countries. The electricity industry has an infrastructural and influential role on the economic growth of any country. In recent years, per capita electricity consumption in Iran has increased significantly in parallel with its production rate. Accordingly, efficiency and determining the appropriate model for companies has always been the concern of managers in the electricity industry. In this article, a real and useful example has been done on Iran’s regional electricity distribution companies and it has used random data that can play an important role as a vision in the electricity industry and solving future challenges.Most studies to find the projection point of each inefficient unit have focused on only one goal. For example, the projection point obtained from the cost model finds an efficient point on efficient frontier with the cost minimization approach. The projection point obtained from the revenue model moves towards the frontier with the revenue maximization approach. Bagheri et al. [18] presented a method that simultaneously tried to meet the projection point of each unit from three perspectives, cost minimization, revenue maximization, and the shortest distance to the efficient frontier as much as possible and it was proved that the efficiency of the model obtained for inefficient units is higher than the unit under evaluation. The results of solving this method for managers to make decisions were based on past information, so we tried to use stochastic programming. In this paper, for each inefficient unit, a projection point is obtained, which is very important in two stages. First, the resulting projection for each unit meets as many as possible the three technical, cost, and revenue concepts. Secondly, random data and stochastic programming in linear form are used in it. Therefore, it will greatly help managers to make better decisions in the future and improve the performance of units.This paper has been organized as follows: technical, cost, and revenue efficiency with deterministic and random data are examined in Sections2 and 3, respectively. The random models of the linearly shaped models have been presented in Section 3. In Section 4, the method proposed by Bagheri et al. [18] to determine the ideal projection point has been implemented this time in stochastic programming. In Section 5, we will provide a real and practical example, and finally in Section 6, the conclusion of the article is presented.


## 2. Deterministic DEA Models

In this section, we briefly review some basic DEA models in deterministic environment. Technical efficiency measure was first introduced by Debreu [19] and Farrell [11]. The traditional DEA-based models CCR model of Charnes et al. [1] and its subsequent extension BCC model of Banker et al. [2] was developed on classical efficiency analysis model of Farrell [11]. The BCC model for evaluating a specific DMUp is formulated in the input-oriented deterministic environment as follows:(1)Minθ∑j=1nλjxij≤θxip,i=1,…,m,∑j=1nλjyrj≥yrp,r=1,…,s,∑j=1nλj=1,λj≥0,j=1,…,n,θpfree variable.In the above model,θ is the input abatement factor. DMUp is said to be technically efficient if and only θ=1 and all slack variables are equal to zero. Corresponding to each inefficient point, a frontier point is determined by reducing inputs. Removing the third convexity constraint ∑j=1nλj=1 leads to CCR model.For each inefficient unit, there is one or more efficient units on the efficient frontier that is called the reference set. The reference set of a specific inefficient DMU (denoted byEo) is defined as follows:(2)Eo=j|λj∗>0j∈1,…,n.In technical efficiency models we have less information on inputs and outputs of DMUs. In other words, all we know the quantitative values of inputs and outputs and there is no information on inputs and outputs prices. More information will be available if we also have the costs of the inputs and the prices of the outputs. In this case, we are motivated to calculate the allocative efficiency of the DMUs. We first introduce cost allocative model. Suppose thatC=c1,c2,…,cmt is the vector of input prices and X¯=x¯1,x¯2,…,x¯mt is cost minimizing vector of input quantities. Farrell (1957) proposed a measure of cost efficiency. The cost efficiency model in VRS environment is formulated as follows:(3)Min∑i=1mcix¯i,∑j=1nλjxij≤x¯i,i=1,…,m,∑j=1nλjyrj≥yrp,r=1,…,s,∑j=1nλj=1λj≥0,j=1,…,n,The cost efficiency ofDMUp is defined as the ratio of the optimal cost to the actual cost as follows:(4)CEp=∑i=1mcix¯i∑i=1mcixip.Clearly,CEp≤1 and when this score is equal to one, then DMUp is called cost efficient, otherwise we say DMUp is cost inefficient.Now suppose that we are given information on the price of the outputs. In this case, the revenue efficiency model in VRS environment is formulated as follows:(5)Max∑r=1spry¯r,∑j=1nλjxij≤xip,i=1,…,m,∑j=1nλjyrj≥y¯r,r=1,…,s,∑j=1nλj=1λj≥0,j=1,…,n.In the above equation,P=p1,p2,…,pst is vector of output prices and Y¯=y¯1,y¯2,…,y¯st is revenue maximizing vector of output quantities.The revenue efficiency ofDMUp is defined as the ratio of optimal revenue to the actual revenue, i.e., REp=∑r=1spry¯r/∑r=1spryrp.It is easy to see thatREp≥1. If this score is equal to one, then DMUp is called revenue efficient, otherwise we say DMUp is revenue inefficient.


## 3. Basic DEA Models in Stochastic Environment

As we stated before, in many real applications, we often encounter uncertainty, so it is necessary to generalize the models to the uncertainty mode, especially the random mode. SDEA models were first proposed by Charnes and Cooper (1959). Land et al. (1993) introduced the CCR model in the stochastic mode considering the inputs and outputs and estimating the probability distribution prevailing over them, after estimating the efficiency in the future, that by considering the convexity condition, the BCC model is obtained randomly in the nature of the input. In this model, we show the random input and output values byx˜ij and y˜rj, respectively. These values show the mathematical expectations of random inputs and outputs. Also, aij and brj are the standard deviations of random inputs and outputs, respectively. Using a random variable error structure, Cooper et al. [20] proposed the following input-oriented DEA model:(6)Minθ∑j=1nλjxij−∅−1∝σpi++pi−≤θxip,i=1,…,m,∑j=1nλjaij−θaip=pi+−pi−,i=1,…,m,∑j=1nλjyrj+∅−1∝σqr++qr−≥yrp,r=1,…,s,∑j=1nλjbrj−θbrp=qr+−qr−,r=1,…,s,∑j=1nλj=1λj,≥0,j=1,…,n,pi+,pi−≥0,i=1,…,m,qr+,qr−≥0,r=1,…,s.The first two constraints in model (4) are related to input variables and the next two constraints are related to output variables. The goal is to minimize the level of inputs by the abatement factor θ. σ has a standard normal distribution and ∅−1 shows the inverse of the standard normal distribution function. The positive deviation variables pi+,pi−, qr+,qr− have been used to transform the model into a linear form.Corresponding to each inefficient DMU, the projection point was obtained from the above stochastic model through the formula∑j=1nλjxij,∑j=1nλjyrj.Producing a certain level of outputs from the least cost of inputs is called cost efficiency, in which information about the price of inputs is considered. In the following, we present the cost efficiency model in stochastic environment. Inspired by the deterministic form of a stochastic model, we assume thatci is the i-th input price whose value is constant for all units. Also, xi¯ is considered to be the ideal input whose value is unknown.Suppose thatxij and yrj are the mathematical expectations of the inputs and outputs of x˜ij and y˜rj, respectively. Moreover, suppose aij and brj are their corresponding standard deviations. The cost efficiency model is formulated as follows:(7)Min∑i=1mcix¯i,∑j=1nλjxij−∅−1∝σpi++pi−≤x¯i,i=1,…,m,∑j=1nλjaij−θaip=pi+−pi−,i=1,…,m,∑j=1nλjyrj+∅−1∝σqr++qr−≥yrp,r=1,…,s,∑j=1nλjbrj−brp=qr+−qr−,r=1,…,s,∑j=1nλj=1,λj,≥0,j=1,…,n,pi+,pi−≥0,i=1,…,m,qr+,qr−≥0,r=1,…,s.In this model, the objective is to minimize the total cost of inputs. This model is formulated under variable returns to scale and the cost efficiency ofDMUp is defined as follows:(8)CEp=∑i=1mcix¯i∑i=1mcixip.Corresponding to each inefficient point, the efficient projection point is obtained from the above stochastic model through the formula (∑j=1nλjxij,∑j=1nλjyrj). It can be easily shown that the obtained projection point is efficient at the desired confidence level. σ has a standard normal distribution and ∅−1 shows the inverse of the standard normal distribution function in the above linear model.Now, we propose the revenue efficiency model. As it is known in DEA literature, producing the highest level of outputs from a certain level of inputs is called revenue efficiency. Assume thatwr is the r-th output price whose value is constant for all units. Also, yr¯ is considered to be the r-th optimal output value whose value is unknown. Suppose that xij and yrj are the mathematical expectations of the inputs and outputs of x˜ij and y˜rj, respectively. Also, aij and brj are the standard deviations of the inputs and outputs x˜ij and y˜rj, respectively. The revenue efficiency model is formulated in a linear form in a random environment as follows:(9)Max∑r=1swryr¯,s.t:∑j=1nλjxij−∅−1∝σpi++pi−≤xip,i=1,…,m,∑j=1nλjaij−θaip=pi+−pi−,i=1,…,m,∑j=1nλjyrj+∅−1∝σqr++qr−≥yr¯,r=1,…,s,∑j=1nλjbrj−brp=qr+−qr−,r=1,…,s,∑j=1nλj=1λj,≥0,j=1,…,n,pi+,pi−≥0,i=1,…,mqr+,qr−≥0,r=1,…,s.The first two constraints in the above model are related to input variables and the next two constraints are related to output variables. In this model, the goal is to maximize the revenue of outputs. This model is formulated under variable returns to scale and its efficiency score is obtained from the following formula.(10)REp=∑r=1swry¯r∑r=1swryrp.The projection point is obtained from the above linearized random revenue model through the formula (∑j=1nλjxij,∑j=1nλjyrj). It can be easily proved that the obtained projection point is efficient at the desired confidence level. σ has a standard normal distribution and ∅−1 shows the inverse of the standard normal distribution function in the above linear model.


## 4. The Best Projection Point

In this section, the best projection point in a random environment is determined. As you know, there is an ideal pattern for any inefficient unit. The ideal pattern means the projection point obtained as each decision unit. Bagheri et al. [18] presented a method in which an ideal model was obtained for each unit, which first dominates the unit under evaluation. Therefore, it will be much more efficient. Secondly, in the convex combination of projection points obtained from cost, revenue, and technical models, this projection point is simultaneously examined from three technical, minimum cost, and maximum revenue perspectives. One of the defects of such models is that they rely on past information, so they cannot provide the desired result. One way to solve this problem is to use random data envelopment analysis. In such models, a random error is added to the model in the form of a random component, in which the efficiency is defined according to the error level (∝.)and the correlation between the input and output variables is taken into account in the model. Toward this end, we proceed as follows:Step 1: First, we obtain the projection points of the technical, cost, and random revenue models in linear form under the variable returns to scale (models (4)–(6)), which were described in the previous section, naming xipT,yrpT, xipC,yrpC, and xipR,yrpR, such that random cost and revenue models in linear form are an innovation of this article.Step 2: Then, using the proposed model, we introduce the point distancex¯i,y¯r from each of the projection points of the BCC, cost, and revenue linearized random models (models (4)–(6)), which are represented by d1,d2,d3, respectively, and minimized under norm 2 in the objective function of the model (8) so that this point of the new projection is located in the convex shell of these points and at the same time overcome the unit under its evaluation. Therefore, the obtained model will be more efficient for each inefficient unit and also more efficient than the unit under evaluation.(11)Mind1+d2+d3∑i=1mx¯i−xipT2+∑r=1sy¯r−yrpT2=d1,j=1,…,n,∑i=1mx¯i−xipC2+∑r=1sy¯r−yrpC2=d2,j=1,…,n,∑i=1mx¯i−xipR2+∑r=1sy¯r−yrpR2=d3,j=1,…,n,λ1xipT+λ2xipC+λ3xipR=x¯i,i=1,…,mj=1,…,n,λ1yrpT+λ2yrpC+λ3yrpR=y¯r,r=1,…,sj=1,…,n,∑k=13λk=1,d1≥0,d2≥0,d3≥0x¯i≥0,y¯r≥0,λk≥0Considering that the expression under the radical of the first three constraints is greater than or equal to zero, therefore, the model can be rewritten as follows:(12)Mind1+d2+d3s.t.∑i=1mx¯i−xipT2+∑r=1sy¯r−yrpT2=d12,j=1,…,n,∑i=1mx¯i−xipC2+∑r=1sy¯r−yrpC2=d22,j=1,…,n,∑i=1mx¯i−xipR2+∑r=1sy¯r−yrpR2=d32,j=1,…,n,λ1xipT+λ2xipC+λ3xipR=x¯ii=1,…,mj=1,…,n,λ1yrpT+λ2yrpC+λ3yrpR=y¯r,r=1,…,sj=1,…,n,∑k=13λk=1,d1≥0,d2≥0,d3≥0,x¯i≥0,y¯r≥0,λk≥0.In the above model, we intend to minimize the total sum of the distance between projection points and ideal points from each of the projection points obtained from the BCC random, cost, and revenue models under norm 2 (L2) in the objective function of the proposed model so that the ideal point are in the convex combination of these points and at the same time dominate the unit under evaluation. This projection point obtained for each unit simultaneously meets the three technical, cost, and revenue objectives as much as possible. Secondly, random data and technical, cost, and revenue stochastic models are used in a linear form that has not been addressed so far, so obtaining such a pattern with these scores is very important.


## 5. Practical Example

The electricity industry has an infrastructural and influential role on the economic growth of any country. Iran Regional Electricity Company is one of the most important companies in the field of electricity industry, whose main task is to produce, distribute, and transmit electricity. In this practical example, fuel consumption and nominal power are used in the power generation process. The nominal power of a propulsion device has been written by the manufacturer on its specifications plate for certain conditions in terms of horsepower or megawatts. In small machines, the nominal power is specified in kilowatts. The power transmission section is done by posts and lines of the transmission network. The substation is part of a network, which has been concentrated in a given location and is used to selectively connect and disconnect electrical circuits within a network. Power distribution is also used to send and sell energy to companies. Iranian electricity distribution companies are responsible for managing and coordinating subsidiary units and supplying the production, transmission and sale of electricity and are subsidiaries of Tavanir specialized parent company.Considering the importance of this industry in the economic growth and development of the country, the model is evaluated based on data related to 16 regional electricity companies in Iran, in a 9-year period (from 2005 to 2014) derived from the statistical yearbook of the Ministry of Energy (https://www.moe.gov.ir) according to Jafari et al. [21]. Therefore, for each company, six inputs and two outputs have been considered, which are defined as follows. The selection of these indicators is based on a study of research conducted in the electricity industry.Input and output variables have been organized in Table2.Table 2
Input and output variables.
Input variablesOutput variables(1) Fuel consumption in terms of million liters(1) Sending energy to neighboring and overseas power companies at the transmission voltage level(2) Nominal power(3) Length of transmission lines(4) Capacity of transfer posts(5) The length of the above distribution lines(2) Selling energy to electricity distribution companies and industrial companies and the above distribution voltage level(6) Capacity of the above distribution postsInput and output indices have a normal distribution, the mean value and standard deviation for each of which have been shown in Tables3 and 4, respectively.Table 3
The mean value inputs and outputs.
DMUsElectricity distribution companiesThe mean value of the first inputThe mean value of the second inputThe mean value of the third inputThe mean value of the fourth inputThe mean value of the fifth inputThe mean value of the sixth inputThe mean value of the first outputThe mean value of the second output1Azerbaijan2921277637235808423445747420122112Esfahan3526220432008406554567006708207523Kermanshah2898236039757767556254467849138184Tehran7994768242142140046791386497513901005Khorasan3570353818834166840069132591164346Khuzestan29062286369112983387487167274281957Zanjan1773241170323123822989227066958Semnan109269100224018591066152726299Sistan and Baluchestan112710023716204629352077411446810The west1609148129584682436933305012618211Fars314434965321109827933834977291645612Kerman1756185843714967430234903399957713Guilan1628143311153267134821016110436814Mazandaran2992221224606449323645235894950115Hormozgan3126244421106555352240103678998316Yazd818930156125141659186818515294Table 4
Standard deviation of inputs.
DMUsElectricity distribution companiesDeviation of the first input criterionDeviation of the second input criterionDeviation of the third input criterionDeviation of the fourth input criterionDeviation of the fifth input criterionDeviation of the sixth input criterionDeviation of the first output criterionDeviation of the second output criterion1Azerbaijan1089.8831237.542215.05581070.038270.2961775.2903831.9791673.5982Esfahan1498.915705.8272319.76711188.98428.8438995.109933.74352313.4443Kermanshah322.66860289.07781296.611561.172693.5971678.10841340.5664Tehran2489.0942504.38983.294661720.807249.07231123.0761623.3693402.2315Khorasan881.8044908.9224431.5762771.0019502.0667926.43561356.3911607.0836Khuzestan297.5013217.1843366.6322131.656276.44891200.9062192.8672713.7547Zanjan237.1139313.7117335.8794695.6687375.4917513.3965681.8732918.47268Semnan134.3652264.0417280.8594654.7603106.6865160.7389812.4838315.3959Sistan and Baluchestan196.5782197.9394496.7202603.1526300.5495327.6523185.0081685.650110The west371.9556476.3024396.28651144.02255.3703425.39391682.792763.474311Fars1379.9461118.915634.97562660.342739.181423.8414040.3262041.42312Kerman224.6264256.3416285.8916713.4466312.9601314.8651626.46471252.8113Guilan632.6855562.7966150.9967497.6073109.0138286.4961697.3034606.642414Mazandaran251.61886.6332553.64699924.1342529.2353893.8775472.11651412.19715Hormozgan118.1397216377.5991084.7041882.73824.2087667.46242207.23416Yazd185.1756195.287315.0286648.906873.67496224.0982903.47881074.818Table5 shows the efficiency number obtained from technical, stochastic cost, and revenue models at three error levels (α = 0.05, 0.3, and 0.5). According to the results obtained from this table, units 1, 3, 4, 6, 7, 8, 9, 10, 13, and 16 are technical efficiencies. Units 1, 3, 4, 8, and 13 are cost efficient. Revenue efficiency includes units 4, 7, 8, and 9. Therefore, among the companies, only two electricity distribution companies (units 4 and 8) in Tehran and Semnan provinces are efficient at all three levels of error in technical, stochastic cost, and revenue models.Table 5
Technical efficiency numbers, cost and random revenue in three error levels.
DMUsElectricity distribution companiesThe stochastic efficiency of companies∝=0.05∝=0.3∝=0.51AzerbaijanTE = 1TE = 1TE = 1CE = 1CE = 1CE = 1RE = 3.29RE = 3.49RE = 3.542EsfahanTE = 0.91TE = 0.88TE = 0.86CE = 0/62CE = 0.59CE = 0.57RE = 3.46RE = 3.55RE = 3.603KermanshahTE = 1TE = 1TE = 1CE = 1CE = 1CE = 1RE = 3.79RE = 4.34RE = 4.604TehranTE = 1TE = 1TE = 1CE = 1CE = 1CE = 1RE = 1RE = 1RE = 15KhorasanTE = 0.78TE = 0.77TE = 0.76CE = 0.39CE = 0.36CE = 0.35RE = 2.08RE = 2.13RE = 2.156KhuzestanTE = 1TE = 1TE = 1CE = 0.81CE = 0.70CE = 0.62RE = 2.49RE = 2.79RE = 2.957ZanjanTE = 1TE = 1TE = 1CE = 0.73CE = 0.70CE = 0.69RE = 1RE = 1RE = 18SemnanTE = 1TE = 1TE = 1CE = 1CE = 1CE = 1RE = 1RE = 1RE = 19Sistan and BaluchestanTE = 1TE = 1TE = 1CE = 0.55CE = 0.53CE = 0.52RE = 1RE = 1RE = 110The westTE = 0.94TE = 0.81TE = 0.76CE = 0.59CE0.55CE = 0.52RE = 3.74RE = 3.97RE = 4.711FarsTE = 1TE = 1TE = 0.90CE = 1CE = 0.89CE = 0.68RE = 5.3RE = 5.58RE = 5.7012KermanTE = 0.63TE = 0.61TE = 0.60CE = 0.44CE = 0.43CE = 0.42RE = 4.15RE = 4.24RE = 4.2713GuilanTE = 1TE = 1TE = 1CE = 1CE = 1CE = 1RE = 1RE = 1.21RE = 1.4014MazandaranTE = 0.80TE = 0.70TE = 0.65CE = 0.56CE = 0.53CE = 0.51RE = 4.84RE = 5.01RE = 5.0315HormozganTE = 0.53TE = 0.53TE = 0.52CE = 0.51CE = 0.45CE = 0.42RE = 5.72RE = 5.93RE = 6.116YazdTE = 1TE = 1TE = 1CE = 0.77CE = 0.75CE = 0.75RE = 1.02RE = 1.07RE = 1.09TE: technical efficiency, CE: cost efficiency, and RE: revenue efficiency.The projections obtained from the technical, stochastic cost, and revenue models can be seen at three levels of error in Tables6–8. For each company, the projection points of the technical, cost, and revenue models have been shown in three rows with TE, CE, and RE, respectively. The results have been set in separate tables for three error levels. These tables have been presented in the appendix.Also, the projection point of the proposed model for each of the electricity distribution companies with three levels of error has been shown in Tables9–11, which is attached to the article. According to the results obtained for most units, the projection point of the proposed dominant model of the unit under evaluation has been obtained. Dominance means having fewer inputs and more outputs than the unit under evaluation. For example, consider the electricity distribution company of Khorasan province (unit 5) at the error level of 0.05. The model point of the model presented in this paper shows that all inputs have decreased and all outputs have increased compared to the unit under evaluation of this unit. These results can also be generalized to units 10, 12, 14, and 15. The projection point obtained from the model presented in this paper for units 1, 3, 4, 7, 8, 9, 11, and 13 equals to the unit under evaluation is obtained at an error level of 0.05. These results can be analyzed for other error levels.In Table12, the ideal cost and revenue are compared with the observed cost and revenue at three levels of error for all companies. The results in this table show that for all companies in all three levels of error, the ideal cost is not more than the observed cost and the ideal revenue is not less than the observed revenue as per unit, which indicates an improvement in performance of all companies. For example, consider the Electricity Distribution Company of Khuzestan Province (DMU 6) at the error level of 0.05. Although in the proposed model, the second and fifth inputs have been increased, but the ideal cost of this unit is less than the observed cost. Also, despite the decrease in the first output of the proposed model for this province, its ideal income has been obtained more than its observed income. This analysis can be generalized to other units.Table 12
Comparing ideal cost and revenue with observed one.
DMUElectricity distribution companiesObserved costObserved revenue∝=0.05∝=0.3∝=0.5Ideal costIdeal revenueIdeal costIdeal revenueIdeal costIdeal revenue1Azerbaijan64869466826486946682648694668264869466822Esfahan825156162849525.0563197.1348718.9562461.347259.561873.593Kermanshah78727511837873151183787315118378731511834Tehran1636778094531636898094531636898094531636898094535Khorasan720174064125976.5742635.9126570.2841505.526915.941074.786Khuzestan1001577821286018.8781868.5282519.3282819.3164110.2782608.587Zanjan30236202003023620200302362020030236202008Semnan1780698391780698391780698391780698399Sistan and Baluchestan355471016935547101693554710169355471016910The west520682740031034.9929389.3628884.5528146.2328863.4127739.0311Fars110172560991101725609996631.43210432.875343.3456320.1612Kerman582372935124871.530382.825958.7929965.926578.3329753.8513Guilan293482706629348270662934827066293482706614Mazandaran601263668435456.7538660.0630563.1137062.4730816.236767.7515Hormozgan598323100025772.5132923.3425486.1531614.5225352.13100016Yazd26018555322159.9616913.1322411.5516627.3422567.8716509.86The input price vector and the output price vector for all units are given a fixed value as follows:(13)c1=2,c2=1,c3=3,c4=4,c5=3,c6=2,p1=3,p2=2.For example, the ideal cost and the observed cost along with the ideal revenue and the observed revenue for DMU11 at the error level of 0.03 are obtained from the following formula, which can be generalized to all units.(14)Observed Cost=c1x1+c2x2+c3x3+c4x4+c5x5+c6x6=2∗3144+3496∗1+5321∗3+10982∗4+7933∗3+8349∗2=110172,Ideal Cost=c1x1∗+c2x2∗+c3x3∗+c4x4∗+c5x5∗+c6x6∗=2∗4.00E+03∗2+3477.71∗1+4025.19∗3+10630.16∗4+5376.55∗3+7213.93∗2=966314,Observed Cost=p1y1+p2y2=7729∗3+16456∗2=56099,Ideal Cost=p1y1∗+p2y2∗=8248.45∗3+92843.7∗2=210433.


## 6. Conclusions

Decision-making units are divided into two categories: efficient and inefficient units. For each inefficient DMU, there is a projection point on the efficient frontier so that the inefficient unit moves towards the efficient frontier with different approaches and becomes efficient.In this paper, in the first step, by technical efficiency models, linear stochastic cost, and revenue under variable returns to scale, we identify projection points such as each DMU. Then, through the presented model, a projection point is obtained that the distance from each of the projection points obtained in the first step under norm 2 is the least. The projection point obtained for each inefficient DMU is very important in two ways. First, the projection obtained for each DMU simultaneously meets the three main aspects of technical efficient, minimal cost, and maximal revenue as much as possible. Secondly, random data and stochastic programming in linear form are used in it. Therefore, it will help managers to make better decisions in the future and will enjoy the benefits of random planning. Considering the importance of the electricity industry in the economic growth and development of each country, a practical example has been implemented on 16 regional electricity companies in Iran in 9 periods. Comparison between cost efficiency and ideal income with observed cost and revenue efficiency showed that the ideal cost of each unit is less than or equal to its observed cost and the ideal revenue of each unit is greater than or equal to the observed revenue of that unit at three levels of error. According to the obtained results, the managers of Iranian electricity distribution companies can plan by identifying inefficient companies, identifying the projection point of the units and setting goals in order to increase their efficiency before the actual performance of the companies. For further studies, the method presented in this paper can be applied to other types of uncertain data such as fuzzy or interval data and can also be focused on other aspects.



---

*Source: 9198737-2022-10-12.xml*