## Sampling on networks: estimating eigenvector centrality on incomplete graphs

Nicolò Ruggeri 1,2, * and Caterina De Bacco 1, †

1 Max Planck Institute for Intelligent Systems, Max-Planck ring 4, Tübingen, 72076, Germany 2 Università di Padova, Department of Mathematics 'Tullio Levi-Civita', via Trieste 63, Padova, 35100, Italy

We develop a new sampling method to estimate eigenvector centrality on incomplete networks. Our goal is to estimate this global centrality measure having at disposal a limited amount of data. This is the case in many real-world scenarios where data collection is expensive, the network is too big for data storage capacity or only partial information is available. The sampling algorithm is theoretically grounded by results derived from spectral approximation theory. We studied the problem on both synthetic and real data and tested the performance comparing with traditional methods, such as random walk and uniform sampling. We show that approximations obtained from such methods are not always reliable and that our algorithm, while preserving computational scalability, improves performance under different error measures.

## I. Introduction

leads to sensible improvements for several network structures, in particular for undirected ones.

One key task in inference in networks is to evaluate quantities when network information is limited as when the available network is only a subsample of the full one. Such evaluations impacts problems like modeling dynamical processes [1, 2], network statistics estimation [3], data compression [4] and survey design [5]. The sampling procedure used to collect network information usually influences these estimates [6-8], it is therefore important to choose a sampling algorithm properly based on what is the quantity one seeks to estimate. The challenge is to provide a reliable estimate while keeping the sampling algorithm scalable. While several sampling strategies can be found in the literature [3], these are usually empiricallydriven and lack of analytical groundings .

Here we take a theoretically motivated approach to approximate in-sample eigenvector centrality [9], a network centrality measure that is widely used across domains and that is highly sensitive to missing information compared to other standard centrality measures [10]. To set up the model, we minimize a proper distance between the centrality vector estimated from the subsample and the ground truth one. This allows to preserve the relative ranking of nodes inside the sample, ordered by their estimated centralities.

To tackle this problem, we develop a novel online sampling algorithm. The idea is that a proper sampling procedure should be capable of selecting nodes in an online fashion. We argue that, in order to achieve scalability, only local information should be used. This is challenging, as eigenvector centrality is considered a global centrality measure, i.e. its calculation requires global network information.

Using results derived from numerical methods for eigenvalues approximation [11], our sampling algorithm provides valid centrality estimates. We derive an effective criterion based on theoretical bounds that we develop and adapt to our case.

Evaluation is performed on both synthetic and real data and performance is measured in terms of recovering the relative ranks of nodes inside the sample. We found that our algorithm

* nicolo.ruggeri@tuebingen.mpg.de

† caterina.debacco@tuebingen.mpg.de

Related work. The problem of reliable network sampling has been treated in the literature mainly using empirical approaches [10, 12] and focusing on recovering standard topological properties like degree distribution, diameter or clustering coefficient [3, 13-18]. To the best of our knowledge, sampling has not been previously applied to estimating eigenvalue centrality, despite the vast amount of literature in this area and the number of applications that use this centrality measure. In general, few theoretical results have been developed in sampling networks, making ours one of the few method that goes beyond heuristics or empirical reasoning. In the case of PageRank score [19] (which can be thought of as the eigenvector centrality for a specific class of adjacency matrices) work has been done towards estimating single nodes' scores or approximating the external information missing for reliable estimation of the scores in the sample [20, 21].

## II. The problem

Consider the adjacency matrix A of a network G = ( /Vcal , /Ecal ) of size |/Vcal | = V , where /Vcal = { 1,... , V } is the set of nodes and /Ecal the set of edges, Aij ∈ /Rbbb ≥ 0 is the weight of the edge from j to i . Notice that, while in practice many real network applications consider only discrete, and possibly binary, entries, our results are valid in general for real-valued non-negative edge weights. We define d G in ( ) i : = ∑ j ∈/Vcal Aij the weighted indegree of node i . The problem of sampling the network under a particular sampling technique can be viewed as selecting a principal submatrix A m ′ of size m ≤ V , induced by selecting a subset of nodes /Ical ⊆ /Vcal , where |/Ical | = m . The subsample network is denoted as Gm =( /Ical , /Ecal m ) , and /Ecal m ⊆/Ecal is the set of edges in the subsample. Ideally, the subnetwork Gm should be selected as to give a good representation of the entire network G . This means that network-related quantities estimated on Gm should be a good approximation for the same quantities estimated on G . In our case we are interesting in estimating the eigenvector centrality, a popular centrality measure that is been used in several contexts [9]. The centrality ui of node i

is formally obtained as:

$$\lambda \mu _ { i } = \sum _ { j \to i } A _ { i j } \mu _ { j } \ \, \quad \ \ ( 1 ) \quad \text{for} \quad,$$

where λ is the maximum eigenvalue of A . Formally, if we denote with µ ∈ /Rbbb V and ˜ µ ∈ /Rbbb m the principal eigenvectors for A and A m ′ respectively and consider µ m ∈ /Rbbb m the restriction to the sample of µ , we want to find a subsample network Gm that minimizes the total distance d ( µ m , ˜ µ ) , under a vectorial distance metric d ( · , · ) . Notice that this means that our goal is to well approximate the centrality of sampled nodes (with ˜ µ ), rather than reconstructing the entire centrality µ . Given that the eigenvector centrality is often used as a measure for ranking nodes, we will consider two error measures relevant in this application, the Spearman and Kendall correlations, see section III.

## A. A theoretically motivated sampling method

Since we are interested in estimating eigenvectors, we use projection methods for spectral approximation. The idea is finding eigenvalues and eigenvectors by computing them on smaller subspaces and then projecting them backwards into the original (bigger) space. As anticipated above, in our case we are interested on better approximating only the eigenvector values relative to nodes in Gm , i.e. the values computed on the projected space.

In the context of projection methods, theoretical bounds for the goodness of approximation of the spectrum of A have been developed [11].

We derive a theorem that specializes these findings to our case as follows. Denote with Pm ∈ /Rbbb V V , the projector onto the subspace Km : = &lt; e i &gt; i ∈/Ical , where e i ∈ /Rbbb V are the vectors of the naive basis of /Rbbb V . In our case we are interested in approximating Pm µ , the eigenvector centrality of the the nodes in the subsample. The zero-padding submatrix of A where the rows and columns belonging to /Ical are selected is defined as Am : = Pm AP m = /parenleftbig4 A m ′ 0 0 0 /parenrightbig4 . Notice that Am ∈ /Rbbb V V , and A m ′ ∈ /Rbbb m m , ; although different in dimension, in practice they represent the same graph.

## Theorem 1. Let A Am , ≥ 0 and irreducible. Then:

$$\sin ( P _ { m } \mu \,, \tilde { \mu } ) \leq \frac { \gamma } { | | A _ { m } - \lambda \i 1 | | _ { 2, K _ { m } \cap \tilde { \mu } ^ { \perp } } } \tan ( \mu, K _ { m } ) \quad ( 2 ) \quad \text{(s: } \text{$\text{$\text{$row:}$} }$$

where:

$$\gamma \coloneqq | | P _ { m } A ( 1 - P _ { m } ) | | _ { 2 } \, \quad \quad ( 3 ) \quad \dot { \ m i 2 }$$

$$t a n ( \mu, K _ { m } ) = \sqrt { \frac { \sum _ { i \notin \mathcal { I } } \mu _ { i } ^ { 2 } } { \sum _ { i \in \mathcal { I } } \mu _ { i } ^ { 2 } } } = \sqrt { \frac { | | \mu | | _ { 2 } ^ { 2 } } { \sum _ { i \in \mathcal { I } } \mu _ { i } ^ { 2 } } - 1 } \quad ( 4 ) \quad \text{of } A. \\ a n d \, | | M | | _ { \dots } \, \text{ is the induced 2 non-m} of M \text{ restricted to the} } \quad \text{read}.$$

