A heuristic rating estimation algorithm for the pairwise comparisons method

The pairwise comparisons method is a powerful tool used for establishing the relative order between different concepts in situations in which it is difficult (or sometimes even impossible) to provide explicit rating. Appropriate ratings are determined by solving eigenvalue problem for the pairwise comparisons matrix. This study presents a new iterative heuristic rating estimation algorithm that tries to deal with the situation when exact estimations for some concepts (stimulus) CK\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_{K}$$\end{document} are a priori known and fixed, whilst the estimates for the others (unknown concepts CU\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_{U}$$\end{document}) need to be computed. The relationship between the local estimation error, understood as the average absolute error E(c)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E(c)$$\end{document} over all direct estimates for the concept c∈CU\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c\in C_{U}$$\end{document} and the pairwise comparisons matrix inconsistency index is shown. The problem of convergence of subsequent intermediate results is discussed and the convergence conditions are given.


Introduction
The pairwise comparisons (PC) method was introduced in its early form by Fechner (1966), then it was popularized and developed by Thurstone (1994). The introduction of hierarchical structures by Saaty (2008) was another important contribution to the PC method, providing the methodology and practical ways to deal with the large amounts of criterion parameters. Initially the PC method was used in the scientific study of psychometrics and psychophysic (Thurstone 1994), however, it then came to be used in other areas of applications, such as complex decision theory (Saaty 2008), economics (Peterson and Brown 1998), voting systems (Tideman 1987) and others. Since the data, which are input to the PC method are the result of human judgment, it is very easy for inaccuracy to occur. Hence, the input data set is frequently ambiguous and does not allow users to draw firm conclusions. There are several indexes of the data inconsistency (Bozóki and Rapcsak 2008), including the best known Saaty's eigenvector method (Saaty 1980), Least Squares Method, Chi Squares Method, Koczkodaj's distance based inconsistency index (Koczkodaj 1993), and others. Using these indexes, the reliability of the data can be assessed, hence, the confidence in the result can be evaluated. The answer to the question of how much input data must be consistent to ensure the result reliability is the subject of empirical research. For instance, according to Saaty's recommendation every occurrence of the consistency ratio greater than or equal to 0.1 should be the subject of re-examination of the pairwise judgments until the inconsistency becomes less than or equal to the desired value (Bozóki and Rapcsak 2008;Triantaphyllou et al. 1990). The main criticism of this approach relates to its separation from the data and lack of localizing the most problematic matrix elements (Bozóki and Rapcsak 2008;Koczkodaj 1993;Triantaphyllou et al. 1990). In contrast, Koczkodaj's inconsistency index has a meaningful interpretation and provides information about the inconsistency location, but it does not provide an exact answer to the question of how good the average data sample is. Inconsistency identified as too high, must be reduced to an acceptable level (ideally to zero). Since the ratio coefficients, which are the input to the PC method, frequently represent experts' judgements, thus a natural way of inconsistency reduction is to call the expert panel once again and ask the professionals gathered to agree on the opinion (Gomes 1993). Because such a solution is usually time-consuming and expensive, heuristic algorithms of inconsistency reduction have been proposed (Koczkodaj and Orłowski 1999;Koczkodaj and Szarek 2010;Gomes 1993;Xu and Wei 1999;Bozóki 2008;Temesi 2006;Cao et al. 2008). The result of these algorithms is a new set of data, which preserves most of the decision maker's original judgment structure and significantly reduces the data inconsistency.
The proposed innovative solution approaches the problem differently. It does not attempt to minimize inconsistency in the data, but rather proposes a way of using the data, which takes into account their inconsistency. Hence, knowing the exact values for a few concepts 1 and some inconsistent set of ratios between them, the data analyst is able to estimate values of all other concepts with errors depending on the degree of data inconsistency. The presented approach is comparable to the inconsistency reduction methods mentioned above, since the set of concepts for which the estimates are known can be easily transformed into a consistent set of data (as addressed in Sect. 8).
The first and second sections of the article focus on the presentation of the necessary facts and definitions concerning the pairwise comparisons methods. Section three formulates formally the problem considered in the work. The fourth and fifth, sections provide the Heuristic Rating Estimation (HRE) algorithm together with a numerical example demonstrating the algorithm application in practice. The sixth section shows the relationship between errors of estimations obtained by using the HRE algorithm and the data inconsistency. The next, seventh section addresses the case of which the subsequent estimation sets converge to some limit. Finally, the last two sections contain the closing discussion and summary.

