p -L APLACIAN B ASED G RAPH N EURAL N ETWORKS
Anonymous authors
Paper under double-blind review
A BSTRACT
Graph neural networks (GNNs) have demonstrated superior performance for semisupervised node classification on graphs, as a result of their ability to exploit node features and topological information simultaneously.
However, most GNNs implicitly assume that the labels of nodes and their neighbors in a graph are the same or consistent, which does not hold in heterophilic graphs, where the labels of linked nodes are likely to differ.
Moreover, when the topology is non-informative for label prediction, ordinary GNNs may work significantly worse than simply applying multi-layer perceptrons (MLPs) on each node.
To tackle the above problem, we propose a new p -Laplacian based GNN model, termed as p GNN, whose message passing mechanism is derived from a discrete regularization framework and could be theoretically explained as an approximation of a polynomial graph filter defined on the spectral domain of p -Laplacians.
The spectral analysis shows that the new message passing mechanism works simultaneously as low-pass and high-pass filters, thus making p GNNs are effective on both homophilic and heterophilic graphs.
Empirical studies on real-world and synthetic datasets validate our findings and demonstrate that p GNNs significantly outperform several stateof-the-art GNN architectures on heterophilic benchmarks while achieving competitive performance on homophilic benchmarks.
Moreover, p GNNs can adaptively learn aggregation weights and are robust to noisy edges.
I NTRODUCTION
In this paper, we explore the usage of graph neural networks (GNNs) for semi-supervised node classification on graphs, especially when the graphs admit strong heterophily or noisy edges.
Semisupervised learning problems on graphs are ubiquitous in a lot of real-world scenarios, such as user classification in social media (Kipf & Welling, 2017), protein classification in biology (Velickovic et al., 2018), molecular property prediction in chemistry (Duvenaud et al., 2015), and many others (Marcheggiani & Titov, 2017; Satorras & Estrach, 2018).
Recently, GNNs are becoming the de facto choice for processing graph structured data.
They can exploit the node features and the graph topology by propagating and transforming the features over the topology in each layer and thereby learn refined node representations.
A series of GNN architectures have been proposed, including graph convolutional networks (Bruna et al., 2014; Henaff et al., 2015; Defferrard et al., 2016; Kipf & Welling, 2017; Wu et al., 2019), graph attention networks (Velickovic et al., 2018; Thekumparampil et al., 2018), and other representatives (Hamilton et al., 2017; Xu et al., 2018; Pei et al., 2020).
Most of the existing GNN architectures work under the homophily assumption, i.e. the labels of nodes and their neighbors in a graph are the same or consistent, which is also commonly used in graph clustering (Bach & Jordan, 2004; von Luxburg, 2007; Liu & Han, 2013) and semi-supervised learning on graphs (Belkin et al., 2004; Hein, 2006; Nadler et al., 2009).
However, recent studies (Zhu et al., 2020; 2021; Chien et al., 2021) show that in contrast to their success on homophilic graphs, most GNNs fail to work well on heterophilic graphs, in which linked nodes are more likely to have distinct labels.
Moreover, GNNs could even fail on graphs where their topology is not helpful for label prediction.
In these cases, propagating and transforming node features over the graph topology could lead to worse performance than simply applying multi-layer perceptrons (MLPs) on each of the nodes independently.
Several recent works were proposed to deal with the heterophily issues of GNNs.
Zhu et al.
(2020) finds that heuristically combining ego-, neighbor, and higher-order embeddings improves GNN performance on heterophilic graphs.
Zhu et al.
(2021) uses a compatibility matrix to model the graph homophily or heterphily level.
Chien et al.
(2021) incorporates the generalized PageRank algorithm with graph convolutions so as to jointly optimize node feature and topological information extraction for both homophilic and heterophilic graphs.
However, the problem of GNNs on graphs with non-informative topologies (or noisy edges) remains open.
Unlike previous works, we tackle the above issues of GNNs by proposing the discrete p -Laplacian based message passing scheme, termed as p -Laplacian message passing.
It is derived from a discrete regularization framework and is theoretically verified as an approximation of a polynomial graph filter defined on the spectral domain of the p -Laplacian.
Spectral analysis of p -Laplacian message passing shows that it works simultaneously as low-pass and high-pass filters 1 and thus is applicable to both homophilic and heterophilic graphs.
Moreover, when p ̸ = 2 , our theoretical results indicate that it can adaptively learn aggregation weights in terms of the variation of node embeddings on edges (measured by the graph gradient (Amghibech, 2003;
Zhou & Sch¨olkopf, 2005; Luo et al., 2010)), and work as low-pass or both low-pass and high-pass filters on a node according to the local variation of node embeddings around the node (measured by the norm of graph gradients).
Based on p -Laplacian message passing, we propose a new GNN architecture, called p GNN, to enable GNNs to work with heterophilic graphs and graphs with non-informative topologies.
Several existing GNN architectures, including SGC (Wu et al., 2019), APPNP (Klicpera et al., 2019) and
GPRGNN (Chien et al., 2021), can be shown to be analogical to p GNN with p
= 2 .
Our empirical studies on real-world benchmark datasets (homophilic and heterophilic datasets) and synthetic datasets (cSBM (Deshpande et al., 2018)) demonstrate that p GNNs obtain the best performance on heterophilic graphs and competitive performance on homophilic graphs against state-of-the-art
GNNs.
Moreover, experimental results on graphs with different levels of noisy edges show that p
GNNs work much more robustly than GNN baselines and even as well as MLPs on graphs with completely random edges.
Additional experiments (reported in Appendix F.5) illustrate that intergrating p GNNs with existing GNN architectures (i.e. GCN (Kipf & Welling, 2017),
JKNet (Xu et al., 2018)) can significantly improve their performance on heterophilic graphs.
In conclusion, our contributions can be summarized as below:
(1) New methodologies.
We propose p -Laplacian message passing and p GNN to adapt GNNs to heterophilic graphs and graphs where the topology is non-informative for label prediction.
(2) Superior performance.
We empirically demonstrate that p GNNs is superior on heterophilic graphs and competitive on homophilic graphs against state-of-the-art GNNs.
Moreover, p GNNs work robustly on graphs with noisy edges or non-informative topologies.
(3) Theoretical justification.
We theoretically demonstrate that p -Laplacian message passing works as both low-pass and high-pass filters and the message passing iteration
is guarantee
to
converge with
proper
settings.
(4) New paradigm
of designing GNN architectures
We bridge
the gap between discrete regularization framework and
GNNs, which could further inspire researchers to develop new graph convolutions or message passing schemes using other regularization techniques with explicit assumptions on graphs.
Due to space limit, we defer the discussions on related work and future work and all proofs to the Appendix.
P RELIMINARIES AND B ACKGROUND
Notation.
Let G =
( V , E , W ) be an undirected graph, where V = { 1 , 2 , . .
, N } is the set of nodes,
E ⊆
V × V is the set of edges, W ∈ R N × N is the adjacency matrix and W i,j = W j,i , W i,j >for [ i, j ] ∈ E , W i,j = 0 , otherwise.
N i = { j } [ i,j ] ∈E denotes the set of neighbors of node i ,
D ∈ R N × N = diag( D 1 , 1 , . . .
D N,N ) denotes the diagonal degree matrix with D
i,i = (cid:80) Nj =1 W i,j , for i = 1 , . . .
, N .
f : V → R and g : E → R are functions defined on the vertices and edges of G , respectively.
F V denotes the Hilbert space of functions endowed with the inner product ⟨ f, ˜ f ⟩ F V :=
(cid:80) i ∈V f ( i )
˜
f
( i ) .
Similarly define F E .
We also denote by [ K ] = { 1 , 2 , . . .
, K } , ∀ K ∈ N
and we use ∥ x ∥ =
∥ x ∥ 2 = ( (cid:80) di =1 x 2 i ) 1 / 2 , ∀ x ∈ R d to denote the Frobenius norm of a vector.
Problem Formulation.
Given a graph G =
( V , E , W ) , each node i ∈ V has a feature vector X i, : which is the i -th row of X and a subset of nodes in G have labels from a label set L = { 1 , . . .
, L } .
The goal of semi-supervised node classification on G is to learn a mapping M : V → L and predict the labels of unlabeled nodes.
Homophily and Heterophily.
The homophily or heterophily of a graph is used to describe the relation of labels between linked nodes in the graphs.
The level of homophily of a graph can be measured by H ( G ) =
E i ∈V (cid:2)(cid:12)(cid:12) { j } j ∈N i ,y i =
y j (cid:12)(cid:12) /
|N
i | (cid:3) (Pei et al., 2020; Chien et al., 2021), where (cid:12)(cid:12) { j } j ∈N i ,y i =
y j (cid:12)(cid:12) denotes the number of neighbors of i ∈ V that share the same label as i and
H ( G ) → 1 corresponds to strong homophily while H ( G ) → 0 indicates strong heterophily.
We say that a graph is a homophilic (heterophilic) graph if it has strong homophily (heterophily).
Graph Gradient.
The graph gradient of an edge [ i, j ] , i, j ∈ V is defined to be a measurement of the variation of a function f 2 : V → R on the edge [ i, j ] .
Definition 1 (Graph Gradient) .
Given a graph G =
( V , E ) and a function f : V → R , the graph gradient is an operator ∇ : F V → F E defined by [Equation] for all [ [Equation]
j ]
[Equation]
E [Equation] For [ i, j ] / ∈ E , ( ∇ f )([ i, j ]) :
= 0 .
The graph gradient of a function f at vertex i is defined to be ∇ f ( i ) :
= (( ∇ f )([ i, 1]) , . . .
, ( ∇ f )([ i, N ])) and its Frobenius norm is given by ∥∇ f ( i ) ∥ 2 :=
( (cid:80) Nj =1 ( ∇ f )
([ i, j ])) 1 / 2 , which measures the variation of f around node i .
We measure the variation of f over the whole graph G by S p ( f ) where it is defined to be
[Equation]
Note that the definition of S p here is different with the p -Dirichlet form in Zhou & Sch ¨ olkopf (2005).
Graph Divergence.
The graph divergence is defined to be the adjoint of the graph gradient:
Definition 2 (Graph Divergence) .
Given a graph G =
( V , E ) , and functions f : V → R , g : E → R , the graph divergence is an operator div :
F E → F V which satisfies
[Equation]
The graph divergence can be computed by
[Equation]
Fig. 4 in Appendix E.1 gives a tiny example of illustration of graph gradient and graph divergence.
Graph p -Laplacian Operator.
By the definitions of graph gradient and graph divergence, we reach the definition of graph p -Laplacian operator as below.
Definition 3 (Graph p -Laplacian 3 ) .
Given a graph G =
( V , E ) and a function f : V → R , the graph p -Laplacian is an operator ∆ p : F V → F V defined by
[Equation]
where ∥ · ∥ p − 2 is element-wise, i.e. ∥∇ f ( i ) ∥
p − 2 = ( ∥ ( ∇ f )([ i, 1]) ∥
p − 2 , . .
∥ ( ∇ f )([ i, N ])
∥ p − 2 ) .
Substituting Eq.
(1) and Eq.
(4) into Eq.
(5), we obtain
[Equation]
The graph p -Laplacian is semi-definite: ⟨ f, ∆ p f ⟩ = S p ( f ) ≥
and we have
[Equation]
When p = 2 , ∆ 2 is refered as the ordinary Laplacian operator and ∆ 2 =
I − D − 1 / 2 WD − 1 /and
when p = 1 , ∆ 1 is refered as the Curvature operator and ∆ 1 f := − 12 div( ∥∇ f ∥ − 1 ∇ f ) .
Note that Laplacian ∆ 2 is a linear operator, while in general for p ̸ = 2 , p -Laplacian is nonlinear since ∆ p ( af ) ̸
= a ∆ p ( f ) for a ∈ R .
p -L APLACIAN BASED G RAPH N EURAL N ETWORKS
In this section, we derive the p -Laplacian message passing scheme from a p -Laplacian regularization framework and present p GNN, a new GNN architecture developed upon the new message passing scheme.
We theoretically characterize how p -Laplacian message passing adaptively learns aggregation weights and profits p GNN for being effective on both homophilic and heterophilic graphs.
3. p -L APLACIAN R EGULARIZATION F RAMEWORK
Given an undirected graph G =
( V , E ) and a signal function with c ( c ∈ N ) channels f : V → R c , let X =
( X ⊤ 1 , : , . . .
, X ⊤ N, : )
⊤ ∈ R N × c with X i, : ∈ R 1 × c , i ∈
[ N ] denoting the node features of G and
F = ( F ⊤ 1 , : , . . .
, F ⊤ N, : )
⊤ ∈ R N × c be a matrix whose i th row vector F i, : ∈ R 1 × c ,
i ∈
[ N ] represents the function value of f at the i -th vertex in G .
We present a p -Laplacian regularization problem whose cost function is defined to be
[Equation]
where µ ∈ (0 , ∞ ) .
The first term of the right-hand side in Eq.
(8) is a measurement of variation of the signal over the graph based on p -Laplacian.
As we will discuss later, different choices of p result in different smoothness constraint on the signals.
The second term is the constraint that the optimal signals F ∗ should not change too much from the input signal X , and µ provides a trade-off between these two constraints.
Regularization with p = 2 .
When p = 2 , the solution of Eq.
(8) satisfies ∆ 2 F ∗ + µ ( F ∗ − X ) =
and we can obtain the closed form (Zhou et al., 2003; Zhou & Sch ¨ olkopf, 2005)
[Equation]
Then, we could use the following iteration algorithm to get an approximation of Eq.
[Equation]
where k represents the iteration index, α = 1 1+ µ and β = µ 1+ µ = 1 − α .
The iteration converges to a closed-form solution as k goes to infinity (Zhou et al., 2003; Zhou & Sch ¨ olkopf, 2005).
We could relate the the result here with the personalized PageRank (PPR) (Page et al., 1999; Klicpera et al., 2019) algorithm
(proof defered to Appendix D.1):
Theorem 1 (Relation to personalized PageRank (Klicpera et al., 2019)) .
µ (∆ 2 + µ I N ) − 1 in the closed form solution of Eq.
(9) is equivalent to the personalized PageRank matrix.
Regularization with p > 1 .
For p > 1 , the solution of Eq.
(8) satisfies p ∆ p F ∗ + 2 µ ( F ∗ − X ) = 0 .
By Eq.
(6) we have that, for all i ∈ [ N ] ,
[Equation]
Based on which we can construct a similar iterative algorithm to obtain a solution (Zhou &
Sch¨olkopf, 2005):
[Equation]
M ( k ) i,j (cid:112) D i,i D j,j F ( k ) j, : + β ( k )
i,i X i, : , for all i ∈ [ N ] , (11) with M ( k ) ∈ R N × N , α ( k ) =
diag( α ( k ) 1 , 1 , . .
α ( k ) N,N ) , β ( k ) =
diag( β ( k ) 1 , 1 , . . .
, β ( k ) N,N ) updated by
M ( k ) i,j = W i,j (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W
i,j D i,i F ( k ) i, : − (cid:115) W i,j
D j,j F ( k ) j, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −, for all i, j ∈ [ N ] , (12) α ( k ) i,i = 1 (cid:14)  
N (cid:88) j =
M ( k ) i,j
D i,i + 2 µ p 
 , β ( k ) i,i = 2 µ p α
i,i , for all i ∈
[ N ] , (13)
Note that in Eq.
(12), when (cid:13)(cid:13)(cid:13)(cid:113) W i,j D i,i F ( k ) i, : − (cid:113)
W
i,j D j,j F ( k ) j, : (cid:13)(cid:13)(cid:13) = 0 , we set M ( k ) i,j = 0 .
It is easy to see that Eq.
(10) is the special cases of Eq.
(14) with p = 2 .
Remark 1 (Discussion on p = 1 ) .
For p = 1 , when f is a real-valued function ( c = 1 ), ∆ 1 f is a step function, which could make the stationary condition of the objective function Eq.
(8) become problematic.
Additionally, ∆ 1 f is not continuous at ∥ ( ∇ f )([ i, j ]) ∥
= 0 .
Therefore, p = 1 is not allowed when f is a real value function.
On the other hand, note that there is a Frobenius norm in ∆ p f .
When f is a vector-valued function ( c > 1 ), the step function in ∆ 1 f only exists on the axes.
The stationary condition will be fine if the node embeddings F are not a matrix of vectors that has only one non-zero element, which is true for many graphs.
p = 1 may work for these graphs.
Overall, we suggest to use p > 1 in practice but p = 1 may work for graphs with multiple channel signals as well.
We conduct experiments for p > 1 (e.g., p = 1 . 5 , 2 , 2 . 5 ) and p = 1 in Sec. 5.
3. p -L APLACIAN M ESSAGE P ASSING AND p GNN A RCHITECTURE
p -Laplacian Message Passing.
Rewrite Eq.
(11) in a matrix form we obtain
[Equation]
Eq.
(14) provides a new message passing mechanism, named p -Laplacian message passing.
Remark 2.
α D − 1 / 2 MD − 1 / 2 in Eq.
(14) can be regarded as the learned aggregation weights at each step for message passing.
It suggests that p -Laplacian message passing could adaptively tune the aggregation weights during the course of learning, which will be demonstrated theoretically and empirically in the sequel of this paper.
β X in Eq.
(14) can be regarded as a residual unit, which helps the model escape from the oversmoothing issue (Chien et al., 2021).
We present the following theorem to show the shrinking property of p -Laplacian message passing.
Theorem 2 (Shrinking Property of p -Laplacian Message Passing) .
Given a graph G =
( V , E , W ) with node features X , β ( k ) , F ( k ) , M ( k ) , α ( k ) are updated accordingly to Equations (11) to (13) for k = 0 , 1 , . . .
, K and F (0)
=
X . Then there exist some positive real value µ > 0 which depends on
X , G , p and p > 1 such that
[Equation]
Proof see Appendix D.2.
Thm.
2 shows that with some proper positive real value µ and p > 1 , the loss of the objective function Eq.
(8) is guaranteed to decline after taking one step p -Laplacian message passing.
Thm. 2 also demonstrates that the iteration Equations (11) to (13) is guaranteed to converge for p > 1 with some proper µ which is chosen depends on the input graph and p .
p GNN Architecture.
We design the architecture of p GNNs using p -Laplacian message passing.
Given node features X ∈ R N × c , the number of node labels L , the number of hidden units h , the maximum number of iterations K , and M , α , and β updated by Equations (12) and (13) respectively, we give the p GNN architecture as following: [Equation]
[Equation]
where Z ∈ R N ×
L is the output propbability matrix with Z i,j is the estimated probability that the label at node i ∈ [ N ] is j ∈
[ L ] given the features X and the graph G , Θ (1) ∈ R c
× h and Θ (2) ∈
R
h
×
L are the first- and the second-layer parameters of the neural network, respectively.
Remark 3 (Connection to existing GNN variants) .
The message passing scheme of p GNNs is different from that of several GNN variants (say, GCN, GAT, and GraphSage), which repeatedly stack message passing layers.
In contrast, it is similar with SGC (Wu et al., 2019), APPNP (Klicpera et al., 2019), and GPRGNN (Chien et al., 2021).
SGC is an approximation to the closed-form in
Eq.
(9) (Fu et al., 2020).
By Thm. 1, it is easy to see that APPNP, which uses PPR to propagate the node embeddings, is analogical to p GNN with p = 2 , termed as 2 .
GNN.
APPNP and 2 .
GNN work analogically and effectively on homophilic graphs.
GNN can also work effectively on heterophilic graphs by letting Θ (2) be negative.
However, APPNP fails on heterophilic graphs as its PPR weights are fixed (Chien et al., 2021).
Unlike APPNP, GPRGNN, which adaptively learn the generalized
PageRank (GPR) weights, works similarly to 2 .
GNN on both homophilic and heterophilic graphs.
However, GPRGNN needs more supervised information in order to learn optimal GPR weights.
On the contrary, p GNNs need less supervised information to obtain similar results because Θ (2) acts like a hyperplane for classification.
p GNNs could work better under weak supervised information.
Our analysis is also verified by the experimental results in Sec. 5.
We also provide an upper-bounding risk of p GNNs by Thm. 4 in Appendix C.1 to study the effect of the hyperparameter µ on the performance of p GNNs.
Thm. 4 shows that the risk of p GNNs is upperbounded by the sum of three terms: the risk of label prediction using only the original node features
X , the norm of p -Laplacian diffusion on X , and the magnitude of the noise in X .
µ controls the trade-off between these three terms.
The smaller µ , the more weights on the p -Laplacian diffusion term and the noise term and the less weights on the the other term and vice versa.
S PECTRAL V IEWS OF p -L APLACIAN M ESSAGE P ASSING
In this section, we theoretically demonstrate that p -Laplacian message passing is an approximation of a polynomial graph filter defined on the spectral domain of p -Laplacian.
We show by spectral analysis that p -Laplacian message passing works simultaneously as low-pass and high-pass filters.
p -Eigenvalues and p -Eigenvectors of the Graph p
-Laplacian.
We first introduce the definitions of p -eigenvalues and p -eigenvectors of p -Laplacian.
Let ϕ p : R → R defined as ϕ p ( u ) =
∥
u ∥ p − 2 u, for u ∈ R , u ̸ = 0 .
Note that ϕ 2 ( u ) = u .
For notational simplicity, we denote by ϕ p ( u ) =
( ϕ p ( u 1 ) , . . . ,
ϕ p ( u N )) ⊤ for u ∈ R N and Φ p ( U ) =
( ϕ p ( U : , 1 ) , . . . , ϕ p ( U : ,N )) for
U ∈ R N × N
and U : ,i ∈ R N is the i th column vector of U .
Definition 4 ( p -Eigenvector and p -Eigenvalue) .
A vector u ∈ R N is a p -eigenvector of ∆ p if it satisfies the equation
[Equation]
where λ ∈ R is a real value referred as a p -eigenvalue of ∆ p associated with the p -eigenvector u .
Definition 5 ( p -Orthogonal
Luo, Huang, Ding, and Nie, 2010)) .
Given two vectors u , v ∈ R N with u , v ̸ = 0 , we call that u and v is p -orthogonal if ϕ p ( u ) ⊤ ϕ p ( v ) =
(cid:80) Ni =1 ϕ p ( u i ) ϕ
p ( v i )
= 0 .
Luo et al.
(2010) demonstrated that the p -eigenvectors of ∆ p are p -orthogonal to each other (see
Thm. 5 in Appendix C.2 for details).
Therefore, the space spanned by the multiple p -eigenvectors of ∆ p is p -orthogonal.
Additionally, we demonstrate that the p -eigen-decomposition of ∆ p is given by: ∆ p =
Φ p ( U ) Λ Φ p ( U ) ⊤ (see Thm. 6 in Appendix C.3 for details), where U is a matrix of p -eigenvectors of ∆ p and Λ is a diagonal matrix in which the diagonal is the p -eigenvalues of ∆ p .
Graph Convolutions based on p -Laplacian.
Based on Thm.
5, the graph Fourier Transform ˆ f of any function f on the vertices of G can be defined as the expansion of f in terms of Φ( U ) where U is the matrix of p -eigenvectors of ∆ p : ˆ f = Φ p ( U )
⊤ f .
Similarly, the inverse graph Fourier transform is then given by: f = Φ p ( U ) ˆ f .
Therefore, a signal X ∈ R N × c being filtered by a spectral filter g θ can be expressed formally as: g θ ⋆ X =
Φ p ( U )
ˆ g θ ( Λ )Φ p ( U ) ⊤ X , where Λ denotes a diagonal matrix in which the diagonal corresponds to the p -eigenvalues { λ l } l =0 ,...,N − 1 of ∆ p and ˆ g θ ( Λ ) denotes a diagonal matrix in which the diagonal corresponds to spectral filter coefficients.
Let ˆ g
θ be a polynomial filter defined as ˆ g θ =
(cid:80) K − 1 k =
θ k λ kl , where the parameter θ =
[ θ 0 , . .
, θ K − 1 ] ⊤ ∈ R K is a vector of polynomial coefficients.
By the p -eigen-decomposition of p -Laplacian, we have
[Equation]
Theorem 3.
The K -step p
-Laplacian message passing is a K -order polynomial approximation to the graph filter given by Eq.
Proof see Appendix D.3.
Thm. 3 indicates that p -Laplacian message passing mechanism is implicitly a polynomial spectral filter defined on the spectral domain of p -Laplacian.
Spectral Analysis of p -Laplacian Message Passing.
Here, we analyze the spectral propecties of p -Laplacian message passing.
We can approximately view p -Laplacian message pasing as a filter of a linear combination of K spectral filters g ( Λ ) (0) , g ( Λ ) (1) , . .
, g ( Λ ) ( K − 1) with each spectral filter defined to be g ( Λ ) ( k ) :
= ( α D − 1 / 2 MD − 1 / 2 ) k where M i,j = W i,j ∥ (cid:113) W
i,j D i,i F i, : − (cid:113) W
i,j D j,j F j, : ∥ p −for i, j ∈ [ N ] and F is the matrix of node embeddings.
We can study the properties of p -Laplacian message passing by studying the spectral properties of α D − 1 / 2 MD − 1 / 2 as given below.
Proposition 1.
Given a connected graph G =
( V , E , W ) with node embeddings F and the pLaplacian ∆ p with its p -eigenvectors { u ( l ) } l =0 , 1 ,...,N − 1 and the p -eigenvalues { λ l } l =0 , 1 ,...,N − 1 .
Let g p ( λ i − 1 ) :
= α
i,i (cid:80) j D − 1 / 2 i,i
M i,j D − 1 / 2 j,j for i ∈
[ N ] be the filters defined on the spectral domain of ∆ p , where M i,j = W i,j ∥∇ f
([ i, j ]) ∥
p − 2 , ( ∇ f )([ i, j ]) is the graph gradient of the edge between node i and j and ∥∇ f ( i ) ∥ is the norm of graph gradient at i .
N i denotes the number of edges connected to i , N min = min { N j } j ∈ [ N ] , and k = arg max j
( {| u ( l ) j | / (cid:112) D l,l } j ∈ [ N ]; l =0 ,.
..,N − 1 ) , then 1.
When p = 2 , g p ( λ i − 1 ) works as both low-pass and high-pass filters.
When p > 2 , if ∥∇ f ( i ) ∥ ≤ 2 ( p − 1) /
( p − 2) , g p ( λ i − 1 ) works as both low-pass and high-pass filters on node i and g p ( λ i − 1 ) works as low-pass filters on i when ∥∇ f ( i ) ∥ ≥ 2 ( p − 1) /
( p − 2) .
When 1 ≤ p < 2 , if 0 ≤ ∥∇ f ( i ) ∥ ≤ 2(2 √ N k ) 1 / ( p − 2) , g p ( λ i − 1 ) works as low-pass filters on node i and g p ( λ i − 1 ) works as both low-pass and high-pass filters on i when ∥∇ f ( i ) ∥ ≥ 2 (cid:0) 2 √ N k (cid:1) 1 / ( p − 2) .
Specifically, when p = 1 , N k can be replaced by N min .
Proof see Appendix D.7.
Proposition 1 shows that when p ̸ = 2 , p -Laplacian message passing adaptively works as low-pass or both low-pass and high-pass filters on node i in terms of the degree of local node embedding variation around i , i.e. the norm of the graph gradient ∥∇ f ( i ) ∥ at node i .
When p = 2 , p -Laplacian message passing works as both low-pass and high-pass filters on node i regardless of the value of ∥∇ f ( i ) ∥ .
When p > 2 , p -Laplacian message passing works as lowpass filters on node i for large ∥∇ f ( i ) ∥ and works as both low-pass and high-pass filters for small ∥∇ f
( i ) ∥ .
Therefore, p GNNs with p > 2 can work very effectively on graphs with strong homophily.
When 1 ≤ p < 2 , p -Laplacian message passing works as low-pass filters for small ∥∇ f
( i ) ∥ and works as both low-pass and high-pass filters for large ∥∇ f
( i ) ∥ .
Thus, p GNNs with 1 ≤ p < 2 can work effectively on graphs with low homophily, i.e. heterophilic graphs.
The results here confirms our analysis of the aggregation weights of p -Laplacian message passing presented in Thm. 2.
E MPIRICAL S TUDIES
In this section, we empirically study the effectiveness of p GNNs for semi-supervised node classification using and real-world benchmark and synthetic datasets with heterophily and strong homophily.
The experimental results are also used to validate our theoretical findings presented previously.
Datasets and Experimental Setup.
We use seven homophilic benchmark datasets: citation graphs
Cora, CiteSeer, PubMed (Sen et al., 2008), Amazon co-purchase graphs Computers, Photo, coauthor graphs CS, Physics (Shchur et al., 2018), and six heterophilic benchmark datasets: Wikipedia graphs
Chameleon, Squirrel (Rozemberczki et al., 2021), the Actor co-occurrence graph, webpage graphs
Wisconsin, Texas, Cornell (Pei et al., 2020).
The node classification tasks are conducted in the transductive setting.
Following Chien et al.
(2021), we use the sparse splitting ( 2 . 5% / 2 .
5% / 95% ) and the dense splitting ( 60% / 20% / 20% ) to randomly split the homophilic and heterophilic graphs into training / validation / testing sets, respectively.
Dataset statistics and their levels of homophily are presented in Appendix E.
Baselines.
We compare p GNN with seven models, including MLP, GCN (Kipf & Welling, 2017),
SGC (Wu et al., 2019), GAT (Velickovic et al., 2018), JKNet (Xu et al., 2018), APPNP (Klicpera et al., 2019), GPRGNN (Chien et al., 2021).
We use the Pytorch Geometric library (Fey & Lenssen, 2019) to implement all baselines except GPRGNN.
For GPRGNN, we use the code released by the authors 4 .
The details of hyperparameter settings are deferred to Appendix E.3.
[Table]
Superior Performance on Real-World Heterophilic Datasets.
The results on homophilic benchmark datasets are deferred to Appendix F.1, which show that p GNNs obtains competitive performance against state-of-the-art GNNs on homophilic datasets.
Table 1 summarizes the results on heterophilic benchmark datasets.
Table 1 shows that p GNNs significantly dominate the baselines and 1 .
GNN obtains the best performance on all heterophilic graphs except the Texas dataset.
For
Texas, 2 .
GNN is the best.
We also observe that MLP works very well and significantly outperforms most GNN baselines, which indicates that the graph topology is not informative for label prediction on these heterophilic graphs.
Therefore, propagating and transforming node features over the graph topology could lead to worse performance than MLP.
Unlike ordinary GNNs, p GNNs can adaptively learn aggregation weights and ignore edges that are not informative for label prediction and thus could work better.
It confirms our theoretical findings presented in previous sections.
Note that GAT can also learn aggregation weights, i.e. the attention weights.
However, the aggregation weights learned by GAT are significantly distinct from that of p GNNs, as we will show following.
[Figure]
Interpretability of the Learned Aggregation Weights of p GNNs.
We showcase the interpretability of the learned aggregation weights α i,i D − 1 / 2 i,i
M i,j D − 1 / 2 j,
[Equation]
of p GNNs
by studying its entropy distribution, along with the attention weights of GAT on real-world datasets.
Denote { A i,j } j ∈N i as the aggregation weights of node i and its neighbors.
For GAT, { A i,j } j ∈N i are referred as the attention weights (in the first layer) and for p GNNs are α i,i D − 1 / 2 i,
i
M i,j D − 1 / 2 j,j .
For any node i , { A i,j } j ∈N i forms a discrete probability distribution over all its neighbors with the entropy given by
H ( { A i,j } j ∈N i ) = − (cid:80) j ∈N i A
i,j log( A i,j ) .
Low entropy means high degree of concentration and vice versa.
An entropy of zero means all aggregation weights or attentions are on one source node.
The uniform distribution has the highest entropy of log ( D i,i ) .
Fig.
1 reports the results on Computers, Wisconsin and we defer more results on other datasets to Appendix F.2 due to space limit.
Fig.
1 shows that the aggregation weight entropy distributions of GAT and p GNNs on Computers (homophily) are both similar to the uniform case.
It indicates the original graph topology of Computers is very helpful for label prediction and therefore GNNs could work very well on Computers.
However, for Wisconsin (heterophily), the entropy distribution of p GNNs is significantly different from that of GAT and the uniform case.
Most entropy of p GNNs is around zero, which means that most aggregation weights are on one source node.
It states that the original graph topology of
Wisconsin is not helpful for label prediction, which explains why MLP works well on Wisconsin.
[Figure]
On the contrary, the entropy distribution of GAT is similar to the uniform case and therefore GAT works similarly to GCN and is significantly worse than p GNNs on Wisconsin.
Similar results can be observed on the experiments on more datasets in Appendix F.
[Figure]
We exam the performance of p GNNs on heterophilic graphs whose topology is informative for label prediction using synthetic graphs generated by cSBM (Deshpande et al., 2018) with ϕ ∈ {− 1 , − 0 .
We use the same settings of cSBM used in Chien et al.
Due to the space limit, we refer the readers to Chien et al.
(2021) for more details of cSBM dataset.
Fig.
2 reports the results on cSBM using sparse splitting (for results on cSBM with dense splitting see Appendix F.3).
Fig.
2 shows that when ϕ ≤ − 0 .
5 (heterophilic graphs), 2 .
GNN obtains the best performance and p GNNs and
GPRGNN significantly dominate the others.
It validates the effectiveness of p GNNs on heterophilic graphs.
Moreover, 2 .
GNN works better than GPRGNN and it again confirms that 2 .
GNN is more superior under weak supervision (2.5% training rate), as stated in Remark 3.
Note that 1 .
GNN and 1 .
GNN are not better than 2 .
GNN, the reason could be the iteration algorithms Eq.
(11) with p = 1 , 1 . 5 are not as stable as the one with p = 2 .
When the graph topology is almost non-informative for label prediction ( ϕ = − 0 .
25 , 0 ), The performance of p GNNs is close to MLP and they outperform the other baselines.
Again, it validates that p GNNs can erase noninformative edges and work as well as MLP and confirms the statements in Thm. 4.
When the graph is homophilic ( ϕ ≥ 0 . 25 ), 1 . 5 GNN is the best on weak homophilic graphs ( ϕ = 0 .
25 , 0 . 5 ) and p
GNNs work competitively with all GNN baselines on strong homophilic graphs ( ϕ ≥ 0 . 75 ).
Results on Datasets with Noisy Edges.
We conduct experiments to evaluate the performance of p
GNNs on graphs with noisy edges by randomly adding edges to the graphs and randomly remove the same number of original edges.
We define the random edge rate as r :=
# random edges # all edges .
The experiments are conducted on 4 homophilic datasets (Computers, Photo, CS, Physics) and 2 heterophilic datasets (Wisconsin, Texas) with r = 0 .
Fig.
3 reports the results on Computers, Wisconsin and we defer more results to Appendix F.4.
Fig.
3 shows that p GNNs significantly outperform all baselines.
Specifically, 1 . 5 GNN obtains the best performance on Computers, and 1 . 5 GNN and 2 .
GNN even work as well as MLP on Computers with completely random edges ( r = 1 ).
For
Wisconsin, 1 .
GNN is the best, and 1 .
GNN and 1 .
GNN significantly dominate the others.
We also observed that APPNP and GPRGNN, whose architectures are analogical to 2 .
GNN, also work better than other GNNs.
Nevertheless, they are significantly outperformed by p GNNs overall.
Similar results can be observed in the experiments conducted on more datasets as presented in Appendix F.4.
C ONCLUSION .
We have addressed the problem of generalizing GNNs to heterophilic graphs and graphs with noisy edges.
To this end, we derived a novel p -Laplacian message passing scheme from a discrete regularization framework and proposed a new p GNN architecture.
We theoretically demonstrate our method works as low-pass and high-pass filters and thereby applicable to both homophilic and heterophilic graphs.
We empirically validate our theoretical results and show the advantages of our methods on heterophilic graphs and graphs with non-informative topologies.
R EPRODUCIBILITY S TATEMENT
In order to ensure reproducibility we have made the efforts in the following respects: (1) Provide a sampled code as the supplementary material; (2) Provide self-contained proofs of the main claims in Appendices C and D; (3) Provide more details on experimental configurations in Appendix E and experimental results in Appendix F.
All the datasets are publicly available as described in the main text.
We will fully release our training and evaluation code, as well as our train/validation/test splits
Appendix C ONTENTS
A Related Work
B Discussions and Future Work
C
Additional Theorems
C.1 Theorem 4 (Upper-Bounding Risk of p GNN) . . . . . . . . . . . . . . . . . . . .
C.2 Theorem 5 ( p -Orthogonal Theorem (Luo et al., 2010)) . . . . . . . . . . . . . . .
C.3 Theorem 6 ( p -Eigen-Decomposition of ∆ p ) . . . . . . . . . . . . . . . . . . . . .
C.4 Theorem 7 (Bounds of p -Eigenvalues) . . . . . . . . . . . . . . . . . . . . . . . .
D Proof of Theorems
D.1 Proof of Theorem 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
D.2 Proof of Theorem 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
D.3 Proof of Theorem 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
D.4 Proof of Theorem 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
D.
Proof of Theorem 6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
D.6 Proof of Theorem 7 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
D.7 Proof of Proposition 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
E Dataset Statistics and Hyperparameters
E.1 Illustration of Graph Gradient and Graph Divergence . . . . . . . . . . . . . . . .
E.2 Dataset Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
E.3 Hyperparameter Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
F Additional Experiments
F.
Experimental Results on Homophilic Benchmark Datasets . . . . . . . . . . . . .
F.
Experimental Results of Aggregation Weight Entropy Distribution . . . . . . . . .
F.
Experimental Results on cSBM . . . . . . . . . . . . . . . . . . . . . . . . . . . .
F.
Experimental Results on Graphs with Noisy Edges . . . . . . . . . . . . . . . . .
F.
Experimental Results of Intergrating p GNNs with GCN and JKNet . . . . . . . . .
F.
Experimental Results of p GNNs on PPI Dataset for Inductive Learning . . . . . . .
F.
Experimental Results of p GNNs on OGBN arXiv Dataset . . . . . . . . . . . . . .
F.
Running Time of p GNNs . .
F.
Results
Units . . . . . . . . .
F.10 Training Curves p = 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
F.11 Visualization Results of
Embeddings . . . . . . . . . . . . . . . . . . . . .
A R ELATED W ORK
Graph Neural Networks.
Graph neural networks (GNNs) are a variant of neural networks for graph-structured data, which can propagate and transform the node features over the graph topology and exploit the information in the graphs.
Graph convolutional networks (GCNs) are one type of GNNs whose graph convolution mechanisms or the message passing schemes were mainly inspired by the field of graph signal processing.
Bruna et al.
(2014) defined a nonparametric graph filter using the Fourier coefficients.
Defferrard et al.
(2016) introduced Chebyshev polynomial to avoid computational expensive eigen-decomposition of Laplacian and obtain localized spectral filters.
GCN (Kipf & Welling, 2017) used the first-order approximation and reparameterized trick to simplify the spectral filters and obtain the layer-wise graph convolution.
SGC (Wu et al., 2019) further simplify GCN by removing non-linear transition functions between each layer.
Chen et al.
(2018) propose importance sampling to design an efficient variant of GCN.
Xu et al.
(2018) explored a jumping knowledge architecture that flexibly leverages different neighborhood ranges for each node to enable better structure-aware representation.
Atwood & Towsley (2016); Liao et al.
Abu-El-Haija et al.
(2019) exploited multi-scale information by diffusing multi-hop neighbor information over the graph topology.
Wang & Leskovec (2020) used label propagation to improve GCNs.
Klicpera et al.
(2019) incorporated personalized PageRank with GCNs.
Liu et al.
(2021) introduced a l 1 norm-based graph smoothing term to enhance the local smoothnesss adaptivity of
GNNs.
Hamilton et al.
; Zeng et al.
(2020) proposed sampling and aggregation frameworks to extent GCNs to inductive learning settings.
Another variant of GNNs is graph attention networks (Velickovic et al., 2018; Thekumparampil et al., 2018; Abu-El-Haija et al., 2018), which use attention mechanisms to adaptively learn aggregation weights based on the nodes features.
There are many other works on GNNs (Pei et al., 2020)
(Ying et al., 2018; Xinyi & Chen, 2019; Velickovic et al., 2019; Zeng et al., 2020), we refer to Zhou et al.
; Battaglia et al.
; Wu et al.
(2021) for a comprehensive review.
Most GNN models implicitly assume that the labels of nodes and their neighbors should be the same or consistent, while it does not hold for heterophilic graphs.
Zhu et al.
(2020) investigated the issues of GNNs on heterophilic graphs and proposed to separately learn the embeddings of ego-node and its neighborhood.
Zhu et al.
(2021) proposed a framework to model the heterophily or homophily levels of graphs.
Chien et al.
(2021) incorporated generalized PageRank with graph convolution to adapt GNNs to heterophilic graphs.
There are also some works on the interpretability of GNNs proposed recently.
Li et al.
; Ying et al.
Fu et al.
(2020) showed that spectral graph convolutions work as conducting Laplacian smoothing on the graph signals and Wu et al.
; NT & Maehara (2019) demonstrated that
GCN, SGC work as low-pass filters.
Gama et al.
(2020) studied the stability properties of GNNs.
Xu et al.
; Oono & Suzuki (2020); Loukas (2020) studied the expressiveness of GNNs.
Verma & Zhang (2019); Garg et al.
(2020) work on the generalization and representation power of GNNs.
Graph based Semi-supervised Learning.
Graph-based semi-supervised learning works under the assumption that the labels of a node and its neighbors shall be the same or consistent.
Many methods have been proposed in the last decade, such as Smola & Kondor (2003);
Zhou et al.
Belkin et al.
(2004) use Laplacian regularization techniques to force the labels of linked nodes to be the same or consistent.
Zhou & Sch¨olkopf (2005) introduce discrete regularization techniques to impose different regularizations on the node features based on p -Laplacian.
Lable propagation (Zhu et al., 2003) recursively propagates the labels of labeled nodes over the graph topology and use the convergence results to make predictions.
To mention but a few, we refer to Zhou & Sch¨olkopf (2005); van Engelen & Hoos (2020) for a more comprehensive review.
B D ISCUSSIONS AND F UTURE W ORK
In this section, we discuss the future work of p GNNs.
Our theoretical results and experimental results could lead to several potential extensions of p GNNs.
New Paradigm of Designing GNN Architectures.
We bridge the gap between discrete regularization framework, graph-based semi-supervised learning, and GNNs, which provides a new paradigm of designing new GNN architectures.
Following the new paradigm, researchers could introduce more regularization techniques, e.g., Laplacian regularization (Smola & Kondor, 2003; Belkin et al., 2004), manifold regularization (Sindhwani et al., 2005; Belkin et al., 2005; Niyogi, 2013), highorder regularization (Zhou & Belkin, 2011), Bayesian regularization (Liu et al., 2014), entropy regularization (Grandvalet & Bengio, 2004), and consider more explicit assumptions on graphs, e.g. the homophily assumption, the low-density region assumption (i.e. the decision boundary is likely to lie in a low data density region), manifold assumption (i.e. the high dimensional data lies on a low-dimensional manifold), to develop new graph convolutions or message passing schemes for graphs with specific properties and generalize GNNs to a much broader range of graphs.
Moreover, the paradigm also enables us to explicitly study the behaviors of the designed graph convolutions or message passing schemes from the theory of regularization (Belkin & Niyogi, 2008; Niyogi, 2013;
Slepcev & Thorpe, 2017).
Applications of p GNNs to learn on graphs with noisy topologies.
The empirical results (as shown in Fig.
3 and Tables 6 and 7) on graphs with noisy edges show that p GNNs are very robust to noisy edges, which suggests the applications of p -Laplacian message passing and p GNNs on the graph learning scenarios where the graph topology could potentially be seriously intervened.
Integrating with existing GNN architectures.
As shown in Table 9, the experimental results on heterophilic benchmark datasets illustrate that integrating GCN, JKNet with p GNNs can significantly improve their performance on heterophilic graphs.
It shows that p GNN could be used as a plug-and-play component to be integrated into existing GNN architectures and improve their performance on real-world applications.
Inductive learning for p GNNs.
p GNNs are shown to be very effective for inductive learning on
PPI datasets as reported in Table 10.
p GNNs even outperforms GAT on PPI, while using much fewer parameters than GAT.
It suggests the promising extensions of p GNNs to inductive learning on graphs.
C A DDITIONAL T HEOREMS C. T HEOREM 4 (U PPER -B OUNDING R ISK OF p GNN)
Theorem 4 (Upper-bounding risks of p GNNs) .
Given a graph G =
( V , E , W ) with N nodes, let
X ∈ R N × c be the node features and y ∈ R N be the node labels and M ( k ) , α ( k ) , β k , F k are updated accordingly by Equations (12) to (14) for k = 0 , 1 , . .
, K − 1 and F (0) =
X , K ∈ N .
Assume that G is d -regular and the ground-truth node features X ∗ =
X + ϵ , where ϵ ∈ R N × c represents the noise in the node features and there exists a L -Lipschitz function σ : R N × c
→ R N such that σ ( X ∗ )
=
y . let ˜ y
( k +1) =
α ( k ) D − 1 / 2 M ( k ) D − 1 / 2 σ ( F ( k ) )
+ β ( k ) σ (cid:0) F (0) (cid:1) , we have
N
N (cid:88) i =
[Equation]
(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
[Equation]
K
[Equation]
i
[Equation]
: +
K
[Equation]
(cid:88) k
=
K − 2 (cid:89) l = k  
N (cid:88) j =
M ( l ) i,j d   ∆ ( k ) p
X i, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
[Equation]
Proof see Appendix D.4.
Thm.
4 shows that the risk of p GNNs is upper-bounded by the sum of three terms: The first term of the r.h.s in the above inequation represents the risk of label prediction using only the original node features X , the second term is the norm of p -Laplacian diffusion on the node features X , and the third term is the magnitude of the noise in the node features.
α i,i and β i,i control the trade-off between these three terms and they are related to the hyperparameter µ in Eq.
The smaller µ , the smaller β i,i and larger α
i,i , thus the more important of the p -Laplacian diffusion term but also the more effect from the noise.
Therefore, for graphs whose topological information is not helpful for label prediction, we could impose more weights on the first term by using a large µ so that p GNNs work more like MLPs which simply learn on node features.
While for graphs whose topological information is helpful for label prediction, we could impose more weights on the second term by using a small µ so that p GNNs can benefit from p -Laplacian smoothing on node features.
In practice, to choose a proper value of µ one may first simply apply MLPs on the node features to have a glance at the helpfulness of the node features.
If MLPs work very well, there is not much space for the graph’s topological information to further improve the prediction performance and we may choose a large µ .
Otherwise, there could be a large chance for the graph’s topological information to further improve the performance and we should choose a small µ .
C.
[Equation]
Theorem 5 ( p -Orthogonal Theorem (Luo et al., 2010)) .
If u ( l ) and u ( r ) are two eigenvectors of p -Laplacian ∆ p associated with two different non-zero eigenvalues λ l and λ r , W is symmetric and p ≥ 1 , then u ( l ) and u ( r ) are p -orthogonal up to the second order Taylor expansion.
Thm. 5 implies that ϕ p ( u ) ( l ) ⊤ ϕ p ( u ( r ) )
≈ 0 , for all l, r = 0 , . .
, N − 1 and λ l ̸ = λ r .
Therefore, the space spanned by the multiple eigenvectors of the graph p -Laplacian is p -orthogonal.
C. T HEOREM 6 ( p -E IGEN -D ECOMPOSITION OF ∆ p )
Theorem 6 ( p -Eigen-Decomposition of ∆ p ) .
Given the p -eigenvalues { λ l ∈ R } l =0 , 1 ,...,N − 1 , and the p -eigenvectors { u ( l ) ∈ R N } l =0 , 1 ,...,N − 1 of p -Laplacian ∆ p and ∥ u ( l ) ∥ p
= ( (cid:80) Ni =
( u ( l ) i ) p ) 1 /p
= 1 , let U be a matrix of p -eigenvectors with U
=
( u (0) , u (1) , . . .
, u ( N − 1) ) and Λ be a diagonal matrix with Λ
= diag( λ 0 , λ 1 , . .
, λ N − 1 ) , then the p -eigen-decomposition of p -Laplacian ∆ p is given by
[Equation]
When p = 2 , it reduces to the standard eigen-decomposition of the Laplacian matrix.
Proof see Appendix D.5.
C. T HEOREM 7 (B OUNDS OF p -E IGENVALUES )
Theorem 7 (Bounds of p -Eigenvalues) .
Given a graph G =
( V , E , W ) , if G is connected and λ is a p -eigenvalue associated with the p -eigenvector u of ∆ p , let N i denotes the number of edges connected to node i , N min = min { N i } i =1 , 2 ,...,N , and k = arg max
( {| u i | / (cid:112) D i,i } i =1 , 2 ,...,N ) , then
for p ≥ 2 , 0 ≤ λ ≤ 2 p − 1 ; 2.
for 1 < p < 2 , 0 ≤ λ ≤ 2 p − 1 √ N k ; 3.
for p = 1 , 0 ≤ λ ≤ √ N min .
Proof see Appendix D.6.
D P ROOF OF T HEOREMS D. P ROOF OF T HEOREM
Proof.
Let i be the one-hot indicator vector whose i -th element is one and the other elements are zero.
Then, we can obtain the personalized PageRank on node i , denoted as π PPR ( i ) , by using the recurrent equation (Klicpera et al., 2019):
[Equation]
where k is the iteration step, 0 < α < 1 and β =
(1 − α ) represents the restart probability.
Without loss of generality, suppose π (0) PPR ( i ) =
i .
Then we have,
[Equation]
Since 0 < α < 1 and the eigenvalues of D − 1 / 2 WD − 1 / 2 in [ − 1 , 1] , we have
[Equation]
and we also have
[Equation]
Therefore,
[Equation]
where we let α = 1 1+ µ and β = µ 1+ µ , µ > 0 .
Then the fully personalized PageRank matrix can be obtained by substituting i with I N : Π PPR = µ (∆ 2 + µ I N ) − 1 .
D. P ROOF OF T HEOREM
Proof.
By the definition of L p ( f ) in Eq.
(8), we have for some positive real value µ, µ >
[Equation]
N (cid:88) i =
N (cid:88)
j = (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115)
[Equation]
j
D
[Equation]
: − (cid:115)
[Equation]
j
D
[Equation]
: (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p + µ
N (cid:88) i = ∥
F i, : − X i, : ∥ 2 .
and by Eq.
[Equation]
(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115)
[Equation]
i
[Equation]
D i
[Equation]
i
[Equation]
( k ) i, : − (cid:115) W
i,j
D j,j F ( k ) j, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −
Then, we have ∂ L p ( F ( k ) )
∂ F ( k ) i, : = p
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13) (cid:13)(cid:13)(cid:13)(cid:13) (cid:115) W i,j
D i,i F ( k ) i, : − (cid:115) W i,j
D j,j F ( k ) j, : (cid:13) (cid:13)(cid:13)(cid:13)(cid:13) p − 2 (cid:32)(cid:115) W i,j
D i,i F ( k )
[Equation]
which indicates that
[Equation]
For all i, j ∈ [ N ] , v ∈ R 1 × c , denote by
[Equation]
i, : ,
M ′ ( k ) i,j :=
W i,j (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W
i,j D i,i ( F ( k ) i, : + v ) − (cid:115) W i,j
D j,j F ( k ) j, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −, α ′ ( k ) i,i :=(cid:44) 
N (cid:88) j =
M ′ ( k ) i,j
D i,i + 2 µ p 
 , β ′ ( k ) i,i :=
2 µ p α ′ ( k ) i,i
F ′ ( k +1) i, : := α ′ ( k ) i,i
N (cid:88) j =
M ′ ( k ) i,j (cid:112) D i,i D j,j F ( k ) j, : + β ′
i,i X i, : .
Then
[Equation]
) i,i ∥ v ∥ + p (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) (cid:13)  
N (cid:88)
j
[Equation]
[Equation] [Equation]
i
[Equation]
j
[Equation]
N (cid:88)
[Equation]
=
[Equation]
( k )
i
[Equation]
j
[Equation]
i
[Equation]
i 
[Equation] [Equation]
−
N (cid:88)
[Equation]
=
[Equation] [Equation]
j (cid:112)
[Equation]
j F ( k )
j
[Equation]
: −
[Equation]
i
[Equation] [Equation]
N (cid:88) j =
M ( k ) i,j (cid:112) D i
[Equation]
i D j
[Equation]
j F ( k ) j, : + 2 µ p X
i, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) (cid:13) =
[Equation]
[Equation]
(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
N (cid:88) j =
M ′ ( k ) i
[Equation]
j
[Equation]
i
[Equation]
i
[Equation]
N (cid:88) j =
[Equation]
( k )
i
[Equation]
j
[Equation]
i,i 
[Equation]
F ( k ) i
[Equation]
: −
N (cid:88) j =
M ′ ( k ) i,j (cid:112)
[Equation]
i
[Equation]
i
[Equation]
j
[Equation]
j F ( k ) j, :
[Equation]
N (cid:88) j =
M ( k ) i,j (cid:112) D i,
i
[Equation]
j
[Equation]
j F ( k ) j, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13) =
[Equation]
(cid:88)
[Equation]
N (cid:88)
[Equation]
Therefore
[Equation]
there exists some real positive value µ ∈ o ( p, v , X , G ) > 0 such that
[Equation]
Let γ =
( γ 1 , . . . , γ N ) ⊤ ∈ R N and η ∈ R N × c .
By Taylor’s theorem
[Equation]
we have:
[Equation]
i ⟨ η i, : , ∂ L p ( F ( k ) i, : ) ⟩ + γ i (cid:90) 0 ⟨ ∂ L p ( F ( k ) i, : + ϵγ i η i, : ) − ∂ L p ( F ( k ) i, : ) , η i, : ⟩ d ϵ ≤ L p ( F ( k ) i, : ) + γ i ⟨ η i, : , ∂ L p ( F ( k ) i, : ) ⟩ + γ i (cid:90)
∥ ∂ L p ( F ( k ) i, : + ϵγ i η i, : ) − ∂ L p ( F ( k ) i, : ) ∥∥
η i, : ∥ d ϵ ≤ L p ( F ( k ) i, : )
+ γ i ⟨ η i, : ,
∂ L p ( F ( k ) i, : ) ⟩ + p 2 α ( k ) i,i γ 2 i ∥
η
i, : ∥
Let η = −∇L
p ( F ( k ) )
and choose some positive real value µ which depends on X , G , p and p > 1 , i.e. µ ∈ o ( p, X , G ) .
By Eq.
(19), we have for all i ∈
[ N ] ,
[Equation]
∥
[Equation]
Then for all i ∈ [ N ] , when 0 ≤ γ i ≤ 2 α ( k ) i,i p , we have S p ( F ( k ) i, : − γ i ∂ S p ( F ( k ) i, : )) ≤ S p ( F ( k ) i, : ) and γ i = α ( k ) i,i p minimizes S p ( F ( k ) i, : − γ i ∂ S p ( F ( k ) i, : )) .
Therefore,
[Equation]
D. P ROOF OF T
HEOREM
Proof.
Without loss of generality, suppose F (0) =
X .
Denote ˜ M
( k ) =
D − 1 / 2 M ( k ) D − 1 / 2 , by
Eq.
(14), we have for K ≥ 2 ,
F ( K ) =
α ( K − 1) D − 1 / 2 M
( K − 1)
D − 1 / 2 F
( K − 1)
+ β ( K − 1)
X
=
α ( K − 1) ˜
M K − 1 F ( K − 1)
+ β ( K − 1)
X
=
α ( K − 1) ˜
M K − 1 (cid:16) α ( K − 2) ˜
M K − 2 F ( K − 2)
+ β ( K − 2)
X (cid:17) + β ( K − 1)
X
=
α ( K − 1) α
( K − 2)
˜ M K − 1 ˜ M K − 2 F ( K − 2)
+ α
( K − 1)
˜ M K − 1 β ( K − 2)
X
+ β ( K − 1) X
= (cid:32) K − 1 (cid:89) k
= α ( k ) (cid:33) (cid:32) K − 1 (cid:89) k
= ˜
M ( k ) (cid:33)
F (0) + K − 1 (cid:88) k
= (cid:32) K − 1 (cid:89) l
=
K − k α ( l ) ˜ M
( l ) (cid:33) β
( K − 1 − k ) X
+ β ( K − 1) X
= (cid:32) K − 1 (cid:89) k
= α ( k ) (cid:33) (cid:32) K − 1 (cid:89) k
= ˜
M ( k ) (cid:33)
X +
K − 1 (cid:88) k
= (cid:32) K − 1 (cid:89) l
=
K − k α ( l ) ˜ M
( l ) (cid:33) β
( K − 1 − k ) X
+ β ( K − 1) X .
Recall Equations (12) and (13), we have ˜ M ( k ) i,j =
[Equation]
i
[Equation]
D
[Equation]
(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W i,j
D i,i F ( k ) i, : − (cid:115) W i,j
D j,j F ( k ) j, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −, for all i, j = 1 , 2 , . . . , N, and
[Equation]
, for all i = 1 , 2 , . . .
, N.
Note that the eigenvalues of ˜ M are not infinity and 0 < α
i,i < 1 for all i = 1 , . . .
, N .
Then we have
[Equation]
and lim
[Equation]
= 0 .
Therefore,
[Equation]
By Equations (6) and (12)
[Equation]
we have
[Equation]
By Eq.
[Equation]
we have
[Equation]
Equations (22) and (23) show that
[Equation]
which indicates
[Equation]
Eq.
(25) shows that α ( k )
˜
M
( k ) is linear w.r.t ∆ p and therefore can be expressed by a linear combination in terms of ∆ p :
[Equation]
where θ ′ =
diag( θ ′ 0 , θ ′ 1 , . .
, θ ′ N − 1 ) are the parameters.
Therefore, we have
[Equation]
where θ ′′ ( k ) =
diag( θ ′′ ( k ), θ ′′ ( k ), . .
, θ ′′ ( k ) N ) defined as θ ′′ (0)
=
β ( K − 1) and θ ′′ ( k ) i =
β ( K − 1 − k )
i,i
K − 1 (cid:89) l = K − k θ ′ ( l ) i , for k = 1 , 2 , . .
, K − 1 .
Let θ =
( θ 0 , θ 1 , . .
, θ K − 1 ) defined as θ k = (cid:80) Ni =1 θ ′′ ( k ) i for all k = 0 , 1 , . . .
, K − 1 , then
[Equation]
Therefore complete the proof
[Equation]
D. P ROOF OF T HEOREM
Proof.
The first-order Taylor expansion with Peano’s form of remainder for σ at X ∗ i, : is given by:
[Equation]
Note that in general the output non-linear layer σ ( · ) is simple.
Here we assume that it can be well approximated by the first-order Taylor expansion and we can ignore the Peano’s form of remainder.
For all i = 1 , . . .
, N , D i,i =
D j,j = d , we have α i,i (cid:80) Nj =1 D − 1 / 2 i,
i
[Equation]
Then
[Equation]
i, : )
[Table] [Table]
(1 − β ( K − 1) i,i ) (cid:13)(cid:13)(cid:13)(cid:13) ∂σ ( X ∗ i, : )
∂ X (cid:13)(cid:13)(cid:13)(cid:13) ∥
ϵ i, : ∥ ≤ β ( K − 1) i,i | y i − σ ( X i, : ) | + α
( K − 1)
i,i
L (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) (cid:13)
N (cid:88) j =
M ( K − 1) i,j d (cid:16) F ( K − 1) j, : − X i, : (cid:17)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13) (cid:13) + (1 − β ( K − 1) i,i )
L ∥
ϵ i, : ∥ =
β ( K − 1) i,i | y
i − σ ( X i, : ) | + α
( K − 1)
i,i
L (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
N (cid:88) j =
M ( K − 1) i,j d (cid:16) F ( K − 1) j, : − F ( K − 1) i, : + F ( K − 1) i, : − X i, : (cid:17)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
+ (1 − β
( K − 1)
i,i )
L ∥
ϵ i, : ∥ =
β ( K − 1) i,i | y
i − σ ( X i, : ) | + α
( K − 1)
i,i
L (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
∆ p F ( K − 1) i, : +
N (cid:88) j =
M ( K − 1) i,j d ·
N (cid:88) j =
M ( K − 2) i,j d (cid:16) F ( K − 2) j, : − X i, : (cid:17)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
+ (1 − β
( K − 1)
i,i )
L ∥
ϵ i, : ∥ =
β ( K − 1) i,i | y
i − σ ( X i, : ) | + α
( K − 1)
i,i
L (cid:13) (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) ∆ ( K − 1) p
F ( K − 1) i, : +
K − 2 (cid:88) k
=
K − 2 (cid:89) l = k  
N (cid:88) j =
M ( l ) i,j d   ∆ ( k ) p
X i, : (cid:13) (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
+ (1 − β
( K − 1)
i,i )
L ∥
ϵ i, : ∥ =
β ( K − 1) i,i | y
i − σ ( X i, : ) | + α
( K − 1)
i,i
L (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13) ∆ ( K − 1) p
F ( K − 1) i, : +
K − 2 (cid:88) k
=
K − 2 (cid:89) l = k  
N (cid:88) j =
M ( l ) i,j d   ∆ ( k ) p
X i, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
[Equation]
Therefore,
N N
[Equation]
(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
[Equation]
p
F
[Equation]
K
[Equation]
1) i, :
[Equation]
K − 2 (cid:88) k
=
K − 2 (cid:89) l = k  
N (cid:88) j =
M ( l ) i,j d   ∆ ( k ) p
X i, : (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13)
[Equation]
D. P ROOF OF T HEOREM
Proof.
Note that ϕ
[Equation]
i | p
[Equation]
∥ u ∥ pp = 1 , then we have
[Equation]
Therefore, ∆ p = Φ p ( U ) Λ Φ p ( U ) ⊤ .
When p = 2 , by Φ 2 ( U ) =
U , we get ∆ 2 =
Φ 2 ( U ) Λ Φ 2 ( U )
⊤ =
UΛU ⊤ .
D. P ROOF OF T HEOREM
Proof.
By the definition of graph p -Laplacian, we have for all i
= 1 , 2 , . .
, N , (∆ p u ) i =
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115)
W i,j
D i,i u i − (cid:115) W i,j
D j,j u j (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p − 2 (cid:32)(cid:115) W i,j
D i,i u i − (cid:115) W i,j
D j,j u j (cid:33) = λϕ p ( u i ) .
Then, for all i = 1 , 2 , . .
. , N , λ
= ϕ p ( u i )
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115)
W i,j
D i,i u i − (cid:115) W i,j
D j,j u j (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p − 2 (cid:32)(cid:115) W i,j
D i,i u i − (cid:115) W i,j
D j,j u j (cid:33) =
∥ u i ∥ p − 2 u
i
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13) (cid:13)(cid:13)(cid:13)(cid:13) (cid:115) W i,j
D i,i u i − (cid:115) W i,j
D j,j u j (cid:13) (cid:13)(cid:13)(cid:13)(cid:13) p − 2 (cid:32)(cid:115) W i,j
D i,i u i − (cid:115) W i,j
D j,j u j (cid:33) =
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13)(cid:13)(cid:13)(cid:113) W
i,j D i,i u i − (cid:113)
W
i,j D j,j u j (cid:13)(cid:13)(cid:13) p −∥ u i ∥ p −(cid:32)(cid:115)
W i,j
D i,i − (cid:115) W i,j
D j,j u j u i (cid:33) =
N (cid:88)
j = (cid:32) W i,j
D i,i −
W i,j (cid:112)
D i,i D j,j u j u i (cid:33)   (cid:13)(cid:13)(cid:13)(cid:113) W
i,j D i,i u i − (cid:113)
W
i,j D j,j u j (cid:13)(cid:13)(cid:13) ∥ u i ∥ 
 p −=
N (cid:88)
j = (cid:32) W i,j
D i,i −
W i,j (cid:112)
D i,i D j,j u j u i (cid:33) (cid:13) (cid:13) (cid:13)(cid:13)(cid:13)(cid:13) (cid:113) W
i,j D i,i u i − (cid:113)
W
i,j D j,j u j u i (cid:13) (cid:13) (cid:13)(cid:13)(cid:13)(cid:13) p −=
N (cid:88)
j = (cid:32) W i,j
D i,i −
W i,j (cid:112)
D i,i D j,j u j u i (cid:33) (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W
i,j
D i,i − (cid:115) W i,j
D j,j u j u
i
cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −
Let l = arg max {∥ u i ∥} i =1 , 2 ,...,N , the above equation holds for all i = 1 , 2 , . .
, N , then
[Equation]
− (cid:115)
[Equation]
j
D
[Equation]
u
[Equation]
u l (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −≥
N (cid:88) j =
[Equation]
D
[Equation]
(cid:112)
D
[Equation]
D
[Equation]
(cid:12)(cid:12)(cid:12)(cid:12)
[Equation]
(cid:12)(cid:12)(cid:12)(cid:12)(cid:33) (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W l,j
D l
[Equation]
l − (cid:115) W l,j
D j,j u j
[Equation]
l (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −≥
[Equation]
D
[Equation]
W
[Equation]
D
[Equation]
(cid:33) (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W l,j
D l
[Equation]
l − (cid:115) W l,j
D j,j u j u l (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −≥ 0 .
When p = 1 , λ =
N (cid:88)
j = (cid:32) W i,j
D i,i −
W i,j (cid:112)
D i,i D j,j u j u i (cid:33) (cid:13)(cid:13)(cid:13) (cid:13) (cid:13) (cid:115) W i,j
D i,i − (cid:115) W i,j
D j,j u j u i (cid:13)(cid:13)(cid:13) (cid:13) (cid:13) −≤
N (cid:88)
j = (cid:115) W i,j
D i,i ≤ (cid:118)(cid:117)(cid:117)(cid:116) N i N (cid:88) j =
W i,j
D i,i = (cid:112) N i , where the last inequality holds by using the Cauchy-Schwarz inequality.
The above inequality holds for all i = 1 , 2 , . .
, N , therefore, λ ≤
N (cid:88)
j = (cid:115) W i,j
D i,i ≤ (cid:112) N min .
When p > 1 , we have for i = 1 , 2 , . .
. , N , λ =
N (cid:88)
j = (cid:32) W i,j
D i,i −
W i,j (cid:112)
D i,i D j,j u j u i (cid:33) (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115) W
i,j
D i,i − (cid:115) W i,j
D j,j u j u
i
cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −≤
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115)
W i,j
D i,i − (cid:115) W i,j
D j,j u j u
i
cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −≤
N (cid:88)
j = (cid:115) W i,j
D i,i (cid:13)(cid:13)(cid:13)(cid:13)(cid:13)(cid:115)
W i,j
D i,i + (cid:115) W i,j
D j,j (cid:12)(cid:12)(cid:12)(cid:12) u
j u i (cid:12)(cid:12)(cid:12)(cid:12)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −=
N (cid:88)
j = (cid:32)(cid:115) W
i,j
D i,i (cid:33) p (cid:13)(cid:13)(cid:13)(cid:13)(cid:13) 1 + (cid:115) D i,
i
D j,j (cid:12)(cid:12)(cid:12)(cid:12) u
j u i (cid:12)(cid:12)(cid:12)(cid:12)(cid:13)(cid:13)(cid:13)(cid:13)(cid:13) p −.
Without loss of generality, let k = arg max( {| u i | / (cid:
) D i,i } i =1 , 2 ,...,N ) .
Because the above inequality holds for all i = 1 , 2 , . .
, N , then we have λ ≤
N (cid:88)
j = (cid:32)(cid:115) W k,j
D k,k (cid:33) p (cid:32) 1 + (cid:115)
D k,k
D j,j (cid:12)(cid:12)(cid:12)(cid:12) u j u k (cid:12)(cid:12)(cid:12)(cid:12)(cid:33) p −≤ 2 p − 1 N (cid:88)
j
= (cid:32)(cid:115) W k,j
D k,k (cid:33) p .
For p ≥ 2 , λ ≤ 2 p − 1 N (cid:88)
j = (cid:32)(cid:115) W k,j
D k,k (cid:33) p ≤ 2 p − 1 N (cid:88)
j = (cid:32)(cid:115) W k,j
D k,k (cid:33)= 2 p − 1 .
For 1 < p < 2 ,
[Equation]
D. P ROOF OF P ROPOSITION
Proof.
We proof Proposition 1 based on the bounds of p -eigenvalues as demonstrated in Thm. 7.
By Eq.
(6) and Eq.
(12), we have
[Equation]
By Eq.
[Equation]
we have
[Equation]
Equations (27) and (28) show that
[Equation]
which indicates
[Equation]
For i = 1 , 2 , . . .
, N , let ˜ α :=
(˜ α 1 , . . . , ˜ α N ) ,
˜ α
i := 1 / (cid:80) Nj =1 M i,j D i,i , then α i,i
[Equation]
Recall the Eq.
(12) that
[Equation]
When p = 2 , for all i = 1 , . .
. , N , ˜ α
i = 1 and 0 ≤ λ i − 1 ≤ 2 , g 2 ( λ i − 1 ) works as both low-pass and high-pass filters.
When
p
>
, by
Thm.
we have for all i = 1 , . . .
, N , 0 ≤ λ i − 1 ≤ 2 p − 1 .
If 0 ≤ ˜ α i ≤ 2 1 − p , then 0 ≤ 1 − ˜ α
i
λ i ≤ 1 , which indicates that g p ( λ i − 1 ) works as a low-pass filter; If
˜ α
i > 2 1 − p , then g p ( λ i − 1 ) works as both low-pass and high-pass filters.
Since
N (cid:88) j =
[Equation]
i
[Equation]
j
[Equation]
which indicates that ˜ α i ≥ ∥∇ f ( i ) ∥
2 − p . 0 ≤ ˜ α i ≤ 2 1 − p directly implies that 0 ≤ ∥∇ f ( i ) ∥ 2 − p ≤ 2 1 − p , i.e. ∥∇ f
( i ) ∥ ≥ 2 ( p − 1) / ( p − 2) and when ∥∇ f ( i ) ∥ 2 − p ≥ 2 1 − p , i.e.
∥∇ f ( i ) ∥ ≤ 2 ( p − 1) / ( p − 2) , ˜ α i ≥ 2 1 − p always holds.
Therefore, if ∥∇ f ( i ) ∥ ≤ 2 ( p − 1) /
( p − 2) , g p ( λ i − 1 ) works as both low-pass and high-pass filters on node i ; If g p ( λ i − 1 ) works as a low-pass filter, ∥∇ f
( i ) ∥ ≥ 2 ( p − 1) /
( p − 2) .
When 1 ≤ p < 2 , by Thm.
7 we have for all i = 1 , . . .
, N , 0 ≤ λ i − 1 ≤ 2 p − 1 √ N k .
If 0 ≤ ˜ α i ≤ 2 1 − p / √ N k , 0 ≤ 1 − ˜ α
i λ i ≤ 1 , which indicates that g p ( λ i − 1 ) work as low-pass filters; If ˜ α i ≥ 2 1 − p / √ N k , g p ( λ i − 1 ) work as both low-pass and high-pass filters.
By
N (cid:88) j =
W i,j
[Equation]
i,
[Equation]
=
[Equation]
we have
[Equation]
˜ α
i ≥ 2 1 − p / √ N k directly implies that ∥∇ f
( i ) ∥
2 − p ≥ 2 1 − p / √ N k , i.e. ∥∇ f
( i ) ∥ ≥ 2(2 √ N k ) 1 / ( p − 2) and when 0 ≤ ∥∇ f ( i ) ∥ 2 − p ≤ 2 1 − p / √ N k , i.e. 0 ≤ ∥∇ f
( i )
∥ ≤ 2(2 √ N k ) 1 / ( p − 2) , 0 ≤ ˜ α i ≤ 2 1 − p / √ N k always holds.
Therefore, if 0 ≤ ∥∇ f ( i ) ∥ ≤ 2(2 √ N k ) 1 / ( p − 2) , g p ( λ i − 1 ) work as low-pass filters; If g p ( λ i − 1 ) work as both low-pass and high-pass filters, ∥∇ f ( i ) ∥ ≥ 2 (cid:0) 2 √ N k (cid:1) 1 / ( p − 2) .
Specifically, when p = 1 , by Thm.
7 we have for all i = 1 , . . .
, N , 0 ≤ λ i − 1 ≤ 2 p − 1 √ N min .
Following the same derivation above we attain if 0 ≤ ∥∇ f ( i ) ∥ ≤ 2(2 √ N min ) 1 / ( p − 2) , g p ( λ i − 1 ) work as low-pass filters; If g p ( λ i − 1 ) work as both low-pass and high-pass filters, ∥∇ f ( i ) ∥ ≥ 2 (cid:0) 2 √ N min (cid:1) 1 / ( p − 2) .
E D ATASET S TATISTICS AND H YPERPARAMETERS E. I LLUSTRATION OF G RAPH G RADIENT AND G RAPH D IVERGENCE
𝑓𝛻𝑓 [1,2] =
𝑊 1,2 𝐷 2,2 𝑓 2 − 𝑊 1,2 𝐷 1,1 𝑓𝑓𝛻𝑓 [1,3] = 𝑊 1,3 𝐷 3,3 𝑓 3 − 𝑊 1,3 𝐷 1,1 𝑓𝑓𝛻𝑓 [2,1] =
𝑊 1,2 𝐷 1,1 𝑓 1 − 𝑊 1,2 𝐷 2,2 𝑓𝛻𝑓 [3,1] =
𝑊 1,3 𝐷 1,1 𝑓 1 − 𝑊 1,3 𝐷 3,3 𝑓(a) Graph gradient.
𝑔([1,3]) 𝑔([2,1]) div𝑔
=
𝑊 1,2 𝐷 1,
𝑔 1,2 −𝑔 2,1 + 𝑊 1,3 𝐷 1,1 𝑔 1,3 −𝑔 3,    𝑔([3,1]) 𝑔([1,2])
(b) Graph divergence.
E. D ATASET S TATISTICS
Table 2 summarizes the dataset statistics and the levels of homophily H ( G ) of all benchmark datasets.
Note that the homophily scores here is different with the scores reported by Chien et al.
There is a bug in their code when computing the homophily scores (doing division with torch integers) which caused their homophily scores to be smaller.
[Table]
E. H YPERPARAMETER S ETTINGS
We set the number of layers as 2, the maximum number of epochs as 1000, the number for early stopping as 200, the weight decay as 0 or 0.0005 for all models.
The other hyperparameters for each model are listed as below:
GNN, 1 . 5 GNN, 2 .
GNN, 2 . 5 GNN: – Number of hidden units:– Learning rate: { 0 . 001 , 0 . 01 , 0 . 05 } – Dropout rate: { 0 , 0 . 5 } – µ : { 0 . 01 , 0 . 1 , 0 . 2 , 1 , 10 } – K : 4 , 6 ,
• MLP: – Number of hidden units:– Learning rate: { 0 . 001 , 0 . 01 } – Dropout rate: { 0 , 0 . 5 } • GCN: – Number of hidden units:– Learning rate: { 0 . 001 , 0 . 01 } – Dropout rate: { 0 , 0 . 5 } • SGC: – Number of hidden units:– Learning rate: { 0 . 2 , 0 . 01 } – Dropout rate: { 0 , 0 . 5 } – K :• GAT: – Number of hidden units:– Number of attention heads:– Learning rate: { 0 . 001 , 0 . 005 } – Dropout rate: { 0 , 0 . 6 } • JKNet: – Number of hidden units:– Learning rate: { 0 . 001 , 0 . 01 } – Dropout rate: { 0 , 0 . 5 } – K :– α : { 0 . 1 , 0 . 5 , 0 .
7 , 1 } – The number of GCN based layers:– The layer aggregation: LSTM with 16 channels and 4 layers • APPNP: – Number of hidden units:– Learning rate: { 0 . 001 , 0 . 01 } – Dropout rate: { 0 , 0 . 5 } – K :– α : { 0 . 1 , 0 . 5 , 0 . 7 , 1 } • GPRGNN: – Number of hidden units:– Learning rate: { 0 . 001 , 0 . 01 , 0 . 05 } – Dropout rate: { 0 , 0 . 5 } – K :– α : { 0 , 0 . 1 , 0 . 2 , 0 . 5 , 0 . 7 , 0 . 9 , 1 } – dprate: { 0 , 0 . 5 , 0 . 7 }
F A DDITIONAL E XPERIMENTS F. E XPERIMENTAL R ESULTS ON H OMOPHILIC B ENCHMARK D ATASETS
Competitive Performance on Real-World Homophilic Datasets.
Table 3 summarizes the averaged accuracy (the micro-F1 score) and standard deviation of semi-supervised node classification on homophilic benchmark datasets.
Table 3 shows that the performance of p GNN is very close to
APPNP, JKNet, GCN on Cora, CiteSeer, PubMed datasets and slightly outperforms all baselines on
Computers, Photo, CS, Physics datasets.
Moreover, we observe that p GNNs outperform GPRGNN on all homophilic datasets, which confirms that p GNNs work better under weak supervised information (2.5% training rate) as discussed in Remark 3.
We also see that all GNN models work significantly better than MLP on all homophilic datasets.
It illustrates that the graph topological information is helpful for the label prediction tasks.
Notably, 1 .
GNN is slightly worse than the other p
GNNs with larger p , which suggests to use p ≈ 2 for homophilic graphs.
Overall, the results of
3 indicates that p GNNs obtain competitive performance against all baselines on
[Table]
F. E XPERIMENTAL R ESULTS OF A GGREGATION W EIGHT E NTROPY D ISTRIBUTION
Here we present the visualization results of the learned aggregation weight entropy distribution of p GNNs and GAT on all benchmark datasets.
Fig.
5 and Fig.
6 show the results obtained on homophilic and heterophilic benchmark datasets, respectively.
We observe from Fig.
5 that the aggregation weight entropy distributions learned by p GNNs and
GAT on homophilic benchmark datasets are similar to the uniform cases, which indicates that aggregating and transforming node features over the original graph topology is very helpful for label prediction.
It explains why p GNNs and GNN baselines obtained similar performance on homophilic benchmark datasets and all GNN models significantly outperform MLP.
Contradict to the results on homophilic graphs shown in Fig.
5, Fig.
6 shows that the aggregation weight entropy distributions of p GNNs on heterophilic benchmark datasets are very different from that of GAT and the uniform cases.
We observe from Fig.
6 that the entropy of most of the aggregation weights learned by p GNNs are around zero, which means that most aggregation weights are on one source node.
It indicates that the graph topological information in these heterophilic benchmark graphs is not helpful for label prediction.
Therefore, propagating and transforming node features over the graph topology could lead to worse performance than MLPs, which validates the results in
Table 3 that the performance of MLP is significantly better most GNN baselines on all heterophilic graphs and closed to p GNNs.
[Figure]
[Figure]
F. E XPERIMENTAL R ESULTS ON C SBM
In this section we present the experimental results on cSBM using sparse splitting and dense splitting, respectively.
We used the same settings in Chien et al.
(2021) in which the number of nodes n = 5000 , the number of features f = 2000 , ϵ = 3 .
25 for all experiments.
Table 4 reports the results on cSBM with sparse splitting setting, which also are presented in Fig.
2 and discussed in
Sec. 5.
Table 5 reports the results on cSBM with dense splitting settings.
[Table] [Table]
Table 5 shows that p GNNs obtain the best performance on weak homophilic graphs ( ϕ = 0 , 0 . 25 ) while competitive performance against GPRGNN on strong heterophilic graphs ( ϕ = − 0 .
75 , − 1 ) and competitive performance with state-of-the-art GNNs on strong homophilic graphs ( ϕ = 0 .
We also observe that GPRGNN is slightly better than p GNNs on weak heterophilic graphs ( ϕ = − 0 .
25 , − 0 .
5 ), which suggests that GPRGNN could work very well using strong supervised information ( 60% training rate and 20% validation rate).
However, as shown in Table 4, p GNNs work better than GPRGNN under weak supervised information ( 2 . 5% training rate and 2 .
5% ) on all heterophilic graphs.
The result is reasonable, as discussed in Remark 3 in Sec. 3.2, GPRGNN can adaptively learn the generalized PageRank (GPR) weights and it works similarly to 2 .
GNN on both homophilic and heterophilic graphs.
However, it needs more supervised information in order to learn optimal GPR weights.
On the contrary, p GNNs need less supervised information to obtain similar results because Θ (2) acts like a hyperplane for classification.
Therefore, p GNNs can work better under weak supervised information.
F. E XPERIMENTAL R ESULTS ON G RAPHS WITH N OISY E DGES
Here we present more experimental results on graph with noisy edges.
Table 6 reports the results on homophilic graphs (Computers, Photo, CS, Physics) and Table 7 reports the results on heterophilic graphs (Wisconsin Texas).
We observe from Tables 6 and 7 that p GNNs dominate all baselines.
Moreover, p GNNs even slightly better than MLP when the graph topologies are completely random, i.e. the noisy edge rate r = 1 .
We also observe that the performance of GCN, SGC, JKNet on homophilic graphs dramatically degrades as the noisy edge rate r increases while they do not change a lot for the cases on heterophilic graphs.
It is reasonable since the original graph topological information is very helpful for label prediction on these homophilic graphs.
Adding noisy edges and remove the same number of original edges could significantly degrade the performance of ordinary
GNNs.
On the other hand, since we find that the original graph topological information in Wisconsin and Texas is not helpful for label prediction.
Therefore, adding noisy edges and removing original edges on these heterophilic graphs would not affect too much their performance.
[Table] [Table]
F. E XPERIMENTAL R ESULTS OF I NTERGRATING p GNN S WITH GCN AND JKN ET
Here we further conduct experiments to study whether p GNNs can be intergrated into existing GNN architectures and improve their performance on heterophilic graphs.
We use two popular GNN architectures: GCN (Kipf & Welling, 2017) and JKNet (Xu et al., 2018)
[Table]
To incorporate p GNNs with GCN
[Table]
we use the p GNN layers as the first layer of the combined models, termed as p GNN + GCN
[Table]
and GCN layer as the second layer.
Specifically, we use the aggregation weights α D − 1 / 2 MD − 1 / 2 learned by the p GNN in the first layer as the input edge weights of GCN layer in the second layer.
To combine p GNN with JKNet, we use the p GNN layer as the GNN layers in the JKNet framework, termed as p GNN + JKNet.
Tables 8 and 9 report the experimental results on homophilic and heterophilic benchmark datasets, respectively.
[Table] [Table]
We observe from Table 8 that intergrating p GNNs with GCN and JKNet does not improve their performance on homophilic graphs.
The performance of GCN slightly degrade after incorporating p
GNNs.
The performance of JKNet also slightly degrade on Cora, CiteSeer, and PubMed but is improved on Computers, Photo, CS, Physics.
It is reasonable since GCN and JKNet can predict well on these homophilic benchmark datasets based on their original graph topology.
However, for heterophilic benchmark datasets, Table 9 shows that there are significant improvements over GCN, and JKNet after intergrating with p GNNs.
Moreover, p GNNs + JKNet obtain advanced performance on all heterophilic benchmark datasets and even better than p GNNs on Squirrel.
The results of Table 9 demonstrate that intergrating p GNNs with GCN and JKNet can sigificantly improve their performance on heterophilic graphs.
F.
E XPERIMENTAL R ESULTS OF p GNN S ON PPI D ATASET FOR I NDUCTIVE L EARNING
Additionally, we conduct comparison experiments of p GNNs against GAT on PPI dataset (Zitnik &
Leskovec, 2017) using the inductive learning settings as in Velickovic et al.
(2018) (20 graphs for training, 2 graphs for validation, 2 graphs for testing).
We use three layers of GAT architecture with 256 hidden units, use 1 attention head for GAT (1 head) and 4 attention heads for GAT (4 heads).
We use three p GNN layers and a MLP layer as the first layer for p GNNs, set µ = 0 .
01 , K =
[Table]
and use 256 hidden units for p GNN-256 and 512 hidden units for p GNN-512.
The experimental results are reported in Table 10.
[Table]
From Appendix F.6 we observe that the results of 2 . 5 GNN on PPI slightly better than GAT withattention heads and other p GNNs are very close to it.
Moreover, all results of p GNNs significantly outperform GAT with one attention head.
The results of p GNNs on PPI is impressive.
p GNNs have much less parameters than GAT with 4 attention heads while obtain very completitive performance on PPI.
When we use more hidden units, 512 hidden units, p GNNs-512 significantly outperform
GAT, while p GNNs-512 still have less parameters.
It illustrates the superior potential of applying p
GNNs to inducting learning on graphs.
F. E XPERIMENTAL R ESULTS OF p GNN S ON OGBN AR X IV D ATASET
[Table]
Here we present the experimental of p GNNs on OGBN arXiv dataset (Hu et al., 2020).
We use the official data split setting of OGBN arXiv.
We use three layers p GNN architecture and 256 hidden units with µ = 0 .
5 , K = 2 .
We also combine p GNNs with correct and smooth model (C&S) (Huang et al., 2021) and introduce residual units.
The results of MLP, GCN, JKNet, DeepGCN (Li et al.,
, GCN with residual units, C&S model are extracted from the leaderboard for OGBN arXiv.
dataset 5 .
Table 11 summaries the results of p GNNs against the baselines.
We observe from Table 11 that p GNNs outperform MLP, GCN, JKNet
[Table]
and
DeepGCN.
The performance of p GNNs can be further improved by combining it with C&S model and residual units and
GNN + residual (6 layers) + C&S obtains the best performance against the baselines.
F. R UNNING T IME OF p GNN S
Tables 12 and 13 report the averaged running time of p GNNs and baselines on homophilic and heterophilic benchmark datasets, respectively.
[Table] [Table]
F. E XPERIMENTAL R ESULTS ON B ENCHMARK D ATASETS FOR 64 H IDDEN U NITS
[Table] [Table]
[Figure]
[Table]
[Figure] [Figure]
[Figure] [Figure]
[Figure] [Figure]
[Figure] [Figure]
[Figure] [Figure]