and || M || 2, V is the induced 2 norm of M restricted to the vectorial subspace V .

Equation (4) is derived by definition of cosine between a vector µ and a vectorial subspace Km . A formal proof can be found in the supplementary (proof 1.1).

The assumptions of having non-negative and irreducible matrices are needed to obtain uniquely defined eigenvalue centralities, as from Perron-Frobenius theorem [22]. In our case, both these hypothesis apply. Except for signed networks, adjacency matrices are always non-negative. In addition, we consider only the strong connected component of the networks considered in our experiments, these have irreducible adjacency matrices. This theoretical bound motivates the main idea behind our algorithm. We aim at reducing the bound on the sin distance between the real and approximated eigenvalue centrality. Specifically we focus at reducing the factor γ . Although the sine distance is not the main metric we consider in performance evaluation, we use it as a proxy for the minimization that we are interested in performing. In practice, we observe that this approach works, see section III. Nonetheless, as said above, the main optimization goal is not on the sine value itself, but on statistics that measure the relative rankings of nodes. In principle, we could reduce the bound in the right-hand side of Eq. (2) by minimizing tan ( µ , Km ) or maximizing || Am -λ /onebbb || 2, Km ∩ ˜ µ ⊥ . However, these two quantities are difficult to treat analytically and we therefore rely on acting on γ , which is instead accessible. We will partially address the minimization of tan ( µ , Km ) using theoretically grounded heuristics below.

From (3), we see that γ depends on the choice of the sample. In other words, if we add one node to the sample, the value of γ will change. This motivates our proposed sampling algorithm: we consider the online addition to the current sample of one node at a time by selecting the best node from the set of non-sampled nodes. The best node is defined as the one whose inclusion in the sample reduces the factor γ the most.

We now explicitly quantify the magnitude of γ as a function of adding one candidate node at a time. Splitting the ad-′

equivalent to the largest singular value s ( B ) of B , i.e. the maximum eigenvalue of the covariance matrix B T B [22]. This implies that minimizing γ is equivalent to minimizing s ( B ) . B contains directed edges from outside to inside the sample. Hence, reducing γ is related to reducing the number of incoming edges obtained from the sampling. C contains outgoing edges, while D contains edges between nodes outside the sample. While factors not related to edges incoming towards the sample from outside do not appear explicitly in the bound, these are included in the tangent factor 4.

$$\text{nodes} \quad \text{sources} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad\text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s}" \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \tilde { s } \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \text{s} \quad \end{cases}$$

It can be shown, see supplementary section 1.2, that minimizing the entries of B T B also reduces s ( B ) , hence we now quantify B T B . Consider a partial sample of size k , and a new candidate node to be chosen j ∈ /Vcal \ /Ical . Without loss of generality, we can think of this node as the k + 1 entry of A . This yields the decomposition B = /parenleftbig4 b 1 U b 2 b T 3 /parenrightbig4 , where b 1 ∈ /Rbbb k -1 are the edges pointing from j to the nodes already in the subsample, b 2 ∈ /Rbbb is the entry corresponding to j , b 3 ∈ /Rbbb n -k + 1 are edges from nodes outside the sample to-

wards j , U ∈ /Rbbb k -1, n -k + 1 are the edges from nodes outside the sample towards nodes in it, j excluded. We assume b 2 = 0 , i.e. no self-loops, but the method can be easily generalized to include them. We then greedily pick the node in the graph that, if chosen, minimizes the entries of

$$B ^ { T } B = \begin{pmatrix} b _ { 1 } ^ { T } b _ { 1 } & b _ { 1 } ^ { T } U \\ U ^ { T } b _ { 1 } & U ^ { T } U + b _ { 3 } b _ { 3 } ^ { T } \end{pmatrix} \ \. \quad \ ( 5 ) \quad \text{and} \quad \text{ing} \ ;$$

This corresponds to maximizing the following quantity:

$$| | b _ { 1 } | | _ { 2 } ^ { 2 } + | | b _ { 1 } ^ { T } U | | _ { 2 } ^ { 2 } - | | b _ { 3 } | | _ { 2 } ^ { 2 } \. \quad \quad ( 6 ) \quad \text{on} \ \text{$\quad$f$}$$

Notice that we consider all the terms in 5 except U T U . This choice is due to the cost of computing || U T U || , since U T U has size V . Give the relatively big size of the matrix, we assume that removing one node will not highly impact its norm, making this part of the optimization negligible compared to the other terms.

We now make few observations. Firstly, from the definition of eigenvalue centrality (1) we aim at minimizing the information missed by not sampling informative incoming neighbors of nodes in the sample. We define the border of a sampled subgraph as the set of all the incoming neighbors not included in the sample itself. In our online algorithm, we choose to select candidate nodes from those in the border. This choice is arbitrary, but motivated by the following observations. While selecting candidate nodes as any arbitrary node in the network might not be doable in all sampling scenarios [23], we assume that accessing neighbors for each node already in the sample is instead possible. This is the case for many real-world applications, e.g. social networks like Twitter and Facebook where follower/followee information is accessible. Second, the effectiveness of this choice is motivated by observing the performance's improvements moving from random nodes sampling (which does not take connecting edges into account) to random walk techniques (which is inherently sampling from neighbors). For comparative results see section III. Finally, notice that external nodes not in the border have b 1 = 0 , thus yielding a negative values for the theoretical criterion 6. This implies that, in practice, they are unlikely the one chosen when maximizing 6.

## B. An extension of the sampling algorithm

Empirical experiments have shown that addressing the minimization of the bound 2 by considering only γ is effective for undirected networks. This doesn't necessarily hold for directed ones. Our guess is that in these cases the term tan ( µ , Km ) 4 also plays an important role in the tightness of the bound in 2. In order to reduce this quantity, we may then select as candidates, nodes with high centrality. This choice stems from eq 4: the bigger the sum of sampled nodes' centralities, the smaller the tangent term. Since µ is not available, we approximate it. A possible choice is based on the results of [24], suggesting µ j ≈ d G in ( j ) , the (weighted) in-degree of node j . We instead consider µ j ≈ d Gk in ( j ) , the (weighted) indegree calculated considering only the incoming edges from nodes that are in the sample.

We incorporate this insight into our original criterion (6) by considering the following convex combination:

$$\mathfrak { h } ( 1 - \alpha ) \left ( | | b _ { 1 } | | _ { 2 } ^ { 2 } + | | b _ { 1 } ^ { T } U | | _ { 2 } ^ { 2 } - | | b _ { 3 } | | _ { 2 } ^ { 2 } \right ) + \alpha \, d _ { i n } ^ { G _ { k } } ( j ) \ \, \ \ ( 7 ) \\ \text{ where } \alpha \text{ is an hyperparameter that can be tuned empirically}$$

where α is an hyperparameter that can be tuned empirically and it allows to interpolate between the two approaches. Using a value of α closer to one means finding nodes with higher degree, therefore aiming to reduce the tangent term. On the other side, setting α = 0 means using the vanilla algorithm, only aiming at decreasing the γ term.

