Next Article in Journal
Study on the Reinforcement Mechanism of High-Energy-Level Dynamic Compaction Based on FDM–DEM Coupling
Previous Article in Journal
Approximate Roots and Properties of Differential Equations for Degenerate q-Special Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiple-Criteria Heuristic Rating Estimation

by
Anna Kędzior
1 and
Konrad Kułakowski
2,*
1
Faculty of Applied Mathematics, AGH University of Krakow, al. Mickiewicza 30, 30-059 Kraków, Poland
2
Faculty of Electrical Engineering, Automatics, Computer Science and Biomedical Engineering, AGH University of Krakow, al. Mickiewicza 30, 30-059 Kraków, Poland
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(13), 2806; https://doi.org/10.3390/math11132806
Submission received: 25 May 2023 / Revised: 18 June 2023 / Accepted: 20 June 2023 / Published: 22 June 2023

Abstract

:
One of the most widespread multi-criteria decision-making methods is the Analytic Hierarchy Process (AHP). AHP successfully combines the pairwise comparisons method and the hierarchical approach. It allows the decision-maker to set priorities for all ranked alternatives. But what if, for some of them, their ranking value is known (e.g., it can be determined differently)? The Heuristic Rating Estimation (HRE) method proposed in 2014 tried to bring the answer to this question. However, the considerations were limited to a model only considering a few criteria. This work analyzes how HRE can be used as part of the AHP hierarchical framework. The theoretical considerations are accompanied by illustrative examples showing HRE as a multiple-criteria decision-making method.

1. Introduction

We constantly face the need to make decisions, which inevitably involve comparisons. The easiest way is to compare one thing to another—objects in a pair. We are accustomed to making such comparisons daily. Take the market, for example, where people try to choose the heavier box of cookies (or the lighter bicycle in a bike store). By comparing the weight of one object in our hand with the weight of another, we can easily estimate the relative “importance” of two similar-looking things. Often we can make the right choice without relying on external tools to indicate weight. In practice, we often have to deal with much more complicated comparisons than assessing the usefulness of products during shopping. There is usually no natural way to make an accurate comparison, and the object does not have a pre-assigned “weight”. Moreover, many different factors often need to be evaluated and compared. The pairwise comparison (PC) method comes to the rescue in such a case. It allows us to use our innate ability to compare alternatives, enabling us to deal with these challenges more effectively. When creating a decision model based on pairwise comparisons of alternatives, one of the problems we face is determining the nature of the objects under study. Some may be tangible, i.e., they can be measured and weighed; others are intangible and difficult to decide on by direct measurement [1]. For example, available methods such as AHP cannot deal with this problem effectively. Hence, in AHP, we prefer to deal with relative expert-based assessments even if we consider distances between cities [2].
The Heuristic Rating Estimation (HRE) method is an attempt to deal with this problem. The HRE was initially defined as two alternative methods (additive and multiplicative) for calculating the ranking [3,4]. In this article, we extend this model to include a hierarchy that allows multiple criteria for evaluating alternatives to be considered. The HRE calculation procedure is based on the same concept as EVM (Eigenvalue Method) [5] and GMM (Geometric Mean Method) [6]. That is, the priority of a given alternative is equal to the weighted average value of the priorities of the other alternatives [7]. This compatibility allows the original HRE idea to be easily combined with these methods within one hierarchical model (Example 2). Thus, naturally, HRE becomes a multi-criteria method that complements AHP in cases where some alternatives have measurable or previously known ratings.
The presented work is a follow-up to [3,4]. The preliminary section introduces the reader to the issues related to pairwise comparisons and priority-deriving methods, such as EVM and GMM. The HRE calculation procedure is described in (Section 4). The next section (Section 5) fits HRE into a multi-criteria hierarchical model used in the context of AHP. The theoretical description is supplemented by practical examples showing the applications of HRE (Example 1) and its integration with AHP (Example 2).
A list of the important symbols used in the work can be found in Table 1 below.

2. Related Works

The pairwise comparisons (PC) method has a centuries-old history. The first considerations of its use come from a thirteenth-century treatise by Ramon Lull, Ars Electionis [8,9]. Later, it was the subject of research by Condorcet, Fechner, Thurston and others [10,11,12,13]. Currently, pairwise comparisons are a valuable source of preferential information in many decision-making methods, including TOPSIS [14], PROMETHEE [15], MACBETH [16], BWM [17] or multiple-criteria sorting [18]. One of the best known methods of multiple-criteria decision making based on pairwise comparisons is the Analytic Hierarchy Process (AHP) [7,19]. The method defined by Saaty in 1977 [5] quickly gained popularity and became widely used in many areas of application, including the performance assessment of employees [20], country competitiveness [21], supplier segmentation [22], portfolio management [23], or sustainability [24].
This method competes with other decision-making methods in many areas, including logistics optimization [25,26], IT management [27,28], military applications [29] or such exotic applications as analyzing the possibility of oil exploration in the Arctic National Wild refuge [19] (pp. 89–100).
More application examples can be found in the reviews [30,31,32]. It also quickly gained opponents who very often rightly pointed out the imperfections of the method [33,34,35]. In addition to a range of applications, AHP has seen numerous extensions and modifications. In particular, various variants of the method have been developed using the representation of uncertain knowledge, including fuzzy sets [36,37,38], intervals [39], and gray numbers [40].

3. Preliminaries

3.1. Pairwise Comparisons

Let A = a 1 , , a n be a set of alternatives and C = [ c i j ] be a pairwise comparisons (PC) matrix, such that c i j R + for i , j { 1 , . . , n } . A single element c i j corresponds to a direct comparison of the i-th and j-th alternatives. A single entry c i j has a quantitative meaning. For example, when an expert decides that the a i alternative is twice more preferred than a j , then c i j takes the value 2. As a result, the elements of the diagonal are ones, i.e., C takes the form
C = 1 c 12 c 1 n c 21 1 c 2 n c n 1 c n 2 1 .
Based on the information stored in C, a priority vector is calculated. For the purpose of the article, we will denote it as [ w ( a 1 ) , , w ( a n ) ] T where w is the function w : A R + . We say that a i is more preferred than a j (denoted as a i a j ) if w ( a i ) > w ( a j ) . Similarly, a i is considered as equally preferred as a j (denoted as a i a j ) when w ( a i ) = w ( a j ) .
In the ideal case c i j = w ( a i ) / w ( a j ) , i.e., the results of comparisons provided by experts perfectly match the computed ranking. In practice, however, we may expect only that c i j w ( a i ) / w ( a j ) . The reason for this inequality is data inconsistency, which is most often the result of expert errors.
Definition 1.
A PC matrix C is called consistent if
c i k = c i j c j k , for every i , j , k = 1 , , n
Every PC matrix is said to be inconsistent if it is not a consistent PC matrix.
Definition 2.
A PC matrix C is called reciprocal if
c i j = 1 c j i for every i , j = 1 , . . . , n
Although the reciprocity condition is natural, the use of non-reciprocal matrices can also be found in the literature [7,41]. It is easily seen that every consistent PC matrix is also reciprocal, but not reversely.
One of the essential issues in the PC methods concerns the question of how the priority vector needs to be determined. The answer is not straightforward because the expert’s judgments may not be consistent (For a consistent matrices, each ranking method will lead to the same priority vector [7] (p. 96)). There are many methods of priority vector estimation in the literature [42]. Two of the most popular prioritization methods are the Eigenvector Method (EVM) [5] and the Geometric Mean Method (GMM) [6].
In EVM, the priorities of alternatives are determined by the real and positive solution of C w = λ max w max , where λ max is the principal eigenvalue (spectral radius) of C, and w max is an appropriate eigenvector. The vector of priorities w is assumed to be w max rescaled so that the sum of its components is 1.
In GMM, the priority of an individual a i is calculated as the geometric mean of the i-th row of C. The sum of components also rescales the created ranking vector. Both EVM and GMM have been repeatedly debated and analyzed, and both methods have their proponents and opponents [43,44].
Depending on how inconsistent the PC matrix is, calculating the priority vector reflects the real preferences of respondents. To determine the degree to which PC matrices are inconsistent, appropriate indices are used. The most popular one comes from Saaty [5] and is defined as:
CI = λ m a x n n 1 ,
where λ m a x denotes the principal eigenvalue of the n × n matrix C. Another interesting inconsistency indicator comes from Koczkodaj [45]. For the n × n PC matrix C, it is defined as:
K ( C ) = df max i , j , k { 1 , . . . , n } i j , j k , k i min 1 c i j c i k c k j , 1 c i k c k j c i j .
In general, the more inconsistent the matrix C, the larger the inconsistency indices. When the value of the inconsistency index is considerable, the ranking results are not credible. Thus, in most cases, it is essential to keep the inconsistency as small as possible. More about different inconsistency indices can be found in [46].

