Network analysis with the aid of the path length matrix

Let a network be represented by a simple graph $\mathcal{G}$ with $n$ vertices. A common approach to investigate properties of a network is to use the adjacency matrix $A=[a_{ij}]_{i,j=1}^n\in\R^{n\times n}$ associated with the graph $\mathcal{G}$, where $a_{ij}>0$ if there is an edge pointing from vertex $v_i$ to vertex $v_j$, and $a_{ij}=0$ otherwise. Both $A$ and its positive integer powers reveal important properties of the graph. This paper proposes to study properties of a graph $\mathcal{G}$ by also using the path length matrix for the graph. The $(ij)^{th}$ entry of the path length matrix is the length of the shortest path from vertex $v_i$ to vertex $v_j$; if there is no path between these vertices, then the value of the entry is $\infty$. Powers of the path length matrix are formed by using min-plus matrix multiplication and are important for exhibiting properties of $\mathcal{G}$. We show how several known measures of communication such as closeness centrality, harmonic centrality, and eccentricity are related to the path length matrix, and we introduce new measures of communication, such as the harmonic $K$-centrality and global $K$-efficiency, where only (short) paths made up of at most $K$ edges are taken into account. The sensitivity of the global $K$-efficiency to changes of the entries of the adjacency matrix also is considered.