Empirically, we observe that the size of the border increases like O m ( ) , with a possibly non-negligible constant factor. Thus, at sampling time, we turn to a Monte Carlo approximation: we calculate the value 7 only for a randomly chosen fraction of size p of the new neighbors inside the border. We store the results in a leaderboard for future sampling. Finally, we impose a maximum size s for the leaderboard, this is for two reasons. First, keeping an ordered leaderboard that allows fast extraction of the maximum element, as well as fast insertion of new ones, is computationally and memory expensive. Secondly, the scores are computed for newly encountered neighbors and stored in the leaderboard for future extraction. The scores given by 7 depend on the current sample. The scores of border nodes that have been computed sometime during past iterations are stored and kept valid for selecting candidate nodes in successive rounds; it is thus advisable to discard old scores by exploiting a small leaderboard size, forcing old nodes to be popped out without being selected when new, potentially higher scoring, nodes are found.

A pseudocode for this is shown in Algorithm 1.

## Algorithm 1 Theoretical criterion sampling (TCEC)

Input: s leaderboard max size, p randomization probability, m final sample size, 1 ≤ k &lt; m random walk initialization size Output: Gm final subgraph

- - store the border of Gk and compute 6 for a random fraction of its nodes
- - sample k nodes via random walk; call the obtained partial sample Gk
- - save these values in a leaderboard
- -while size ( Gk ) &lt; m :

pick v , highest value neighbor from leaderboard

Gk

+

1

Gk

←

+

v

for a fraction p of the neighbors of v , compute 6 and store in leaderboard

## C. Complexity analysis

We compare the complexity of our sampling algorithm with the one of random walk methods. At every new node choice, random walk stochastically picks a node among the out-connections of the currently visited one. Random walk has a step cost of O ( 1 ) .

Call dout , din the average (weighted) out and in degree of the nodes. Our algorithm performs the following steps: i) extract the best node from the leaderboard; ii) find the new din neighbors; iii) compute the theoretical criterion for a fraction p · din

of them; iv) store this in the leaderboard. Since a single computation of the theoretical criterion has cost O din ( + dout ) , each step costs O p din ( · · ( din + dout )) . Implementing the leaderboard with a heap with size limit s allows for a cost of O ( 1 ) for the extraction of the maximum and a cost of O din ( · log ( )) s for the insertion of all new scores. As by design we decide to keep s small, these leaderboard costs are negligible. Finally, notice that with the increase in sample size, the effective number of new neighbors is less than p · din , since many of the neighbors have already been sampled. The final cost is then empirically lower than what expected, and in any case constant with the sample size k .

## III. Empirical studies

We apply the sampling algorithms to real and synthetic datasets. We describe these datasets below, and we provide a summary of the datasets in Table I.

Synthetic data. We consider Erd˝ os-Rényi networks [25], a popular topology when the goal is testing performance on random networks.

Real data. We use four real-world open datasets. These datasets include both directed and undirected and both binary and weighted topologies.

Condensed Matter collaboration network 2 [27]: this dataset contains all the authors in the e-print arXiv Condensed Matter Physics section; an undirected edge is present when two authors collaborated in a publication.

Internet topology 1 [26]: this dataset consists of a a single daily snapshot of a part of the AS (autonomous systems) Internet layer. An undirected binary edge is present between two nodes if a contact between the two has happened during the day.

Epinions 3 [28]: this dataset is a type of who-trusts-whom directed network extracted from the review site Epionions.com, a consumer reviews website where users rate items and mark users as trustworthy.

Stanford Web Network 4 [29]: this dataset represents the pages in the Stanford.edu domain. A directed binary edge represents a hyperlink from one page to another.

Experimental setup. For each dataset, we sample subnetworks of different sizes using several sampling algorithms. In addition to ours, we consider mainstream sampling algorithms as random walk [30], Metropolis-Hastings random walk [31], a degree-weighted random walk and uniform sampling on nodes (UniSam). In the results' tables we report only the best performing results among the various random walk sampling techniques listed above and denote it as RWbest ; we denote our method as TCEC. We considered also snowball sampling [32] and breadth first search (BFS) for comparison. As these provided poorer results then the other methods listed above, we

1 http://irl.cs.ucla.edu/topology/

2 https://snap.stanford.edu/data/ca-CondMat.html

3 https://snap.stanford.edu/data/soc-Epinions1.html

4 https://snap.stanford.edu/data/web-Stanford.html

omit them here for the sake of space. Finally, we also include some modern sampling, techniques, such as forest fire [3] and expansions sampling [33].

Implementation details. Empirically, we find that for undirected graphs the vanilla version (6) of the algorithm, corresponding to α = 0 , leads to good performances; for directed ones we find good performance with α = 0.5 . All experiments presented here are run using these values. Similarly, the leaderboard maximum size is always set to s = 100 , which we find generally good for most applications. We also initialize every sampling round of our algorithm with a random walk exploration of size 1 5 / of the required finaTCECl sample size. This is needed to allow for a degree of reliability in computing the theoretical criterion 7 that requires a minimum sample size to be really discriminative for discerning important nodes. Finally, in order to obtain faster convergence, we set the randomization level p at different values ranging from 0.2 to 0.8 depending on the dataset. In practice, p can be chosen after the random walk initialization, when information about the average degree is available.

Performance evaluation. We compare the ground-truth global vector of eigenvalue centrality restricted to the entries relative to sampled nodes Pm µ , with the approximate centrality on the sampled network ˜ µ . As it is the case for many applications, we are interested in testing whether the relative ranking among nodes is preserved, rather than recovering the absolute values of the centrality scores. For this reason we compare the two vectors using the Spearman and Kendall rank correlation coefficients [34]. As we empirically see that these two metrics give similar results, we report results only in terms of Kendall and leave the ones with Spearman in the supplementary (section 1.3). One should notice that it is usually relatively easy for any algorithm to obtain high scores on a global scale. This is not surprising, as it is generally easy to tell very different nodes apart, obtaining, for example, a good comparison score in the majority of the possible /parenleftbig1 m 2 /parenrightbig1 pairs using Kendall. Things are more difficult on a finer scale, i.e. when selecting nodes nearby in the ranking; in this case nodes are often mistaken in terms of approximate ranking. For this reason we decided to evaluate the two statistics also on a moving window of 10% of the sampled nodes, ordered by global score. We report these results in the supplementary as well (section 1.3).

## A. Results on synthetic data

As we can see in tables II and S1, S2, S3, on an Erd˝ osRényi graph our method performs worse both on a global and local scale. This is explained by the fact that the fundamental assumption underlying our algorithm, i.e. that some nodes are more important than others for the goodness of estimation, is fundamentally wrong in this case. In fact, in an Erd˝s-Rényi o graph all nodes share the same statistical properties; this is shown also by the equivalent performances of random walk algorithms and uniform random sampling on nodes in this case. However, it is very unlikely to observe this statistical symmetry in real networks. This is also the only case in which

Table I. Datasets Description

| name                  |   #nodes | #edges     |   Average Degree | Directed   | Edge Type      |
|-----------------------|----------|------------|------------------|------------|----------------|
| Erd˝s-Rényi o         |    30000 | 4.5 · 10 6 |            300   | No         | Binary         |
| Internet topology     |    49179 | 49179      |              3.8 | No         | Binary         |
| Collaboration Network |    21363 | 91342      |              4.3 | No         | Discrete Count |
| Epinions Network      |    32223 | 443506     |             27.5 | Yes        | Binary         |
| Stanford Web Network  |   150532 | 1576314    |             20.9 | Yes        | Binary         |

