國科會經濟學門效率與生產力分析研習營

國科會經濟學門效率與生產力分析研習營 時間: 2008年5月29日(週四) 地點: 長榮大學國際企業學系主講人: 黃台心(國立政治大學金融系教授) 第二部份 DEA方法 (1): 技術效率之估計

Reference: Coelli, T., D.S.P. Rao, and G.E. Battese (1999), An Introduction to Efficiency and Productivity Analysis, Kluwer Academic Publishers. Chapter 6

Farrell (1957) suggested the use of a non-parametric piece-wise linear convex isoquant, constructed such that no observed point lies to the left or below it. See Figure 6.2.

S 0 Figure 6.2Piece-wise Linear Convex Unit Isoquant

The piece-wise-linear convex hull approach to frontier estimation, proposed by Farrell (1957), was considered by only a few authors in the two decades following Farrell’s paper. Boles (1966) and Afriat (1972) suggested mathematical programming methods which could achieve the task, but the method did not receive wide attention until the paper by Charnes, Cooper and Rhodes (1978, EJOR), in which the term data envelopment analysis (DEA) was first used. Since then there has been a large number of papers which have extended and applied the DEA methodology.

Charnes, Cooper and Rhodes (1978) proposed a model which had an input orientation and assumed constant returns to scale (CRS). Subsequent papers have considered alternative sets of assumptions, such as Banker, Charnes and Cooper (1984, Management Science), in which a variable returns to scale (VRS) model is proposed.

Figure 6.3Input- and Output-Orientated Technical Efficiency Measures and Returns to Scale (a) DRTS (b) CRTS y y f(x) D f(x) D B B A P A P 0 x 0 x C C

Notations:Assume there are data on K inputs and M outputs for each of N firms. For the i-th firm these are represented by the column vectors and , respectively. Theｋｘｎ matrix, X, and the output matrix, Y, represent the data for all N firms.

An intuitive way to introduce DEA is via the ratio form. For each firm, we would like to obtain a measure of the ratio of all outputs over all inputs, such as , where u is an vector of output weights and v is a vector of input weights. The optimal weights are obtained by solving the mathematical programming problem: st (6.8)

This involves finding values for u and v, such that the efficiency measure for the i-th firm is maximized, subject to the constraints that all efficiency measures must be less than or equal to one. One problem with this particular ratio formulation is that it has an infinite number of solutions.[10] To avoid this, one can impose the constraint , which provides: st (6.9) where the change of notation from u and v to and is used to stress that this is a different linear programming problem. The form in equation 6.9 is known as the multiplier form of the DEA linear programming problem. [10] That is, if is a solution, then is another solution, etc.

Using the duality in linear programming, one can derive an equivalent envelopment form of this problem: st (6.10) where is a scalar and is a vector of constants. This envelopment form involves fewer constraints than the multiplier form ( K+M <N+1 ), and hence is generally the preferred form to solve.

The value of obtained will be the efficiency score for the i-th firm. It will satisfy: , with a value of 1 indicating a point on the frontier and hence a technically efficient firm, according to the Farrell (1957) definition. Note that the linear programming problem must be solved N times, once for each firm in the sample. A value of is then obtained for each firm.

Slacks The piece-wise linear form of the non-parametric frontier in DEA can cause a few difficulties in efficiency measurement. The problem arises because of the sections of the piece-wise linear frontier which run parallel to the axes (refer the Figure 6.2) which do not occur in most parametric functions. To illustrate the problem, refer to Figure 6.5 where the firms using input combinations C and D are the two efficient firms which define the frontier, and firms A and B are inefficient firms.

The Farrell (1957) measure of technical efficiency gives the efficiency of firms A and B as and , respectively. However, it is questionable as to whether the point is an efficient point since one could reduce the amount of input used (by the amount ) and still produce the same output. This is known asinput slack in the literature. Once one considers a case involving more inputs and/or multiple outputs, the diagrams are no longer as simple, and the possibility of the related concept of output slack also occurs.

S A B C D 0 Figure 6.5 Efficiency Measurement and Input Slacks

DEA Calculations using the Computer DEA calculations can be conducted using a number of different computer programs. If you are familiar with linear programming (LP), then all you need is access to computer software that can conduct LP. For example, one could use spreadsheet software such as Excel or statistical software such as SAS or SHAZAM. (Lingo, Gams) There are also a few specialist DEA computer packages available, such as IDEAS, Frontier Analyst, WDEA and DEAP. We use DEAP Version 2.1 in this book. A brief description of the DEAP computer program is provided in the Appendix.

To calculate the above simple numerical example using DEAP the user must construct a data file and an instruction file. Note that all data, instruction and output files are text files. The data file for this example, EG1.DTA, (refer to Table 6.3a) contains five observations on one output and two inputs. The output quantities are listed in the first column and the inputs in the next two columns. This data is identical to that listed in Table 6.1.

The instruction file, EG1.INS, is listed in Table 6.3b. The purpose of the majority of entries in the file should be self explanatory, due to the comments on the right-hand side of the file. The first two lines of the file contain the name of the data file (EG1.DTA) and an output file name (here we have used EG1.OUT).

Then on the next four lines we specify the number of firms (5); number of time periods (1); [17] number of outputs (1); and number of inputs (2). On the last three lines, we specify a ‘‘0’’ to indicate CRS; a ‘‘0’’ to indicate an input orientation; and a ‘‘0’’ to indicate that we wish to estimate a standard DEA model. [18] [17]Note that the number of time periods must be equal to 1 unless the Malmquist DEA option is selected. [18] Note that by specifying ‘‘0’’ on the final line we are asking that slacks be calculated using the multi-stage method. If we wished the 1-stage or 2-stage methods to be used we would have used a ‘‘3’’ or a ‘‘4’’, respectively. These different methods for the calculation of slacks are discussed further in the following chapter.

Table 6.3a Listing of Data File, EG1.DTA 1 2 5 2 2 4 3 6 6 1 3 2 2 6 2 Finally, we execute DEAP and type in the name of the instruction file (EG1.INS). The program sends the output to the file (EG1.OUT). This file is reproduced in Table 6.3c. The information presented in this output file should be self explanatory given the preceding discussion. Note that the results are identical to those presented in Table 6.2.

6 1 5 4 3 2 3 4 2 FRONTIER 5 1 0 1 2 3 4 5 6 Figure 6.6 CRS Input-Orientated DEA Example

Table 6.3b Listing of Instruction File, EG1.INS eg1.dta DATA FILE NAME eg1.out OUTPUT FILE NAME 5 NUMBER OF FIRMS 1 NUMBER OF TIME PERIODS 1 NUMBER OF OUTPUTS 2 NUMBER OF INPUTS 0 0=INPUT AND 1=OUTPUT ORIENTATED 0 0=CRS AND 1=VRS 0 0=DEA (MULTI-STAGE), 1=COST-DEA, 2=MALMQUIST-DEA, 3=DEA (1-STAGE), 4=DEA (2-STAGE)

Table 6.3c Listing of Output File, EG1.OUT Results from DEAP Version 2.1 Instruction file = eg1.ins Data file = eg1.dta Input orientated DEA Scale assumption: CRS Slacks calculated using multi-stage method EFFICIENCY SUMMARY: firm te 1 0.500 2 1.000 3 0.833 4 0.714 5 1.000 mean 0.810

Table 6.3c Listing of Output File, EG1.OUT (continuous) SUMMARY OF OUTPUT SLACKS: firm output: 1 1 0.000 2 0.000 3 0.000 4 0.000 5 0.000 mean 0.000 SUMMARY OF INPUT SLACKS: firm input: 1 2 1 0.000 0.500 2 0.000 0.000 3 0.000 0.000 4 0.000 0.000 5 0.000 0.000 mean 0.000 0.100

Table 6.3c Listing of Output File, EG1.OUT (continuous) SUMMARY OF PEERS: firm peers: 1 2 2 2 3 5 2 4 5 2 5 5 SUMMARY OF PEER WEIGHTS: (in same order as above) firm peer weights: 1 0.500 2 1.000 3 0.500 1.000 4 0.286 0.214 5 1.000

Table 6.3c Listing of Output File, EG1.OUT (continuous) PEER COUNT SUMMARY: (i.e., no. times each firm is a peer for another) firm peer count: 1 0 2 3 3 0 4 0 5 2 SUMMARY OF OUTPUT TARGETS: firm output: 1 1 1.000 2 2.000 3 3.000 4 1.000 5 2.000

Table 6.3c Listing of Output File, EG1.OUT (continuous) SUMMARY OF INPUT TARGETS: firm input: 1 2 1 1.000 2.000 2 2.000 4.000 3 5.000 5.000 4 2.143 1.429 5 6.000 2.000