A pairwise comparisons method
. . , n} a PC matrix that expresses some quantitative relation R over the finite set of concepts . . , n}} where C is a non empty universe of concepts and R(c i , c j ) = m i j , R(c j , c i ) = m ji . Traditionally, these concepts are interpreted as subjective stimuli (Thurstone 1994), whilst the values m i j and m ji are considered as the relative importance indicators (stimulus intensities), so that according to the best knowledge of an expert the significance of c i equals m i j c j .
Since the knowledge stored in the PC matrix usually comes from different professionals in the field of relation R, it often results in inaccuracy. In such a case, reasoning using the data gathered in M may give ambiguous results. This observation gave rise to the research on the concept of data consistency. In the ideal case, M is consistent, and there is no doubt as regards the value assigned to the concept c j if the value assigned to c i and m ji is known. Unfortunately, in practice the knowledge in M is inconsistent, and professionals using the PC method have to deal with this incoherence. Thus, it is important to answer the question of how inconsistent the knowledge in the PC matrix is. There are a number of inconsistency indexes, including Eigenvecor Method (Saaty 2008), Least Squares Method, Chi Squares Method (Bozóki and Rapcsak 2008), Koczkodaj's distance based inconsistency index (Koczkodaj 1993) and others. For the purpose of this article, the Koczkodaj's distance based inconsistency index has been adopted since it is the only localizing index amongst the above mentioned indexes.
Definition 2 Koczkodaj's distance based inconsistency index K of n ×n and (n > 2) reciprocal matrix M is equal to Intuitively speaking, since in an "ideal" matrix ∀i, j, k ∈ {1, . . . , n} : m i j · m jk · m ki = 1 the Koczkodaj's index localizes the worse triad (by the Euclidean distance) from this ideal.
PC matrices may be over-complete, i.e. there is more than one estimation describing one and the same relation between v i and v j , but they can be also incomplete, i.e. not all values m ji are defined. While the first situation can be addressed in many ways, for instance the estimates related to the same pair of concepts can be averaged, the other is not desirable. It indicates that the model is lacking of knowledge, although the missing estimates can be somehow compensated (for example, by employing the properties of reciprocity, consistency or transitivity) . Hence, for further consideration, it will be assumed that the n × n PC matrix M is complete in the sense specified below.
Definition 3 A matrix M (n × n) is said to be complete if ∀i, j = 1, . . . , n : m i j is defined.
Due to the further consideration it is useful to define a graph structure over the set C and matrix M.

Definition 4 Let a pairwise comparisons graph
is a function of experts' assessments Wherever it does not raise doubts, instead of M(e) ∧ e = (u, v) the function M will be written with two arguments i.e. M (u, v).