Table II. Global Kendall on Erd˝s-Rényi Synthetic Graph o

Table III. Global Kendall on Internet topology Network

| S ratio   | UniSam            | RW Best                | TCEC          |
|-----------|-------------------|------------------------|---------------|
| 5%        | 0 . 145 ± 0 . 015 | 0.142 ± 0.016 (DW)     | 0.139 ± 0.020 |
| 10%       | 0 . 208 ± 0 . 013 | 0.205 ± 0.014 (DW)     | 0.165 ± 0.011 |
| 15%       | 0.251 ± 0.005     | 0 . 253 ± 0 . 009 (RW) | 0.190 ± 0.009 |
| 20%       | 0 . 294 ± 0 . 007 | 0.293 ± 0.009 (DW)     | 0.222 ± 0.007 |
| 25%       | 0.331 ± 0.004     | 0 . 334 ± 0 . 007 (RW) | 0.246 ± 0.007 |
| 30%       | 0.369 ± 0.004     | 0 . 373 ± 0 . 006 (DW) | 0.276 ± 0.007 |
| 40%       | 0.435 ± 0.004     | 0 . 436 ± 0 . 003 (DW) | 0.345 ± 0.004 |

uniform random sampling on nodes performs as well as other methods (see the next sections).

## B. Results on real data

In this case, we can see from Tables III IV V VI, that our algorithm outperforms all other sampling strategies both on a global and local scale in the majority of the dataset. For the Epinions network, while the relative trend among other sampling methods is unchanged, ours yields lower performances on smaller sample sizes. Only from sample sizes around 20% our method brings statistical improvements with respect to random walk methods. Our guess is that small sample sizes do not allow a good statistical assessment for the computation of the theoretical criterion 7, also because in the directed case we expect a heavy tailed distribution of centrality score, which means that detecting important nodes becomes harder due to their small number. Results are presented in table V (and in supplementary tables S10, S11, S12). The Stanford Web Network is affected by a high level of noise in the results VI (and in supplementary tables S13, S14, S15), signaled by poor performance and negative correlations. This applies to all sampling methods for small sample sizes. Only around 25% sampling size, results start to get more stable and performance improves. On a global scale, statistics grow significantly higher from that point. On the local evaluation window, values are as low as a few percentile points; however, our method obtains statistically significant improvements of around 100% with respect to the best of all random walk sampling algorithms.

## IV. Conclusions

While handling network data a relevant problem arises when estimating quantities if one can only access a limited amount of information. If one has the ability to decide what

| S ratio   | UniSam        | RW Best             | TCEC               |
|-----------|---------------|---------------------|--------------------|
| 5%        | 0.357 ± 0.681 | 0.977 ± 0.0022 (DW) | 0 . 982 ± 0 . 0014 |
| 10%       | 0.303 ± 0.044 | 0.984 ± 0.0006 (MH) | 0 . 984 ± 0 . 0014 |
| 15%       | 0.315 ± 0.002 | 0.986 ± 0.0003 (MH) | 0 . 987 ± 0 . 0011 |
| 20%       | 0.328 ± 0.003 | 0.987 ± 0.0004 (MH) | 0 . 990 ± 0 . 0002 |
| 25%       | 0.355 ± 0.003 | 0.989 ± 0.0003 (DW) | 0 . 993 ± 0 . 0003 |
| 30%       | 0.397 ± 0.001 | 0.990 ± 0.0002 (RW) | 0 . 995 ± 0 . 0001 |
| 40%       | 0.441 ± 0.000 | 0.992 ± 0.0001 (MH) | 0 . 997 ± 0 . 0001 |

| Table IV. Global Kendall on Collaboration Network   | Table IV. Global Kendall on Collaboration Network   | Table IV. Global Kendall on Collaboration Network   | Table IV. Global Kendall on Collaboration Network   |
|-----------------------------------------------------|-----------------------------------------------------|-----------------------------------------------------|-----------------------------------------------------|
| S ratio                                             | UniSam                                              | RW Best                                             | TCEC                                                |
| 5%                                                  | 0.245 ± 0.046                                       | 0.729 ± 0.064 (DW) 0                                | . 935 ± 0 . 003                                     |
| 10%                                                 | 0.339 ± 0.029                                       | 0.819 ± 0.032 (RW) 0 .                              | 958 ± 0 . 001                                       |
| 15%                                                 | 0.380 ± 0.022                                       | 0.874 ± 0.007 (RW) 0 .                              | 966 ± 0 . 001                                       |
| 20%                                                 | 0.403 ± 0.015                                       | 0.900 ± 0.006 (DW) 0 .                              | 970 ± 0 . 002                                       |
| 25%                                                 | 0.451 ± 0.015                                       | 0.921 ± 0.005 (DW) 0 .                              | 974 ± 0 . 001                                       |
| 30%                                                 | 0.486 ± 0.014                                       | 0.938 ± 0.005 (MH) 0 .                              | 977 ± 0 . 001                                       |
| 40%                                                 | 0.545 ± 0.019                                       | 0.960 ± 0.003 (MH) 0                                | . 982 ± 0 . 001                                     |

fraction of the dataset to collect or observe, it is then important to choose the sample wisely in order to recover good estimates of the relevant quantities. We presented here a novel sampling algorithm for estimating a popular centrality measure on networks, eigenvector centrality. The algorithm is theoretically motivated using methods for eigenvalue approximation. It is computationally scalable and valid for both directed and undirected weighted non-negative adjacency matrices.

The model relies on lowering the bound on the sin distance between the ground truth and in-sample estimated eigenvector centrality. This allows to well preserve the relative ranking of nodes inside the sample ordered by their estimated eigenvector centralities. It allows to interpolate between quantities related to the network topology which contribute to dif-

|         | Table V. Global Kendall on Epinions   | Table V. Global Kendall on Epinions   | Network           |
|---------|---------------------------------------|---------------------------------------|-------------------|
| S ratio | UniSam                                | RW Best                               | TCEC              |
| 5%      | 0.358 ± 0.260                         | 0 . 634 ± 0.010 (RW)                  | 0.498 ± 0.009     |
| 10%     | 0.378 ± 0.012                         | 0 . 649 ± 0 . 005 (RW)                | 0.601 ± 0.005     |
| 15%     | 0.389 ± 0.010                         | 0 . 660 ± 0 . 003 (MH)                | 0.658 ± 0.004     |
| 20%     | 0.409 ± 0.014                         | 0.665 ± 0.004 (MH)                    | 0 . 681 ± 0 . 004 |
| 25%     | 0.429 ± 0.017                         | 0.665 ± 0.002 (RW)                    | 0 . 686 ± 0 . 002 |
| 30%     | 0.425 ± 0.007                         | 0.664 ± 0.003 (DW)                    | 0 . 688 ± 0 . 002 |
| 40%     | 0.461 ± 0.011                         | 0.657 ± 0.002 (DW)                    | 0 . 676 ± 0 . 002 |

Table VI. Global Kendall on Stanford Web Network