Introduction
An important characteristic of a network is how well communication can flow in it, i.e., how easy or difficult it is to reach one part of the network from another part by following edges.How well information flows through the whole network can be measured by the diameter of the graph that represents the network or by its global efficiency.Both these measures can be determined with the aid of the path length matrix associated with the network.We will discuss these connections and introduce new measures of communication based on the path length matrix.
Let us introduce some notation and definitions that will be used throughout this paper.A network is represented by a graph.A weighted graph G = V, E, W consists of a set of nodes or vertices V = {v 1 , v 2 , . . ., v n }, a set of edges E = {e 1 , e 2 , . . ., e m } that connect the vertices, and a set of weights W = {a ij } n i,j=1 ; the weights a ij are the entries of the adjacency matrix A = [a ij ] n i,j=1 associated with the graph G; see below.An edge is said to be directed if it starts at a vertex v i and ends at a vertex v j , and is denoted by e(v i → v j ).An edge between the vertices v i and v j is said to be undirected when the pair of vertices is unordered and the weights a ij and a ji are positive and equal.An undirected egde between the vertices v i and v j is denoted by e(v i ↔ v j ).A graph with only undirected edges is said to be undirected ; otherwise the graph is directed.A simple graph is a graph without multiple edges or self-loops.In particular, this implies that the diagonal entries of the adjacency matrix for the graph vanish.This work considers simple graphs.
The adjacency matrix A = [a ij ] n i,j=1 ∈ R n×n for a weighted graph G is determined by the weights a ij of the graph with a ij > 0 if there is an edge e(v i → v j ) in G.If there is no edge e(v i → v j ) in G, then a ij = 0.For an unweighted graph, all positive entries a ij of A equal one.A sequence of k edges (not necessarily distinct) such that {e(v 1 → v 2 ), e(v 2 → v 3 ), . . ., e(v k → v k+1 )} form a walk.If v k+1 = v 1 , then the walk is said to be closed.A sequence of distinct edges such that {e(v 1 → v 2 ), e(v 2 → v 3 ), . . ., e(v k → v k+1 )} form a path.The length of a path is given by the sum of all weights of the edges in the path.(In the unweighted case, the sum of all weights of the edges in a path of length k is k.)For further discussions on networks and graphs; see [8,14].
To construct the path length matrix associated with the network, we will make use of min-plus matrix multiplication, i.e., matrix multiplication in the tropical algebra [12]: where A, B, C ∈ R n×n .We denote by A 1,⋆ = [a ] n i,j=1 ∈ R n×n the matrix obtained by setting to ∞ the vanishing off-diagonal entries of the adjacency matrix A associated with the graph G under consideration.For k > 1, the k th min-plus power of A 1,⋆ is given by ] n i,j=1 ∈ R n×n : a hj }.
Notice that the matrix A k,⋆ gives vertex distances using paths of at most k edges.In detail, the entry a , with i = j, represents the length of the shortest path from v i to v j made up of at most k edges.The diagonal entries of A k,⋆ are zero by definition.One has a The diameter of a graph G is the maximal length d G of the shortest path between any distinct vertices of the graph and provides a measure of how easy it is for the vertices of the graph to communicate.One has Indeed, the entry a ] n i,j=1 yields the length of the shortest path from v i to v j .We will refer to A n−1,⋆ as the path length matrix.Note that the triangle inequality holds for the entries of a path length matrix.Specifically, Consider an unweighted connected graph G with associated adjacency matrix A ∈ R n×n .Recall that the diameter of G is the maximal number of edges in the shortest path between all pairs of distinct vertices of the graph.Given the vertices v i and v j , there is an integer k, 1 ≤ k < n, such that since the graph is connected, whereas for 1 ≤ h < k, one has a (h,⋆) ij = ∞.Thus, information provided by the path length matrix A n−1,⋆ includes information about all powers A k,⋆ for 1 ≤ k < n − 1.
Let v i and v j be distinct vertices in a weighted graph.Then there is an integer k, 1 ≤ k < n − 1, such that a Thus, as in the unweighted case, information provided by the path length matrix A n−1,⋆ refines information given by the powers A k,⋆ for 1 ≤ k < n − 1.However, the information of the minimal number of steps required to reach vertex v j from vertex v i is lost.
As mentioned above, the path length matrix may be constructed by evaluating the min-plus powers of A 1,⋆ n − 2 times; here A 1,⋆ is obtained from the adjacency matrix A by setting all zero off-diagonal entries to ∞.The following MATLAB function, with the adjacency matrix A for a graph and the level = n − 1 as input arguments, returns the path length matrix associated with the adjacency matrix A. The function implements the dynamic programming Bellmann-Ford algorithm for solving the well-known "all-pairs shortest path problem".The algorithm requires O(n 2 m) arithmetic floating point operations (flops), where n is the number of vertices and m is the number of edges of the graph; if the graph is undirected, then the cost of the algorithm is halved.Notice that the function can be applied to determine shortest paths in a weighted graph having positive or negative weights.
n = size(A, 1); 3: for i = 1 : n do 4: for j = 1 : n do 5: if A(i, j) = 0 and i = j then 6: A(i, j) = inf; 7: end if 8: end for 9: end for 10: for count = 1 : level − 1 do 13: for i = 1 : n do 14: for j = 1 : n do 15: if j = i then 16: for k = 1 : n do 17: if B(k, j) = inf and k = j then 18: In line 2 of the function Pathlength matrix, the function call size(A, 1) yields the order of the matrix A ∈ R n×n ; in line 6, the matrix entry A(i, j) is set to ∞.Similarly, in line 11, C is defined as an n × n matrix with all entries equal to ∞.If the graph is unweighted and connected, then the above MATLAB function can be modified by introducing a break before the last end for when there is no entry ∞.The diameter of the graph then is count+1.
For both weighted and unweighted graphs, also when the graph is not connected, the above MATLAB function with argument level = K, where 1 < K < n, computes A K,⋆ , i.e., the matrix of the distances between any distinct vertices of the graph using paths with at most K edges.We note that the triangle inequality might not hold for the entries of A K,⋆ and some entries of this matrix may have the value ∞.
We will see how the matrix A n−1,⋆ associated with a connected graph G sheds light on the communication within the network determined by the graph.In fact, as a measure of the ease of communication in the graph, we like to use the average inverse geodesic length of G (i.e., its global efficiency, cf.Section 3) instead of the maximum geodesic length of G (i.e., its diameter).To this end, we introduce the reciprocal path length matrix ] n i,j=1 obtained by replacing the off-diagonal entries of the path length matrix by their reciprocals, i.e., a where 1/∞ is identified with 0.
We are interested in determining the shortest paths that use at most K edges.We therefore also consider the reciprocal K-path . Thus, the entry Note that the same would happen if K = n − 1, in case the graph G that determines the adjacency matrix A is not connected.The matrix A K,⋆,−1 allows us to define the global K-efficiency of G; see Section 3.
In order to enhance communication using paths with at most K edges, with 1 < K < n, i.e., to increase the global K-efficiency of the graph associated with the adjacency matrix A, we select edge weights by analyzing centrality properties of the vertices of the graph and, if computationally feasible, the spectral properties of the reciprocal K-path length matrix A K,⋆,−1 .In detail, if K ≪ n and the (sparse) non-negative reciprocal K-path length matrix A K,⋆,−1 is irreducible, then we apply the Perron-Frobenius theory by following the approach in [4,7].In our context, the choice of weights is dictated by the analysis of the sensitivity to perturbations in A.
Applications of our approach include city planning and information transmission.As for disease propagation, a recent research study of the Zhejiang City Planning Center (China) pointed out a strong connection between the spread of the Covid-19 epidemic and the shape of the city.Cities with a radial structure (such as Milan) have good internal connections due to the capillarity of the public transport system -buses, trams, subways, and trains.This dynamic made citizens of such cities more vulnerable to the arrival of the Covid-19 virus: the incidence of infections compared to the number of inhabitants was generally larger in cities with a radial urban structure than in cities without this structure, due to the good communication of the people in cities with radial urban structure.This paper is organized as follows: Section 2 analyzes differences and similarities of powers and tropical powers of the adjacency matrix for undirected and unweighted graphs.Section 3 reviews well-known measures that can be easily computed by means of the path length matrix and introduces novel ones.In Section 4 we present two algorithms that determine which edge-weight should be changed in order to boost global efficiency.Changing the edge weights may entail widening streets or increasing the number of trams on a route, decreasing travel times on a highway by increasing the travel speed, or decreasing the waiting time for trams on a route.Finally, numerical tests are reported in Section 5 and concluding remarks can be found in Section 6.

