# A Data Envelopment-Based Clustering Approach for Public Sugar Factories in Privatizing Process

Authors: Ezgi A. Demirtas
Journal: Mathematical Problems in Engineering
 (2011)

---

## Abstract

Turkish Sugar Inc., a public enterprise including 25 factories, is the first corporation of Turkish industry. According to the government policy, public sugar factories (PSFs) will be privatized as geography-based 6 portfolio groups in two years. As performance measures of PSF affect government, sugar producers, and several unions in privatizing process, a systematic approach is necessary to measure efficiencies and grouping factories. This paper uses a new DEA- (Data Envelopment Analysis-) based clustering approach for measuring efficiency scores of PSF and grouping them instead of geography- based portfolio groups. This new approach can help decision makers in privatizing process. At the same time, target values obtained by dual model can be used to eliminate inefficiencies of some PSFs.


---

## Body

## 1. Introduction

Sugar factories are the first corporations of Turkish industry. The first sugar factory was established by the direction of Kemal Ataturk in Alpullu in 1926. Annual sugar demand of Turkey which is supplied by three different sugar producers is 2.3 million ton. These producers are Turkish Sugar Inc., a public enterprise including 25 factories, Pankobirlik (beet producers union) with 6 factories, and starch-based sugar producers with 5 factories. The market share of these producers is 70%, 20%, and 10%, respectively. Turkish sugar Inc. and Pankobirlik use beet to produce sugar instead of starch.According to the government policy, sugar factories of Turkish Sugar Inc. will be privatized as geography-based 6 portfolio groups in two years: A (Kars, Ercis, Agri, Mus, and Erzurum), B (Elazig, Malatya, Erzincan, and Elbistan), C (Kastamonu, Kirsehir, Turhal, Yozgat, Corum, and Carsamba), D (Bor, Eregli, and Ilgin), and E (Usak, Alpullu, Burdur, and Afyon). Some Turkish and foreign corporations aspire to buy sugar factories like Pankobirlik, Keskinkilic, Torunlar (Turkish), Sudzucker (German), Saint Louys Sucre (French), and British Sugar.According to the privatizing supporters, average production season length is shorter in TR (64 days) than that in EU (120 days) and the average number of personnel per factory in TR and EU are 500 and 200, respectively. The number of average personnel per factory is too large in contrast with Europe. On the other hand, usage of noneconomic beets with low polar sugar value brings high costs. According to the opponents, the reasons of high prices are high starch-based sugar quota (15% in TR and 2% in EU) and outdated technology. If the starch based-sugar quota is decreased and new technology is used, PSFs become profitable. At the same time, to prevent external dependency and rural-urban migration and to protect sector, state of council should stop enforcement decision about privatizing. In the final situation, state of council stopped enforcement decision because of the following reasons: (i) the specification contains 5-year production obligation and 50-million-dollar assurance, (ii) it does not guarantee supply-demand balance and stability, and (iii) it does not guarantee production sustainability, and it creates external dependency.Because of the reasons mentioned above, performance measurement of PSF is an important task which affects government, sugar producers, and several unions in privatizing process. A systematic approach is necessary to measure efficiencies of PSF. This paper is the first real-life application of DEA-based clustering approach developed by Po et al. [1] for measuring efficiency scores of PSF and grouping them instead of geography-based portfolio groups. This new approach can help decision makers in privatizing process of PSF. At the same time, target values obtained by DEA model can be used to eliminate inefficiencies and to make inefficient factories profitable. Additionally, to the best of  the author’s knowledge, there is no scientific study for efficiency measurement of public sugar factories in our country or elsewhere.The rest of this paper is organized as follows: Section2 discusses DEA and DEA-based clustering approach which is developed by Po et al. [1]. In this section, the focus is why and how piecewise production functions drawn from DEA models are employed to cluster data. Section 3 illustrates the proposed DEA-based clustering approach for measuring efficiency scores of PSF and grouping them to help decision makers in privatizing process. The results obtained by DEA-based clustering approach are compared with geographic based portfolio groups and target values obtained by CCR model are given to eliminate the inefficiencies of some PSF. Finally, conclusions are stated in Section 4.