| S ratio   | UniSam            | RW Best                | TCEC              |
|-----------|-------------------|------------------------|-------------------|
| 5%        | - 0.350 ± 0.825   | - 0.110 ± 0.630 (MH)   | 0 . 157 ± 0 . 108 |
| 10%       | 0.064 ± 0.008     | 0 . 110 ± 0 . 010 (DW) | - 0.300 ± 0.854   |
| 15%       | 0 . 059 ± 0 . 009 | - 0.060 ± 0.648 (MH)   | - 0.101 ± 0.634   |
| 20%       | 0.055 ± 0.008     | 0 . 164 ± 0 . 023 (MH) | 0.092 ± 0.088     |
| 25%       | - 0.146 ± 0.618   | 0.164 ± 0.007 (RW)     | 0 . 170 ± 0 . 092 |
| 30%       | - 0.144 ± 0.619   | 0.161 ± 0.009 (RW)     | 0 . 258 ± 0 . 070 |
| 40%       | 0.078 ± 0.013     | 0.143 ± 0.007 (RW)     | 0 . 253 ± 0 . 018 |

ferent extent to the quality of the approximation via hyperparameter tuning. We empirically find that on real networks the vanilla model performs well on undirected networks compared to other standard sampling algorithms; for directed ones

- [1] Munmun De Choudhury, Yu-Ru Lin, Hari Sundaram, Kasim Selcuk Candan, Lexing Xie, and Aisling Kelliher. How does the data sampling strategy impact the discovery of information diffusion in social media? In Fourth International AAAI Conference on Weblogs and Social Media , 2010.

[2] Eldar Sadikov, Montserrat Medina, Jure Leskovec, and Hector Garcia-Molina. Correcting for missing data in information cascades. In Proceedings of the fourth ACM international conference on Web search and data mining , pages 55-64. ACM, 2011.

[3] Jure Leskovec and Christos Faloutsos. Sampling from large graphs. In Proceedings of the 12th ACM SIGKDD international conference on Knowledge discovery and data mining , pages 631-636. ACM, 2006.

[4] Micah Adler and Michael Mitzenmacher. Towards compressing web graphs. In Proceedings DCC 2001. Data Compression Conference , pages 203-212. IEEE, 2001.

[5] Ove Frank. Network sampling and model fitting. Models and methods in social network analysis , pages 31-56, 2005.

[6] Jing-Dong J Han, Denis Dupuy, Nicolas Bertin, Michael E Cusick, and Marc Vidal. Effect of sampling on topology predictions of protein-protein interaction networks. Nature biotechnology , 23(7):839, 2005.

[7] Sang Hoon Lee, Pan-Jun Kim, and Hawoong Jeong. Statistical properties of sampled networks. Physical Review E , 73(1):016102, 2006.

[8] Gueorgi Kossinets. Effects of missing data in social networks. Social networks , 28(3):247-268, 2006.

[9] Phillip Bonacich. Factoring and weighting approaches to status scores and clique identification. Journal of mathematical sociology , 2(1):113-120, 1972.

[10] Elizabeth Costenbader and Thomas W Valente. The stability of centrality measures when networks are sampled. Social networks , 25(4):283-307, 2003.

[11] Y. Saad. Numerical Methods for Large Eigenvalues Problems . 2011.

[12] Neli Blagus, Lovro Šubelj, and Marko Bajec. Empirical comparison of network sampling: How to choose the most appropriate method? Physica A: Statistical Mechanics and its Applications , 477:136-148, 2017.

5 https://github.com/cdebacco/tcec

performance is improved when using the full model.

Our model assumes that the network topology is static in time. It would be interesting to investigate extensions of the algorithm to the dynamical case, when node or edge removal/addition are possible. The method presented in this paper can in principle be adapted for estimating other spectral centrality measures, for instance PageRank. We focused here in optimizing the estimate of a specific centrality measure. An open question is whether our algorithm, while optimized to preserve this global network property, is indirectly preserving also other network quantities. While we leave these extensions for future work, we provide an open source implementation of the algorithm . 5

- [13] Fred Morstatter, Jürgen Pfeffer, Huan Liu, and Kathleen M Carley. Is the sample good enough? comparing data from twitter's streaming api with twitter's firehose. In Seventh international AAAI conference on weblogs and social media , 2013.

[14] Daniel Stutzbach, Reza Rejaie, Nick Duffield, Subhabrata Sen, and Walter Willinger. On unbiased sampling for unstructured peer-to-peer networks. IEEE/ACMTransactions on Networking (TON) , 17(2):377-390, 2009.

[15] Christian Hübler, Hans-Peter Kriegel, Karsten Borgwardt, and Zoubin Ghahramani. Metropolis algorithms for representative subgraph sampling. In 2008 Eighth IEEE International Conference on Data Mining , pages 283-292. IEEE, 2008.

[16] Michael PH Stumpf and Carsten Wiuf. Sampling properties of random graphs: the degree distribution. Physical Review E , 72(3):036118, 2005.

[17] Apratim Ganguly and Eric D Kolaczyk. Estimation of vertex degrees in a sampled network. In 2017 51st Asilomar Conference on Signals, Systems, and Computers , pages 967-974. IEEE, 2018.

[18] Nelson Antunes, Shankar Bhamidi, Tianjian Guo, Vladas Pipiras, and Bang Wang. Sampling-based estimation of in-degree distribution with applications to directed complex networks. arXiv preprint arXiv:1810.01300 , 2018.

[19] Sergey Brin and Lawrence Page. The anatomy of a large-scale hypertextual web search engine. Computer networks and ISDN systems , 30(1-7):107-117, 1998.

[20] Yuta Sakakura, Yuto Yamaguchi, Toshiyuki Amagasa, and Hiroyuki Kitagawa. An improved method for efficient pagerank estimation. In International Conference on Database and Expert Systems Applications , pages 208-222. Springer, 2014.

[21] Yen-Yu Chen, Qingqing Gan, and Torsten Suel. Local methods for estimating pagerank values. In Proceedings of the thirteenth ACM international conference on Information and knowledge management , pages 381-389. ACM, 2004.

[22] Gene H Golub and Charles F Van Loan. Matrix computations , volume 3. JHU press, 2012.

[23] Minas Gjoka, Maciej Kurant, Carter T Butts, and Athina Markopoulou. Walking in facebook: A case study of unbiased sampling of osns. In 2010 Proceedings IEEE Infocom , pages 1-9. Ieee, 2010.

[24] Miguel Romance. Local estimates for eigenvector-like centralities of complex networks. Journal of Computational and Applied Mathematics , 235(7):1868-1874, 2011.

- [25] Paul Erdös and Alfréd Rényi. On random graphs, i. Publicationes Mathematicae (Debrecen) , 6:290-297, 1959.
- [26] Ricardo Oliveira, Walter Willinger, Beichuan Zhang, et al. Quantifying the completeness of the observed internet as-level structure. work , 11(15):13-17, 2008.
- [27] Jure Leskovec, Jon Kleinberg, and Christos Faloutsos. Graph evolution: Densification and shrinking diameters. ACM Transactions on Knowledge Discovery from Data (TKDD) , 1(1):2, 2007.
- [28] Matthew Richardson, Rakesh Agrawal, and Pedro Domingos. Trust management for the semantic web. In International semantic Web conference , pages 351-368. Springer, 2003.
- [29] Jure Leskovec, Kevin J Lang, Anirban Dasgupta, and Michael W Mahoney. Community structure in large networks: Natural cluster sizes and the absence of large well-defined clusters. Internet Mathematics , 6(1):29-123, 2009.
- [30] László Lovász et al. Random walks on graphs: A survey. Combinatorics, Paul erdos is eighty , 2(1):1-46, 1993.
- [31] Nicholas Metropolis, Arianna W Rosenbluth, Marshall N Rosenbluth, Augusta H Teller, and Edward Teller. Equation of state calculations by fast computing machines. The journal of chemical physics , 21(6):1087-1092, 1953.
- [32] Leo A Goodman. Snowball sampling. The annals of mathematical statistics , pages 148-170, 1961.
- [33] Arun S Maiya and Tanya Y Berger-Wolf. Sampling community structure. In Proceedings of the 19th international conference on World wide web , pages 701-710. ACM, 2010.
- [34] Gregory W Corder and Dale I Foreman. Nonparametric statistics: A step-by-step approach . John Wiley &amp; Sons, 2014.