Powers versus tropical powers
Consider an undirected and unweighted simple graph G with adjacency matrix A ∈ R n×n .Then the entry a ∈ R n×n counts the number of walks of length k between the vertices v i and v j .A matrix function based on the powers A k that is analytic at the origin, and vanishes there, can be defined by a formal Maclaurin series where we for the moment ignore the convergence properties of this series.Usually long walks are considered less important than short walks, because information flows more easily through short walks than through long ones.Therefore matrix functions applied in network analysis generally have the property that 0 ≤ c k+1 ≤ c k for all k ≥ 1.The most common matrix function used in network analysis is the matrix exponential; see [5,6,8,9,10] for discussions and illustrations.We prefer to use the the modified matrix exponential where I denotes the identity matrix, because the first term in the Maclaurin series of exp(A) has no natural interpretation in the context of network modeling.For the modified matrix exponential, we have c k = 1/k!, and the series (2) converges for any adjacency matrix A.
The communicability between distinct vertices v i and v j , i = j, is defined by see [8] for the analogous definition based on exp(A).The communicability accounts for all possible routes of communication between the vertices v i and v j in the network defined by the adjacency matrix A, and assigns a larger weight to shorter walks than to longer ones.The larger the value of [exp 0 (A)] ij , the better is the communicability between the vertices v i and v j .
Remark 1 Notice that even if there exists an integer k, 1 ≤ k < n − 1, such that a ij > 0, one may have a ( k+1) ij = 0. Information provided by A n−1 does not include information provided by all A h for 1 ≤ h < n − 1.This is one of the reasons for the interest in the matrix functions exp 0 (A) and exp(A).
It is straightforward to show the following result.
Then k is the length of the shortest path that connects v i and v j , i.e., a and p is the number of shortest paths that connect v i and v j .
Example 1 In view of Proposition 1, the information provided by both the path length matrix and a suitable power of the adjacency matrix is of interest.Consider the undirected and unweighted graphs G 1 and G 2 depicted in Figures 1-2.The adjacency matrices of these graphs are For both graphs, the shortest path between the vertices v 1 and v 2 has length 2. However, in G 1 there are three shortest paths that connect v 1 and v 2 , while there is only one in G 2 .Thus, there is surely better communication between these vertices in G 1 than in G 2 , even though this information is not provided by the path length matrices for these graphs.The path length matrices are On the other hand, the second powers of the above adjacency matrices are Remark 2 For weighted graphs, the interpretation of the entry a (k) ij of the matrix A k has to be modified.Indeed, a (k) ij yields the sum of all products of all weights of the edges in the walks of length k between the vertices v i and v j .(In the unweighted case, any product of all weights of the edges in a walk is 1.) Remark 1 holds true, while Proposition 1 does not.

Measures that depend on the path length matrix
Let for now the graph G be undirected.

Closeness centrality
Let the graph G be connected.Then the reciprocal of the sum of all lengths of the shortest paths starting from vertex v i , where 1 i ∈ R n denotes the vector with all zero entries except for the i th entry, which is one, is referred to as the closeness centrality of v i in G; see, e.g., [2].This measure gives a large value to vertices that have small shortest path distances to the other vertices of the graph.