Problem formulation
Since the concepts are linked to each other by a quantitative relation R, then assuming that the exact values of some concepts are known, the values of the others should be proposed. Thus, let μ : C R + be a partial function that assigns to some concepts from C ⊂ C positive values from R + . Hence, the concepts for which the actual value μ is known are denoted by C K ⊂ C and called known concepts, whilst concepts for which μ need to be determined are denoted C U = C\C K and called unknown concepts. The relation between different concepts in terms of the function μ is represented in the form of the PC matrix M, so that m ji μ(v i ) = μ(v j ). Since m ji usually aims to express how much greater or smaller v i is than v j with respect to μ by convention it is assumed that elements of the n × n PC matrix are positive real numbers, i.e. m ji ∈ R + where i, j ∈ 1, . . . , n.
The presented method aims to provide an iterative heuristic estimation algorithm that for all v ∈ C U proposes the appropriate value of μ(v). In this approach The algorithm is iterative and sets the new expected values based on the ones previously determined. It stops either when it reaches a fixed number of iterations or (if convergent) when calculations reach the desired accuracy. Although during the course of the algorithm the new values for μ (v) and v ∈ C U are calculated, M remains unchanged and serves as reference data.
Listing 1: Heuristic ratio estimation algorithm its samples, and then choosing from among all the calculated estimations the one which is optimal. The idea of the procedure comes from the BFS algorithm (Cormen et al. 2009), which, layer by layer, traverses the graph of interest. In the presented approach, the estimates for the next layer are computed on the base of the previous layer, assuming that vertex repetition in different steps is allowed. The algorithm stops when either the appropriate layer is reached or if the algorithm converges for the given G, then it stops when the distance in the sense of the chosen metrics (see Eq. 47) between the subsequent estimates for elements in C U is smaller than some desired . If the algorithm does not converge (when G is fixed) then the number of steps that need to be made is explicitly set at the beginning of the estimation procedure. In the case of standard complete PC matrices the graph is a directed clique, where every single vertex is connected to each other. In such a case it appears that even in the first step all the vertices are visited, whilst in the second step the computed estimations may take into account all possible values 2 gathered in M except those describing ratios between elements in C K . Thus, when the algorithm is not convergent, a few iterations of the HRE procedure seem to be a reasonable choice.
The main principle of the algorithm (Listing 1) seems to be quite straightforward. It starts from assigning followers (in sense E) of all elements of C K to the set L (Listing: 1, line: 2). For each u, for which μ(u) is unknown, all its predecessors in E are scanned. If μ(v) is already known then v becomes an element taken into account during the mean computation (Listing: 1, line: 10). Then, within the two loops while the current level is traversed (the outer loop Listing: 1, lines: 6-20), and appropriate estimates are computed (the inner loop Listing: 1, lines: 8-12). The predicate stop becomes false when either the auxiliary variable l reaches the assumed number of levels or it holds that l > 1 and ρ(x l , x l+1 ) < , where x l = (μ old (c u 1 ), . . . , μ old (c u k )) and x l+1 = (μ(c u 1 ), . . . , μ(c u k )) where the concepts c u i ∈ C U and ρ is one of the metrics defined in (Eq. 47). Otherwise it is true, and the outer loop continues traversing G. As long as L is not empty the inner loop proceeds as follows: removes one element from the set L and assigns it to the variable u (Listing: 1, line: 9), forms the auxiliary set T containing input edges for u so that their beginnings are predecessors of u in E for which μ is already known, and for all u ∈ L adds to the auxiliary mapping tmp : C U → R + a pair (u, μ(v)M(v, u)) (Listing: 1, line: 11). When all the elements of L are processed then the inner loop ends. Next the outer loop rewrites the auxiliary mapping tmp to the result map μ (Listing: 1, lines: 13-14). Since for all the concepts that were previously in L the function μ is known, the next level is devoted to μ calculation for their followers in L, hence the set L is filled back (Listing: 1, lines: 15-16). At the end of the outer loop the step counter is incremented and the auxiliary mapping tmp is emptied. When the outer loop completes its operation, the set Est contains a sequence of subsequent sets of estimates for concepts from C U . Then, at the end of the procedure (Listing: 1, line: 21) an optimal set of estimates needs to be chosen. For the purposes of this work it was assumed that the optimal set of estimates is one for which the average of the absolute mean estimation errors is minimal (Eq. 4). Let us define the absolute mean error e(u) as an error indicator for some concept u and mapping μ as follows: where u ∈ C U , (v i , u) ∈ E and μ(v i ) is already defined. Then the average error for all unknown concepts with respect to μ is defined as: In every iteration of the algorithm every weight M(v, u) (Listing: 1, line: 11) is considered exactly once. Hence, the running time of the outer loop is where r is the number of iterations, and n is the size of C. Similarly, the calculation of the average error (Eq. 3) requires the consideration of every M(v, u) exactly once, hence computing formulae (Eq. 4) also requires at most r · n 2 steps. Thus, the overall running time of the HRE algorithm is also O(r · n 2 ).