## S1. Additional Theorems and Proofs

## A. Proof of theorem 1

We follow the same notation of the proof in [11, p.103/104]. What we derive here is a specialization of their proof to our

Proof. case. Call φ the angle between µ and Km and ω the angle between Pm µ and ˜ µ . Therefore we can decompose

$$\mu = \nu \cos \phi + w \sin \phi$$

where w v , are the normalized projections of µ respectively on Km and K m ⊥ . Multiplying both sides by ( A -λ /onebbb ) , since A µ = λµ , we obtain

$$( A - \lambda 1 ) v \cos \phi + ( A - \lambda 1 ) w \sin \phi = 0$$

By multiplying for P m and taking the norm we obtain

$$| | P _ { m } ( A - \lambda 1 ) v | | \cos \phi = | | P _ { m } ( A - \lambda 1 ) w | | \sin \phi \leq \gamma \sin \phi$$

where the second passage derives from the fact that w =( /onebbb -Pm ) w therefore

$$| | P _ { m } ( A - \lambda 1 ) w | | & = | | P _ { m } ( A - \lambda 1 ) ( 1 - P _ { m } ) w | | \\ & = | | P _ { m } A ( 1 - P _ { m } ) w - \lambda P _ { m } ( 1 - P _ { m } ) w | | \\ & = | | P _ { m } A ( 1 - P _ { m } ) w | | \\ & \leq | | P _ { m } A ( 1 - P _ { m } ) | | = \gamma$$

Now, consider the vector z ∈ Km , z ⊥ ˜ µ such that

$$v = \tilde { \mu } \cos \omega + z \sin \omega$$

Then applying the left multiplication Pm ( A -λ /onebbb ) we obtain

$$P _ { m } ( A - \lambda 1 ) v = ( \tilde { \lambda } - \lambda ) \tilde { \mu } \cos \omega + P _ { m } ( A - \lambda 1 ) z \sin \omega$$

which implies

Then using (S1) we obtain

$$| | ( A _ { m } - \lambda 1 ) | | _ { 2, K _ { m } \cap \tilde { \mu } ^ { \perp } } \sin \omega \cos \phi \leq \gamma \sin \phi$$

which concludes the proof since sin ω = sin ( Pm µ µ , ˜ ) and sin φ = sin ( µ , Km )

## B. Additional Theorems

Lemma 1. Let 0 ≤ B ∈ /Rbbb s t , (where 0 ≤ B is intended elementwise, i.e. 0 ≤ Bi j , ∀ i , j ). Then the largest singular value of B s B , ( ) , satisifies the following

$$s ( B ) & = \max _ { \| x \| = 1 } x ^ { T } B ^ { T } B x = \max _ { \| x \| = 1 } \| B x \| _ { 2 } ^ { 2 } \\ & = \max _ { \| x \| = 1, x \geq 0 } \| B x \| _ { 2 } ^ { 2 } = \max _ { \| x \| = 1, x \geq 0 } x ^ { T } B ^ { T } B x$$

where the first equality derives from Courant-Fischer theorem [11, p24], as s ( B ) is the biggest eigenvalue of B T B , and the third one from the fact that B ≥ 0

$$| | P _ { m } ( A - \lambda 1 ) v | | & \geq | | P _ { m } ( A - \lambda 1 ) z | | \sin \omega \\ & = | | ( P _ { m } A P _ { m } - \lambda 1 ) z | | \sin \omega \\ & = | | ( A _ { m } - \lambda 1 ) | | _ { 2, K _ { m } \cap \tilde { \mu } ^ { \perp } } \sin \omega$$

Theorem 2. Let 0 ≤ B C , ∈ /Rbbb s t , , B T B ≤ C T C . Then s ( B ) ≤ s ( C ) . (This theorem holds the same if we only suppose 0 ≤ B ≤ C , but we use this version for developing the algorithm)

Proof. Using the previous lemma

$$s ( B ) = \max _ { \| x \| = 1, x \geq 0 } x ^ { T } B ^ { T } B x \leq \max _ { \| x \| = 1, x \geq 0 } x ^ { T } C ^ { T } C x = s ( C )$$

(to prove in the case B ≤ C , just use s ( B ) = max || x || = 1, x ≥ 0 || Bx || 2 2 ≤ max || x || = 1, x ≥ 0 || Cx || 2 2 = ( s C ) )

## S2. Supplemental Tables

In this section we provide the results that are not included in the original paper. These are made of the Spearman score both on the global score vectors and on a moving window of size 10% of the entire sample. The moving window scores for the Kendall Tau statistic are reported here as well. Notice that in some of the columns corresponding to the moving window score for the uniform sampling, some NA entries are present. These are placeholders to indicate where it hasn't been possible to perform the computations due numerical reasons.

| Table S1. Global Spearman on Erd˝s-Rényi o Synthetic Graph   | Table S1. Global Spearman on Erd˝s-Rényi o Synthetic Graph   | Table S1. Global Spearman on Erd˝s-Rényi o Synthetic Graph   | Table S1. Global Spearman on Erd˝s-Rényi o Synthetic Graph   |
|--------------------------------------------------------------|--------------------------------------------------------------|--------------------------------------------------------------|--------------------------------------------------------------|
| S ratio                                                      | UniSam                                                       | RW Best                                                      | TCEC                                                         |
| 5%                                                           | 0.103 ± 0.023                                                | 0 . 211 ± 0 . 024 (DW)                                       | 0.206 ± 0.029                                                |
| 10%                                                          | 0.169 ± 0.018                                                | 0 . 304 ± 0 . 020 (DW)                                       | 0.245 ± 0.017                                                |
| 15%                                                          | 0.251 ± 0.007                                                | 0 . 371 ± 0 . 012 (RW)                                       | 0.282 ± 0.014                                                |
| 20%                                                          | 0.265 ± 0.009                                                | 0.293 ± 0.013 (DW)                                           | 0 . 327 ± 0 . 010                                            |
| 25%                                                          | 0.331 ± 0.006                                                | 0 . 484 ± 0 . 010 (RW)                                       | 0.361 ± 0.011                                                |
| 30%                                                          | 0.369 ± 0.005                                                | 0 . 529 ± 0 . 007 (RW)                                       | 0.404 ± 0.010                                                |
| 40%                                                          | 0.435 ± 0.005                                                | 0 . 612 ± 0 . 006 (MH)                                       | 0.500 ± 0.005                                                |

Table S2. Moving Window Spearman on Erd˝ os-Rényi Synthetic Graph