## 2. DEA-Based Clustering Approach

Conventionally, most clustering algorithms are procedures that minimize total dissimilarity; examples of such algorithms are given in the paper of Po et al. [1].A general clustering method is to findc cluster centers z1,z2,…,zc so that the total dissimilarity measure Js(z) with Js(z)=∑i=1c∑j=1naijf(d(xj,zi)) is minimized. Js(z) is usually defined as a distance-based function, and the problem here is to select a useful and reasonable distance measure d(xj,zi).On the other hand, the stated clustering approaches can be seen as a feature analysis technique. An assumption of the underlying feature analysis is to regard the feature itemsA1,A2,…,As as multiple features so that the minimization of Js(z) presents the closer of data among their features and makes it more possible for these DMUs to be classified into the same cluster. However, the clustering results derived from the minimization of the total feature dissimilarity Js(z) may not be helpful in some cases of clustering DMUs, especially in production units. In these cases, we use their production data to cluster them. Suppose that the production data have feature items A1,A2,…,Ak,Ak+1,…,As with A1 to Ak being input items and Ak+1 to As being output items. The clustering information obtained from the conventional clustering approaches can only reveal that DMUs are more similar to another one. However, the more important information we want to know is the production feature (functions) implied from the production data of all DMUs. That is, fj(A1,A2,…,Ak,Ak+1,Ak+2,…,As)=0. From these derived production functions, f1,f2,…, all DMUs are classified into different clusters (production functions). Therefore, each DMU knows not only the cluster that it belongs to but also knows the production function type that it confronts. Each DMU can compare its production feature with the other production functions so that the combination of its input resources or the combination of inputs and outputs can be readjusted. That is, for the case of data feature with input and output items, the cluster derived from production functions is more valuable than that derived from feature dissimilarity measures.The idea of Po et al.’s study [1] is to employ the production functions to cluster production data. The method supporting this idea is DEA, as initiated and developed by Charnes et al. [2]. The DEA is a data-oriented method for evaluating the relative efficiency of DMUs where each DMU is an entity responsible for converting multiple inputs into multiple outputs. Since the fundamental of DEA uses the nonparametric mathematical programming approach to estimate piecewise frontiers and envelop the DMU data sets, in this study, each piecewise frontier is regarded as one cluster of production functions. Therefore, we use all piecewise frontiers as a base to cluster production data. That is, they give up traditional clustering approaches of feature dissimilarity and propose a new approach by adopting the production functions revealed by the observation data to cluster all DMUs.DEA is a nonparametric method for the estimation of production frontiers. It is a useful tool for evaluating the relative efficiency for a group of DMUs. Up to now, DEA has been widely studied and applied in various areas for 30 years since Charnes et al. [2] first proposed the DEA method with the CCR model. Among them, the main forms of DEA models and their extensions include those of BCC model [3], the additive model, [4] and the imprecise DEA models [5, 6]. Modifications and extensions are the assurance region models [7, 8], superefficiency models [9, 10], cone ratio models [11, 12]. Stochastic and chance-constrained extensions are considered by some authors [13–17]. Taxonomy and general model frameworks for DEA can be found in [18, 19]. The CCR is the original model of DEA (see the M1 model) and is used in this study to explain the DEA-based clustering approach.The DEA model generalizes the usual input/output ratio measure of efficiency for a given unit in terms of a fractional linear program formulation. According to the economic notion of Pareto optimality, the DEA method states that a DMU is considered inefficient if some other DMUs or some combinations of other DMUs produce at least the same amount of output with less of the same resources input and not more of any other resources. Conversely, a DMU is considered Pareto efficient if the above is not possible. Suppose that there aren DMUs to be evaluated, xij is the noted amount of the ith input for the jth DMU and yrj is the noted amount of the rth output for the jth DMU. Output multipliers are u1,u2,…,us (one for each item of output) and input multipliers are v1,v2,…,vm (one for each item of input). The mathematical formulation of the method is summarized next, where the relative efficiency of the DMUk is determined [20]. See the M1 model.M1 Model:(2.1)Max⁡u,vEffk=∑r=1suryrk∑i=1mvixik,Subject to∑r=1suryrk∑i=1mvixik≤1,ur≥ε,vi≥ε∀r,i.
The DEA model is essentially a fractional programming problem with a ratio of a weighted sum of outputs to a weighted sum of inputs where the weights for both inputs and outputs are to be selected in a manner that calculates the efficiency of the evaluated unit. Therefore, the original form of the DEA model is both nonlinear and nonconvex problem. Charnes et al. [21] proved that fractional programming problem can be transformed into linear programming formulations. The first formulation is ‘‘input based,” constraining the weighted sum of outputs to be unity and minimizes the inputs that can then be obtained. The second formulation is “output based,” constraining the weighted sum of inputs to be unity and maximizes the outputs that can then be obtained (see the M2 model). Given constant returns to scale assumption, the result from the input-based model is the reciprocal of that from output-based model. If variable returns to scale are assumed, there is no direct relation can be found between these two models.For the clustering approach used in this study, the results can be different for those PSFs which are not on the production frontier according to the way that input-based or output based model is applied. The choice of using an input-based or output-based model depends on the production process characterizing the firm (i.e., minimize the use of inputs to produce a given output or maximize the output with given levels of inputs). The objective of this study is to find the set of coefficients associated with each output and input that will give the PSF being evaluated the highest possible efficiency by using the M2 model. Then, target values are calculated by using this model to eliminate the inefficiencies of some PSFs.M2 Model:(2.2)Max⁡u,vEffk=∑r=1suryrk,Subjectto∑r=1suryrj-∑i=1mvixij≤0j=1,…,n,∑i=1mvixik=1ur≥ε,vi≥ε∀r,i.DEA differs from the production theory of economics in that it is nonparametric. In economics, the production function is a function that summarizes the process of converting multiple inputs into a single output. Thus, a general mathematical form for the production function in economics can be expressed asy=f(x1,x2,x3,…,xn), where y is a quantity of output and x1,x2,x3,…,xn are quantities of inputs. However, DEA is a nonlinear programming model for evaluating a process converting multiple inputs into multiple outputs, that is, g(y1,y2,y3,…,yn)=f(x1,x2,x3,…,xn). Most previous studies had mentioned and discussed the properties of production function that are hidden in DEA methods [8–10, 14, 15, 17, 22].Since the number of DMUs is usually much larger than the number of inputs, we prefer to express the linear programming in its duality form. Further, the duality form can interpret the geometric meaning of DEA and provide information about conservation of resources or expansion of outputs to have DMUs from inefficiency to efficiency.IfEffk* is the optimal value of Effk, the DMUk is said to be efficient if and only if  Effk*=1. If Effk* is less than 1, DMUk is inefficient. According to the efficiency ratio, DMUs may be grouped as good (Effk*=1) and poor (Effk*<1) performers or clustered by assigning different efficiency ratio grades [23–27]. Although clustering by efficiency ratio gives some information about the rationality of output/input, it does not reveal the intrinsic relationship between the input and output production features. Therefore, this study adopts piecewise production functions derived from the DEA method to cluster data.In M2 model, it is obvious that the constraint∑r=1suryrj-∑i=1mvixij≤0 is an inequality formula of production functions. Solving M2 model yields the virtual multipliers ur* and vi*. Thus, ∑r=1sur*yrj-∑i=1mvi*xij=0 is derived. Running M2 model for k=1 to n gives all production functions. Then, all DMUs are classified into different clusters by these piecewise production functions. Thus, a clustering method using production functions via the DEA method is implemented. Po et al. [1] find that there is less consideration in using these production functions as a reference to classify evaluated DMUs, and they propose a clustering approach according to the properties of DEA and its production possibility set such that they can use these production functions as a reference to classify evaluated DMUs. The details about the algorithm used in which the DEA-based clustering method is applied are given in their paper.