Numerical example
Let us illustrate the algorithm defined above by a simple numerical example. Some state agency supports five innovative projects u 1 , . . . , u 5 . After a while the two of them u 4 and u 5 come to an end and their actual cost becomes known μ(u 4 ) = 6, μ(u 5 ) = 2. As both projects exceeded the initial budget, the agency wants to reestimate the expected costs of the other projects using the already acquired knowledge. For this purpose the agency organizes a panel of experts under which a PC matrix M is formed (Eq. 5) reflecting the predicted cost relations between all the five projects.
Since experts were assigned to the pairs (u i , u j ) randomly and did not know each other's estimates, the matrix M is neither reciprocal nor consistent. Based on this data the agency must calculate adequate estimates of the projects u 1 , u 2 and u 3 and take appropriate decisions regarding the future of these projects. The reasoning presented below shows how to calculate the estimated cost of the projects using a HRE procedure (Listing: 1).
The matrix M is used for the generation graph G, the set of known concepts C K is formed by the projects u 4 and u 5 , and finally the unknown concept set is C U = {u 1 , u 2 , u 3 }. Since all concepts are reachable in the first step, and during the second step the μ mapping is computed for every c ∈ C U using all other concepts in C, then for the practical demonstration of the algorithm the number of levels traversed is limited to 2. In fact, for the matrix M as given (Eq. 5) the subsequent estimations calculated by HRE converge. Thus, instead of limiting the number of steps, an appropriately small can be chosen. This case will be discussed later, after the convergence criterion for HRE is defined.
Let us assume that the first vertex considered during the first step of the algorithm is u 1 . Then, according to the presented procedure μ(u 1 ) is computed as follows: hence, and further, in the same way: The first turn of the outer loop (traversing the first level of G) brings estimations for all vertices in C U . Thus, although the estimation process can be terminated, the agency concludes that it is better to take into account more data (for instance, due to the excluding of accidental errors of some experts), then decides to perform the second iteration of the algorithm. The new subsequent values μ(u 1 ), μ(u 2 ) and μ(u 3 ) are calculated below.
Recognizing the achieved result as optimal, the agency finishes the algorithm. It is worth noting that in the matrix M (Eq. 5) the two estimates do not come from the experts and have been introduced just for the matrix completeness. Namely, since the values μ(v 4 ) and μ(v 5 ) were previously known, the values m 4,5 and m 5,4 were calculated as ratios μ(v 4 ) /μ(v 5 ) and μ(v 5 ) /μ(v 4 ) respectively.

Data consistency and estimation error
For any algorithm, the result of which being inaccurate, a natural questions concern the accuracy of the resulting solution, and what should be done to improve the accuracy. The presented procedure is based on the fundamental sample mean estimation equation (Walpole 2012), where for the purpose of the algorithm the product μ(u i ) · M(u i , u j ) is treated as a sample, whilst μ(u j ) denotes the expected value inferred from samples. In this case the natural measure of the algorithm output accuracy is an estimation error understood as the distance between sample and the mean. Of course, the smaller the error, the more accurate the result. According to the popular adage "garbage in, garbage out" it is expected that even the best algorithm is not able to provide good output if the input data are bad. Hence, it might be expected that the estimation errors of HRE depend on data consistency, and are smaller in correlation with data inconsistency. The reasoning below supports this assertion.

Theorem 1 For a complete PC matrix M, and PC graph G over M it holds that
where K (M) is Koczkodaj's distance based inconsistency index for M, e(u) is the mean absolute estimation error for u ∈ C U , and μ is any estimation provided by the HRE procedure. Proof Let us consider some element u ∈ C representing an initially unknown concept i.e. u ∈ C U . Let v 1 , . . . , v n be its predecessors in G so that (v i , u) ∈ E ∧i = 1, . . . , n, for which the value μ(v i ) is known. Since the matrix M is complete there must exist Moreover, let us denote edges between v i 's and u as a = M(v 1 , u), and c 2 = M(v 2 , u), . . . , c n = M(v n , u) (see Fig. 1). Following the Eq. 1, Koczkodaj's distance inconsistency index K (M), in short K , means that the maximal local inconsistence for some maximal triad of three vertices t 1 , t 2 , t 3 ∈ C is K . Thus, in the case of triads composed of the concepts v 1 , v i , u it must hold that: for all i = 2, . . . , n. This implies that one of the following two statements is true: Let us denote α d f = 1 − K then the statements above can be written in the form: Combining these two expressions (Eqs. 17 and 18) we obtain: Let us denote β 1 Since α ≤ 1 (see Eqs. 15 and 16) thus the statement (Eq. 20) implies: and of course During the first step of the algorithm the μ is defined only for known concepts v ∈ C K . Thus, all the concepts v 1 , . . . , v n are in C K since only such elements are taken into account when calculating the value of μ(u). It is assumed that the value of the ratio . Therefore the update equation for μ(u) (Listing: 1, line: 11) can be written in the form: Since every component β i is bounded (Eq. 21) their mean must also be within the same bounds, which leads to the conclusion that: The absolute estimation error for some u with respect to μ at the end of the first step of the algorithm (Listing: 1) can be written and bounded from above as follows: where k ∈ {1, . . . , n}. Because both components of the absolute difference on the right side of the expression 25 have the same lower and upper bounds (Eqs. 22 and 24), then the maximal possible distance between them is limited by the difference between their upper and lower bounds. Thus, then the absolute mean error for the purpose of traversing the first level of G is upper bounded by: The second level is a bit more complicated. Under the conditions of the algorithm there is u ∈ C U , and M is complete. Thus, there exists at least one known concept v ∈ C K which precedes u in E i.e. (v, u) ∈ E. Let us put v 1 = v. This means that during the first step of the algorithm either v i (for i = 2, . . . , n) was in C K thus obviously can be bounded using Eq. 24 (note that in order to use Eq. 24 a needs to be replaced by b i ). This leads to the following inequality: thus, Since the Eq. 21 is valid for each triad, hence and then, which provides the estimation for μ(u) for the purpose of the second step of the algorithm: Once again, the absolute error with respect to μ at the end of the second step can be written and upper bounded as follows: where k ∈ {1, . . . , n}. Then, similarly as in the first step, both components have the same upper and lower bounds (Eqs. 30 and 32), thus the distance between them must not be greater than the distance between these limits. So, which means that the absolute mean error in the second step of the algorithm is bounded from above as follows: Let us consider the r 'th step of the algorithm. Similarly as in the second step we Fig. 1). Let us assume by induction that every v i for i = 2, . . . , n is bounded as follows: (compare with Eq. 28). Then by repeating the same reasoning as for the second step (Eqs. 29-35) we come to the conclusion that: Due to the principle of induction the above inequality holds for eve r ∈ N + . Since every e r (u) is a sum of absolute values then it cannot be negative i.e.
0 ≤ e r (u) Moreover, the definition of α implies Thus, due to the ( 1 /α r − α r ) component on the right side of the inequality 37, when α → 1 then the whole right side of Eq. 37 also approaches to 0. Therefore, it holds that for every step r > 0 it is true that: The above assertion, in the light of the arbitrary choice of u ∈ C U , satisfies the thesis of the theorem.