3.2. Analytic Hierarchy Process

A single PC matrix allows us to compare alternatives with respect to only one criterion at a time. In practice, however, there is often a need to take more than one criterion into account. For example, when buying a car, we do not just choose a “better” vehicle, but also pay attention to its price, terms of warranty, and so on. A framework that allows the use of pairwise comparisons as a multi-criteria method is provided by the Analytic Hierarchy Process (AHP). Saaty proposed the method in 1977 [5] following J. R. Miller idea published in 1966 [47]. using PC matrices hierarchy for processing a large number of criteria. AHP is able to handle multiple criteria by defining a hierarchy (a tree) of alternatives and criteria. Hence, at the bottom of the hierarchical tree defining the decision model, alternatives are compared with respect to the lowest level criteria, then those criteria are compared against each other with respect to the higher level of criteria, and so on. Although the AHP hierarchy can be freely expanded, the most popular models consist of three layers: alternatives A = { a 1 , . . . , a n } , criteria S = { s 1 , . . . , s m } and the root node of the hierarchy called “goal” or “purpose”.
As an example of a hierarchical AHP decision model, let us consider the problem of selecting a company manager. Suppose there are three candidates for this position: Andrew, Benjamin and Christopher. When considering these alternatives, we take into account their experience in managerial positions, education and interpersonal skills, as well as the ability to work under pressure. We compare the candidates in terms of these criteria. The one with the best results in our comparison will become the CEO of the company.
Going into details, according to the AHP method, we create a separate PC matrix for each criterion. Thus, C ( s ) is a PC matrix corresponding to the mutual comparisons of all three candidates with regard to the criterion s. There are four criteria for assessing candidates: experience ( ex ) , education ( ed ) , interpersonal skills ( is ) , and stress resistance ( sr ) . Hence, for example, C ( ex ) corresponds to the comparisons of our candidates in terms of their experience. We assume that the following four matrices were prepared as a result of the expert assessment.
C ( ex ) = 1 2 4 1 2 1 2 1 4 1 2 1 , C ( ed ) = 1 1 2 1 8 2 1 1 4 8 4 1 ,
C ( is ) = 1 1 2 3 2 1 6 1 3 1 6 1 , C ( sr ) = 1 1 2 1 2 1 2 1 1 2 1
For each considered matrix C ( s ) , we calculate, using EVM [5], the priority vector. These are
w ( ex ) = 4 7 , 2 7 , 1 7 T , w ( ed ) = 1 11 , 2 11 , 8 11 T ,
w ( is ) = 3 10 , 6 10 , 1 10 T , w ( sr ) = 1 4 , 2 4 , 1 4 T .
When calculating the final priority vector, we must take into account the fact that the criteria may contribute to a different extent to achieving the goal. For this reason, we compare them in pairs. These comparisons form one more 4 × 4 PC matrix C ^ = c ^ i j where c ^ i j denotes the result of individual comparisons between the criteria s i and s j , where s 1 , , s 4 mean ( ex ) , ( ed ) , ( is ) , and ( sr ) , correspondingly.
C ^ = 1 4 2 8 1 4 1 1 2 2 1 2 2 1 4 1 8 1 2 1 4 1 .
The priority vector w ^ for C ^ is
w ^ = 8 15 , 2 15 , 4 15 , 1 15 T .
Thus, according to the expert opinion, the first criterion; experience, has the priority 8 / 15 , education 2 / 15 , interpersonal skills 4/15, and stress resistance– 1 / 15 . The priorities of criteria determine the degree to which the comparisons of alternatives affect the overall score. Thus, they become weights scaling the results of direct comparisons between alternatives (the final priority vector is a linear combination of vectors w ( s ) where the scaling factors come from w ^ ). The overall ranking w is given as:
w = j = 1 m w ^ j w ( s j ) .
In our case, we have:
w = w ^ 1 w ( e x ) + w ^ 2 w ( e d ) + w ^ 3 w ( c h ) + w ^ 4 w ( a g ) =
= 8 15 4 7 , 2 7 , 1 7 T + 2 15 1 11 , 2 11 , 8 11 T + 4 15 3 10 , 6 10 , 1 10 T + 1 15 1 4 , 2 4 , 1 4 T
0.414 , 0.370 , 0.216 T .
The above priority vector translates to the observation that Andrew is slightly more preferred to be a company manager than Benjamin, and both are more preferred than Christopher.

4. Heuristic Rating Estimation