## 3. DEA-Based Clustering of PSF

In this study, we have an efficiency evaluation problem with 25 PSFs (DMU), each PSF with three inputs and one output obtained by 2009-2010 annual activity reports. Actually processed beet quantity (PBQ), fuel consumption (FC), number of total personnel (TP), sugar production (SP), and molasses production (MP) data are placed in annual activity reports of PSF, and all of them are real and correct. PBQ, FC, and TP are considered as inputs. Only SP is selected as output because it is correlated with MP.The simplified production data of PSF are shown in Table1. This table shows the required quantity of inputs to produce one unit of (one metric ton) sugar. For example, PSF 22 uses 9.96 ton beet, 0.419 ton fuel, and 0.0213 personnel to produce one ton sugar according to Table 1.Table 1
Production data of PSF.
PSF no.PSFPBQFCTPSP(ton)(ton)(ton)1Afyon6.6020.33840.005812Agri7.76040.34060.033513Alpullu7.60270.37080.021614Ankara7.43870.28820.021715Bor7.30620.30630.012416Burdur6.59250.31740.00817Carsamba9.61620.51490.03218Corum6.67250.26340.008219Elazig7.61710.35850.0312110Elbistan7.27890.34770.0207111Ercis6.67690.22890.0367112Eregli7.00980.37730.0042113Erzincan6.8070.26920.018114Erzurum6.16560.30180.018115Eskisehir6.82390.27630.0076116Ilgin6.80720.34840.0061117Kars7.3360.33410.0455118Kastamonu7.16470.33450.015119Kirsehir6.84180.36230.0095120Malatya8.20430.33560.019121Mus6.69490.25250.0172122Susurluk9.960.41910.0213123Turhal7.47560.36090.0104124Usak6.81360.3080.0201125Yozgat6.40540.28180.0121The objective is to find the set of coefficientsu’s associated with each output and v’s associated with each input that will give the PSF being evaluated the highest possible efficiency. By using the M2 model for each PSF its efficiency ratio Effk and the solution of virtual multipliers, v1*,v2*,v3*,u1*are obtained. The multipliers v1*,v2*,v3*are measure of the relative increase in efficiency with each unitreduction of input value, whereu1* is a measure of the relative decrease in efficiency with each unit reduction of output value.The analytical results are shown in Table2.Table 2
Analytical results derived from M2 model for PSF.
PSF no.PSF Namev1*v2*v3*u1*Evaluated by the frontier of1Afyon0.1420.0945.3610.142x1+0.094x2+5.36x32***Agri0.0990.69100.8160.109x1+1.018x2+0.978x30.8163Alpullu0.1150.0764.3580.8130.142x1+0.094x2+5.36x34Ankara0.091.0491.2560.8880.101x1+1.181x2+1.415x35Bor0.0911.0541.2620.8920.101x1+1.181x2+1.415x36Burdur0.14070.0935.3120.9910.142x1+0.094x2+5.36x37*Carsamba0.104000.6410.142x1+0.094x2+5.36x30.6318Corum0.1190.6195.3610.119x1+0.619x2+5.36x39***Elazig0.0990.69200.8180.109x1+1.018x2+0.978x30.81610Elbistan0.0930.870.8360.8540.109x1+1.018x2+0.978x311**Ercis04.369010.101x1+1.181x2+1.415x3112*Eregli0.13905.360.9520.142x1+0.094x2+5.36x30.94913Erzincan0.0981.1421.3670.9670.101x1+1.181x2+1.415x314***Erzurum0.1210.847010.109x1+1.018x2+0.978x3115Eskisehir0.1160.6055.240.9780.119x1+0.619x2+5.36x316Ilgin0.1380.0915.1950.9690.142x1+0.094x2+5.36x317***Kars0.1030.72500.8560.109x1+1.018x2+0.978x30.84418Kastamonu0.1260.0834.7480.8860.142x1+0.094x2+5.36x319*Kirsehir0.13904.9810.9480.142x1+0.094x2+5.36x30.94720Malatya0.0810.941.1260.7960.101x1+1.181x2+1.415x321**Mus03.6594.42810.101x1+1.181x2+1.415x3122Susurluk0.0660.7690.9210.6510.101x1+1.181x2+1.415x323Turhal0.1230.0824.6570.8690.142x1+0.094x2+5.36x324Usak0.1010.9440.9070.9270.109x1+1.018x2+0.978x325Yozgat0.1420.0945.3610.142x1+0.094x2+5.36x3*They are reevaluated byy=0.142x1+0.094x2+5.36x3, because of zero virtual multipliers.**They are reevaluated byy=0.101x1+1.181x2+1.415x3, because of zero virtual multipliers.***They are reevaluated byy=0.109x1+1.018x2+0.978x3. because of zero virtual multipliers.By selecting the set of virtual multipliersv1*,v2*,v3*,u1*to be all nonzero, four frontiers of production functions (PFs) are found:PF(1)y=0.142x1+0.094x2+5.36x3,PF(2)y=0.101x1+1.181x2+1.415x3,PF(3)y=0.119x1+0.619x2+5.36x3,PF(4)y=0.109x1+1.018x2+0.978x3.PSFs with (*), (**), (***) in Table2 confront the degenerative frontier. Po et al. [1] suggest that they should be reclassified into the nearest effective frontier (the frontier with nonzero virtual multipliers). In this application, it is observed that PSFs with (*) confront the nearest effective frontier y=0.142x1+0.09x2+5.36x3, thus their efficiency ratio will be reevaluated by this frontier. However, in complicated applications (with more data items of input and output), it is impossible to judge the nearest effective frontier by observation. Hence, for PSF 7 (Carsamba), we follow the procedure of Po et al. [1], taking (x1,x2,x3,y)=(9.6162t(i),0.5149t(i),0.032t(i),1),   i=1,2,3,4 into (PF1), (PF2), (PF3), and (PF4), respectively. The t(i) value is calculated, giving t(1)=0.6307, t(2)=0.6155, t(3)=0.6118 and t(4)=0.6236. By taking the maximal value, the efficiency ratio for PSF 7 is re-evaluated as 0.6307. In addition, PSF 7 is classified into the cluster determined by the corresponding envelope y=0.142x1+0.09x2+5.36x3 (PF1).In this study, some PSFs achieve 100 percent efficiency and are referred to as the relatively efficient units, whereas other units with efficiency ratings of less than 100 percent are referred to as inefficient units. According to the results of Table2, there are six efficient (PSF1, PSF8, PSF11, PSF14, PSF21, and PSF25) and 19 inefficient PSFs. The 5 PSFs out of 19 have greater than or equal to 0.95 efficiency ratio. Additionally, the net revenues of PSF are supported by the DEA results. According to four different production functions, 25 PSFs are classified into four clusters. Clustering results are shown in Table 3.Table 3
DEA-based PSF clusters.
Cluster (PF) no.PSF statusPSF number(PF1)Efficient1, 25Inefficient3, 6, 7, 12, 15, 16, 18, 19, 23(PF2)Efficient11, 21Inefficient4, 5, 13, 20, 22(PF3)Efficient8Inefficient15(PF4)Efficient14Inefficient2, 9, 10, 17, 24By considering (PF1), (PF2), and (PF3), TP is the most critical input for the PSF in clusters 1, 2, and 3. The multiplier of TP has the biggest value for these clusters. The relative increase in efficiency is 5.36 with each number reduction of TP for the inefficient PSF placed in clusters 1 and 3. Similarly, the relative increase in efficiency is 1.415 for the inefficient PSF placed in cluster 2 (see Table 3). The order of multipliers for other inputs changes. For example, for the PSF in cluster 2, the multipliers of FC and TP are similar and higher than PBQ. On the other hand, FC is the most critical input for the PSF in cluster 4. The relative increase in efficiency is 1.018 with each ton reduction of FC for the inefficient PSF placed in cluster 4 (see Table 3). DEA and geography-based clustering results are compared in Table 4.Table 4
DEA-/Geography-based PSF clusters.
PSF no.PSFDEA-based cluster no.Geography-based cluster no.1Afyon1E2Agri4A3Alpullu1E4Ankara2*5Bor2D6Burdur1E7Carsamba1C8Corum3C9Elazig4B10Elbistan4B11Ercis2A12Eregli1D13Erzincan2B14Erzurum4A15Eskisehir3*16Ilgin1D17Kars4A18Kastamonu1C19Kirsehir1C20Malatya2B21Mus2A22Susurluk2*23Turhal1C24Usak4E25Yozgat1C*These PSFs are out of geography-based portfolio according to government policy.As you can see in Table4, DEA-based clusters contain different geography-based portfolio groups (1-E, C, D), (2-A, B, D), (3-C), and (4-A, B, E). Moreover, the clustering results derived from geography-based portfolio may not be helpful in cases of clustering PSF. From the derived production functions (PF1, PF2, PF3, and PF4), all PSFs are classified into four different clusters (production functions). Therefore, each PSF knows the PF type that it confronts. Additionally, each PSF can compare its production feature with the other production functions so that the combination of its input resources or the combination of inputs and outputs can be readjusted. That is, for the case of data feature with input and output items, the cluster derived from production functions is more valuable than that derived from geography-based portfolio groups. It is possible to eliminate inefficiencies by considering DEA-based clustering. For example, inefficient PSF placed in clusters (1), (2), and (3) should give priority to decrease the number of total personnel. In the same manner, inefficient PSF placed in cluster (4) should decrease fuel consumption at first. It is meaningful to support privatizing decisions by DEA-based clustering results than geography-based portfolio groups.At the same time, target values of inputs are calculated by using slack variables of M2 model and illustrated in Table5 for the inefficient PSF. Target values can help decision makers to eliminate the inefficiencies. For example, 9.96 ton beet, 0.419 ton fuel, and 0.0213 personnel are used to produce one ton sugar in PSF 22. When 6.477 ton beet, 0.273 ton fuel, and 0.0014 personnel are used to produce one ton sugar, PSF 22 becomes efficient.Table 5
Input target values of PSF.
Input targetsPSF no.PSF namePBQFCTP2Agri6.3330.2780.0243Alpullu6.1820.3020.0184Ankara6.6050.2560.01935Bor6.5180.2730.0116Burdur6.5160.3140.0097Carsamba6.6020.3380.00589Elazig6.2270.2930.02010Elbistan6.2080.2970.01812Eregli6.6020.3380.00613Erzincan6.5870.2610.01715Eskisehir6.6110.26780.00916Ilgin6.5920.3380.00617Kars6.2770.2860.02218Kastamonu6.3470.2960.01319Kirsehir6.4870.3290.00920Malatya6.5250.2670.01522Susurluk6.4770.2730.01423Turhal6.5030.3140.00924Usak6.3130.2850.019