6.4 The Variable Returns to Scale (VRS) Model and Scale Efficiencies The CRS assumption is only appropriate when all firms are operating at an optimal scale. Imperfect competition, constraints on finance, etc., may cause a firm to be not operating at optimal scale. Banker, Charnes and Cooper (1984) suggested an extension of the CRS DEA model to account for variable returns to scale (VRS) situations. The use of the CRS specification when not all firms are operating at the optimal scale, results in measures of TE which are confounded by scale efficiencies (SE). The use of the VRS specification permits the calculation of TE devoid of these SE effects.

The CRS linear programming problem can be easily modified to account for VRS by adding the convexity constraint: to equation 6.10 to provide : st (6.12) where N1 is an vector of ones. This approach forms a convex hull of intersecting planes which envelope the data points more tightly than the CRS conical hull and thus provides technical efficiency scores which are greater than or equal to those obtained using the CRS model. The VRS specification has been the most commonly used specification in the 1990s.

Note that the convexity constraint ( ), essentially ensures that an inefficient firm is only “benchmarked” against firms of a similar size. That is, the projected point (for that firm) on the DEA frontier will be a convex combination of observed firms. This convexity restriction is not imposed in the CRS case. Hence, in a CRS DEA, a firm may be benchmarked against firms which are substantially larger (smaller) than it. In this instance the -weights will sum to a value greater than (less than) one.

Calculation of Scale Efficiencies Given that one believes that the technology is VRS, one may then obtain a scale efficiency measure for each firm. This is done by conducting both a CRS and a VRS DEA. One then decomposes the TE scores obtained from the CRS DEA into two components, one due to scale inefficiency and one due to “pure” technical inefficiency. If there is a difference in the CRS and VRS TE scores for a particular firm, then this indicates that the firm has scale inefficiency, and that the scale inefficiency can be calculated from the difference between the VRS and CRS TE scores.

In Figure 6.7, we illustrate scale inefficiency using a one-input, one-output example. The CRS and VRS DEA frontiers are indicated in the figure. Under CRS, the input-orientated technical inefficiency of the point P is the distance . However, under VRS, the technical inefficiency would only be . The difference between these two TE measures, , is due to scale inefficiency.

These concepts can be expressed in ratio efficiency measures as: where all of these measures are bounded by zero and one. We also note that because

Thus, the CRS technical efficiency measure is decomposed into “pure” technical efficiency and scale efficiency. This scale efficiency measure can be roughly interpreted as the ratio of the average product of a firm operating at the point to the average product of the point operating at a point of (technically) optimal scale (point R).

One shortcoming of this measure of scale efficiency is that the value does not indicate whether the firm is operating in an area of increasing or decreasing returns to scale. This latter issue can be determined by running an additional DEA problem with non-increasing returns to scale (NIRS) imposed. This is done by altering the DEA model in equation 6.12 by substituting the restriction with , to provide: st (6.13)

CRS Frontier NIRS Frontier y Q R A P VRS Frontier 0 x Figure 6.7 Calculation of Scale Economies in DEA

The NIRS DEA frontier is also plotted in Figure 6.7. The nature of the scale inefficiencies (i.e., due to increasing or decreasing returns to scale) for a particular firm can be determined by seeing whether the NIRS TE score is equal to the VRS TE score. If they are unequal (as is the case for the point P in Figure 6.7) then increasing returns to scale exist for that firm. If they are equal (as is the case for the point Q in Figure 6.7) then decreasing returns to scale apply. [19] [19]Note that if TECRS= TEVRS, then by definition, the firm is operating under CRS.

Examples of this approach applied to international airlines and electricity are provided in Bureau of Industry Economics (1994) and Färe, Grosskopf and Logan (1985, J. of Public Economics), respectively. Note that the constraint: , ensures that the i-th firm will not be ‘‘benchmarked’’ against firms which are substantially larger than it, but may be compared with firms smaller than it.

We now illustrate these calculations using DEAP computer program. The data file for this example, EG2.DTA (refer to Table 6.6a), contains five observations on one output and one input. The output quantities are listed in the first column and the inputs in the second column. These data are identical to those listed in Table 6.4.