The radius and center of a network
The maximum length over all shortest paths starting from vertex v i in a weighted or unweighted connected graph G, that is to say, , is commonly referred to as the eccentricity of the vertex v i .The diameter may be seen as the maximum eccentricity among the vertices of the network.The radius r G of G is the minimum eccentricity of a vertex.One has A vertex is said to be central if its eccentricity is equal to the radius of the graph.The center of the graph is the set of all central vertices; see, e.g., [8].

Average shortest path length
The average shortest path length of a weighted or unweighted connected graph G computed over all possible pairs of vertices in the network [2] is given by where 1 ∈ R n denotes the vector with all entries one.Let the graph G be unweighted and be formed by a path of n vertices.Then the largest average shortest path length is a G = (n + 1)/3.If a G scales logarithmically with n, then G displays the "small-world phenomenon"; see, e.g., [8].

Harmonic centrality and global efficiency
The efficiency of a path between any two vertices of a weighted or unweighted graph G is defined as the inverse of the length of the path.The sum of the inverses of the length of all shortest paths starting from vertex v i , i.e., the sum of the efficiencies of all shortest paths starting from v i , is referred to as the harmonic centrality of v i ; see, e.g., [2].The latter measure gives a large centrality to vertices v i that have small shortest path distances to the other vertices of the graph.Harmonic centrality will play a central role in our analysis.We define the h-center of a graph as the set of all vertices with the largest harmonic centrality.
If the graph is connected, then the average shortest path efficiency over all possible pairs of vertices, also known as the average inverse geodesic length, is referred to as the global efficiency of the graph [2]: Finally, we remark that in the context of molecular chemistry, the sum of reciprocals of distances between all pairs of vertices of an undirected and unweighted connected graph is known as the Harary index and the reciprocal path length matrix A n−1,⋆,−1 is referred to as the Harary matrix [17].

Remark 3
The measures (4) and ( 5) can be especially useful when the network has more than one connected component, because infinite distances do not contribute to these "harmonic" averages.
Example 2 As an illustration of the above measures, consider again the graphs G 1 and G 2 of Example 1. Table 1 reports the diameter, radius, average shortest path length, and global efficiency of these graphs.Notice that the center of G 1 is given by the set of all vertices and the center of G 2 is made up of the vertex v 3 , only.The h-center of G 1 is formed by the vertices v 1 and v 2 , and the h-center of G 2 is given by the vertex v 3 , only.Table 2 shows the eccentricity, harmonic centrality, and closeness centrality of all the vertices of G 1 and G 2 .All the measures in this example are computed by using the path length matrices (3).

Harmonic K-centrality and global K-efficiency
When one considers shortest paths that are made up of at most K edges, the matrix A K,⋆,−1 takes the role of the reciprocal path length matrix A n−1,⋆,−1 .We are in a position to introduce the harmonic K-centrality of the vertex v i .It is given by The global K-efficiency of a graph G is defined by The set of all vertices of a graph G with the largest harmonic K-centrality is referred to as the h K -center of the graph.

Out-centrality versus in-centrality
Let the graph G be directed.Then the above centrality measures (closeness, eccentricity, harmonic centrality, and harmonic K-centrality) determine the importance of a vertex v j by taking into account the paths that start at v j .These measures therefore may be considered measures of out-centrality.One also may be interested in measuring the importance of a vertex v j by considering the paths that end at it, that is to say by measuring the in-centrality of v j .This can be achieved by replacing the path length matrix A n−1,⋆ in the measures mentioned by its transpose.This allows us to introduce the measures closeness in-centrality, in-eccentricity, harmonic in-centrality, and harmonic K in -centrality.These measures are defined as We also define the in-radius of G, Finally, the vertex v j is said to be in-central if its in-eccentricity equals the in-radius of the graph.The h Kin -center of a graph G is the set of all vertices with the largest harmonic K in -centrality.The notions of h Kin -center and h Kout -center will be of interest in the sequel.

Enhancing network communication
The diameter of a weighted or unweighted graph provides a measure of how easy it is for the vertices of the graph to communicate.When graphs are used as models for communication networks, the diameter plays an important role in the performance analysis and cost optimization.A simple way to decrease the diameter of a graph so that information can be transmitted more easily between vertices of the graph [4] is to decrease the weight of an edge that belongs to all maximal shortest paths, if feasible.
Example 3 Consider the graph G 2 in Example 1.The edge e(v 2 ↔ v 3 ) belongs to the maximal shortest path and is a bridge, that is its removal would make the vertex v 2 unreachable and the perturbed graph G2 disconnected.In detail, the edge e(v 2 ↔ v 3 ) belongs to the shortest path between the vertices v 1 and v 2 , which is the only shortest path of length d G .Thus, decreasing the weight of the edge e(v 2 ↔ v 3 ) decreases both the average shortest path length and the diameter of the graph.Specifically, if one decreases the weights a 23 and a 32 from 1 to 0.5, one obtains the matrices The perturbed graph G2 has diameter d G2 = 1.5, radius r G2 = 1, average shortest path length a G2 = 1, and global efficiency e G2 = 1.22.The eccentricity, harmonic centrality, and closeness centrality of the vertices of the graph G2 are reported in Table 3.

