When is Network Lasso Accurate?

TL;DR

This paper establishes precise conditions on network structure and sampling for the accuracy of network Lasso in recovering graph signals, leveraging compressed sensing concepts.

Contribution

It derives exact topological and sampling conditions that guarantee the accuracy of network Lasso for graph signal estimation.

Findings

Conditions on network topology for accurate recovery

Quantification of error based on node connectivity

Extension of compressed sensing principles to network Lasso

Abstract

The "least absolute shrinkage and selection operator" (Lasso) method has been adapted recently for networkstructured datasets. In particular, this network Lasso method allows to learn graph signals from a small number of noisy signal samples by using the total variation of a graph signal for regularization. While efficient and scalable implementations of the network Lasso are available, only little is known about the conditions on the underlying network structure which ensure network Lasso to be accurate. By leveraging concepts of compressed sensing, we address this gap and derive precise conditions on the underlying network topology and sampling set which guarantee the network Lasso for a particular loss function to deliver an accurate estimate of the entire underlying graph signal. We also quantify the error incurred by network Lasso in terms of two constants which reflect the connectivity of the sampled nodes.

Figures

Equations

$$\Vert x[\cdot ]{\Vert }_{\mathrm{TV}}:=\sum _{\{i,j\}\in \mathcal{E}}{W}_{i,j}|x[j]-x[i]|.$$

$$\Vert x[\cdot ]{\Vert }_{\mathcal{S}}:=\sum _{\{i,j\}\in \mathcal{S}}{W}_{i,j}|x[j]-x[i]|.$$