The DEAP instruction file, EG2.INS, is listed in Table 6.6b. The only changes relative to EG1.INS are that: ◆ the input and output file names are different; ◆ the number of inputs is reduced to 1; and ◆ there is a ‘‘1’’ entered on the 2nd last line to indicate that VRS is required. The output file,EG1.OUT, is reproduced in Table 6.6c. These results are identical to those presented in Table 6.5. Note that when the VRS option is specified, the DEAP program conducts VRS, CRS and NIRS DEA and calculates scale efficiencies as well as technical efficiencies.

Table 6.6a Listing of Data File, EG2.DTA • 1 2 • 2 4 • 3 • 5 • 6

Table 6.6b Listing of Instruction File, EG2.INS eg2.dta DATA FILE NAME eg2.out OUTPUT FILE NAME 5 NUMBER OF FIRMS 1 NUMBER OF TIME PERIODS 1 NUMBER OF OUTPUTS 1 NUMBER OF INPUTS 0 0=INPUT AND 1=OUTPUT ORIENTATED 1 0=CRS AND 1=VRS 0 0=DEA (MULTI-STAGE), 1=COST-DEA, 2=MALMQUIST-DEA, 3=DEA (1-STAGE), 4=DEA(2-STAGE)

Table 6.6c Listing of Output File, EG2.OUT Results from DEAP Version 2.1 Instruction file = eg2.ins Data file = eg2.dta Input orientated DEA Scale assumption: VRS Slacks calculated using multi-stage method EFFICIENCY SUMMARY: firm crste vrste scale 1 0.500 1.000 0.500 irs 2 0.500 0.625 0.800 irs 3 1.000 1.000 1.000 - 4 0.800 0.900 0.889 drs 5 0.833 1.000 0.833 drs mean 0.727 0.905 0.804

Table 6.6c Listing of Output File, EG2.OUT (continuous) Note: crste = technical efficiency from CRS DEA vrste = technical efficiency from VRS DEA scale = scale efficiency = crste / vrste Note also that all subsequent tables refer to VRS results

Table 6.6c Listing of Output File, EG2.OUT (continuous) SUMMARY OF OUTPUT SLACKS: firm Output: 1 1 0.000 2 0.000 3 0.000 4 0.000 5 0.000 mean 0.000 SUMMARY OF INPUT SLACKS: firm input: 1 1 0.000 2 0.000 3 0.000 4 0.000 5 0.000 mean 0.000

Table 6.6c Listing of Output File, EG2.OUT (continuous) SUMMARY OF PEERS: firm peers: 1 1 2 1 3 3 3 4 3 5 5 5 SUMMARY OF PEER WEIGHTS: (in same order as above) firm peer weights: 1 1.000 2 0.500 0.500 3 1.000 4 0.500 0.500 5 1.000

Table 6.6c Listing of Output File, EG2.OUT (continuous) PEER COUNT SUMMARY: (i.e., no. times each firm is a peer for another) firm peer count: 1 1 2 0 3 2 4 0 5 1 SUMMARY OF OUTPUT TARGETS: firm output: 1 1 1.000 2 2.000 3 3.000 4 4.000 5 5.000

Table 6.6c Listing of Output File, EG2.OUT (continuous) SUMMARY OF INPUT TARGETS: firm input: 1 1 2.000 2 2.500 3 3.000 4 4.500 5 6.000

6.5 Input and Output Orientations In the preceding input-orientated models, discussed in Section 6.3 and 6.4, the method sought to identify technical inefficiency as a proportional reduction in input usage, with output levels held constant. This corresponds to Farrell’s input-based measure of technical inefficiency. As discussed in Section 6.2, it is also possible to measure technical inefficiency as a proportional increase in output production, with input levels held fixed. The two measures provide the same value under CRS but are unequal when VRS is assumed (see Figure 6.3). Given that linear programming does not suffer from such statistical problems as simultaneous equation bias, the choice of an appropriate orientation is not as crucial as it is in the case of econometric estimation.

In a number of studies, analysts have tended to select input-orientated models because many firms have particular orders to fill (e.g., as in electricity generation) and hence the input quantities appear to be the primary decision variables, although this argument may not be as strong in all industries. In some industries, the firms may be given a fixed quantity of resources and asked to produce as much output as possible. In this case, an output orientation would be more appropriate. Essentially, one should select the orientation according to which quantities (inputs or outputs) the managers have most control over. Furthermore, in many instances, the choice of orientation has only a minor influence upon the scores obtained.

國科會經濟學門 效率與生產力分析研習營