ELECTRE and AHP MCDM methods versus CP method and the official choice applied to high-speed railway layout alternative election JOSE M. ANTONa, JUAN B. GRAUb, DIEGO ANGUINAC Technical University of Madrid E.T.S. de Ing. Agrónomos. Av. Complutense s/n, Ciudad Universitaria. 28040 Madrid. SPAIN Abstract: - The application of different multicriteria methods for the election of alternatives for the layout of high-speed railway in one specific corridor is considered in this paper using the outranking MCDM methods: ELECTRE-I and AHP, in comparison with Compromise Programming (C.P.) methods. Starting from two reported cost-effectiveness studies for Madrid-Valencia link, the feasibility may be assessed, considering only a small set of criteria well elected and studied and containing relevant variables. Using outranking methods the solution “North” was preferable, in coincidence with the official solution. C.P. methods for these applications present several inconveniences. Thus, we propose outranking methods for such Decision-Making procedures. Key-Words: Transportation systems, Multicriteria Methods, Decision Aiding. TAV net was approved as a “Public-Works_20002007_Plan for TAV rail infrastructure”. The case in this paper refers to layout of the link from Madrid to Valencia at the Mediterranean coast at East, which is an essential and decisive part of the approved TAV SE Sub-net. The officially adopted solution communicates almost at best Valencia with Madrid, and includes the northern city and region of Cuenca in the area of the new facility. For that link MadridValencia the authors have taken for their academic case-study for this paper three basic alternatives for decision extracted from the six public official alternatives for the whole TAV SE Sub-net that were described in the document PS. These three alternatives will be designed as (North, Centre and South) in the following, and they were designed as local optimums because of topography and use of territory. There are, indexed by i or k: a. The North solution (i = 1) is shorter for MadridValencia, at 350 km at first, through Cuenca, which is also a main historical town. The cost is higher than for the other solutions because the terrain is a higher plateau. b. The Centre or central solution (i = 2) follows in the middle section the actual new MadridValencia expressway, is shorter with a length of 370 km roughly. c. The South alternative (i = 3) follows roughly the actual rail main line laid in the 19th century to communicate the whole South and Southwest. As a result of the flatter terrain and of the use of older portions of the present line it is cheaper 1. Introduction. The authors, involved in Decision Making in U.P.M. and having some older experience in real management, set first a prior soft academic casestudy using the methods ELECTRE and AHP applied to the problem of high-speed railway layout alternatives election. In that moment election of alternatives for different corridors are part of a general long term planning to upgrade the Spanish rail network.. Then authors examined the relations of different methods with the big and innovative real project, and present in this paper case their revised models and commentaries about their aptitude to represent a real case. 2. Alternative layouts and criteria for the TAV Madrid-Valencia. The real decision was adopted (see http://www.mfom.es) the 8th of January 2002 by the administrations concerned of the sector SE of Spain with centre in Madrid. The official decision procedure involved the elaboration of a technical study referred later as PS (from Public Survey, ref. [7]), with plans at scale 1/25000) that was exposed in a public information procedure. At present date (winter of 2003) there is a trend for high-speed trains in Spain coexistent with other diverse public work plans. The future of trains require to upgrade totally some main lines to passenger rail standards for speeds of 300-350 km/h, designed here as TAV from "Tren de Alta Velocidad" in Spanish, and a new 1 from Madrid to Valencia, although it has a much bigger length of 470 km. As said an ELECTRE method and an Expert Choice AHP method will be used for the present case-study, by experts using at best the previously referred information, at least to make acceptable trade-offs between variables that will take the form of weights. Other methods could also be considered, ( [2] is a seminal reference, [5] contains a review of MCDM methods). At the present academic level four aggregated criteria that seem to be reasonably independent and complete, are used and are represented by j-criteria variables as follows: 1. Cost (j = 1), is the value in million euros of the initial costs, it is a "more is worse” variable, as the total project budget from the survey PS at state level to have the operating TAV new line from present state. 2. Trip Duration (j = 2). That variable will be in minutes, as travelling time between MadridValencia, and is a "more is worse” variable, represents the important related benefits and is roughly proportional to the distance. 3. Population index (j = 3). That aggregated variable, of type "more is better", will include the population served by the train, as disposed to use it, that is habitants and their tendency to use the train, and with a lower weight the population which will have indirect benefits from the train as said. At present level, as the effects are of different nature, indirect, and not easily predictable, a simple figure as population is used. The North alternative will be given a somehow higher value because it would integrate a larger part of territory as said. 4. Environment index (j = 4). Some main environment effects of that new train will be positive and the proposed alternatives have already avoided some critical sensitive areas, as said. What remains are some unfavourable environmental impacts to populations or to natural areas, and what matters are the differences between alternatives that are not major. The “environmental impact” studies include methods that are often ELECTRE ([9],[11],[12]), which result in indexes obtained from a conventional Multicriteria aggregation, that are “more is worse” index. The natural environment is more sensitive at North to some degree, and that can be saw in the PS study, and that will result in growing indexes. The cost and trip duration can now be well predicted, and are of great relevance because the costs are high and different, and because the trip durations are short and also different. The trip duration may be reduced to 90mn with the TAV new project to compare to almost 4 hours actually by train or bus, and the cost will be of the order of 4.000 million euros. The population index is itself an aggregation, and could be decomposed, maybe in a population variable and in a surface variable. As it is, it must include how important is the region beneficiated by the new transport system, direct or indirectly, at long term. Note that the units for criteria are different, they will be compared specifically later and the intervention of experts will matter then. 3. Method in an ELECTRE format. The ELECTRE-I method from Roy, see Roy in [12] or in [14], also in [9] and in [4], has been applied using Microsoft Office-97 EXCEL, an ELECTRE software made from Roy may be purchased instead. The ELECTRE – IV, reported by Roy in [13] with an academic model case “ex post” for the enlargement of the METRO of Paris was found not appropriate for the case herein, because the goal is to choose one alternative instead of putting in order different ones, and because the variables are simple and because they do not correspond to what he calls “pseudocriteria“”. To start, for each i-alternative a value Iij has been assessed to each variable valuating the j-criterion or attribute, to get the Information Table or Initial Decision Matrix, Table 1, with a row of j being 1 if the j-criterion is “more is better” and –1 if it is “more is worse”, and a row of the values wj of the weights vector wT = ( wj ). Table 1, Information Table or Initial Decision Matrix I = ( Iij ). TRIP POTEN- ENVIRONDURATIAL MENTAL AlternaCOST TION USERS IMPACT tives: NORTH 4207,0847 85 3654286 3 CENTRE 3606,0726 90 3000000 2 SOUTH 3005,0605 125 3000000 1 wj 0,3 0,4 0,2 0,1 j index -1 -1 1 -1 The element Ii1 is the criteria variable Cost in millions of Spanish pesetas (1 euro = 166.386 pts exactly). The Ii2 is the Trip Duration in minutes. In the context of these project and ranges of variation of the corresponding variables, 1202.0242 and 40 respectively, the Trip Duration seems to the authors somewhat more important than the Cost, and so has 2 where got a weight somehow higher, w2 = 0.4 versus w1 = 0.3. That contains a comprehensive trade-off between, estimating that an increase of 1202.242/0.3 in Cost is to be compensated by a decrease in Trip Duration of 40/0.4 . The assignment of the relative values of the weights is essential for the results and it must be done by experts knowing the method and the real case, that is available prior information, documents as the PS [7] with “Demand prevision” and “Profitability Studies”, and moreover some subjective estimation about long term consequences, because the Costs are considerable, but the benefits of short Trip Duration are various ( e.g. [8] is a survey of social local effects, etc... ) and will last many years. The element Ii3 of the Potential users column corresponds to the “Population index” criteria as defined above in a wide sense, and so a higher value has been taken for the North solution, with about 20% increase, by subjective rational belief in fact. The Centre and South alternatives got a similar value 3000000 that represents the habitants in the corridor and in the Valencia extremity beneficiated by the solution. The weight w3 = 0.2 has been taken lower than for the other variables but not too much, because the variable includes long-term benefits for organisation and development of territory, etc…. The element Ii4 of the Environment impact column corresponds to a differential Environment impact index, and it has got a lower weight, w4 = 0.1 because as said in that case the difference among alternatives is not very relevant. The weight values ought to be normalised so that they add 1 and are somewhat imprecise, as a result the values in the tableau are adopted. A sensibility analysis for the weights wj values has proved that results do not change for small changes in weights. In this Information Table no alternative dominates any other for the four criteria, in a “Pareto-like” sense, and to go further ELECTRE method outranking relationships R are obtained. By definition the alternative Ei outranks the Ek , i. e. Ei R E k ,as meaning that “Ei is at least so good as Ek’”, i. if the concordance C(Ei ,Ek ) of Ei to Ek is greater than the concordance threshold, CT , ii. and if the discordance D(Ei ,Ek ) of Ei to Ek is not greater than the discordance threshold DT. First a Concordance Index Matrix C = (Cik) , Table 2, was calculated, as defined for i k as Cik = C(Ei ,Ek ) = jik w j , (1) jik is equal to 1 if j I ij I kj 0 , ½ if I ij I kj , 0 if j I ij I kj 0 . Table 2, Concordance Index Matrix C = (Cik) . NORTH CENTRE SOUTH NORTH 0,6 0,6 CENTRE 0,4 0,5 0,49 SOUTH 0,4 0,5 For CT a first value is to be taken as the average of these 6 elements, getting CT = 0.5. It is convenient to change it slightly to CT = 0.49 after a posterior sensitivity analysis, to get a better structure of results. Then a Standardized Decision Matrix S = (Sij) is obtained as S ij I ij SupI rj Inf I rj , r 1,2,3 r 1,2,3 (2) and to follow the columns of S , that have an equal range of variation equal to 1, are weighted to get the Standardized and Weighted Matrix T = (Tij), see Table 3, in the form Tij = Sij * wj . (3) Table 3. Standardized and Weighted Matrix T = (Tij). TRIP DURACOST TION NORTH 1,05 CENTRE 0,9 SOUTH 0,75 0,85 0,9 1,25 POTENTIAL USERS ENVIRONMENTAL IMPACT 1,117030167 0,917030167 0,917030167 0,15 0,1 0,05 Let the Discordance Index Matrix D = (Dik) be defined as, for (i,k = 1,2,3 ) and (i k) , Max SupTkj Tij j ,0 Dik = D(E ,Ek) = j 1, 2,3, 4 , (4) Max Tij Tkj j 1, 2 , 3, 4 giving the Table 4. Table 4. Discordance Index Matrix D = ( Dik ). NORTH CENTRE NORTH 0,75 CENTRE 1 SOUTH 1 1 SOUTH 0,750 0,428571429 For DT a first value is to be taken as the average of these 6 elements, getting DT = 0.8214328571. That value was conserved after a sensibility analysis, as the results were insensitive to small changes of it. j1, 2 , 3, 4 3 To obtain the final outranking relationships from precedent definition for (i,k = 1,2,3 ) and (i k): i. An Outranking Concordance Matrix U was obtained from the Concordance Index Matrix C with the rule: if Cik > CT = 0.50 0.49 (sensitivity analysis) then Uik = 1, else Uik = 0. ii. Then an Outranking Discordance Matrix W was obtained from the Discordance Index Matrix D with the rule: if Dik < DT = 0.8214328571 then Wik = 1, else Wik = 0. iii. Finally the Aggregated Matrix A ,Table 5, was obtained from U and D with the rule : Aik = Uik * Wik . (5) That was considered appropriated for the case, in view of the data and context, and the resultant graph of Figure 1 was adopted as ELECTRE result. It appears that the logic of judgement between criteria was not the same as in a IRR-NPV method, that many factors and data have intervened and that the decision was in the same direction as was taken by the group of real decision-makers. 4. Using Expert Choice for A.H.P. . The “Analytic Hierarchy Process” was originated as known by Saaty ( see [15] and [16]). It appears as a close school of decision, proclaimed in [17] as a distinct “theory of measurement” applied in “decision making”, as a “theory of decision different and independent from utility theory”. An Expert Choice Inc. Company decision support AHP software (http://www.expertchoice.com/) for PC has been used, obtained in a CD including sets of manuals rich on managerial procedures but hiding the exact mathematical formulas behind. Following [9] for AHP, an assumed decision problem, the goal, is structured as a Hierarchy by experts including alternatives and criteria that are called objectives; sub-criteria levels will be added to criteria if convenient, but that has not been made for the case studied herein. To rank the i-alternatives a Table 5. Aggregated Matrix A = ( A i k ). NORTH NORTH CENTER SOUTH CENTER 1 0 0 SOUTH 1 01 0 In ELECTRE-I as said the alternative Ei outranks the alternative Ek if and only if Aik = 1, and that means that “Ei is at least so good as Ek’”. Note that before the sensitivity analysis CT was 0.5 and hence A23 was 0 , and that afterwards the value 0.49 has been adopted CT obtaining the result A23 = 1 that is preferred. To have a graphical expression of the resulting partial ordering of alternatives an ELECTRE graph was drawn. In it an arc (i,k) was drawn if Ei outranks Ek . The non-dominated alternatives form the kern, in that case containing only the North alternative, which is to be considered more favourable than the others. The resulting graph for CT = 0.50 is in Figure1. set of m values vi adding 1, m i 1 i v 1 , is to be Figure 1.ELECTRE graph and kern obtained through an “Analytic Hierarchy Process” from successive analogous previous valuations quantifying the opinion of experts in view of the real problem, maybe involving successive meetings of experts. Starting from the inferior level, for each jth objective experts have to assess sets of weights wij adding 1, wij 1 , to valuate the ith alternative North respectively to the objective jth. Then to aggregate the j-objectives they have to assess a set of n values i or weights uj adding 1, n j 1 u j 1 . If there is only one level of objectives the final values are obtained wij u j or similarly, but if there were a as vi Centre j South “Hierarchy” of them, experts would have to made a “Process” of similar assessments. For the case study the weights wij were internally assessed for each jth A sensibility analysis was then made to check the results were sensitive to small changes of the values CT and DT that are somehow conventional. It was observed that a small decrease in CT from 0.5 to 0.49 changes U23 from 0 to 1, meaning that the alternative Centre outranks the South alternative. objective by means of pair-wise comparisons of the three i-alternatives made by the authors as experts. Then the weights uj were assessed as result of a similar process of pair-wise comparisons of the four j-objectives. 4 Table 7. Standardized and weighted distances from From the wij and u j , and using internal Ideal, pT 0,75 procedures in two modes, Ideal and Distributive, it obtained the Resulting estimated weights v i that follow in Table 6: COST NORTH CENTER SOUTH Table 6. AHP global weights, v i . Alternative vi NORTH CENTRE 0,484749 0,301688 SOUTH 0,213563 0,3 0,15 0 0,85 1,11703027 TRIP DURA- POTENTIAL TION USERS 0 0,05 0,4 0 0,2 0,2 0,05 ENVIRONMENTAL IMPACT 0,1 0,05 0 the ideal point, P = ( Pij ) . For the ith alternative an “h-disutility” is defined for h 1 as h-distance from the ideal point, These weight values are indexes of “priority” for the alternatives, the first one is clearly higher, and thus they point that the North alternative is the best, and that the Centre one is sensibly better than the South. The EC results contain several possible additional presentations, including 4 charts of sensitivity analysis, Performance, Dynamic, Gradient, Head-to-head. The first chart Performance of Figure 3 contains all the results. 1h h Dh,i Pij for 1 h < , j D,i Inf Pij for h = . (6) j To conclude the alternatives are ranked in opposite sense than the disutilities. For the present case-study the disutilities were calculated for h = 1, 1.5, 2 and are shown in the Table 8. Figure 3. Chart with EC results. Table 8. h-Distance (Minkowski) from the ideal Altrntv. D1 D2 D D1,5 North 0,4 0,31622777 0,3 0,33735053 Centre 0,45 0,25980762 0,2 0,30675555 South 0,6 0,4472136 0,4 0,48945216 point, P = ( Pij ) for distance Dh. For h = 1 then D1,i P ij , and the obtained j ranking for this case the same as in ELECTRE, as in EXPERT CHOICE and as in the official election, all that being satisfactory. In fact there is only a comparison of added rows, and it can be seen that the result is independent on the ideal point that does not represent a feasible solution, it depends on values of variables, weights, in a more primitive form that in ELECTRE. But for the other higher values of h the Centre solution is preferred without rational justification, and that result is unsatisfactory. The ideal point has then an influence that is not desirable in that case, it has the lower cost of the South solution and that enhances the influence of the cost criteria over trip duration and population ones. The introduction of other feasible non-optimal alternatives, as is the South alternative in fact, may thus change the ideal point and the results. In most cases in literature about CP the feasible set of solutions in the criteria variables space is convex and an intermediate alternative is chosen. It is clear that in the present 5. Comparison with a Compromise Programming method. Starting from the table precedent in chapter III “Standardised and weighted matrix” T and from the j , a Compromise Programming method CP (see [9], [10] or [3], or [1] for a public utility case with several methods) is applied now easily using the precedent tables in EXCEL. First an “ideal point or alternative” having the best rating for each criteria was calculated, resulting in the row pT, that has been subtracted from the rows of the matrix T so as to obtain the matrix P = (Pij) of the Table 7. 5 Valencia-Madrid TAV case-study the feasible set is not convex. That is at first because the number of solutions is discrete, and that is probably because the intermediate solutions are not feasible. For there is not a feasible intermediate layout between the North and the South solution, because between them are the great Alarcón dam reservoir in a river Júcar deep valley and a gentle rift that belongs to the separation between Mediterranean and Atlantic rivers. This non-convexity of sets of solutions happens often in Public Works because of natural restrictions. Selection in Madrid Metropolitan Area. Journal of the Operations Research Society of Japan, 46, March 2003. [2] Arrow K. J. Social choice and individual values (2nd edn). Wiley, New York, 1963. [3] Ballestero E. and Romero C. Multiple Criteria Decision Making and its Application to economic problems, Kluwer Academic Publishers, Boston 1998. [4] Barba-Romero S. & Pomerol J.Ch. Decisiones Multicriterio. Servicio de Publicaciones Universidad de Alcalá, Alcalá de Henares, 1997. [5] Gal T. & Stewart T..J. & Hanne T. Multicriteria Decision Making, Advances in MCDM Models, Algorithms, Theory and Applications. Kluwer Academic Publishers, Boston 1999. [6] Keeney R. L. and Raiffa H. Decisions with Multiple Objectives, Preferences and Value Tradeoffs. John Wiley & Sons, New York, 1976. [7] Ministerio de Fomento D. G. de Ferrocarriles 6. Comparison with the Official Choice, Discussion and Conclusions. The Official Choice, that in fact seems appropriated to the real problem, contains the North solution that was also preferred in the precedent methods herein. That was because of lower Trip Duration and higher long-term indirect benefits such as better future Organisation of Territory, as including the Cuenca corridor in the SE TAVSubnet, in spite of higher initial costs. That is satisfactory for the authors, but what does it means in view of practical use? One may think that the academic case-study of this paper completed ex post may be too strictly influenced by the official decision. The authors think that their result is honest, even having to incorporate the reality in a very simple format. It seems to the authors that Methods like ELECTRE I and Expert-Choice are then of use to take into account inhomogeneous factors, social opinion and the knowledge of different kinds of experts. But care must be taken to include the relevant decision the pertinent relevant criteria and to include at least as preliminary the results of some quantitative studies such as the Profitability studies. Only with good technical and decisional information will experts made useful subjective assessments of variables, weights or pair wise comparisons to use methods like ELECTRE or AHP. In the case studied herein the Cost and Trip Duration criteria were relevant, but also were more imprecise factors about the future of territory that were been incorporated and aggregated in a Potential Users variable. The differences on Environment Impact were here less relevant, but they may be decisive in other cases, and in fact the E.I. had a previous great influence on the details of each alternative. Estudio Informativo del Proyecto, línea de Alta Velocidad “Madrid, Castilla La_Mancha, Comunidad Valenciana, Region_de_Murcia”. Official Survey made in collaboration with S.E.N.E.R., Madrid, September 2000. [8] Ribalaygua C. B. et alt., Efectos Territoriales de la alta velocidad ferroviaria. Estrategias para el planeamiento supramunicipal. O.P. Ingeniería y Territorio, 60, Madrid 2002. [9] Romero C. Teoría de la decisión multicriteria: Conceptos,técnicas y aplicaciones. Alianza Universidad Textos, Madrid, 1993. [10] Romero C. and González-Pachón. Distancebased consensus methods: a goal programming approach. Omega Int. Mgmt. Sci., 1999, 27, 341-347. [11] Romero C. and González-Pachón. Aggregation of partial ordinal ranking: an interval goal programming approach. Computers & Operational Research, 2001, 28, 827-834. [12] Roy B., Bouyssou D. Aide Multicritère à la Décision: Méthodes et cas. Economica, Paris, 1993. [13] Roy B., Hugonnard J-Chr., Ranking of suburban line extension projects on the Paris metro system by multicriteria method. Transportation Research 1982, 16, A, 301-312. [14] Roy B. Méthodologie Multicritère d’Aide à la Décision. Economica, Paris, 1985. [15] Saaty T. The Analytic Hierarchy Process. MCGraw-Hill, New York, 1980-1988 [16] Saaty T. Multicriteria Decision Making: The Analytic Hierarchy Process (extended edition Vol. 1, AHP Series), RWS Publications, 1996. [17] Saaty T. Decision Making for Leaders (extended edition Vol. 2, AHP Series), RWS Publications, 1996. References: [1] Anton J., Ballestero E. and Bielza C.. Compromise Based-Approach to Road Project 6
© Copyright 2025 Paperzz