On a Maximum Eigenvalue of Third-Order Pairwise Comparison Matrix in Analytic Hierarchy Process and Convergence of Newton’s Method

Nowadays, the analytic hierarchy process is an established method of multiple criteria decision making in the field of Operations Research. Pairwise comparison matrix plays a crucial role in the analytic hierarchy process. The principal (maximum magnitude) eigenvalue of the pairwise comparison matrix can be utilized for measuring the consistency of the decision maker’s judgment. The simple transformation of the maximum magnitude eigenvalue is known to be Saaty’s consistency index. In this short note, we shed light on the characteristic polynomial of a pairwise comparison matrix of third order. We will show that the only real-number root of the characteristic equation is the maximum magnitude eigenvalue of the third-order pairwise comparison matrix. The unique real-number root appears in the area where it is greater than 3, which is equal to the order of the matrix. By applying usual Newton’s method to the characteristic polynomial of the third-order pairwise comparison matrix, we see that the sequence generated from the initial value of 3 always converges to the maximum magnitude eigenvalue.


Introduction
Analytic hierarchy process (AHP) is a scaling method for priorities in the hierarchical structure [1,2]. When people (decision makers) make decisions, they should often treat multiple alternatives. Decision makers' main task is to define the priorities of alternatives.
Pairwise comparison is one of the most tangible ways to prioritize them. Let C i , i = 1, … , n be alternatives to be prioritized. One defines the value a ij in pairs between C i and C j . The value represents the priority of C i as compared to C j . In the context of AHP, one usually uses its value by Saaty's discrete scale from 1 to 9 and their reciprocals. It is based on the verbal expressions in Table 1 [1].
For instance, suppose that one wishes to determine his/her preference for fruits. In comparison of apples ( C i ) and oranges ( C j ), if he/she prefers apples to oranges moderately, a ij is assigned to 3 and a ji to 1/3. Logically the value a ji in pairs between C j and C i should be 1∕a ij . After all pairwise comparisons, one obtains the following so-called pairwise comparison matrix: In AHP, employing a pairwise comparison matrix enables the treatment of subjective judgment.
In the context of AHP, evaluating the principal (maximum magnitude) eigenvalues of pairwise comparison matrices is a key ingredient. In AHP, solving the following linear system is called by the eigenvector method: where max is the maximum magnitude eigenvalue and w is the associated eigenvector. The maximum magnitude eigenvalue defines the first and most popular consistency index of pairwise comparison matrices [2,3]. The associated eigenvector is utilized as a priority vector (weight vector) [1,2,4,5]. From the viewpoint of mathematical programming, the theoretical foundation of the eigenvector method was established by [6,7]. In this note, we focus on third-order pairwise comparison matrix. The third-order pairwise comparison matrix has importance in AHP. The reason is that it is a source of inconsistency. There are several research articles devoted to third-order pairwise comparison matrix; see [8][9][10] and references therein. In this note, we treat inconsistent third-order pairwise comparison matrix A. Inconsistency is characterized by det A > 0. It is clear that the existence of the real-number root of the characteristic equation because the polynomial to be considered is odd order. In Sect. 2, we will specify the area where the real-number root exists. The root is shown to be the maximum magnitude eigenvalue among all eigenvalues of the pairwise comparison matrix.
The analytic form of the maximum magnitude eigenvalue is known; see [9] and references therein. Nevertheless, in this note, we shed light on Newton's method to obtain the maximum magnitude eigenvalue. We will see mathematical nature of Newton's method for third-order comparison matrix stands out. In Sect. 3, we will construct a sequence by the usual Newton's method. We will take an initial value to be 3, which is equal to the order of the matrix. We rely on the favorable properties of inconsistent third-order pairwise comparison matrix; uniqueness of the real-number root and det A > 0 . We will theoretically prove that the sequence generated by Newton's method converges to the maximum magnitude eigenvalue. We also confirm the property of convergence by Newton's method by computational experiment in Sect. 4. Utilization of Newton's method is expected to extend the result to possible application for fourth or more order pairwise comparison matrix.