$${\mathcal{I}}_{\mathcal{C}}[i]:=\{\begin{array}{ll}1& \text{mbox}ifi\in \mathcal{C}\\ 0& \text{mbox}else.\end{array}$$

$$x[i]=\sum _{l=1}^{|\mathcal{F}|}{a}_l{\mathcal{I}}_{{\mathcal{C}}_l}[i].$$

$$\Vert x[\cdot ]{\Vert }_{\mathrm{TV}}\le 2\underset{l\in \{1,\dots ,|\mathcal{F}|\}}{\max }|a_l|\sum _{\{i,j\}\in \partial \mathcal{F}}{W}_{i,j}.$$

$$\sum _{\{i,j\}\in \mathcal{E}}{W}_{i,j}(x[j]-x[i]{)}^2={\mathbf{x}}^T\mathbf{Lx}.$$

$$L_{i,j}=\{\begin{array}{ll}{\sum }_{k\in \mathcal{V}}{W}_{i,k}& \text{mbox}ifi=j\\ -W_{i,j}& \text{mbox}otherwise.\end{array}$$

$$\tilde{x}[l]:={\mathbf{u}}_l^T\mathbf{x}\text{mbox},forl=1,\dots ,N,$$

$$y[i]=x[i]+e[i]$$

$$\mathcal{M}:=\{i_1,\dots ,i_M\}\subseteq \mathcal{V}.$$

$$\hat{E}(\hat{x}[\cdot ]):=\sum _{i\in \mathcal{M}}|\hat{x}[i]-y[i]|.$$

$$\hat{x}[\cdot ]$$

$$\sum _{j\in \mathcal{N}(i)}h(j,i)-h(i,j)=d[i]\text{mbox},foranyi\in \mathcal{V}$$

$$|h(i,j)|\le W_{i,j}\text{mbox}foranyedge(i,j)\in \mathcal{E}\setminus \partial \mathcal{F}.$$

$$h(i,j)=b_{i,j}\cdot L\cdot W_{i,j}\text{mbox},h(j,i)={\bar{b}}_{i,j}\cdot L\cdot W_{i,j}$$

$$|d[i]|\le K\text{mbox}fori\in \mathcal{M}\text{mbox},andd[i]=0\text{mbox}fori\in \mathcal{V}\setminus \mathcal{M}.$$

$$K\ge \underset{\{i,j\}\in \partial \mathcal{F}}{\max }L{W}_{i,j}.$$

$$K=L\cdot \underset{i\in \mathcal{V}}{\max }|\mathcal{N}(i)\cap \partial \mathcal{F}|.$$

$$\Vert \hat{x}[\cdot ]-x[\cdot ]{\Vert }_{\mathrm{TV}}\le (K+4/(L-1))\sum _{i\in \mathcal{M}}|e[i]|.$$

$$(K+4/(L-1))\sum _{i\in \mathcal{M}}|e[i]|<\eta /2,$$

$$\sum _{i\in \mathcal{V}}(\tilde{x}[i]-x[i]{)}^2\le N{\epsilon }^2.$$

$$K\sum _{i\in \mathcal{M}}|z[i]|+\Vert z[\cdot ]{\Vert }_{\mathcal{E}\setminus \partial \mathcal{F}}\ge L\Vert z[\cdot ]{\Vert }_{\partial \mathcal{F}}$$

$$\sum _{j\in \mathcal{N}(i)}h(j,i)-\sum _{j\in \mathcal{N}(i)}h(i,j)=0\text{mbox}foralli\notin \mathcal{M}$$

$$\displaystyle\bigg{|}\sum_{j\in\mathcal{N}(i)}h(j,i)-\sum_{j\in\mathcal{N}(i)}h(i,j)\bigg{|}\!\leq\!K\mbox{ for all }i\in\mathcal{M}$$

$$|h(i,j)|\le W_{i,j}\text{mbox}for(i,j)\notin \partial \mathcal{F}$$

$$h(i,j)\cdot h(j,i)=0\text{mbox}for\{i,j\}\in \partial \mathcal{F}$$

$$(h(i,j)-h(j,i))(z[i]-z[j])=L{W}_{i,j}|z[i]-z[j]|$$

$$L\Vert z[\cdot ]{\Vert }_{\partial \mathcal{F}}$$

$$\stackrel{\text{eqref}eq{u}_{c}rucia{l}_{i}dent}{=}\sum _{(i,j)\in \partial \mathcal{F}}(z[i]-z[j])h(i,j).$$

$$L\Vert z[\cdot ]{\Vert }_{\partial \mathcal{F}}$$

$$=\sum _{(i,j)\in \mathcal{E}}(z[i]-z[j])h(i,j)-\sum _{(i,j)\in \mathcal{E}\setminus \partial \mathcal{F}}(z[i]-z[j])h(i,j)$$

$$=\sum _{i\in \mathcal{V}}z[i]\sum _{j\in \mathcal{N}(i)}(h(j,i)-h(i,j))$$

$$-\sum _{(i,j)\in \mathcal{E}\setminus \partial \mathcal{F}}(z[i]-z[j])h(i,j)$$

$$\stackrel{\text{eqref}eq{u}_{c}onditio{n}_{f}low}{\le }K\sum _{i\in \mathcal{M}}|z[i]|+\sum _{\{i,j\}\in \mathcal{E}\setminus \partial \mathcal{F}}|z[i]-z[j]|W_{i,j}$$

Full text

I Introduction

In many applications ranging from image processing, social networks to bioinformatics, the observed datasets carry an intrinsic network structure. Such datasets can be represented conveniently by signals defined over a “data graph” which models the network structure inherent to the dataset [30, 12]. The nodes of this data graph represent individual data points which are labeled by some quantity of interest, e.g., the class membership in a classification problem. We represent this label information as a graph signal whose value for a particular node is given by its label [6, 11, 41, 30, 16, 2, 5]. This graph signal representation of datasets allows to apply efficient methods from graph signal processing (GSP) which are obtained, in turn, by extending established methods (e.g., fast filtering and transforms) from discrete time signal processing (over chain graphs) to arbitrary graphs [31, 34, 27].

The resulting graph signals are typically clustered, i.e., these signals are nearly constant over well connected subset of nodes (clusters) in the data graph. Exploiting this clustering property enables the accurate recovery of graph signals from few noisy samples. In particular, using the total variation to measure how well a graph signal conforms with the underlying cluster structure, the authors of [19] obtain the network Lasso (nLasso) by adapting the well-known Lasso estimator which is widely used for learning sparse models [20, 21]. The nLasso can be interpreted as an instance of the regularized empirical risk minimization principle, using total variation of a graph signal for the regularization. Some applications where the use of nLasso based methods has proven beneficial include housing price prediction and personalised medicine [19, 40]

A scalable implementation of the nLasso has been obtained via the alternating direction method of multipliers (ADMM) [7]. However, the authors of [7] do not discuss conditions on the underlying network structure which ensure success of the network Lasso. We close this gap in the understanding of the performance of network Lasso, by deriving sufficient conditions on the data graph (cluster) structure and sampling set such that nLasso is accurate. To this end, we introduce a simple model for clustered graph signals which are constant over well connected groups or clusters of nodes. We then define the notion of resolving sampling sets, which relates the cluster structure of the data graph to the sampling set. Our main contribution is an upper bound on the estimation error obtained from nLasso when applied to resolving sampling sets. This upper bound depends on two numerical parameters which quantify the connectivity between sampled nodes and cluster boundaries.

Much of the existing work on recovery conditions and methods for graph signal recovery (e.g., [29, 36, 13, 14, 32, 38]), relies on spectral properties of the data graph Laplacian matrix. In contrast, our approach is based directly on the connectivity properties of the underlying network structure. The closest to our work is [33, 39], which provide sufficient conditions such that a special case of the nLasso (referred to as the “edge Lasso”) accurately recovers piece-wise constant (or clustered) graph signals from noisy observations. However, these works require access to fully labeled datasets, while we consider datasets which are only partially labeled (as it is typical for machine learning applications where label information is costly).

Outline. The problem setting considered is formalized in Section II. In particular, we show how to formulate the problem of learning a clustered graph signal from a small amount of signal samples as a convex optimization problem, which is underlying the nLasso method. Our main result, i.e., an upper bound on the estimation error of nLasso is stated in Section III. Numerical experiments which illustrate our theoretical findings are discussed in Section IV.

Notation. We will conform to standard notation of linear algebra as used, e.g., in [17]. For a binary variable $b$, we denote its negation as $\bar{b}$.

II Problem Formulation

We consider datasets which are represented by a network model, i.e., a data graph $\mathcal{G}\!=\!(\mathcal{V},\mathcal{E},\mathbf{W})$ with node set $\mathcal{V}=\{1,\ldots,N\}$, edge set $\mathcal{E}$ and weight matrix $\mathbf{W}\in\mathbb{R}_{+}^{N\times N}$. The nodes $\mathcal{V}$ of the data graph represent individual data points. For example, the node $i\!\in\!\mathcal{V}$ might represent a (super-)pixel in image processing, a neuron of a neural network [18] or a social network user profile [15].

Many applications naturally suggest a notion of similarity between individual data points, e.g., the profiles of befriended social network users or greyscale values of neighbouring image pixels. These domain-specific notions of similarity are represented by the edges of the data graph $\mathcal{G}$, i.e., the nodes $i,j\!\in\!\mathcal{V}$ representing similar data points are connected by an undirected edge $\{i,j\}\!\in\!\mathcal{E}$. We denote the neighbourhood of the node $i\in\mathcal{V}$ by $\mathcal{N}(i):=\{j\in\mathcal{V}:\{i,j\}\in\mathcal{E}\}$. It will be convenient to associate with each undirected edge $\{i,j\}$ a pair of directed edges, i.e., $(i,j)$ and $(j,i)$. With slight abuse of notation we will treat the elements of the edge set $\mathcal{E}$ either as undirected edges $\{i,j\}$ or as pairs of two directed edges $(i,j)$ and $(j,i)$.

In some applications it is possible to quantify the extent to which data points are similar, e.g., via the physical distance between neighbouring sensors in a wireless sensor network application [43]. Given two similar data points $i,j\!\in\!\mathcal{V}$, which are connected by an edge $\{i,j\}\in\mathcal{E}$, we will quantify the strength of their connection by the edge weight $W_{i,j}\!>\!0$ which we collect in the symmetric weight matrix $\mathbf{W}\in\mathbb{R}_{+}^{N\times N}$. The absence of an edge between nodes $i,j\in\mathcal{V}$ is encoded by a zero weight $W_{i,j}\!=\!0$. Thus the edge structure of the data graph $\mathcal{G}$ is fully specified by the support (locations of the non-zero entries) of the weight matrix $\mathbf{W}$.

II-A Graph Signals

Beside the network structure, encoded in the data graph $\mathcal{G}$, datasets typically also contain additional labeling information. We represent this additional label information by a graph signal defined over $\mathcal{G}$. A graph signal $x[\cdot]$ is a mapping $\mathcal{V}\rightarrow\mathbb{R}$, which associates every node $i\!\in\!\mathcal{V}$ with the signal value $x[{i}]\!\in\!\mathbb{R}$ (which might representing a label characterizing the data point). We denote the set of all graph signals defined over a graph $\mathcal{G}=(\mathcal{V},\mathcal{E},\mathbf{W})$ by $\mathbb{R}^{\mathcal{V}}$.

Many machine learning methods for network structured data rely on a “cluster hypothesis” [11]. In particular, we assume the graph signals $x[\cdot]$ representing the label information of a dataset conforms with the cluster structure of the underlying data graph. Thus, any two nodes $i,j\in\mathcal{V}$ out of a well-connected region (“cluster”) of the data graph tend to have similar signal values, i.e., $x[{i}]\approx x[{j}]$. Two important application domains where this cluster hypothesis has been applied successfully are digital signal processing where time samples at adjacent time instants are strongly correlated for sufficiently high sampling rate (cf. Fig. 1-(a)) as well as processing of natural images whose close-by pixels tend to be coloured likely (cf. Fig. 1-(b)). The cluster hypothesis is verified also often in social networks where the clusters are cliques of individuals having similar properties (cf. Fig. 1-(c) and [28, Chap. 3]).

In what follows, we quantify the extend to which a graph signal $x[\cdot]\in\mathbb{R}^{\mathcal{V}}$ conforms with the clustering structure of the data graph $\mathcal{G}=(\mathcal{V},\mathcal{E},\mathbf{W})$ using its total variation (TV)

[TABLE]

For a subset of edges $\mathcal{S}\subseteq\mathcal{E}$, we use the shorthand

[TABLE]

For a supervised machine learning application, the signal values $x[{i}]$ might represent class membership in a classification problem or the target (output) value in a regression problem. For the house price example considered in [19], the vector-valued graph signal $\mathbf{x}[i]$ corresponds to a regression weight vector for a local pricing model (used for the house market in a limited geographical area represented by the node $i$).

Consider a partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ of the data graph $\mathcal{G}$ into $|\mathcal{F}|$ disjoint subsets $\mathcal{C}_{l}$ of nodes (“clusters”) such that $\mathcal{V}=\bigcup_{l=1}^{|\mathcal{F}|}\mathcal{C}_{l}$. We associate a subset $\mathcal{C}\subseteq\mathcal{V}$ of nodes with a particular “indicator” graph signal

[TABLE]

A simple model of clustered graph signals is then obtained by piece-wise constant or clustered graph signals of the form

[TABLE]

In Figure 2, we depict a clustered graph signal for a chain graph with $10$ nodes which are partitioned into two clusters: $\mathcal{C}_{1}$ and $\mathcal{C}_{2}$.

It will be convenient to define, for a given partition $\mathcal{F}$, its boundary $\partial\mathcal{F}\subseteq\mathcal{E}$ as the set of edges $\{i,j\}\in\mathcal{E}$ which connect nodes $i\!\in\!\mathcal{C}_{a}$ and $j\!\in\!\mathcal{C}_{b}$ from different clusters, i.e., with $\mathcal{C}_{a}\!\neq\!\mathcal{C}_{b}$. With a slight abuse of notation, we will use the same symbol $\partial\mathcal{F}$ also to denote the set of nodes which are connected to a node from another cluster.

The TV of a clustered graph signal of the form (4) can be upper bounded as

[TABLE]

Thus, for a partition $\mathcal{F}$ with small weighted boundary $\sum_{\{i,j\}\in\partial\mathcal{F}}W_{i,j}$, the associated clustered graph signals (4) have small TV $\|x[\cdot]\|_{\rm TV}$ due to (5).

The signal model (4), which also has been used in [33, 39], is closely related to the stochastic block model (SBM) [26]. Indeed, the SBM is obtained from (4) by choosing the coefficients $a_{\mathcal{C}}$ uniquely for each cluster, i.e., $a_{\mathcal{C}}\in\{1,\ldots,|\mathcal{F}|\}$. Moreover, the SBM provides a generative (stochastic) model for the edges within and between the clusters $\mathcal{C}_{l}$.

We highlight that the clustered signal model (4) is somewhat dual to the model of band-limited graph signals [29, 29, 13, 11, 41, 30, 16, 2]. The model of band-limited graph signals is obtained by the subspaces spanned by the eigenvectors of the graph Laplacian corresponding to the smallest (in magnitude) eigenvalues, i.e., the low-frequency components. Such band-limited graph signals are smooth in the sense of small values of the Laplacian quadratic form [3]

[TABLE]

Here, we used the vector representation $\mathbf{x}=(x[{1}],\ldots,x[{N}])^{T}$ of the graph signal $x[\cdot]$ and the graph Laplacian matrix $\mathbf{L}\in\mathbb{R}^{N\times N}$ defined element-wise as

[TABLE]

A band-limited graph signal $x[\cdot]$ is characterized by a clustering (within a small bandwidth) of their graph Fourier transform (GFT) coefficients [38]

[TABLE]

with the orthonormal eigenvectors $\{\mathbf{u}_{l}\}_{l=1}^{N}$ of the graph Laplacian matrix $\mathbf{L}$. In particular, by the spectral decomposition of the psd graph Laplacian matrix $\mathbf{L}$ (cf. (7)), we have $\mathbf{L}=\mathbf{U}{\bm{\Lambda}}\mathbf{U}^{H}$ with \mathbf{U}=\big{(}\mathbf{u}_{1},\ldots,\mathbf{u}_{N}\big{)} and the diagonal matrix ${\bm{\Lambda}}$ having (in decreasing order) the non-negative eigenvalues $\lambda_{l}$ of $\mathbf{L}$ on its diagonal.

In contrast to band-limited graph signals, a clustered graph signal of the form (4) will typically have GFT coefficients which are spread out over the entire (graph) frequency range. Moreover, while band-limited graph signals are characterized by having a sparse GFT, a clustered graph signal of the form (4) has a dense (non-sparse) GFT in general. On the other hand, while a clustered graph signal of the form (4) has sparse signal differences $\{x[i]\!-\!x[j]\}_{\{i,j\}\in\mathcal{E}}$, the signal differences of a band-limited graph signal are dense (non-sparse).

Let us illustrate the duality between the clustered graph signal model (4) and the model of band-limited graph signals (cf. [29, 29, 2]) by considering a dataset representing a finite length segment of a time series. The data graph $\mathcal{G}_{0}$ underlying this time series data is chosen as a chain graph (cf. Fig. 2), consisting of $N=100$ nodes which represent the individual time samples. The time series is partitioned into two clusters $\mathcal{C}_{1},\mathcal{C}_{2}$, each cluster consisting of $50$ consecutive nodes (time samples). We model the correlations between successive time samples using edge weight $W_{i,j}=1$ for data points $i,j$ belonging to the same cluster and a smaller weight $W_{i,j}=1/2$ for the single edge $\{i,j\}$ connecting the two clusters $\mathcal{C}_{1}$ and $\mathcal{C}_{2}$.

A clustered graph signal (time series) $x_{0}[i]=a_{1}\mathcal{I}_{\mathcal{C}_{1}}[i]+a_{2}\mathcal{I}_{\mathcal{C}_{2}}[i]$ (cf. (4) ) defined over $\mathcal{G}_{0}$ is characterized by very sparse signal differences $\{x_{0}[i]\!-\!x_{0}[j]\}_{\{i,j\}\in\mathcal{E}}$. Indeed the signal difference $x_{0}[i]\!-\!\mathbf{x}_{0}[j]$ of the clustered graph signal $\mathbf{x}_{0}[\cdot]$ is non-zero only for the single edge $\{i,j\}$ which connects $\mathcal{C}_{1}$ and $\mathcal{C}_{2}$. In stark contrast, the GFT of $x_{0}[\cdot]$ is spread out over the entire (graph) frequency range (cf. Fig. 3), i.e., the graph signal $x_{0}[\cdot]$ does not conform with the band-limited signal model.

On the other hand, we illustrate in Fig. 4 a graph signal $x_{\rm BL}[\cdot]$ with GFT coefficients $\tilde{x}_{\rm BL}[l]=1$ (cf. (8)) for $l=1,2$ and $\tilde{x}_{\rm BL}[l]=0$ otherwise. Thus, the graph signal is clearly band-limited (it has only two non-zero GFT coefficients) but the signal differences $x_{\rm BL}[i]-x_{\rm BL}[j]$ across the edges $\{i,j\}\in\mathcal{E}$ are clearly non-sparse.

II-B Recovery via nLasso

Given a dataset with data graph $\mathcal{G}=(\mathcal{V},\mathcal{E},\mathbf{W})$, we aim at recovering a graph signal $x[\cdot]\in\mathbb{R}^{\mathcal{V}}$ from its noisy values

[TABLE]

provided on a (small) sampling set

[TABLE]

Typically $M\ll N$, i.e., the sampling set is a small subset of all nodes in the data graph $\mathcal{G}$.

The recovered graph signal $\hat{x}[\cdot]$ should incur only a small empirical (or training) error

[TABLE]

Note that the definition (11) of the empirical error involves the $\ell_{1}$-norm of the deviation $\hat{x}[\cdot]{i}\!-\!y[i]$ between recovered and measured signal samples. This is different from the error criterion used in the ordinary Lasso, i.e., the squared-error loss $\sum_{i\in\mathcal{M}}(\hat{x}[{i}]\!-\!y[i])^{2}$ [9]. The definition (11) is beneficial for applications with measurement errors $e_{i}$ (cf. (9)) having mainly small values except for a few large outliers [36, 10]. However, by contrast to plain Lasso, the error function in (11) does not satisfy a restricted strong convexity property [1], which might be detrimental for the convergence speed of the resulting recovery methods (cf. Section IV).

In order to recover a clustered graph signal with a small TV $\|\hat{x}[\cdot]\|_{\rm TV}$ (cf. (5)) from the noisy signal samples $\{y[i]\}_{i\in\mathcal{M}}$ it is sensible to consider the recovery problem

[TABLE]

This recovery problem amounts to a convex optimization problem [8], which, as the notation already indicates, might have multiple solutions $\hat{x}[\cdot]$ (which form a convex set). In what follows, we will derive conditions on the sampling set $\mathcal{M}$ such that any solution $\hat{x}[\cdot]$ of (12) allows to accurately recover clustered a graph signal $x[\cdot]$ of the form (4).

Any graph signal obtained from (12) balances the empirical error $\widehat{E}(\hat{x}[\cdot])$ with the TV $\|\hat{x}[\cdot]\|_{\rm TV}$ in an optimal manner. The parameter $\lambda$ in (12) allows to trade off a small empirical error against the amount to which the resulting signal is clustered, i.e., having a small TV. In particular, choosing a small value for $\lambda$ enforces the solutions of (12) to yield a small empirical error, whereas choosing a large value for $\lambda$ enforces the solutions of (12) to have small TV. Our analysis in Section III provides a selection criterion for the parameter $\lambda$ which is based on the location of the sampling set $\mathcal{M}$ (cf. (10)) and the partition $\mathcal{F}$ underlying the clustered graph signal model (4). Alternatively, for sufficiently large sampling sets one might choose $\lambda$ using a cross-validation procedure [20].

Note that the recovery problem (12) is a particular instance of the generic nLasso problem studied in [19]. There exist efficient convex optimization methods for solving the nLasso problem (12) (cf. [44] and the references therein). In particular, the alternating method of multipliers (ADMM) has been applied to the nLasso problem in [19] to obtain a scalable learning algorithm which can cope with massive heterogeneous datasets.

III When is Network Lasso Accurate?

The accuracy of graph signal recovery methods based on the nLasso problem (12), depends on how close the solutions $\hat{x}[\cdot]$ of (12) are to the true underlying graph signal $x[\cdot]\in\mathbb{R}^{\mathcal{V}}$. In what follows, we present a condition which guarantees any solution $\hat{x}[\cdot]$ of (12) to be close to the underlying graph signal $x[\cdot]$ if it is clustered of the form (4).

A main contribution of this paper is the insight that the accuracy of nLasso methods, aiming at solving (12), depends on the topology of the underlying data graph via the existence of certain flows with demands [23]. Given a data graph $\mathcal{G}$, we define a flow on it as a mapping $h[\cdot]:\mathcal{V}\times\mathcal{V}\rightarrow\mathbb{R}$ which assigns each directed edge $(i,j)$ the value $h[(i,j)]$, which can be interpreted as the amount of some quantity flowing through the edge $(i,j)$ [23]. A flow with demands has to satisfy the conservation law

[TABLE]

with a prescribed demand $d[i]$ for each node $i\in\mathcal{V}$. Moreover, we require flows to satisfy the capacity constraints

[TABLE]

Note that the capacity constraint (14) applies only to intra-cluster edges and does not involve the boundary edges $\partial\mathcal{F}$. The flow values $h(i,j)$ at the boundary edges $(i,j)\in\partial\mathcal{F}$ take a special role in the following definition of the notion of resolving sampling sets.

Definition 1**.**

Consider a dataset with data graph $\mathcal{G}=(\mathcal{V},\mathcal{E},\mathbf{W})$ which contains the sampling set $\mathcal{M}\subseteq\mathcal{V}$. The sampling set $\mathcal{M}$ resolves a partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ with constants $K$ and $L$ if, for any $b_{i,j}\in\{0,1\}$ with $\{i,j\}\!\in\!\partial\mathcal{F}$, there exists a flow $h[\cdot]$ on $\mathcal{G}$ (cf. (13), (14)) with

[TABLE]

for every boundary edge $\{i,j\}\in\partial\mathcal{F}$ and demands (cf. (13)) satisfying

[TABLE]

This definition requires nodes of a resolving sampling set to be sufficiently well connected with every boundary edge $\{i,j\}\in\partial\mathcal{F}$. In particular, we could think of injecting (absorbing) certain amounts of flow into (from) the data graph at the sampled nodes. At each sampled node $i\in\mathcal{M}$, we can inject (absorb) a flow of level at most $K$ (cf. (16)). The injected (absorbed) flow has to be routed from the sampled nodes via the intra-cluster edges to each boundary edge such that it carries a flow value $L\cdot W_{i,j}$. Clearly, this is only possible if there are paths of sufficient capacity between sampled nodes and boundary edges available.

The definition of resolving sampling sets is quantitive as it involves the numerical constants $K$ and $L$. Our main result stated below is an upper bound on the estimation error of nLasso methods which depends on the value of these constants. It will turn out that resolving sampling sets with a small values of $K$ and large values of $L$ are beneficial for the ability of nLasso to recover the entire graph signal from noisy samples observed on the sampling set. However, the constants $K$ and $L$ are coupled via the flow $h[\cdot]$ used in Definition 1. E.g., the constant $K$ always has to satisfy

[TABLE]

Thus, the minimum possible value for $K$ depends on the values of the edge weights $W_{i,j}$ of the data graph. Moreover, the minimum possible value for $L$ depends on the precise connectivity of sampled nodes with the boundary edges $\partial\mathcal{F}$. Indeed, Definition 1 requires to route (by satisfying the capacity constraints (14)), an amount of flow given by $LW_{i,j}$ from a boundary edge $\{i,j\}\in\partial\mathcal{F}$ to the sampled nodes in $\mathcal{M}$.

In order to make (the somewhat abstract) Definition 1 more transparent, let us state an easy-to-check sufficient condition for a sampling set $\mathcal{M}$ such that it resolves a given partition $\mathcal{F}$.

Lemma 2**.**

Consider a partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ of the data graph $\mathcal{G}$ which contains the sampling set $\mathcal{M}\subseteq\mathcal{V}$. If each boundary edge $\{i,j\}\in\partial\mathcal{F}$ with $i\!\in\!\mathcal{C}_{a}$, $j\!\in\!\mathcal{C}_{b}$ is connected to sampled nodes, i.e., $\{m,i\}\!\in\!\mathcal{E}$ and $\{n,j\}\!\in\!\mathcal{E}$ with $m\!\in\!\mathcal{M}\!\cap\!\mathcal{C}_{a}$, $n\!\in\!\mathcal{M}\!\cap\!\mathcal{C}_{b}$, and weights $W_{m,i},W_{n,j}\geq LW_{i,j}$, then the sampling set $\mathcal{M}$ resolves the partition $\mathcal{F}$ with constants $L$ and

[TABLE]

In Fig. 1-(c) we depict a data graph consisting of two clusters $\mathcal{F}=\{\mathcal{C}_{1},\mathcal{C}_{2}\}$. The data graph contains the sampling set $\mathcal{M}=\{m,n\}$ which resolves the partition $\mathcal{F}$ with constants $K=L=4$ according to Lemma 2.

The sufficient condition provided by Lemma 2 can be used to guide the choice for the sampling set $\mathcal{M}$. In particular Lemma 2 suggests to sample more densely near the boundary edges $\partial\mathcal{F}$ which connect different clusters. This rationale allows to cope with applications where the underlying partition $\mathcal{F}$ is unknown. In particular, we could use highly scalable local clustering methods (cf. [35]) to find the cluster boundaries $\partial\mathcal{F}$ and then select the sampled nodes in their vicinity. Another approach to cope with lack of information about $\mathcal{F}$ is based on using random walks to identify the subset of nodes with a large boundary which are sampled more densely [4].

We now state our main result which is that solutions of the nLasso problem (12) allow to accurately recover the true underlying clustered graph signal $x[\cdot]$ (conforming with the partition $\mathcal{F}$ (cf. (4)) from the noisy measurements (9) whenever the sampling set $\mathcal{M}$ resolves the partition $\mathcal{F}$.

Theorem 3**.**

Consider a clustered graph signal $x[\cdot]$ of the form (4), with underlying partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ of the data graph into disjoint clusters $\mathcal{C}_{l}$. We observe the noisy signal values $y[i]$ at the samples nodes $\mathcal{M}\subseteq\mathcal{V}$ (cf. (9)). If the sampling set $\mathcal{M}$ resolves the partition $\mathcal{F}$ with parameters $K>0,L>1$, any solution $\hat{x}[\cdot]$ of the nLasso problem (12) with $\lambda\!:=\!1/K$ satisfies

[TABLE]

Thus, if the sampling set $\mathcal{M}$ is chosen such that it resolves the partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ (cf. Definition 1), nLasso methods (cf. (12)) recover a clustered graph signal $\mathbf{x}[\cdot]$ (cf. (4)) with an accuracy which is determined by the level of the measurement noise $e[i]$ (cf. (9)).

Let us highlight that the knowledge of the partition $\mathcal{F}$ underlying the clustered graph signal model (4) is only needed for the analysis of nLasso methods leading to Theorem 3. In contrast, the actual implementation methods of nLasso methods based on (12) does not require any knowledge of the underlying partition. What is more, if the true underlying graph signal $\mathbf{x}[\cdot]$ is clustered according to (4) with different signal values $a_{l}$ for different clusters $\mathcal{C}_{l}$, the solutions of the nLasso (12) could be used for determining the clusters $\mathcal{C}_{l}$ which constitute the partition $\mathcal{F}$.

We also note that the bound (19) characterizes the recovery error in terms of the semi-norm $\|\hat{x}[\cdot]\!-\!x[\cdot]\|_{\rm TV}$ which is agnostic towards a constant offset in the recovered graph signal $\hat{x}[\cdot]$. In particular, having a small value of $\|\hat{x}[\cdot]\!-\!x[\cdot]\|_{\rm TV}$ does in general not imply a small squared error $\sum_{i\in\mathcal{V}}(\hat{x}[i]\!-\!x[{i}])^{2}$ as there might be an arbitrarily large constant offset contained in the nLasso solution $\hat{x}[\cdot]$.

However, if the error $\|\hat{x}[\cdot]\!-\!x[\cdot]\|_{\rm TV}$ is sufficiently small, we might be able to identify the boundary edges $\{i,j\}\in\partial\mathcal{F}$ of the partition $\mathcal{F}$ underlying a clustered graph signal of the form (4).

Indeed, for a clustered graph signal of the form (4), the signal difference $x[{i}]-x[{j}]$ across edges is non-zero only for boundary edges $\{i,j\}\in\partial\mathcal{F}$. Lets assume the signal differences of $x[\cdot]$ across boundary edges $\{i,j\}\in\mathcal{F}$ are lower bounded by some positive constant $\eta>0$ and the nLasso error satisfies $\|\hat{x}[\cdot]\!-\!x[\cdot]\|_{\rm TV}<\eta/2$. As can be verified easily, we can then perfectly recover the boundary $\partial\mathcal{F}$ of the partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ as precisely those edges $\{i,j\}\in\mathcal{E}$ for which $|\hat{x}[{i}]-\hat{x}[{j}]|\geq\eta/2$. Given the boundary $\partial\mathcal{F}$, we can recover the partition $\mathcal{F}$ and, in turn, average the noisy observations $y[i]$ over all sampled nodes $i\in\mathcal{M}$ belonging to the same cluster. This simple post-processing of the nLasso estimate $\hat{x}[i]$ is summarized in Algorithm 1.

Lemma 4**.**

Consider the setting of Theorem 3 involving a clustered graph signal $x[\cdot]$ of the form (4) with coefficients $a_{l}$ satisfying $|a_{l}-a_{l^{\prime}}|>\eta$ for $l\neq l^{\prime}$ with a known positive threshold $\eta>0$. We observe noisy signal samples $y[i]$ (cf. (9)) over the sampling set $\mathcal{M}$ with a bounded error $e[i]\leq\epsilon$. If the sampling set $\mathcal{M}$ resolves the partition $\mathcal{F}$ with parameters $K>0,L>1$ such that

[TABLE]

then the signal $\tilde{x}[\cdot]$ delivered by Algorithm 1 satisfies

[TABLE]

IV Numerical Experiments

In order to illustrate the theoretical findings of Section III we report the results of some illustrative numerical experiments involving the recovery of clustered graph signals of the form (4) from a small number of noisy measurements (9). To this end, we implemented the iterative method ADMM [7] to solve the nLasso (12) problem. We applied the resulting semi-supervised learning algorithm to two synthetically generated data sets. The first data set represents a time series, which can be represented as a graph signal over a chain graph. The nodes of the chain graph, which represent the discrete time instants are partitioned evenly into clusters of consecutive nodes. A second experiment is based on data sets generated using a recently proposed generative model for complex networks.

IV-A Chain Graph

Our first experiment, is based on a graph signal defined over a chain graph $\mathcal{G}_{\rm chain}$ (cf. Fig. 2) with $N=10^{5}$ nodes $\mathcal{V}=\{1,2,\dots,N\}$, connected by $N-1$ undirected edges. The nodes of the data graph $\mathcal{G}_{\rm chain}$ are partitioned into $N/10$ equal-sized clusters $\mathcal{C}_{l}$, $l=1,\dots,N/10$, each constituted by $10$ consecutive nodes. The intrinsic clustering structure of the chain graph $\mathcal{G}_{\rm chain}$ matches the partition $\mathcal{F}_{\rm chain}=\{\mathcal{C}_{l}\}_{l=1}^{N/10}$ via the edge weights $W_{i,j}$. In particular, the weights of the edges connecting nodes within the same cluster are chosen i.i.d. according to $W_{i,j}\sim|\mathcal{N}(2,1/4)|$ (i.e., the absolute value of a Gaussian random variable with mean 2 and variance 1/4). The weights of the edges connecting nodes from different clusters are chosen i.i.d. according to $W_{i,j}\sim|\mathcal{N}(1,1/4)|$.

We then generate a clustered graph signal $x[\cdot]$ of the form (4) with coefficients $a_{l}\in\{1,5\}$, where the coefficients $a_{l}$ and $a_{l^{\prime}}$ of consecutive clusters $\mathcal{C}_{l}$ and $\mathcal{C}_{l^{\prime}}$ are different. The graph signal $x[\cdot]$ is observed via noisy samples $y[i]$ (cf. (9) with $e[i]\sim\mathcal{N}(0,1/4)$) obtained for the nodes $i\in\mathcal{V}$ belonging to a sampling set $\mathcal{M}$. We consider two different choices for the sampling set, i.e., $\mathcal{M}=\mathcal{M}_{1}$ and $\mathcal{M}=\mathcal{M}_{2}$. Both choices contain the same number of nodes, i.e., $|\mathcal{M}_{1}|=|\mathcal{M}_{2}|=2\cdot 10^{4}$. The sampling set $\mathcal{M}_{1}$ contains neighbours of cluster boundaries $\partial\mathcal{F}_{\rm chain}$ and conforms to Lemma 2 with constants $K=5.39$ and $L=2$ (which have been determined numerically). In contrast, the sampling set $\mathcal{M}_{2}$ is obtained by selecting nodes uniformly at random from $\mathcal{V}$ and thereby completely ignoring the cluster structure $\mathcal{F}_{\rm chain}$ of $\mathcal{G}_{\rm chain}$.

The noisy measurements $y[i]$ are then input to an ADMM implementation for solving the nLasso problem (12) with $\lambda=1/K$. We run ADMM for a fixed number of $300$ iterations and using ADMM-parameter $\rho=0.01$ [7]. In Fig. 5 we illustrate the recovered graph signals (over the first $100$ nodes of the chain graph) $\hat{x}[\cdot]$, obtained from noisy signal samples over either sampling set $\mathcal{M}_{1}$ or $\mathcal{M}_{2}$.

As evident from Fig. 5, the recovered signal obtained when using the sampling set $\mathcal{M}_{1}$, which takes the partition $\mathcal{F}_{\rm chain}$ into account, better resembles the original graph signal $x[\cdot]$ than when using the randomly selected sampling set $\mathcal{M}_{2}$. The favourable performance of $\mathcal{M}_{1}$ is also reflected in the empirical normalized mean squared errors (NMSE) between the real and recovered graph signals, which are ${\rm NMSE}_{\mathcal{M}_{1}}\!=\!3.3\cdot 10^{-2}$ and ${\rm NMSE}_{\mathcal{M}_{2}}\!=\!2.192\cdot 10^{-1}$, respectively.

We have repeated the above experiment with the same parameters but considering noiseless initial samples $y[i]$ for both sampling sets $\mathcal{M}_{1}$ and $\mathcal{M}_{2}$. The recovered graph signals $\hat{x}[\cdot]$ for the first $100$ nodes of the chain are presented in Fig. 6. It can be observed that the recovery starting from the sampling set $\mathcal{M}_{1}$ (conforming to the partition $\mathcal{F}_{\rm chain}$) perfectly resembles the original graph signal $x[\cdot]$, as expected according to our upper bound in (19). The NMSE obtained after running ADMM for $300$ iterations for solving the nLasso problem (12) are ${\rm NMSE}_{\mathcal{M}_{1}}\!=\!7.5\cdot 10^{-6}$ and ${\rm NMSE}_{\mathcal{M}_{2}}\!=\!1.475\cdot 10^{-1}$, respectively.

IV-B Complex Network

In this second experiment, we generate a data graph $\mathcal{G}_{\rm lfr}$ using the generative model introduced by Lancichinetti, Fortunato, and Radicchi [24], in what follows referred to as LFR model. The LFR model aims at imitating some key characteristics of real-world networks such as power law distributions of node degrees and community sizes. The data graph $\mathcal{G}_{\rm lfr}$ contains a total of $N=10^{5}$ nodes which are partitioned into $1399$ clusters, $\mathcal{F}_{\rm lfr}=\{\mathcal{C}_{1},\dots,\mathcal{C}_{1399}\}$. The nodes $\mathcal{V}$ of $\mathcal{G}_{\rm lfr}$ are connected by a total of $9.45\cdot 10^{5}$ undirected edges $\mathcal{E}$.

The edge weights $W_{i,j}$, which are also provided by the LFR model, conform to the cluster structure of $\mathcal{G}_{\rm lfr}$, i.e., inter-cluster edges $\{i,j\}\in\mathcal{E}$ with $i,j\in\mathcal{C}_{l}$ have larger weights compared to intra-cluster edges $\{i,j\}\in\mathcal{E}$ with $i\in\mathcal{C}_{l}$ and $j\in\mathcal{C}_{l^{\prime}}$. Given the data graph $\mathcal{G}_{\rm lfr}$ and partition $\mathcal{F}_{\rm lfr}$ we generate a clustered graph signal according to (4) as $x[{i}]=\sum_{j=1}^{1399}a_{j}\mathcal{I}_{\mathcal{C}_{j}}[i]$ with coefficients $a_{j}$ randomly chosen i.i.d. according to a uniform distribution $\mathcal{U}(1,50)$.

We then try to recover the entire graph signal $x[\cdot]$ by solving the nLasso problem (12) using noisy measurements $y[i]$, according to (9) with i.i.d. measurement noise $e[i]\sim\mathcal{N}(0,1/4)$, obtained at the nodes in a sampling set $\mathcal{M}$. As in Section IV-A, we consider two different choices $\mathcal{M}_{1}$ and $\mathcal{M}_{2}$ for the sampling set which both contain the same number of nodes, i.e., $|\mathcal{M}_{1}|=|\mathcal{M}_{2}|=10^{4}$. The nodes in sampling set $\mathcal{M}_{1}$ are selected according to Lemma 2, i.e., by choosing nodes which are well connected (close) to boundary edges $\partial\mathcal{F}_{\rm lfr}$ which connect different clusters of the partition $\mathcal{F}_{\rm lfr}$. In contrast, the sampling set $\mathcal{M}_{2}$ is constructed by selecting nodes uniformly at random, i.e., the partition $\mathcal{F}_{\rm lfr}$ is not taken into account.

In order to construct the sampling set $\mathcal{M}_{1}$, we first sorted the edges $\{i,j\}\in\mathcal{E}$ of the data graph $\mathcal{G}_{\rm lfr}$ in ascending order according to their edge weight $W_{i,j}$. We then iterate over the the edges according to the list, starting with the edge having smallest weight, and for each edge $\{i,j\}\in\mathcal{E}$ we select the neighbouring nodes of $i$ and $j$ with highest degree and add them to $\mathcal{M}_{1}$, if they are not already included there. This process continues until the sampling set $\mathcal{M}_{1}$ has reached the prescribed size of $10^{4}$. Using Lemma 2, we then verified numerically that the sampling set $\mathcal{M}_{1}$ resolves $\mathcal{F}_{\rm lfr}$ with constants $K=142.6$ and $L=2$ (cf. Definition 1).

The measurements $y[i]$ collected for each sampling sets $\mathcal{M}_{1}$ and $\mathcal{M}_{2}$ are fed into the ADMM algorithm (using parameters $\rho=1/100$) for solving the nLasso problem (12) with $\lambda=1/K$. The evolution of the NMSE achieved by the ADMM output for an increasing number the iterations is shown in Fig. 7. According to Fig. 7 the signal recovered from the sampling set $\mathcal{M}_{1}$ approximates the true graph signal $x[\cdot]$ more closely compared to when using the sampling set $\mathcal{M}_{2}$. The NMSE achieved after $300$ iterations of ADMM is ${\rm NMSE}_{\mathcal{M}_{1}}\!=\!1.56\cdot 10^{-2}$ and ${\rm NMSE}_{\mathcal{M}_{2}}\!=\!4.25\cdot 10^{-2}$, respectively.

Finally, we compare the recovery accuracy of nLasso to that of plain label propagation (LP) [42], which relies on a band-limited signal model (cf. Section II-A). In particular, LP quantifies signal smoothness by the Laplacian quadratic form (6) instead of the total variation (1), which underlies nLasso (12). The signals recovered after running the LP algorithm for $300$ iterations for the two sampling sets $\mathcal{M}_{1}$ and $\mathcal{M}_{2}$ incur an NMSE of ${\rm NMSE}_{\mathcal{M}_{1}}\!=\!3.1\cdot 10^{-2}$ and ${\rm NMSE}_{\mathcal{M}_{2}}\!=\!7.43\cdot 10^{-2}$, respectively. Thus, the signals recovered using nLasso are more accurate compared to LP, as illustrated in Fig. 8. However, our results indicate that LP also benefits by using the sampling set $\mathcal{M}_{1}$ whose construction is guided by our theoretical findings (cf. Lemma 2).

V Proofs

The high-level idea behind the proof of Theorem 3 is to adapt the concept of compatibility conditions for Lasso type estimators [9]. This concept has been championed for analyzing Lasso type methods [9]. Our main technical contribution is to verify the compatibility condition for a sampling set $\mathcal{M}$ which resolves the partition $\mathcal{F}$ underlying the signal model (4) (cf. Lemma 6 below).

V-A The Network Compatibility Condition

As an intermediate step towards proving Theorem 3, we adopt the compatibility condition [37], which has been introduced to analyze Lasso methods for learning sparse signals, to the clustered graph signal model (4). In particular, we define the network compatibility condition for sampling graph signals with small total variation (cf. (1)).

Definition 5**.**

Consider a data graph $\mathcal{G}=(\mathcal{V},\mathcal{E},\mathbf{W})$ whose nodes $\mathcal{V}$ are partitioned into disjoint clusters $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$. A sampling set $\mathcal{M}\subseteq\mathcal{V}$ is said to satisfy the network compatibility condition, with constants $K,L>0$, if

[TABLE]

for any graph signal $z[\cdot]\in\mathbb{R}^{\mathcal{V}}$.

It turns out that any sampling set $\mathcal{M}$ which resolves the partition $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$ with constants $K$ and $L$ (cf. Definition 1) also satisfies the network compatibility condition (22) with the same constants.

Lemma 6**.**

Any sampling set $\mathcal{M}$ which resolves the partition $\mathcal{F}$ with parameters $K,L>0$ satisfies the network compatibility condition with parameters $K,L$.

Proof.

Let us consider an arbitrary but fixed graph signal $z[\cdot]\in\mathbb{R}^{\mathcal{V}}$. Since the sampling set $\mathcal{M}$ resolves the partition $\mathcal{F}$ there exists a flow $h[e]$ on $\mathcal{G}$ with (cf. Definition 1)

[TABLE]

Moreover, due to (15), we have the important identity

[TABLE]

which holds for all boundary edges $\{i,j\}\!\in\!\partial\mathcal{F}$. This yields, in turn,

[TABLE]

Since \mathcal{E}\!=\!\partial\mathcal{F}\!\cup\!\big{(}\mathcal{E}\!\setminus\!\partial\mathcal{F}\big{)}, we can develop (V-A) as

[TABLE]

which verifies (22). ∎

The next result shows that if the sampling set satisfies the network compatibility condition, any solution of the nLasso (12) allows to accurately recover a clustered graph signal (cf. (4)).

Lemma 7**.**

Consider a clustered graph signal $x[\cdot]$ of the form (4) defined on the data graph $\mathcal{G}=(\mathcal{V},\mathcal{E},\mathbf{W})$ whose nodes $\mathcal{V}$ are partitioned into the clusters $\mathcal{F}=\{\mathcal{C}_{1},\ldots,\mathcal{C}_{|\mathcal{F}|}\}$. We observe the noisy signal values $y[i]$ at the sampled nodes $\mathcal{M}\!\subseteq\!\mathcal{V}$ (cf. (9)). If the sampling set $\mathcal{M}$ satisfies the network compatibility condition with constants $L>1,K>0$, then any solution of the nLasso problem (12), for the choice $\lambda:=1/K$, satisfies

[TABLE]

Proof.

Consider a solution $\hat{x}[\cdot]$ of the nLasso problem (12) which is different from the true underlying clustered signal $x[\cdot]$ (cf. (4)). We must have (cf. (9))

[TABLE]

since otherwise the true underlying signal $x[\cdot]$ would achieve a smaller objective value in (12) which, in turn, would contradict the premise that $\hat{x}[\cdot]$ is optimal for the problem (12).

Let us denote the difference between the solution $\hat{x}[\cdot]$ of (12) and the true underlying clustered signal $x[\cdot]$ by $\tilde{x}[\cdot]:=\hat{x}[\cdot]-x[\cdot]$. Since $x[\cdot]$ satisfies (4),

[TABLE]

Applying the decomposition property of the semi-norm $\|\cdot\|_{\rm TV}$ to (28) yields

[TABLE]

Therefore, using (29) and the triangle inequality,

[TABLE]

Since $\sum_{i\in\mathcal{M}}|\hat{x}[i]-y[i]|\geq 0$, (V-A) yields

[TABLE]

i.e., for sufficiently small measurement noise $e[i]$, the signal differences of the recovery error $\tilde{x}[\cdot]\!=\!\hat{x}[\cdot]\!-\!x[\cdot]$ cannot be concentrated across the edges within the clusters $\mathcal{C}_{l}$. Moreover, using

[TABLE]

the inequality (V-A) becomes

[TABLE]

Thus, since the sampling set $\mathcal{M}$ satisfies the network compatibility condition, we can apply (22) to $\tilde{\mathbf{x}}[\cdot]$ yielding

[TABLE]

Inserting (35) into (34), with $\lambda=1/K$, yields

[TABLE]

Combining (32) and (36) yields

[TABLE]

∎

V-B Proof of Theorem 3

Combine Lemma 6 with Lemma 7.

VI Conclusions

Given a known cluster structure of the data graph, we introduced the notion of resolving sampling sets. A sampling set resolves a cluster structure if there exists a sufficiently large network flow between the sampled nodes, with prescribed flow values over boundary edges which connect different clusters. Loosely speaking, this requires to choose the sampling set mainly in the boundary regions between different clusters in the data graph. Thus, we can leverage efficient clustering methods for identifying the cluster boundary regions in order to find sampling sets which resolve the intrinsic cluster structure of the network structure underlying a dataset.

The verification if a particular sampling set resolves a given partition requires to consider all possible sign patterns for the boundary edges, which is intractable for large graphs. An important avenue for follow-up work is the investigation if resolving sampling sets can be characterized easily using probabilistic models for the underlying network structure and sampling sets. Moreover, we plan to extend our analysis to nLasso methods using other loss functions, e.g., the squared error loss and also the logistic loss function in the context of classification problems.

Acknowledgement

The authors are grateful to Madelon Hulsebos for a careful proof-reading of an early manuscript. Moreover, the constructive comments of anonymous reviewers are appreciated sincerely. This manuscript is available as a pre-print at the following address: https://arxiv.org/abs/1704.02107. Copyright of this pre-print version rests with the authors