Increasing the global K-efficiency
We propose two approaches to increase the global efficiency of a network.

The function eKG1
Let for now the graph G be directed.The first approach is based on the observation that the most important vertices with respect to the global K-efficiency live in the vertex subsets h Kout -center and h Kin -center of the graph.These vertices may be interpreted as important intermediaries, that quickly collect information from many vertices and quickly broadcast it to many others vertices.Indeed, strengthening an existing connection from a vertex of the h Kin -center to a vertex of the h Kout -center is likely to strengthen their communicability by having new shorter paths with at most K steps that exploit these connections.This is likely to increase the global K-efficiency more than strengthening an existing connection between vertices with lower harmonic K in -and K out -centrality.Here we consider graphs whose edge weights represent travel times or waiting times.Hence, strengthening is achieved by decreasing appropriate weights.
We construct the perturbed adjacency matrix If the graph is undirected, the above approach simplifies, because the h Kout -and h Kin -centers coincide.Hence, the idea is to strengthen the connection between vertices with the largest harmonic K-centrality.The perturbed adjacency matrix Ã will be The MATLAB function eKG1 describes the necessary computations.The operator == in line 5 of the function eKG1 stands for logical equal to, and the symbol ./ in line 6 denotes element-wise division.The function call sum(M ) for a matrix M ∈ R n×n computes a row vector m ∈ R n , whose j th component is the sum of the entries of the j th column; the function call sum(M, 2) for a matrix M ∈ R n×n computes a column vector, whose i th entry is the sum of the elements of row i of M .The blip in line 8 denotes transposition, and the operator .* in line 11 stands for vector-vector element-wise product.

The function eKG2
Assume for now that the graph G is directed.If computing the Perron root ρ K and the unique positive left and right eigenvectors of unit norm (the Perron vectors) of the reciprocal K-path length matrix A K,⋆,−1 is not computationally feasible or if this matrix is not irreducible, then one can use the function eKG1.However, if A K,⋆,−1 is irreducible and its left and right Perron vectors x K = (x K,i ) and y K = (y K,i ) can be computed, then these vectors determine the Wilkinson perturbation W K = y K x T K ; see [20,Section 2].Following [7], to induce the maximal perturbation in ρ K , one chooses the indices (h 1 , h 2 ) such that W K (h 1 , h 2 ) is the largest entry of W K and A(h 1 , h 2 ) > 0, i.e., the indices of the largest entry of the Wilkinson perturbation "projected" onto the zero-structure of A; see, e.g.[15].Thus, As in function eKG1, one strengthens the edge e(v h1 → v h2 ) by halving its weight.The perturbed adjacency matrix then is given by (6).
When the graph G is undirected, the left and right Perron vectors coincide and the perturbed adjacency matrix is constructed as in (7).The outlined approach is implemented by the MATLAB function eKG2.We recall that the function abs(v) of a vector v, used in the function eKG2, returns a vector, whose components are the absolute value of the components of v.The expression sum(sum(P r)) on line 7 sums all entries of the matrix P r.The operator .* on line 20 denotes the Hadamard product of two matrices; the entries of the matrix A > 0 are one if the corresponding entry of A is positive; they are zero otherwise.The function ind2sub determines the equivalent subscript values corresponding to a given single index into an array.
We expect the global K-efficiency to increase the most when decreasing the edge-weight that makes the Perron root ρ K change the most.The perturbation of the Perron root generated by the function eKG2 typically is larger than the perturbation determined by the function eKG1.The difference in these perturbations is analogous to the difference between considering the most important vertex in a graph the one with the largest degree and the one with maximal eigenvector centrality.

Numerical tests
The numerical tests reported in this section have been carried out using MATLAB R2022b on a 3.2 GHz Intel Core i7 6 core iMac.The Perron root and left and right Perron vectors for small to moderately sized graphs can easily be evaluated by using the MATLAB function eig.For large-scale graphs these quantities can be computed by the MATLAB function eigs or by the two-sided Arnoldi algorithm, introduced by Ruhe [18] and improved by Zwaan and Hochstenbach [21].
Example 5 Consider the adjacency matrix for the network Air500 in [1].This data set describes flight connections for the top 500 airports worldwide based on total passenger volume.The flight connections between airports are for the year from 1 July 2007 to 30 June 2008.The network is represented by a directed unweighted connected graph with 500 vertices and 24009 directed edges.The vertices of the network are the airports and the edges represent direct flight routes between two airports.
The path length matrix A 5,⋆ = A 499,⋆ yields the diameter and the radius of the graph 5 and 3, respectively.The information provided by the vector of the harmonic centralities and the Perron vector for the reciprocal path length matrix A 5,⋆,−1 is the same as the one given by the vector of harmonic K-centralities and the Perron vector for A K,⋆,−1 with K = 2; cf.Table 4. Therefore, the perturbation that increases the global K-efficiency the most also will enhance the global efficiency the most.The information provided by Table 4 suggests that the number of flights from the Frankfurt FRA Airport (vertex v 161 ) to the JFK Airport in New York (vertex v 224 ) should be doubled in order to half the wait time between these flights.Doubling the number of flights corresponds to halving the weight for the corresponding edge.Example 6 This example considers an undirected unweighted connected graph G that represents the German highway system network Autobahn.The graph is available at [1].Its 1168 vertices are German locations and its 1243 edges represent highway segments that connect them.
Let A be the adjacency matrix associated with G.The path length matrix A 62,⋆ = A 1167,⋆ shows that the diameter and the radius of G are 62 and 34, respectively, whereas its global efficiency equals 6.7175 • 10 −2 .One notices that there is only one shortest path of length 62, which connects the vertices v 116 and v 1154 .The diameter of the graph can be decreased by halving the weight of the edges a 120,116 and a 116,120 (since the graph is undirected), because this is the unique edge that connects the vertex v 116 to the other vertices of the graph.
Let Â denote the perturbed adjacency matrix and Ĝ the corresponding graph.The global efficiency of Ĝ is 6.7177 • 10 −2 and the diameter is 61.5.We turn to the application of the functions eKG1 and eKG2 to increasing the global efficiency.Table 5 reports the global K-efficiency (for several values of K) of the graph G associated with the adjacency matrix Ã obtained by halving both entries a h1,h2 and a h2,h2 of A, computed by the function eKG1.Table 6 shows the global K-efficiency (for several values of K) of the graph G associated with the adjacency matrix Ã, computed by the function eKG2.Also in this example, the information provided by the vector of harmonic centralities and the Perron vector for the reciprocal path length matrix A 62,⋆,−1 is exactly the same information that is provided by the vector of harmonic K-centralities and the Perron vector for A K,⋆,−1 with K ≥ 4. We note that the latter vectors are less expensive to determine than the former.The information provided by both tables suggests that one should double the width of the highway that connects the cities of Duisburg (vertex v 219 ) and Krefeld (vertex v 565 ) to half the travel time.These cities are 10 miles apart.Doubling the width corresponds to halving the weight associated with the corresponding edge.

Concluding remarks
The adjacency matrix of a graph is a well-known tool for studying properties of a network defined by the graph.The path length matrix associated with a graph also sheds light on properties of the network, but so far has not received much attention.A review of measures that can be defined in terms of the path length matrix is provided, and new such measures are introduced.The sensitivity of the transmission of information to perturbations of the entries of the adjacency matrix is investigated.

2
even though other choices of existing edges are equally valid.The global 2-efficiency e 2 G1 of Ã1 is 0.95; compare with the global 2-efficiency e 2 G1 = 0.80 of A 1 .We remark that a perturbation of another existing edge would have led to the same increase of the global 2-efficiency.Consider the graph G 2 in Example 1.One has A

Table 1
Example 2. Diameter, radius, average shortest path length, and global efficiency of the graphs G 1 and G 2 depicted in Figures1-2.Note that the average shortest path length of G 2 is maximal; indeed, (n + 1)/3 = 1.33.

Table 2 Example
2. Eccentricity, harmonic centrality, and closeness centrality of the vertices of G 1 (left table) and of G 2 (right table).

Table 3
Example 3. Eccentricity, harmonic centrality, and closeness centrality of the vertices of G2 .

Table 4
Example 5. Indices chosen by the functions eKG1 and eKG2 and the global K-efficiency of both the given graph, G, and the perturbed graph, G, for K = 2, 3, 4, 5.