Existence of the Root
We consider the third-order pairwise comparison matrix: where a ij > 0 . Its characteristic polynomial has the following form [11]: We can calculate the determinant of the pairwise comparison matrix as follows.
The last inequality follows from the relationship between the arithmetic mean and the geometric mean. The equality holds if and only if a 12 a 23 a 13 = a 13 a 12 a 23 , which is equivalent to a 12 a 23 = a 13 . We say the matrix A is consistent if a 12 a 23 = a 13 . For the thirdorder pairwise comparison matrix, as we mentioned above, it is obvious that the consistency of A and det A = 0 is equivalent. If A is consistent, the characteristic a 12 a 13 1∕a 12 1 a 23 1∕a 13 1∕a 23 1 a 12 a 13 1∕a 12 1 a 23 1∕a 13 1∕a 23 1 . In this case, the roots of the characteristic equation are = 3, 0 (0 is a multiple root). Hereafter, we treat the case when A is inconsistent (i.e., det A > 0 ). Taking the derivatives of P A ( ) , we have P � A ( ) = 3 ( − 2) . So P A ( ) takes a local maximum at the point = 0 and it takes a local minimum at the point = 2 . Moreover, for > 2 , P � A ( ) = 3 ( − 2) > 0 , so P A ( ) is monotone increasing in this area. The vertical axis section is − det A < 0 . We indicate the shapes of the graph of P A ( ) in Fig. 1.
An easy calculation shows the following: We immediately have the following existence theorem of a root for the characteristic equation.  (3, 3 + det A) . This root is the maximum magnitude eigenvalue.

Theorem 1 For the inconsistent pairwise comparison matrix A, the characteristic equation has a unique real-number root in the interval
Proof By the intermediate value theorem (p.206 in [12]), the existence of a realnumber root in the interval (3, 3 + det A) is clear. From the shape of the graph, this root is the unique real-number root. We denote it by * . In general, other two conjugate complex-number roots exist, say a ± bi , b ≠ 0.
From the general theory of eigenvalues and the trace of the matrix [2], we have * + (a + bi) + (a − bi) = 3 . So we have

Convergence of Newton's Method
Here we use the usual Newton's method (Chapter 8 in [13]). We give an elementary proof of convergence. We set the initial value to be 0 = 3 and generate a sequence by the following iteration.
Denote the maximum magnitude eigenvalue guaranteed by Theorem 1 by max . So P A ( max ) = 0 and 3 < max < 3 + det A.
Lemma 1 For all n ≥ 1 , we have n > max .
Proof We prove it by induction. Set n = 1 .
So we have 1 > max . Assume n > max is true. Then from this assumption, we can take max ≤ < n . Since P � A ( ) is monotone increasing for > 1 (because P �� A ( ) = 6( − 1) ) and max > 3 , we have Thus we have So we conclude Taking account into P A ( max ) = 0 , we have Lemma 2 The sequence { n } is monotone decreasing for n ≥ 1.
Proof From Lemma 1, we have 3 < max < n . Since P A ( ) is monotone increasing for > 3 , we have P A ( n ) > P A ( max ) = 0 . Obviously P � A ( n ) > 0 and P A ( n ) > 0 , so we have From Lemmas 1 and 2, the sequence is monotone decreasing and bounded below. So it converges, see p.183 in [12].

Theorem 2
The sequence { n } converge to max .
Proof From Lemmas 1 and 2, the sequence is monotone decreasing and bounded below. So it converges to, say ̂≥ max . Taking the limit of the iteration , for n ≥ 0.
we have by the continuity of P A ( ) and P � A ( ) . Since ̂≥ max > 3 , the inequality P � A (̂) > 0 holds. So we have P A (̂) = 0 . Because the root of P A ( ) = 0 is unique, ̂ is identical to max . Now we state the rate of convergence.

Theorem 3
The sequence { n } converge to max quadratically.

Computational Experiment
We have a computational experiment to compare Newton's method with other rootfinding methods. The experiment consists of the following steps.
Step 1. Generate an inconsistent pairwise comparison matrix A of order 3, at random.
Step 2. Iterate the following three methods until the differences between max and the obtained values become below the threshold and store the numbers of iterations.
-Newton's method with the initial value of 0 = 3 . ) have the same sign, for n = 0, 1, 2. …. Figure 2 shows the distributions of the numbers of iterations until convergence by each method, after repeating these steps 5,000 times. We set the threshold of convergence to 10 −12 here.
The summary statistics of the 5,000 numbers of iterations until convergence by each method are shown in Table 2.
For example, for the pairwise comparison matrix of  Fig. 3.
This experiment confirms the speed of convergence of Newton's method empirically.

Conclusion
In this note, we entirely use the favorable properties from which A is third order. About other properties concerning third-order pairwise comparison matrices, see [9]. The key contribution of this note is findings of the possible area where the root of the characteristic equation and to show the convergence of Newton's method from the initial value 0 = 3 . We also confirm the superiority of Newton's method compared to other root-finding methods-the secant method and the bisection method. This is considered as a direct consequence of the quadratic convergence of Newton's method. For the matrices for fourth and more order dimensions, expanding this note is left to future research. -For fourth-order pairwise comparison matrix, there usually exist two roots of the characteristic equation. So, if we set the initial value 0 = 4 , Newton's method does not generally converge to max . -We need to study the appropriate initial value and the convergence property.