| S ratio   | UniSam            | RW Best                | TCEC              |
|-----------|-------------------|------------------------|-------------------|
| 5%        | 0 . 026 ± 0 . 003 | 0.025 ± 0.006 (MH)     | 0 . 026 ± 0 . 004 |
| 10%       | 0 . 036 ± 0 . 003 | 0 . 036 ± 0 . 002 (MH) | 0.029 ± 0.003     |
| 15%       | 0.044 ± 0.002     | 0 . 045 ± 0 . 003 (RW) | 0.035 ± 0.003     |
| 20%       | 0 . 053 ± 0 . 002 | 0 . 053 ± 0 . 003 (DW) | 0.039 ± 0.002     |
| 25%       | 0.061 ± 0.002     | 0 . 062 ± 0 . 002 (RW) | 0.045 ± 0.001     |
| 30%       | 0.069 ± 0.002     | 0 . 070 ± 0 . 002 (DW) | 0.050 ± 0.001     |
| 40%       | 0 . 087 ± 0 . 002 | 0 . 087 ± 0 . 001 (DW) | 0.066 ± 0.001     |

Table S3. Moving Window Kendall on Erd˝ os-Rényi Synthetic Graph

| S ratio   | UniSam            | RW Best                | TCEC          |
|-----------|-------------------|------------------------|---------------|
| 5%        | 0 . 018 ± 0 . 002 | 0.017 ± 0.004 (MH)     | 0.017 ± 0.002 |
| 10%       | 0 . 024 ± 0 . 002 | 0 . 024 ± 0 . 001 (MH) | 0.019 ± 0.002 |
| 15%       | 0.029 ± 0.001     | 0 . 030 ± 0 . 002 (RW) | 0.023 ± 0.002 |
| 20%       | 0 . 036 ± 0 . 001 | 0.035 ± 0.002 (DW)     | 0.026 ± 0.001 |
| 25%       | 0 . 041 ± 0 . 001 | 0 . 041 ± 0 . 001 (MH) | 0.030 ± 0.001 |
| 30%       | 0.046 ± 0.001     | 0 . 047 ± 0 . 001 (DW) | 0.033 ± 0.001 |
| 40%       | 0 . 058 ± 0 . 001 | 0 . 058 ± 0 . 000 (DW) | 0.044 ± 0.001 |

| Table S4. Global Spearman on Internet Topology Network   | Table S4. Global Spearman on Internet Topology Network   | Table S4. Global Spearman on Internet Topology Network   | Table S4. Global Spearman on Internet Topology Network   |
|----------------------------------------------------------|----------------------------------------------------------|----------------------------------------------------------|----------------------------------------------------------|
| S ratio                                                  | UniSam                                                   | RW Best                                                  | TCEC                                                     |
| 5%                                                       | 0.095 ± 0.702                                            | 0 . 999 ± 0 . 000 (DW)                                   | 0 . 999 ± 0 . 000                                        |
| 10%                                                      | 0.388 ± 0.057                                            | 0 . 999 ± 0 . 000 (MH)                                   | 0 . 999 ± 0 . 000                                        |
| 15%                                                      | 0.402 ± 0.072                                            | 0 . 999 ± 0 . 000 (MH)                                   | 0 . 999 ± 0 . 000                                        |
| 20%                                                      | 0.422 ± 0.109                                            | 0 . 999 ± 0 . 000 (DW)                                   | 0 . 999 ± 0 . 000                                        |
| 25%                                                      | 0.454 ± 0.057                                            | 0 . 999 ± 0 . 000 (DW)                                   | 0 . 999 ± 0 . 000                                        |
| 30%                                                      | 0.502 ± 0.077                                            | 0 . 999 ± 0 . 000 (DW)                                   | 0 . 999 ± 0 . 000                                        |
| 40%                                                      | 0.552 ± 0.094                                            | 0 . 999 ± 0 . 000 (MH)                                   | 0 . 999 ± 0 . 000                                        |

Table S5. Moving Window Spearman on Internet Topology Network

| S ratio   | UniSam   | RW Best                | TCEC            |
|-----------|----------|------------------------|-----------------|
| 5%        | NA       | 0.928 ± 0.011 (RW) 0 . | 950 ± 0 . 006   |
| 10%       | NA       | 0.955 ± 0.003 (MH) 0 . | 958 ± 0 . 006   |
| 15%       | NA       | 0.965 ± 0.002 (DW) 0 . | 969 ± 0 . 004   |
| 20%       | NA       | 0.971 ± 0.002 (DW) 0 . | 980 ± 0 . 001   |
| 25%       | NA       | 0.976 ± 0.001 (DW) 0 . | 989 ± 0 . 001   |
| 30%       | NA       | 0.980 ± 0.001 (DW) 0 . | 993 ± 0 . 001   |
| 40%       | NA       | 0.987 ± 0.000 (RW) 0   | . 997 ± 0 . 000 |

Table S6. Moving Window Kendall on Internet topology Network

| S ratio   | UniSam   | RW Best                | TCEC          |
|-----------|----------|------------------------|---------------|
| 5%        | NA       | 0.810 ± 0.014 (RW) 0 . | 839 ± 0 . 012 |
| 10%       | NA       | 0.852 ± 0.005 (MH) 0 . | 854 ± 0 . 011 |
| 15%       | NA       | 0.869 ± 0.004 (DW) 0 . | 880 ± 0 . 010 |
| 20%       | NA       | 0.880 ± 0.003 (DW) 0 . | 907 ± 0 . 002 |
| 25%       | NA       | 0.892 ± 0.002 (DW) 0 . | 934 ± 0 . 004 |
| 30%       | NA       | 0.903 ± 0.002 (RW) 0 . | 952 ± 0 . 002 |
| 40%       | NA       | 0.923 ± 0.001 (MH) 0 . | 974 ± 0 . 001 |

Table S7. Global Spearman on Collaboration Network

| S ratio   | UniSam   | RW Best                | TCEC              |
|-----------|----------|------------------------|-------------------|
| 5%        | NA       | 0.893 ± 0.050 (DW) 0 . | 994 ± 0 . 001     |
| 10%       | NA       | 0.950 ± 0.016 (RW) 0 . | 997 ± 0 . 000     |
| 15%       | NA       | 0.974 ± 0.003 (RW) 0   | . 998 ± 0 . 000   |
| 20%       | NA       | 0.982 ± 0.002 (DW) 0   | . 998 ± 0 . 000   |
| 25%       | NA       | 0.988 ± 0.002 (MH)     | 0 . 999 ± 0 . 000 |
| 30%       | NA       | 0.992 ± 0.001 (MH)     | 0 . 999 ± 0 . 000 |
| 40%       | NA       | 0.996 ± 0.001 (MH)     | 0 . 999 ± 0 . 000 |

Table S8. Moving Window Spearman on Collaboration Network

| S ratio   | UniSam   | RW Best                    | TCEC              |
|-----------|----------|----------------------------|-------------------|
| 5%        | NA       | 0.250 ± 0.063 (DW) 0 . 692 | ± 0 . 019         |
| 10%       | NA       | 0.370 ± 0.055 (RW) 0 .     | 826 ± 0 . 007     |
| 15%       | NA       | 0.503 ± 0.024 (RW) 0 .     | 875 ± 0 . 008     |
| 20%       | NA       | 0.593 ± 0.018 (RW) 0 .     | 899 ± 0 . 007     |
| 25%       | NA       | 0.677 ± 0.025 (DW) 0 .     | 920 ± 0 . 006     |
| 30%       | NA       | 0.752 ± 0.015 (DW) 0 .     | 932 ± 0 . 004     |
| 40%       | NA       | 0.855 ± 0.013 (MH)         | 0 . 950 ± 0 . 002 |

Table S9. Moving Window Kendall on Collaboration Network

| S ratio   | UniSam   | RW Best                | TCEC          |
|-----------|----------|------------------------|---------------|
| 5%        | NA       | 0.729 ± 0.064 (DW) 0 . | 935 ± 0 . 003 |
| 10%       | NA       | 0.819 ± 0.032 (RW) 0 . | 958 ± 0 . 001 |
| 15%       | NA       | 0.874 ± 0.007 (RW) 0 . | 966 ± 0 . 001 |
| 20%       | NA       | 0.900 ± 0.006 (DW) 0 . | 970 ± 0 . 002 |
| 25%       | NA       | 0.921 ± 0.005 (DW) 0 . | 974 ± 0 . 001 |
| 30%       | NA       | 0.938 ± 0.005 (MH) 0 . | 977 ± 0 . 001 |
| 40%       | NA       | 0.960 ± 0.003 (MH) 0 . | 982 ± 0 . 001 |

Table S10. Global Spearman on Epinions Network

| S ratio   | UniSam        | RW Best                | TCEC              |
|-----------|---------------|------------------------|-------------------|
| 5%        | 0.455 ± 0.033 | 0 . 822 ± 0 . 009 (RW) | 0.681 ± 0.011     |
| 10%       | 0.490 ± 0.015 | 0 . 840 ± 0 . 004 (RW) | 0.792 ± 0.005     |
| 15%       | 0.510 ± 0.013 | 0 . 852 ± 0 . 003 (MH) | 0.847 ± 0.004     |
| 20%       | 0.541 ± 0.018 | 0.856 ± 0.003 (MH)     | 0 . 868 ± 0 . 003 |
| 25%       | 0.571 ± 0.023 | 0.856 ± 0.002 (RW)     | 0 . 873 ± 0 . 002 |
| 30%       | 0.571 ± 0.006 | 0.854 ± 0.001 (MH)     | 0 . 874 ± 0 . 002 |
| 40%       | 0.619 ± 0.015 | 0.856 ± 0.001 (RW)     | 0 . 862 ± 0 . 002 |

Table S11. Moving Window Spearman on Epinions Network

| S ratio   | UniSam   | RW Best                | TCEC              |
|-----------|----------|------------------------|-------------------|
| 5%        | NA       | 0 . 168 ± 0 . 007 (RW) | 0.119 ± 0.003     |
| 10%       | NA       | 0 . 170 ± 0 . 002 (DW) | 0.162 ± 0.004     |
| 15%       | NA       | 0.173 ± 0.003 (DW)     | 0 . 181 ± 0 . 003 |
| 20%       | NA       | 0.174 ± 0.001 (MH)     | 0 . 188 ± 0 . 002 |
| 25%       | NA       | 0.174 ± 0.002 (RW)     | 0 . 188 ± 0 . 002 |
| 30%       | NA       | 0.175 ± 0.002 (DW)     | 0 . 188 ± 0 . 001 |
| 40%       | NA       | 0.176 ± 0.002 (DW)     | 0 . 185 ± 0 . 001 |

Table S12. Moving Window Kendall on Epinions Network

| S ratio   | UniSam   | RW Best                | TCEC              |
|-----------|----------|------------------------|-------------------|
| 5%        | NA       | 0 . 114 ± 0.005 (RW)   | 0.081 ± 0.002     |
| 10%       | NA       | 0 . 115 ± 0 . 001 (DW) | 0.109 ± 0.003     |
| 15%       | NA       | 0 . 116 ± 0 . 002 (DW) | 0.122 ± 0.002     |
| 20%       | NA       | 0.117 ± 0.002 (MH)     | 0 . 127 ± 0 . 002 |
| 25%       | NA       | 0.117 ± 0.001 (RW)     | 0 . 126 ± 0 . 001 |
| 30%       | NA       | 0.118 ± 0.001 (DW) 0   | . 127 ± 0 . 001   |
| 40%       | NA       | 0.118 ± 0.001 (DW) 0 . | 124 ± 0 . 001     |

Table S13. Global Spearman on Stanford Web Network

| S ratio   | UniSam            | RW Best                  | TCEC              |
|-----------|-------------------|--------------------------|-------------------|
| 5%        | - 0.336 ± 0.832   | - 0 . 064 ± 0 . 646 (MH) | 0.228 ± 0.156     |
| 10%       | 0.083 ± 0.011     | 0 . 167 ± 0 . 015 (DW)   | - 0.254 ± 0.881   |
| 15%       | 0 . 079 ± 0 . 012 | 0.006 ± 0.671 (MH)       | - 0.055 ± 0.654   |
| 20%       | 0.075 ± 0.010     | 0 . 240 ± 0 . 031 (MH)   | 0.127 ± 0.125     |
| 25%       | - 0.125 ± 0.625   | 0 . 239 ± 0 . 010 (RW)   | 0.234 ± 0.132     |
| 30%       | - 0.122 ± 0.626   | 0.235 ± 0.012 (RW)       | 0 . 365 ± 0 . 097 |
| 40%       | 0.110 ± 0.019     | 0.209 ± 0.010 (RW)       | 0 . 355 ± 0 . 025 |

Table S14. Moving Window Spearman on Stanford Web Network

| S ratio   | UniSam   | RW Best                | TCEC                |
|-----------|----------|------------------------|---------------------|
| 5%        | NA       | - 0.018 ± 0.606 (MH)   | 0 . 005 ± 0 . 030   |
| 10%       | NA       | 0 . 020 ± 0 . 003 (DW) | - 0.372 ± 0.815     |
| 15%       | NA       | - 0.170 ± 0.610 (MH)   | - 0 . 167 ± 0 . 611 |
| 20%       | NA       | 0.035 ± 0.006 (MH)     | 0 . 054 ± 0 . 028   |
| 25%       | NA       | 0.037 ± 0.004 (MH)     | 0 . 064 ± 0 . 021   |
| 30%       | NA       | 0.037 ± 0.002 (RW)     | 0 . 067 ± 0 . 016   |
| 40%       | NA       | 0.032 ± 0.001 (RW)     | 0 . 062 ± 0 . 005   |

Table S15. Moving window Kendall on Stanford Web Network

| S ratio   | UniSam   | RW Best                  | TCEC                |
|-----------|----------|--------------------------|---------------------|
| 5%        | NA       | - 0.187 ± 0.604 (MH)     | 0 . 035 ± 0 . 020   |
| 10%       | NA       | 0 . 014 ± 0 . 002 (DW) - | 0.381 ± 0.810       |
| 15%       | NA       | - 0.179 ± 0.607 (MH)     | - 0 . 178 ± 0 . 608 |
| 20%       | NA       | 0.024 ± 0.004 (MH)       | 0 . 037 ± 0 . 019   |
| 25%       | NA       | 0.025 ± 0.003 (MH)       | 0 . 044 ± 0 . 015   |
| 30%       | NA       | 0.025 ± 0.001 (RW)       | 0 . 046 ± 0 . 011   |
| 40%       | NA       | 0.022 ± 0.001 (RW)       | 0 . 043 ± 0 . 004   |
Nicolò Ruggeri Caterina De Bacco
Sampling on networks: estimating eigenvector centrality on incomplete graphs
We develop a new sampling method to estimate eigenvector centrality on
incomplete networks. Our goal is to estimate this global centrality measure
having at disposal a limited amount of data. This is the case in many
real-world scenarios where data collection is expensive, the network is too big
for data storage capacity or only partial information is available. The
sampling algorithm is theoretically grounded by results derived from spectral
approximation theory. We studied the problem on both synthetic and real data
and tested the performance comparing with traditional methods, such as random
walk and uniform sampling. We show that approximations obtained from such
methods are not always reliable and that our algorithm, while preserving
computational scalability, improves performance under different error measures.