Convergence of solution
One of the immediate questions to come up is about the optimal number of iterations. The answer is not obvious, since the result of the algorithm depends on the input data.
In particular, it is easy to construct simple graphs over the inconsistent PC matrix M in cases where every further step of the algorithm significantly increases the absolute error of estimation. For instance, the graph G such that G = ( or all are less than or greater than 1. In the first case the subsequent estimates μ(v 2 ) and μ(v 3 ) tend to 0, and reversely, if all the values are greater than 1, then μ(v 2 ) and μ(v 3 ) tends to ∞. In the general case, the inequality 37 (Th. 1) suggests that with every subsequent step errors may increase. This leads to the conclusion that the optimal strategy is "the fewer steps the better". Therefore, if the algorithm is not convergent, traversing at most one or two levels seems like a good idea. If there are some edges with weights below 1 and some edges with weights above 1, the behavior of the algorithm is not so obvious. During the conducted experiments, it turned out that very often the subsequent iterations produce series of estimations convergent to some fixed positive results. To explain this phenomenon let us write the algorithm in the form of an appropriate system of equations, describing the onstep calculations (for the second step and the following ones). For this purpose, let us assume that the first iteration of the procedure was performed, hence the values μ(v) are assigned to all the concepts v ∈ C. For simplicity, let us assume that C U = {c 1 , . . . , c k }, C K = {c k+1 , . . . , c n } and denote b i for all i = 1, . . . , k as Thus, during the second and subsequent iterations the algorithm calculates the new estimation value μ(c i ) for each unknown concepts c i ∈ C U according to one of the following equations: Let us denote: Then, the equation system takes the form: Let us define the operator T : R k → R k (see Bronstein et al. 2005) as follows: The considered equation system can be written in the form of a fixed point problem in the metric space R k : Let us assume one of the following metrics: (47) and x = (ξ 1 , . . . , ξ k ), y = (η 1 , . . . , η k ). It holds that if one of the following conditions Q 1 , Q 2 or Q 3 is less than 1, then T turns out to be a contracting operator (Bronstein et al. 2005). According to the Banach fixed point theorem, there exists only one such point. Hence, in such a case there is only one set of values μ(c 1 ), . . . , μ(c k ) which is the limit of the sequence of HRE estimations. Due to the fraction 1 /n−1 the smaller a i j is the larger C is and the smaller M(c j , c i ) is. Moreover, the conditions (Eq. 48) can be more easily met with fewer concepts in C U (that is because the summations embedded in the conditions (Eq. 48) include fewer elements). In other words, the estimation procedure has a high chance to be convergent if 1. The set C U is relatively small (C K is relatively large), 2. The estimated values μ(v) are similar.
Both of these conditions are quite intuitive and, in practice, are likely to be satisfied. The first of them reflects the natural desire to provide the experts with rather more than the lower number of known, reference concepts. The second corresponds to the common-sense observation that all the considered concepts should be similar to each other, because then, it is easy to compare them. In other words, the expert estimates are more reliable when the compared projects are more similar.
The convergence of the algorithm implies the convergence of the estimation error. Unfortunately, the limit towards which the estimation error tends is not necessarily the smallest possible error value. Hence, in order to minimize the error, the user needs to choose the best estimation from all the estimations generated during the course of the algorithm (Listing: 1, line: 21).

Discussion
The idea underlying the HRE algorithm is the assumption that experts work independently and try to do their job as best as they can. Thus, they may be wrong in a random manner, and if so, it makes sense to accept their estimates as a part of samples, and the expected value of the sample (the arithmetic mean) as the alleged value of the estimated μ for the unknown concept. The absolute average estimation error may indicate how good such an estimation is. Depending on the data, the algorithm may or may not converge. If the algorithm does not converge (subsequent estimates are getting bigger) or converges to 0, selecting from among the first few estimates the one with the smallest error seems to be the best choice. If the algorithm converges to a non 0 value, it is very often useful to generate subsequent estimates until the current estimate does not differ from the limit (fixed point) with some , then to choose a set of estimations (among those generated) for which the average estimation error is minimal. Of course, this is not the only possible strategy. For instance, if some subset of C U is particularly important, only errors for its elements may be taken into account while determining the optimal set of estimations.
However, the problem considered in the work relies on determining values of μ for elements from C U = {c 1 , . . . , c k } on the basis of concepts C K = {c k+1 , . . . , c n } using the matrix M, it can be reformulated as computing a consistent approximation of a matrix M where certain elements are fixed. Indeed, at the end of the presented algorithm all the concepts c ∈ C have assigned some values μ(c). Thus, defining In particular it holds that m 45 = m 45 = 3 and m 54 = m 54 = 1 3 . Although the presented considerations in (Sect. 6 and 7) assume that the matrix M is complete (see Def. 3), the HRE algorithm seems to work without this assumption. For incomplete matrices, the value μ(u) for u ∈ C U can be determined as long as a path exists in the pairwise comparisons graph G over M (Def. 4) between u and some element v ∈ C K . Hence, there is a chance that the HRE algorithm may support different pairwise comparisons techniques also when the input data are incomplete. In particular, it might be useful for the AHP (Saaty 1977) approach. The properties of the presented algorithm for the incomplete sets of input data will be the subject of future research.

Summary
The HRE algorithm presented here for computing estimations of initially unknown concepts C U using information about known concepts C K and the matrix M proposes a new approach to the pairwise comparisons method. It defines the intuitive algorithm of using the pairwise comparisons matrix M for determining the most probable values of unknown concepts (original stimulus) on the basis of the known concepts. The presented procedure iteratively generates sets of estimations, then chooses the set which has the smallest average absolute mean estimation error. According to the proven theorem, the size of the estimation error depends on K -Koczkodaj's distance based inconsistency index shows that the lower the inconsistency index, the lower the estimation errors. The theorem can be particularly useful when the number of iterations of the HRE algorithm is small. In such cases, it may in practice be used to estimate the size of the estimation errors.
For some input data sets, the subsequent estimations produced by the HRE algorithm converges. If this happens, the estimation errors also converge to some limit, thus the number of estimation sets produced by the HRE algorithm does not need to be limited to a few elements. In such a case the sets of estimations can be generated until the limit towards which the HRE results converge will not be close enough. Such a situation has also been addressed in the paper. The given conditions of convergence have an intuitive explanation and in many practical situations are likely to be met.
The HRE algorithm is suitable for any matrix with positive elements, i.e., even in those situations where the applicability of the classical eigenvector method can be limited (finding the largest absolute eigenvalue of a nonreciprocal matrix may be difficult). The presented algorithm remains open for much more data than can be stored in a single pairwise comparisons matrix. Due to the graph representation of the problem, multiple values defining ratios between the same pairs of concepts can be easily encompassed within the algorithm as multiple arcs between the same pairs of vertices.