## 4. Conclusions

This study develops a DEA-based clustering approach for the evaluation of PSF. The proposed approach employs the piecewise production functions derived from the DEA method to cluster the data with input and output items. Compared with geography-based clustering that only considers geographical location of PSF, our proposed approach reveals the input-output relationships hidden in the data items of input and output. Thus, for each evaluated PSF, we know not only the cluster that it belongs to but also the production function type that it confronts. It is very important for managerial decision making where decision makers are interested in knowing the changes required in combining input resources so that it can be reclassified into a different and desired cluster/class in privatizing process.The focus of this paper is to examine the CCR model of DEA and then establish the DEA-based clustering. Without loss of generality, while this approach has been carried out for the CCR model, the proposed approach can be easily extended to other DEA models. The clustering results drawn from the DEA-based clustering are unit invariant, meaning that they are not affected by the scale of data.The DEA-based clustering approach is suitable for most clustering problems, where there are inputs-and-outputs or cause-and-effect relationships between the features. For example, we use the proposed approach in the analysis of industry classification, sorting of PSF by input-output data.In summary, in view of the advantages of the DEA-based clustering approach, it is uniquely poised for clustering problems. We believe that future researches are necessary to unleash the full potential of this DEA-based clustering approach. It thus has tremendous potential to be used for various clustering problems. DEA-based clustering algorithm developed by Po et al. [1] is robust to a slight change in the input and output data sets, but not to outliers. Future researches will consider developing a robust-type DEA-based clustering algorithm.



---

*Source: 946549-2011-08-10.xml*