When comparing alternatives, sometimes, we already know priorities for some of them. In such a case, we can use them and not ask experts for unnecessary comparisons, saving time and money. This observation gave rise to the Heuristic Rating Estimation (HRE) method [3,4]. In HRE, the set of alternatives A is composed of two disjoint subsets: A K -alternatives for which the final priorities are known, and A U -alternatives for which priority weights need to be estimated (also referred as to unknown alternatives). Knowledge about the priority values of elements in A K combined together with pairwise comparisons of all the elements from A = A K A U allow us to calculate priorities for A U . Hence, A K can be considered as the set of references that provides a benchmark for new alternatives in A U . For the sake of simplicity, in this study we will denote A U = a 1 , . . . , a k and A K = a k + 1 , . . . , a n . If a i , a j A K , then we will assume that w ( a i ) and w ( a j ) are known. Therefore, the value c i j is also known and equals c i j = w ( a i ) / w ( a j ) , so querying experts about c i j is also not necessary. The comparisons matrix C in HRE takes the form:
C = 1 c 1 , k c 1 k + 1 c 1 , n c k , 1 1 c k , k + 1 c k , n c k + 1 , 1 c k + 1 , k 1 w ( a k + 1 ) / w ( a n ) c n , 1 c n , k w ( a n ) / w ( a k + 1 ) 1 .
The above matrix and vector of reference values w ( a k + 1 ) , . . . , w ( a n ) T allow us to calculate the complete ranking for A (Figure 1).
The HRE method has been designed as an extension of AHP, where part of the ranking is already known. Therefore, it is based on the same premise that
w ( a i ) c i j w ( a j ) .
Thus, as in AHP, we can request that the priority of the i-th alternative w ( a i ) is a weighted average of priorities of all the other alternatives [7], i.e.,
w ( a i ) = 1 n 1 j = 1 j i n c i j w ( a j ) , for i = 1 , , k .
This leads to a matrix equation
M w = b ,
where M is given as
M = 1 1 n 1 c 1 , 2 1 n 1 c 1 , k 1 n 1 c 2 , 1 1 n 1 c 2 , k 1 n 1 c k , 1 1 n 1 c k , 2 1 ,
the constant term vector,
b = 1 n 1 c 1 , k + 1 w ( a k + 1 ) + . . . + 1 n 1 c 1 , n w ( a n ) 1 n 1 c 2 , k + 1 w ( a k + 1 ) + . . . + 1 n 1 c 2 , n w ( a n ) 1 n 1 c k , k + 1 w ( a k + 1 ) + . . . + 1 n 1 c k , n w ( a n )
and the ranking vector for unknown alternatives is given as:
w = w ( a 1 ) , . . . , w ( a k ) T .
The solution of (7)-if it exists and is acceptable, forms the desired vector of priorities. It has been shown [48] that for a small value of K ( C ) and sufficiently large set of the reference alternatives A K the feasible solution of (7) always exists.
In HRE, as in AHP, many methods for calculating the ranking can be defined. The GMM equivalent in AHP is geometric HRE [4]. Thus, starting from (5), we may request
w ( a i ) = j = 1 j i n c i j w ( a j ) 1 n 1 for i = 1 , , n .
After raising both sides of (8) to the power n 1 , we obtain the non-linear equation system
w n 1 ( a 1 ) = c 1 , 2 w ( a 2 ) · c 1 , 3 w ( a 3 ) · . . . . . . . . . . · c 1 , n w ( a n ) w n 1 ( a 2 ) = c 2 , 1 w ( a 1 ) · c 2 , 3 w ( a 3 ) · . . . . . . . . . . · c 2 , n w ( a n ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . w n 1 ( a k ) = c k , 1 w ( a 1 ) · c k , 2 w ( a 2 ) · . . . · c k , n 1 w ( a n 1 ) .
Thanks to the logarithmic transformation, the above is equivalent to the linear equation system
N μ = d ,
where N is a k × k auxiliary matrix in the form
N = ( n 1 ) 1 1 1 1 1 1 ( n 1 ) ,
d is a constant term vector
d = lg c 1 , 2 c 1 , 3 · . . . · c 1 , k c 1 , k + 1 w ( a k + 1 ) · . . . · c 1 , n w ( a n ) lg c 2 , 1 c 2 , 3 · . . . · c 2 , k · c 2 , k + 1 w ( a k + 1 ) · . . . · c 2 , n w ( a n ) lg c k , 1 c k , 2 · . . . · c k , k 1 · c k , k + 1 w ( a k + 1 ) · . . . · c k , n w ( a n ) ,
and
μ = μ 1 μ k = lg w ( a 1 ) lg w ( a k ) .
The solution μ of (10) induces the weight vector w—a solution of the original non-linear problem (9)
w = [ exp μ 1 , . . . , exp μ k ] T .
The feasible solution of (10) always exists and is optimal [4].

5. Multiple-Criteria Heuristic Rating Estimation Method

Sometimes when we want to compare a few objects or select several items from many others, we follow different criteria with varying degrees of importance. Moreover, we may have preliminary information about specific things. For example, when choosing a school for a child, parents pay attention to the level of education and the percentage of pupils entering the next school, the educational system, location, social group from which the pupils come, or additional activities offered by the school. Specific information about some of the considered schools can be known, although the set of known criteria may differ in each school. As another example, consider a company planning to expand its business by offering some new goods for sale in addition to existing ones. Before introducing new products or services, it is necessary to conduct market research and then select the most-profitable products/services. The company has data on the currently offered goods (i.e., sales volume, profitability, and popularity), only the data on the goods considered for introduction to the market are unknown. During pairwise comparisons, the different criteria could be taken into account.
The cases described above have two distinctive features: some of the alternatives considered have initially known priorities, and all the alternatives are compared with respect to more than one criterion. This observation leads us to the idea of the hierarchical HRE method. Such an approach would reduce the multi-criteria problem to a series of simple pairwise comparisons and calculate the weights of concepts related to individual criteria. Finally, the priorities of the alternatives are determined. Moreover, based on previously known estimations, we would be able to estimate actual values of different decision options. Hence, the meaning of the priority is not only relative but real, expressing the current value or unit. In general, using a hierarchy is nothing new. One of the first hierarchical methods —AHP—was introduced by Saaty and is still one of the most popular decision-making methods. The structure of the hierarchical HRE is shown in Figure 2.
The primary data for the hierarchical HRE method are two sets: a set of criteria Q = { q 1 , . . , q s } and a set of alternatives A = { a 1 , . . . , a n } . For each criterion q t ( t 1 , 2 , , s ) we construct a PC matrix C ( q t ) = [ c i j ( t ) ] containing comparisons of elements from A. For this purpose, we establish a set of reference alternatives (with known priorities) A K ( q t ) , and a set of unknown alternatives (i.e., with unknown priorities) A U ( q t ) . Obviously, those sets for each t 1 , 2 , , s , must satisfy the following conditions A = A K ( q t ) A U ( q t ) , A K ( q t ) A U ( q t ) = O . It is worth noting that the set of initially known alternatives and set of unknown alternatives may differ among criteria.
We construct the PC matrix C ( q t ) based on pairwise comparisons of alternatives, where at least one is unknown. Then, we apply the HRE prioritization method for the matrix C ( q t ) . More specifically, we construct the following linear equation:
M ( q t ) w ( q t ) = b ( q t ) ,
where w ( q t ) is a priority vector of alternatives a i A U ( q t ) , M ( q t ) is a k t × k t matrix, where k t means the number of known alternatives ( k t = A K ( q t ) ) in the form
M ( q t ) = 1 1 n 1 c 1 , 2 ( t ) 1 n 1 c 1 , k t ( t ) 1 n 1 c 2 , 1 ( t ) 1 n 1 c 1 , k t ( t ) 1 n 1 c k t , 1 ( t ) 1 n 1 c k t , 2 ( t ) 1
and the constant term vector is given as
b ( q t ) = 1 n 1 c 1 , k t + 1 ( t ) w K ( q t ) ( a k t + 1 ( t ) ) + + 1 n 1 c 1 , n ( t ) w K ( q t ) ( a n ( t ) ) 1 n 1 c 2 , k t + 1 ( t ) w K ( q t ) ( a k t + 1 ( t ) ) + + 1 n 1 c 2 , n ( t ) w K ( q t ) ( a n ( t ) ) 1 n 1 c k t , k t + 1 ( t ) w K ( q t ) ( a k t + 1 ( t ) ) + + 1 n 1 c k t , n ( t ) w K ( q t ) ( a n ( t ) ) ,
where w K ( q t ) is a priority vector for alternatives from A K ( q t ) . We calculate the priority vector w ( q t ) by solving the Equation (11).
In the geometric version of the HRE approach [4] for each criterion q t and the appropriate PC matrix C ( q t ) , we construct the following linear equation:
N ( q t ) μ ( q t ) = d ( q t ) ,
where
N ( q t ) = ( n 1 ) 1 1 1 1 1 1 ( n 1 )
is a matrix of dimensions k t × k t , ( k t = A K ( q t ) ) and
d ( q t ) = l o g c 1 , 2 ( t ) c 1 , 3 ( t ) · · c 1 , k t ( t ) c 1 , k t + 1 ( t ) w ( a k t + 1 ( t ) ) · · c 1 , n ( t ) w ( a n ( t ) ) l o g c 2 , 1 ( t ) c 2 , 3 ( t ) · · c 2 , k t ( t ) c 2 , k t + 1 ( t ) w ( a k t + 1 ( t ) ) · · c 2 , n ( t ) w ( a n ( t ) ) l o g c k t , 1 ( t ) c k t , 2 ( t ) · · c k t , k t 1 ( t ) c k t , k t + 1 ( t ) w ( a k t + 1 ( t ) ) · · c k t , n ( t ) w ( a n ( t ) ) ,
and μ ( q t ) is a logarithmized priority vector of alternatives a i A U ( q t ) . We compute the vector μ ( q t ) by solving (13) and then we calculate the vector w U ( q t ) = w U ( q t ) ( a 1 ( t ) ) , ,   w U ( q t ) ( a k t ( t ) ) using the dependency
μ ( q t ) = log w U ( q t ) ( a 1 ( t ) ) , , log w U ( q t ) ( a k t ( t ) ) for every a i A U ( q t ) .
Based on the priorities given in w U ( q t ) , we create an n dimensional vector w ( q t ) , being a catenation of a priori known and just calculated priorities:
w ( q t ) = w U ( q t ) , w K ( q t ) T .
Similarly, as in AHP, when calculating the final priority vector, we must take into account the fact that the criteria q 1 , . . . , q s may contribute at different extents to the goal. Thus, we construct a PC matrix C ^ , in which every single entry c ^ k l corresponds to the result of comparison q k against q l . Then (using the HRE approach), we calculate the priority vector
w ^ = w ^ ( q 1 ) , , w ^ ( q s ) T .
The final vector of weights is a linear combination of vectors w ( q t ) where scaling factors come from w ^ , i.e.,
w = t = 1 s w ^ ( q t ) w ( q t ) .
Note that if an additive approach is used, the existence of a solution to each HRE Equation (7) used must be verified.

6. Numerical Examples

The HRE method can be used independently or with other strategies for calculating priorities for pairwise comparison matrices. In the first example, we will show the use of HRE in a model with three criteria and seven alternatives. In the second, the HRE approach is combined with the EVM method in a single hierarchical multi-criteria model.
Example 1.
A company, managed by Mr. Smith, runs a sports facility consisting of a sports swimming pool, gym and fitness club. The company is developing dynamically. Therefore, Mr. Smith intends to expand the facility with two further investments. He considers a bowling alley, a professional massage parlor, a trampoline fitness point, or a recreational pool. He would spend a similar amount of money on each project. However, to increase profitability, Mr. Smith intends to conduct a market survey first. Market research will be carried out based on several criteria: Q = { p r — average monthly income (profitability), d u — durability of the equipment (how many months it will serve), and p o p — possible increase in the popularity of a given form of spending free time (also as a result of advertising campaigns)}. A group of management professionals (from now on referred to as experts) will apply the HRE hierarchical approach to the analysis of the obtained data. Let A = { a 1 — a bowling alley, a 2 — a professional massage salon, a 3 — a trampoline fitness point, a 4 — a recreational pool, a 5 — a sports pool, a 6 − gym, and a 7 — fitness club } . The set of reference alternatives is the same for each criterion A K = { a 5 , a 6 , a 7 } .
The first considered criterion is profitability ( p r ) . For objects from the reference set, the values w ( p r ) ( a 5 ) = 20 thous . , w ( p r ) ( a 6 ) = 12 thous . , w ( p r ) ( a 7 ) = 9 thous . represent the average monthly income from the previous year. Considering profitability, experts create a PC matrix C ( p r ) . Each c i j in the matrix C ( p r ) corresponds to the relative profitability of a i with respect to a j . We calculate the values c i j for i , j { 5 , 6 , 7 } , i j as w ( a i ) w ( a j ) , therefore
C ( p r ) = 1 2 3 2 1 2 1 2 1 3 2 3 2 1 2 2 3 1 2 1 2 1 2 1 2 1 1 3 1 4 1 2 3 2 3 2 3 1 2 3 3 2 1 2 2 4 2 3 1 20 12 20 9 1 1 1 3 2 12 20 1 20 9 2 3 1 2 3 2 1 9 20 9 20 1 .
In order to estimate the vector w ( p r ) , we use the HRE method, so according to (11) we construct the matrix M ( p r ) and vector b ( p r )
M ( p r ) = 1 1 n 1 c 1 , 2 1 n 1 c 1 , 3 1 n 1 c 1 , 4 1 n 1 c 2 , 1 1 1 n 1 c 2 , 3 1 n 1 c 2 , 4 1 n 1 c 3 , 1 1 n 1 c 3 , 2 1 1 n 1 c 3 , 4 1 n 1 c 4 , 1 1 n 1 c 4 , 2 1 n 1 c 4 , 3 1
b ( p r ) = 1 n 1 c 1 , 5 w ( p r ) ( a 5 ) + 1 n 1 c 1 , 6 w ( p r ) ( a 6 ) + 1 n 1 c 1 , 7 w ( p r ) ( a 7 ) 1 n 1 c 2 , 5 w ( p r ) ( a 5 ) + 1 n 1 c 2 , 6 w ( p r ) ( a 6 ) + 1 n 1 c 2 , 7 w ( p r ) ( a 7 ) 1 n 1 c 3 , 5 w ( p r ) ( a 5 ) + 1 n 1 c 3 , 6 w ( p r ) ( a 6 ) + 1 n 1 c 3 , 7 w ( p r ) ( a 7 ) 1 n 1 c 4 , 5 w ( p r ) ( a 5 ) + 1 n 1 c 4 , 6 w ( p r ) ( a 6 ) + 1 n 1 c 4 , 7 w ( p r ) ( a 7 ) .
The Equation (11) M ( p r ) w ( p r ) = b ( p r ) , written numerically, takes the form
1 0.111 0.333 0.083 0.25 1 0.333 0.111 0.083 0.083 1 0.056 0.333 0.25 0.5 1 w ( p r ) ( a 1 ) w ( p r ) ( a 2 ) w ( p r ) ( a 3 ) w ( p r ) ( a 4 ) = 5.917 6.667 3.833 6.722
After solving the equation, we obtain the following weights w ( p r ) ( a 1 ) = 11.164 ,   w ( p r ) ( a 2 ) = 13.667 , w ( p r ) ( a 3 ) = 6.863 , w ( p r ) ( a 4 ) = 17.292 , the values of which represent the estimated average monthly income for the objects a 1 , a 2 , a 3 , a 4 . After normalization, the resulting weights are
w ( p r ) = [ 0.124 0.152 0.076 0.192 0.222 0.133 0.1 ] T .
According to the second considered criterion ( d u ) , we compare objects in terms of the durability of the equipment. The values c i j in the matrix C ( d u ) indicate how many times the durability of the equipment of the object a i is better than the durability of the equipment of the object a j . For reference alternatives, the values w ( d u ) ( a 5 ) = 72 , w ( d u ) ( a 6 ) = 24 , and w ( d u ) ( a 7 ) = 36 represent the average number of months the appliance has been used in the last few years. The PC matrix C ( d u ) prepared by the experts has the form
C ( d u ) = 1 3 2 1 1 2 2 3 2 1 3 1 3 2 1 3 1 4 1 2 1 2 1 2 2 3 1 1 3 1 4 2 3 1 2 1 3 3 1 4 5 2 2 2 4 4 5 4 1 72 24 72 36 1 2 2 3 2 1 2 24 72 1 24 36 2 3 2 2 1 2 36 72 36 24 1 .
Using HRE for the matrix C ( d u ) , similarly to the first criterion (see (11) and (12)), we obtain the matrix M ( d u ) and vector b ( d u ) :
M ( d u ) = 1 0.5 0.333 0.167 0.056 1 0.25 0.056 0.083 0.111 1 0.056 0.167 0.5 0.5 1 ,
b ( d u ) = 23 8 8.667 29.6 .
Solving the equation M ( d u ) w ( d u ) = b ( d u ) gives values that indicate the estimated average durability in months:
w ( d u ) ( a 1 ) = 47.183 , w ( d u ) ( a 2 ) = 18.119 , w ( d u ) ( a 3 ) = 17.688 , w ( d u ) ( a 4 ) = 55.367 ,
Scaling the solution vector gives
w ( d u ) = [ 0.174 0.067 0.065 0.205 0.266 0.089 0.133 ] T .
Analysis of the third criterion—the possibility of increasing the popularity of activities offered by individual facilities ( p o p ) —the values c i j in the matrix C ( p o p ) mean how many times the popularity of the object a i can increase in relation to the increase in the popularity of the object a j . The reference alternatives have the values w ( p o p ) ( a 5 ) = 5 , w ( p o p ) ( a 6 ) = 20 , w ( p o p ) ( a 7 ) = 25 , and represent the average percentage increase in popularity in recent years. In this case, the PC-matrix C ( p o p ) proposed by the experts is
C ( p o p ) = 1 1 1 5 3 3 2 2 1 1 1 3 3 3 2 2 5 3 1 4 5 3 3 1 3 1 3 1 4 1 3 2 1 2 1 2 1 3 1 3 1 5 2 3 1 5 20 5 25 1 2 1 2 1 3 2 20 5 1 20 25 1 2 1 2 1 3 2 25 5 25 20 1 .
The auxiliary matrix M ( p o p ) and the vector b ( p o p ) derived from the matrix C ( p o p ) are
M ( p o p ) = 1 0.167 0.033 0.5 0.167 1 0.056 0.5 0.833 0.5 1 0.667 0.056 0.056 0.042 1 ,
b ( p o p ) = 17.5 17.5 26.667 5 .
Using the HRE method to solve the equation M ( p o p ) w ( p o p ) = b ( p o p ) , we obtain the following values:
w ( p o p ) ( a 1 ) = 31.459 , w ( p o p ) ( a 2 ) = 32.93 , w ( p o p ) ( a 3 ) = 77.21 , w ( p o p ) ( a 4 ) = 11.794 ,
These numbers represent the estimated average growth in popularity of objects a 1 , a 2 , a 3 , a 4 . After normalization, the solution vector of the weights is equal:
w ( d u ) = [ 0.155 0.162 0.38 0.058 0.025 0.098 0.123 ] T .
Since in the considered case the weights of the criteria are known (in the past, the company has encountered comparable cases and has used similar weights to address this set of criteria) and equal
w ^ = [ 0.5 0.2 0.3 ] T .
the final priorities of individual objects are calculated according to the schema (15):
w ( a i ) = 0.5 · w ( p r ) ( a i ) + 0.2 · w ( d u ) ( a i ) + 0.3 · w ( p o p ) ( a i ) .
So the final weight vector is equal
w = [ 0.143 0.138 0.165 0.154 0.172 0.114 0.114 ] T ,
which means that investing in a trampoline fitness point and a recreational pool may be slightly more profitable than any of the other two new facilities.
Example 2.
Mr. Kowalski became interestedin collecting porcelain. He paid particular attention to the cups. Mr. Kowalski visited a local antique store and found five exciting items. He would like to buy two of them, which are the most valuable. So, he asked our art appraiser’s office for help in this matter. Unfortunately, there is no universal algorithm to value porcelain step by step. Moreover, the valuation of porcelain is complex and unmeasurable due to the peculiarities of different markets. Identical cups can have different prices depending on the region where they are offered. For example, Polish collectors are fond of German products, which undervalue Polish porcelain. At the same time, the English are proud of their wares and are eager to buy them, and they also value Chinese porcelain. Therefore, when analyzing trends in the porcelain market, we should also consider the place of purchase and ask an expert to evaluate the various criteria by which they make their purchasing decision. Let us assume that Mr. Kowalski lives in Poland, so we will take the values of the different criteria according to the Polish market and help him make the right decision. Our market research shows that the final value of mugs here depends on several factors.
The most straightforward and obvious rule is that the more difficult it is to obtain products from a specific series, period, or manufacturer, the higher the price it receives. Due to the fragility of porcelain, whole sets are more expensive than single pieces. Some units have special status: they are almost unattainable on the market, are valued as “unlimited”, and are usually auctioned at exorbitant prices. The products from the first period of operation of the factory in Meissen, Germany, run by the “father of European porcelain”, Friedrich Böttger, are a good example. Moreover, products of a particular manufacturer may be more in demand at one time and less at another.
The second criterion is the reputation of the manufacturer. It is closely influenced by the history of a given factory, the quality of the porcelain produced there, or its uniqueness. It is also worth remembering that reputation is tied to the current fashion—it depends on time and location. For example: in Germany, Meissen porcelain is the most valued, while in Poland, collectors like Rosenthal more.
Traces, called signatures, can often be found on porcelain. They define its origin, age, model name, and sometimes the project’s author. Some signatures even specify the type of porcelain the item is made of. Thanks to the signature, we can estimate the value of the product. Usually, the older it is, the more valuable it is. The first characters were discovered on Chinese porcelain from the 14th century. The well-known collector’s rarities include, among others, unique items from the Ming Dynasty. There are a maximum of a dozen of them in museum collections worldwide. They are almost nonexistent at auctions; they reach exorbitant prices if they appear. In Europe, the first porcelain products were made in Saxony in the 18th century during the reign of Augustus II and were marked with letters and numbers. Soon, other manufacturers started putting their signatures on the bottom of dishes. Depending on the creation time, the signatures differ in color and how they were made (painting, printing, decal). For example, the markings also changed within one factory with a change of ownership. On their basis, the experienced collector can accurately determine the creation time.
Ancient porcelain is exceptionally fragile, which, of course, increases its value. However, in the event of damage, even porcelain from an excellent manufacturer is worthless. Therefore, before buying, you should evaluate the condition of a single piece of porcelain. An ideal copy will be a product with no acquired or primary defects. By primary defects, we mean any flaws visible in the product that are traces of an ineffective porcelain forming process - stains, painting, or decals errors. Reputable manufacturers sold only first-class products. For example, all defective products are considered broken at Rosenthal factories, so they do not end up on the market. Unlike acquired defects, primary defects usually do not influence the product’s price. Acquired defects affecting porcelain price include nicks, sticking together, painting with other paint, or visible spider web shape cracks.
The quality of porcelain is primarily determined by the appropriate selection of raw materials, especially a large amount of kaolin. Good porcelain is hard, durable, and scratch-resistant. It requires manual work as machines do not provide a top-quality finish. Therefore, each item is processed precisely and accurately by hand in good factories. It guarantees the highest quality and timeless beauty of the final product. Ornaments on porcelain also determine its quality. The most valuable is hand-decorated porcelain. Often, ornaments are made of noble materials—real gold or, for example, platinum.
Based on the above information, we will compare the value of considered porcelain cups against five criteria: factory, period of creation, uniqueness, state of preservation and quality. Let
Q = { man , per , un , st , qua }
denote a set of significant criteria, where man , per , un , st , qua mean factory, period of creation, uniqueness, state of preservation and quality, respectively.
Let A = { a 1 , a 2 , a 3 , a 4 , a 5 } denote a considered collection of cups. First, we compare the cups in terms of the first criterion. Depending on the reputation of the factory, each cup is assigned a value from 1 to 10, where 1 denotes the lowest possible score. All cups come from different factories. We know three of them (for the alternatives a 3 , a 4 and a 5 ) and we can assign values to them based on our research (e.g., the indicated values may be the average list price of cups produced by the indicated manufacturers expressed in tens of dollars). They are w ( man ) ( a 3 ) = 8.7 , w ( man ) ( a 4 ) = 4.2 ,   w ( man ) ( a 5 ) = 7.2 , respectively. As it is the first time that we are dealing with the products of manufacturers from which the cups a 1 and a 2 come from, we will apply the HRE procedure. As a result, we obtain the matrix:
C ( m a n ) = 1 2 1 3 3 1 1 2 1 1 5 1 1 2 3 5 1 8.7 4.2 8.7 7.2 1 3 1 4.2 8.7 1 4.2 7.2 1 2 7.2 8.7 7.2 4.2 1 .
Since the weight of cups a 3 , a 4 and a 5 are known, the pairs ( a 3 , a 4 ) , ( a 3 , a 5 ) and ( a 4 , a 5 ) are not evaluated, but c i j ( m a n ) is calculated as w ( man ) ( i ) w ( man ) ( j )   ( 3 i < j 5 ) . To estimate the expected values of factories A and B, we follow the HRE procedure. Since A U ( ( man ) = { a 1 , a 2 } and A K ( man ) = { a 3 , a 4 , a 5 } , the matrix M ( man ) and the constant term vector b ( man ) have the form (see (11))
M ( man ) = 1 1 n 1 c 1 , 2 1 n 1 c 2 , 1 1 ,
b ( m a n ) = 1 n 1 c 1 , 3 w ( man ) ( a 3 ) + 1 n 1 c 1 , 4 w ( man ) ( a 4 ) + 1 n 1 c 1 , 5 w ( man ) ( a 5 ) 1 n 1 c 2 , 3 w ( man ) ( a 3 ) + 1 n 1 c 2 , 4 w ( man ) ( a 4 ) + 1 n 1 c 2 , 5 w ( man ) ( a 5 ) .
Hence, numerically,
M ( man ) = 1 0.476 0.119 1 , b ( m a n ) = 5.675 2.385 .
As M ( man ) w ( man ) = b ( man ) has an admissible solution (i.e., real and positive), the ranking vector w ( m a n ) takes the form
w ( m a n ) = [ 7.22 3.244 8.7 4.2 7.2 ] T .
Hence, according to experts’ judgments w ( m a n ) ( a 1 ) = 7.22 and w ( m a n ) ( a 2 ) = 3.244 . After rescaling
w ( m a n ) = [ 0.236 0.106 0.284 0.137 0.235 ] T .
The second considered criterion applies to the date the cup was manufactured. Based on the signatures placed on the cups a 1 , a 3 and a 5 , we can read the years of production and assign the appropriate number of points according to age. The older the cup, the more points it becomes. Finally, we receive the results w ( per ) ( c 1 ) = 8 ,   w ( per ) ( c 3 ) = 10 , w ( per ) ( c 5 ) = 5 . The cup a 2 has no signature and the signature of the cup a 4 does not include the year, so in order to establish the period of creation for these two cups we ask an expert for help. In the case of this criterion A U ( per ) = { a 2 , a 4 } and A K ( per ) = { a 1 , a 3 , a 5 } . As a result of the expert’s assessment, we obtain the matrix:
C ( p e r ) = 1 2 8 10 4 8 5 1 2 1 1 3 2 1 10 8 3 1 5 10 5 1 4 1 2 1 5 1 1 2 5 8 1 5 10 2 1 .
In order to obtain the values of the priority vector w ( per ) for cups a 2 and a 4 , it is enough to create the matrix M ( per ) and vector b ( p e r ) using the HRE method and then solve the equation M ( per ) w ( per ) = b ( per ) . The matrix M ( per ) can be obtained from the matrix C ( per ) by removing the columns and rows corresponding to known alternatives from the set A K ( per ) and multiplying values out of a diagonal by 1 n 1 . The values removed only from the columns (but not from the rows) are used to form the vector b ( per ) . Thus,
M ( per ) = 1 1 n 1 c 2 , 4 1 n 1 c 4 , 2 1 ,
b ( per ) = 1 n 1 c 2 , 1 w ( per ) ( a 1 ) + 1 n 1 c 2 , 3 w ( per ) ( a 3 ) + 1 n 1 c 2 , 5 w ( per ) ( a 5 ) 1 n 1 c 4 , 1 w ( per ) ( a 1 ) + 1 n 1 c 4 , 3 w ( per ) ( a 3 ) + 1 n 1 c 4 , 5 w ( per ) ( a 5 ) ,
and the equation M ( per ) w ( per ) = b ( per ) takes the shape
1 0.498 0.124 1 w ( per ) ( a 2 ) w ( per ) ( a 4 ) = 37 12 13 8 .
The solution of the above equation is the following vector:
w ( per ) = [ 8 4.15 10 2.14 5 ] T .
Thus, in the category “period of creation”, the expected number of points on a scale from 1 to 10 for cups a 2 and a 4 are 4.15 and 2.14 , respectively. After normalization,
w ( p e r ) = [ 0.273 0.141 0.341 0.073 0.171 ] T .
The third considered criterion concerns uniqueness, although this time the set of reference alternatives is empty. Thus, once again, we ask an expert to compare in pairs the profitability of purchasing a given cup in terms of its uniqueness on the market. His assessments form the following PC matrix:
C ( u n ) = 1 1 2 1 2 1 2 1 3 2 1 2 1 1 2 2 1 2 1 1 1 2 2 1 1 1 1 2 3 2 2 2 1 .
In such a case, we can use EVM (Section 3.1), i.e., solve the equation C ( un ) w = λ max w . Hence, after appropriate rescaling of the vector w, we obtain
w ( u n ) = [ 0.097 0.214 0.161 0.182 0.345 ] T .
The condition of the cup depends mainly on two factors: primary damage and acquired damage. Therefore, the experts assess all our alternatives against two sub-criteria: level of primary damage and level of acquired damage. As a result of comparing the degree of primary damage, we obtain the following matrix:
C ( pd ) = 1 9 9 8 3 5 3 1 9 1 4 9 1 9 1 8 1 9 9 4 1 1 9 1 8 3 8 9 9 1 2 3 5 8 8 1 2 1 .
The comparisons of acquired damage are as follows:
C ( ad ) = 1 1 9 1 9 4 5 1 9 9 1 3 4 9 6 5 9 4 3 1 9 9 4 5 4 1 9 1 9 1 1 6 9 5 6 4 9 6 1 .
The priority vectors for the C ( pd ) and C ( ad ) matrices calculated using EVM are as follows:
w ( pd ) = 0.416 0.029 0.04 0.289 0.223 T ,
and
w ( ad ) = 0.033 0.301 0.39 0.038 0.235 T .
Comparing the significance of both types of damage, the expert found that the acquired damage is two times more important than the primary damage. Hence, the ranking of both damage significance (ds) sub-criteria is as follows:
w ( ds ) = w ( ds ) ( pd ) w ( ds ) ( ad ) = 0.333 0.666 .
The state of preservation for the i-th cup is calculated as
w ( st ) ( a i ) = w ( ds ) ( pd ) · w ( pd ) ( a i ) + w ( ds ) ( ad ) · w ( ad ) ( a i ) ,
for i = 1 , , 5 , hence, finally the assessment of the state of the cups’ preservation is as follows:
w ( st ) = 0.16 0.211 0.27 0.12 0.23 T .
When assessing the quality of porcelain, the expert considered its resistance to scratching. In the case of two cups a 1 and a 2 , it was known and calculated as the average number of significant scratches among a similar class of products appearing at online auctions in the last five years. Thus, A U ( qua ) = { a 3 , a 4 , a 5 } and A K ( qua ) = { a 1 , a 2 } , where w ( qua ) ( a 1 ) = 5.7 and w ( qua ) ( a 2 ) = 2.4 . As a result of comparing the quality of alternatives by an expert, we obtain the matrix
C ( qua ) = 1 5.7 2.4 9 3 7 5 2.4 5.7 1 9 9 5 8 3 1 9 1 9 1 1 9 1 9 1 3 5 9 9 1 6 5 5 7 3 8 9 5 6 1 .
Based on the HRE approach, the auxiliary matrix and the constant term vector are
M ( qua ) = 1 0.0259 0.0259 2.104 1 0.28 2.104 0.194 1 , b ( qua ) = 0.225 0.808 1.242 .
Solving the HRE equation M ( qua ) w = b ( qua ) reveals that the expected number of scratches for a 3 , a 4 and a 5 are 0.344 , 2.2 and 2.39 . This allows us to form the ranking vector
w ^ ( qua ) = [ 5.7 2.4 0.344 2.2 2.39 ] T .
However, since the ranking values represent the expected number of scratches, it means that the higher the value, the worse the product. To change this, we will raise each of the elements of w ^ ( qua ) to the −1 power and then normalize it. As a result, we obtain the vector:
w ( qua ) = [ 0.04 0.0954 0.665 0.104 0.0955 ] T .
It is worth noting that, thanks to the last transformation, if the i-th alternative is two times worse than the j-th one according to w ^ ( qua ) , then according to w ( qua ) , it is two times better than the j-th alternative, etc. The obtained vector is a ranking of the cups against the porcelain quality criterion.
In order to establish the weights of criteria ( q 1 = m a n , q 2 = p e r , q 3 = u n , q 4 = s t , q 5 = q u a ) the expert constructed a PC-matrix, in which c ^ i j denotes how many times the criterion q i is more preferred than the criterion q j .
C ^ = 1 1 3 2 3 3 1 1 3 2 3 3 2 3 2 3 1 2 2 1 3 1 3 1 2 1 1 1 3 1 3 1 2 1 1
Then, using EVM again, we obtain a priority vector for the criteria:
w ^ = [ 0.3 0.3 0.2 0.1 0.1 ] T .
The final weights of individual cups a 1 , . . . , a 5 are calculated according to the well-known schema:
w ( a i ) = 0.3 w ( m a n ) ( a i ) + 0.3 w ( p e r ) ( a i ) + 0.2 w ( u n ) ( a i ) + 0.1 w ( s t ) ( a i ) + 0.1 w ( q u a ) ( a i ) .
Thus, we obtain
w = [ 0.1922 0.1478 0.3138 0.1222 0.2237 ] T .
The recommendation presented to Mr. Kowalski puts in the first place the cup a 3 with the rank 0.3138 , followed by a 5 with the rank 0.2237 , then a 1 , a 2 and a 4 with the ranking values 0.1922 , 0.1478 and 0.1222 .

7. Discussion and Summary

We demonstrate the hierarchical HRE as a multiple-criteria decision-making method based on the same assumptions found in EVM or GMM. Unlike AHP, HRE allows the use of actual measured values. HRE fills a gap by combining tangible and intangible criteria within a single hierarchical model. The ranking obtained in HRE maintains the same scales on the input and output of the priority deriving procedure, e.g., if the reference values are expressed in a particular currency, the vector of weights in the output contains the “valuation” of the alternatives for that specific currency. This property makes HRE an exciting proposition for anyone who needs a well-defined scale in the output of a ranking procedure. Combining HRE with other methods can relativize the outcome and make the loss of a “unit” at the output of the ranking procedure. Future research will explore uncertain knowledge representations and hybrid multi-criteria models, providing practical and theoretical insights. Guidelines will be developed to facilitate decision-making models using HRE.

Author Contributions

Conceptualization, A.K.; Methodology, K.K.; Investigation, A.K.; Supervision, K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are also grateful to anonymous reviewers for their insightful observations and comments. Special thanks are due to Ian Corkill for his editorial help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vasconcelos, G.R.; Maria De Miranda Mota, C. Exploring Multicriteria Elicitation Model Based on Pairwise Comparisons: Building an Interactive Preference Adjustment Algorithm. Math. Probl. Eng. 2019, 2019, 2125740. [Google Scholar] [CrossRef]
  2. Harker, P.T.; Vargas, L.G. The Theory of Ratio Scale Estimation: Saaty’s Analytic Hierarchy Process. Manag. Sci. 1987, 33, 1383–1403. [Google Scholar] [CrossRef]
  3. Kułakowski, K. Heuristic Rating Estimation Approach to The Pairwise Comparisons Method. Fundam. Inform. 2014, 133, 367–386. [Google Scholar] [CrossRef] [Green Version]
  4. Kułakowski, K.; Grobler-Dębska, K.; Wąs, J. Heuristic rating estimation: Geometric approach. J. Glob. Optim. 2015, 62, 529–543. [Google Scholar] [CrossRef] [Green Version]
  5. Saaty, T.L. A scaling method for priorities in hierarchical structures. J. Math. Psychol. 1977, 15, 234–281. [Google Scholar] [CrossRef]
  6. Crawford, R.; Williams, C. A note on the analysis of subjective judgement matrices. J. Math. Psychol. 1985, 29, 387–405. [Google Scholar] [CrossRef]
  7. Kułakowski, K. Understanding the Analytic Hierarchy Process; Chapman and Hall/CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar] [CrossRef]
  8. Colomer, J.M. Ramon Llull: From ‘Ars electionis’ to social choice theory. Soc. Choice Welf. 2011, 40, 317–328. [Google Scholar] [CrossRef]
  9. Hägele, G.; Pukelsheim, F. Llull’s writings on electoral systems. Stud. Lul. 2001, 41, 3–38. [Google Scholar]
  10. Condorcet, M. Essai sur L’application de L’analyse à la Probabilité des Décisions Rendues à la Pluralité des Voix; Imprimerie Royale: Paris, France, 1785. [Google Scholar]
  11. Dzhafarov, E.; Colonius, H. The Fechnerian idea. Am. J. Psychol. 2011, 124, 127–140. [Google Scholar] [CrossRef] [Green Version]
  12. Thurstone, L.L. A Law of Comparative Judgment, reprint of an original work published in 1927. Psychol. Rev. 1994, 101, 266–270. [Google Scholar] [CrossRef]
  13. Ramik, J. Pairwise Comparisons Method: Theory and Applications in Decision Making; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  14. Yuen, K.K.F. Combining compound linguistic ordinal scale and cognitive pairwise comparison in the rectified fuzzy TOPSIS method for group decision making. Fuzzy Optim. Decis. Mak. 2014, 13, 105–130. [Google Scholar] [CrossRef]
  15. Brans, J.; Mareschal, B. Promethee Methods. In Multiple Criteria Decision Analysis: State of the Art Surveys; Figueira, J., Greco, S., Ehrgott, M., Eds.; Springer: Boston, MA, USA; Dordrecht, The Netherlands; London, UK, 2016; pp. 187–219. [Google Scholar]
  16. Bana e Costa, C.; De Corte, J.M.; Vansnick, J. On the Mathematical Foundation of MACBETH. In Multiple Criteria Decision Analysis: State of the Art Surveys; Figueira, J., Greco, S., Ehrgott, M., Eds.; Springer: Boston, MA, USA; Dordrecht, The Netherlands; London, UK, 2016; pp. 421–463. [Google Scholar]
  17. Rezaei, J. Best-worst multi-criteria decision-making method. Omega 2015, 53, 49–57. [Google Scholar] [CrossRef]
  18. Kadziński, M.; Ciomek, K.; Słowiński, R. Modeling assignment-based pairwise comparisons within integrated framework for value-driven multiple criteria sorting. Eur. J. Oper. Res. 2015, 241, 830–841. [Google Scholar] [CrossRef]
  19. Saaty, T.L.; Vargas, L.G. Decision Making with the Analytic Network Process: Economic, Political, Social and Technological Applications with Benefits, Opportunities, Costs and Risks; Springer Science and Business Media: Boston, MA, USA, 2013; Volume 195. [Google Scholar]
  20. Lidinska, L.; Jablonsky, J. AHP model for performance evaluation of employees in a Czech management consulting company. Cent. Eur. J. Oper. Res. 2018, 26, 239–258. [Google Scholar] [CrossRef]
  21. Kramulová, J.; Jablonský, J. AHP model for competitiveness analysis of selected countries. Cent. Eur. J. Oper. Res. 2016, 24, 335–351. [Google Scholar] [CrossRef]
  22. Rezaei, J.; Ortt, R. Multi-criteria supplier segmentation using a fuzzy preference relations based AHP. Eur. J. Oper. Res. 2013, 225, 75–84. [Google Scholar] [CrossRef]
  23. Charouz, J.; Ramík, J. A multicriteria decision making at portfolio management. Ekon. Manag. 2010, 13, 44–52. [Google Scholar]
  24. Szabo, Z.K.; Szádoczki, Z.; Bozóki, S.; Stanciulescu, G.C.; Szabo, D. An analytic hierarchy process approach for prioritisation of strategic objectives of sustainable development. Sustainability 2021, 13, 2254. [Google Scholar] [CrossRef]
  25. Jamshidi, A.; Jamshidi, F.; Ait-Kadi, D.; Ramudhin, A. A review of priority criteria and decision-making methods applied in selection of sustainable city logistics initiatives and collaboration partners. Int. J. Prod. Res. 2019, 57, 5175–5193. [Google Scholar] [CrossRef] [Green Version]
  26. Lu, F.; Bi, H.; Feng, W.; Hu, Y.; Wang, S.; Zhang, X. A Two-Stage Auction Mechanism for 3PL Supplier Selection under Risk Aversion. Sustainability 2021, 13, 9745. [Google Scholar] [CrossRef]
  27. Lu, F.; Hu, Y.; Bi, H.; Huang, M.; Zhao, M. An Auction Approach for Cost and Schedule Management of IT Outsourcing Project. Asia-Pac. J. Oper. Res. 2018, 35, 1850034. [Google Scholar] [CrossRef]
  28. Al-Subhi Al-Harbi, K.M. Application of the AHP in project management. Int. J. Proj. Manag. 2001, 19, 19–27. [Google Scholar] [CrossRef]
  29. Hyde, R.A.; Davis, K., A. Military applications of the analytic hierarchy process. Int. J. Multicriteria Decis. Mak. 2012, 2, 267. [Google Scholar] [CrossRef]
  30. Kou, G.; Ergu, D.; Lin, C.S.; Chen, Y. Pairwise comparison matrix in multiple criteria decision making. Technol. Econ. Dev. Econ. 2016, 22, 738–765. [Google Scholar] [CrossRef] [Green Version]
  31. Ho, W.; Ma, X. The state-of-the-art integrations and applications of the analytic hierarchy process. Eur. J. Oper. Res. 2018, 267, 399–414. [Google Scholar] [CrossRef]
  32. Darko, A.; Chan, A.P.C.; Ameyaw, E.E.; Owusu, E.K.; Pärn, E.; Edwards, D.J. Review of application of analytic hierarchy process (AHP) in construction. Int. J. Constr. Manag. 2019, 19, 436–452. [Google Scholar] [CrossRef]
  33. Bana e Costa, C.A.; Vansnick, J. A critical analysis of the eigenvalue method used to derive priorities in AHP. Eur. J. Oper. Res. 2008, 187, 1422–1428. [Google Scholar] [CrossRef]
  34. Barzilai, J.; Golany, B. AHP rank reversal, normalization and aggregation rules. Inf. Syst. Oper. Res. 1994, 32, 57–64. [Google Scholar] [CrossRef]
  35. Dyer, J.S. Remarks on the analytic hierarchy process. Manag. Sci. 1990, 36, 249–258. [Google Scholar] [CrossRef]
  36. Simon, D.; Zhao, M.; Qin, S.; Li, Q.; Lu, F.; Shen, Z. The Likelihood Ranking Methods for Interval Type-2 Fuzzy Sets Considering Risk Preferences. Math. Probl. Eng. 2015, 2015, 680635. [Google Scholar] [CrossRef] [Green Version]
  37. Sun, F.; Lu, F.; Bi, H.; Yu, C. FAHP based evaluation of IT Outsourcing risk. In Proceedings of the Proceeding of the 11th World Congress on Intelligent Control and Automation, Shenyang, China, 29 June–4 July 2014; pp. 602–604. [Google Scholar] [CrossRef]
  38. Ramik, J.; Korviny, P. Inconsistency of pair-wise comparison matrix with fuzzy elements based on geometric mean. Fuzzy Sets Syst. 2010, 161, 1604–1613. [Google Scholar] [CrossRef]
  39. Faramondi, L.; Oliva, G.; Setola, R.; Bozóki, S. Robustness to rank reversal in pairwise comparison matrices based on uncertainty bounds. Eur. J. Oper. Res. 2023, 304, 676–688. [Google Scholar] [CrossRef]
  40. Duleba, S.; Çelikbilek, Y.; Moslem, S.; Esztergár-Kiss, D. Application of grey analytic hierarchy process to estimate mode choice alternatives: A case study from Budapest. Transp. Res. Interdiscip. Perspect. 2022, 13, 100560. [Google Scholar] [CrossRef]
  41. Hovanov, N.V.; Kolari, J.W.; Sokolov, M.V. Deriving weights from general pairwise comparison matrices. Math. Soc. Sci. 2008, 55, 205–220. [Google Scholar] [CrossRef]
  42. Ishizaka, A.; Lusti, M. How to derive priorities in AHP: A comparative study. Cent. Eur. J. Oper. Res. 2006, 14, 387–400. [Google Scholar] [CrossRef] [Green Version]
  43. Barzilai, J. Deriving weights from pairwise comparison matrices. J. Oper. Res. Soc. 1997, 48, 1226–1232. [Google Scholar] [CrossRef]
  44. Saaty, T.L. Decision-making with the AHP: Why is the principal eigenvector necessary. Eur. J. Oper. Res. 2003, 145, 85–91. [Google Scholar] [CrossRef]
  45. Koczkodaj, W.W. A new definition of consistency of pairwise comparisons. Math. Comput. Model. 1993, 18, 79–84. [Google Scholar] [CrossRef]
  46. Brunelli, M. A survey of inconsistency indices for pairwise comparisons. Int. J. Gen. Syst. 2018, 47, 751–771. [Google Scholar] [CrossRef] [Green Version]
  47. Miller, J.R. The Assessment of Worth: A Systematic Procedure and Its Experimental Validation. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 1966. [Google Scholar]
  48. Kułakowski, K. Notes on the existence of a solution in the pairwise comparisons method using the heuristic rating estimation approach. Ann. Math. Artif. Intell. 2016, 77, 105–121. [Google Scholar] [CrossRef] [Green Version]
Figure 1. How does the HRE priority-deriving method work?
Figure 1. How does the HRE priority-deriving method work?
Mathematics 11 02806 g001
Figure 2. The structure of hierarchical HRE.
Figure 2. The structure of hierarchical HRE.
Mathematics 11 02806 g002
Table 1. A list of the important symbols.
Table 1. A list of the important symbols.
SymbolDescription
Aset of alternatives
A K set of alternatives with a priori known priorities (in HRE)
A U set of alternatives with a priori unknown priorities (in HRE)
a i i-th alternative
Cpairwise comparisons matrix
c i j comparison of the i-th and j-th alternatives
w ( a i ) weight (priority) of the i-th alternative
wvector of priority weights
CI Saaty’s consistency index
K Koczkodaj’s inconsistency index
Mauxiliary matrix (in the additive HRE method)
bconstant terms vector (in the additive HRE method)
Nauxiliary matrix (in the geometric HRE method)
dconstant terms vector (in the geometric HRE method)
μ logarithmized priority vector (in the geometric HRE method)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kędzior, A.; Kułakowski, K. Multiple-Criteria Heuristic Rating Estimation. Mathematics 2023, 11, 2806. https://doi.org/10.3390/math11132806

AMA Style

Kędzior A, Kułakowski K. Multiple-Criteria Heuristic Rating Estimation. Mathematics. 2023; 11(13):2806. https://doi.org/10.3390/math11132806

Chicago/Turabian Style

Kędzior, Anna, and Konrad Kułakowski. 2023. "Multiple-Criteria Heuristic Rating Estimation" Mathematics 11, no. 13: 2806. https://doi.org/10.3390/math11132806

APA Style

Kędzior, A., & Kułakowski, K. (2023). Multiple-Criteria Heuristic Rating Estimation. Mathematics, 11(13), 2806. https://doi.org/10.3390/math11132806

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop