ADAGCN: ADABOOSTING GRAPH CONVOLUTIONAL NETWORKS INTO DEEP MODELS

Abstract

The design of deep graph models still remains to be investigated and the crucial part is how to explore and exploit the knowledge from different hops of neighbors in an efficient way. In this paper, we propose a novel RNN-like deep graph neural network architecture by incorporating AdaBoost into the computation of network; and the proposed graph convolutional network called AdaGCN (Adaboosting Graph Convolutional Network) has the ability to efficiently extract knowledge from high-order neighbors of current nodes and then integrates knowledge from different hops of neighbors into the network in an Adaboost way. Different from other graph neural networks that directly stack many graph convolution layers, AdaGCN shares the same base neural network architecture among all "layers" and is recursively optimized, which is similar to an RNN. Besides, We also theoretically established the connection between AdaGCN and existing graph convolutional methods, presenting the benefits of our proposal. Finally, extensive experiments demonstrate the consistent state-of-the-art prediction performance on graphs across different label rates and the computational advantage of our approach AdaGCN 1 .

1. INTRODUCTION

Recently, research related to learning on graph structural data has gained considerable attention in machine learning community. Graph neural networks (Gori et al., 2005; Hamilton et al., 2017; Veličković et al., 2018) , particularly graph convolutional networks (Kipf & Welling, 2017; Defferrard et al., 2016; Bruna et al., 2014) have demonstrated their remarkable ability on node classification (Kipf & Welling, 2017) , link prediction (Zhu et al., 2016) and clustering tasks (Fortunato, 2010) . Despite their enormous success, almost all of these models have shallow model architectures with only two or three layers. The shallow design of GCN appears counterintuitive as deep versions of these models, in principle, have access to more information, but perform worse. Oversmoothing (Li et al., 2018) has been proposed to explain why deep GCN fails, showing that by repeatedly applying Laplacian smoothing, GCN may mix the node features from different clusters and makes them indistinguishable. This also indicates that by stacking too many graph convolutional layers, the embedding of each node in GCN is inclined to converge to certain value (Li et al., 2018) , making it harder for classification. These shallow model architectures restricted by oversmoothing issue limit their ability to extract the knowledge from high-order neighbors, i.e., features from remote hops of neighbors for current nodes. Therefore, it is crucial to design deep graph models such that high-order information can be aggregated in an effective way for better predictions. There are some works (Xu et al., 2018b; Liao et al., 2019; Klicpera et al., 2018; Li et al., 2019; Liu et al., 2020) that tried to address this issue partially, and the discussion can refer to Appendix A.1. By contrast, we argue that a key direction of constructing deep graph models lies in the efficient exploration and effective combination of information from different orders of neighbors. Due to the apparent sequential relationship between different orders of neighbors, it is a natural choice to incorporate boosting algorithm into the design of deep graph models. As an important realization of boosting theory, AdaBoost (Freund et al., 1999) is extremely easy to implement and keeps competitive in terms of both practical performance and computational cost (Hastie et al., 2009) . Moreover, boosting theory has been used to analyze the success of ResNets in computer vision (Huang et al., 2018) and AdaGAN (Tolstikhin et al., 2017) has already successfully incorporated boosting algorithm into the training of GAN (Goodfellow et al., 2014) . In this work, we focus on incorporating AdaBoost into the design of deep graph convolutional networks in a non-trivial way. Firstly, in pursuit of the introduction of AdaBoost framework, we refine the type of graph convolutions and thus obtain a novel RNN-like GCN architecture called AdaGCN. Our approach can efficiently extract knowledge from different orders of neighbors and then combine these information in an AdaBoost manner with iterative updating of the node weights. Also, we compare our AdaGCN with existing methods from the perspective of both architectural difference and feature representation power to show the benefits of our method. Finally, we conduct extensive experiments to demonstrate the consistent state-of-the-art performance of our approach across different label rates and computational advantage over other alternatives.

2. OUR APPROACH: ADAGCN

2.1 ESTABLISHMENT OF ADAGCN Consider an undirected graph G = (V, E) with N nodes v i ∈ V, edges (v i , v j ) ∈ E. A ∈ R N ×N is the adjacency matrix with corresponding degree matrix D ii = j A ij . In the vanilla GCN model (Kipf & Welling, 2017) for semi-supervised node classification, the graph embedding of nodes with two convolutional layers is formulated as: Z = Â ReLU( ÂXW (0) )W (1) (1) where Z ∈ R N ×K is the final embedding matrix (output logits) of nodes before softmax and K is the number of classes. X ∈ R N ×C denotes the feature matrix where C is the input dimension. Â = D-1 2 Ã D-1 2 where Ã = A + I and D is the degree matrix of Ã. In addition, W (0) ∈ R C×H is the input-to-hidden weight matrix for a hidden layer with H feature maps and W (1) ∈ R H×K is the hidden-to-output weight matrix. Our key motivation of constructing deep graph models is to efficiently explore information of highorder neighbors and then combine these messages from different orders of neighbors in an AdaBoost way. Nevertheless, if we naively extract information from high-order neighbors based on GCN, we are faced with stacking l layers' parameter matrix W (i) , i = 0, ..., l -1, which is definitely costly in computation. Besides, Multi-Scale Deep Graph Convolutional Networks (Luan et al., 2019) also theoretically demonstrated that the output can only contain the stationary information of graph structure and loses all the local information in nodes for being smoothed if we simply deepen GCN. Intuitively, the desirable representation of node features does not necessarily need too many nonlinear transformation f applied on them. This is simply due to the fact that the feature of each node is normally one-dimensional sparse vector rather than multi-dimensional data structures, e.g., images, that intuitively need deep convolution network to extract high-level representation for vision tasks. This insight has been empirically demonstrated in many recent works (Wu et al., 2019; Klicpera et al., 2018; Xu et al., 2018a) , showing that a two-layer fully-connected neural networks is a better choice in the implementation. Similarly, our AdaGCN also follows this direction by choosing an appropriate f in each layer rather than directly deepen GCN layers. Thus, we propose to remove ReLU to avoid the expensive joint optimization of multiple parameter matrices. Similarly, Simplified Graph Convolution (SGC) (Wu et al., 2019 ) also adopted this prac- θ sharing the same neural network architecture f θ . w l and θ l denote node weights and parameters computed after the l-th base classifier, respectively. tice, arguing that nonlinearity between GCN layers is not crucial and the majority of the benefits arises from local weighting of neighboring features. Then the simplified graph convolution is: Z = Âl XW (0) W (1) • • • W (l-1) = Âl X W , where we collapse 1) as W and Âl denotes Â to the l-th power. In particular, one crucial impact of ReLU in GCN is to accelerate the convergence of matrix multiplication since the ReLU is a contraction mapping intuitively. Thus, the removal of ReLU operation could also alleviate the oversmoothing issue, i.e. slowering the convergence of node embedding to indistinguishable ones (Li et al., 2018) . Additionally, without ReLU this simplified graph convolution is also able to avoid the aforementioned joint optimization over multiple parameter matrices, resulting in computational benefits. Nevertheless, we find that this type of stacked linear transformation from graph convolution has insufficient power in representing information of high-order neighbors, which is revealed in our experiment described in Appendix A.2. Therefore, we propose to utilize an appropriate nonlinear function f θ , e.g., a two-layer fully-connected neural network, to replace the linear transformation W in Eq. 2 and enhance the representation ability of each base classifier in AdaGCN as follows: W (0) W (1) • • • W (l- Z (l) = f θ ( Âl X), where Z (l) represents the final embedding matrix (output logits before Softmax) after the l-th base classifier in AdaGCN. This formulation also implies that the l-th base classifier in AdaGCN is extracting knowledge from features of current nodes and their l-th hop of neighbors. Due to the fact that the function of l-th base classifier in AdaGCN is similar to that of the l-th layer in other traditional GCN-based methods that directly stack many graph convolutional layers, we regard the whole part of l-th base classifier as the l-th layers in AdaGCN. As for the realization of Multi-class AdaBoost, we apply SAMME (Stagewise Additive Modeling using a Multi-class Exponential loss function) algorithm (Hastie et al., 2009) , a natural and clean multi-class extension of the two-class AdaBoost adaptively combining weak classifiers. As illustrated in Figure 1 , we apply base classifier f (l) θ to extract knowledge from current node feature and l-th hop of neighbors by minimizing current weighted loss. Then we directly compute the weighted error rate err (l) and corresponding weight α (l) of current base classifier f (l) θ as follows: err (l) = n i=1 w i I c i = f (l) θ (x i ) / n i=1 w i α (l) = log 1 -err (l) err (l) + log(K -1), where w i denotes the weight of i-th node and c i represents the category of current i-th node. To attain a positive α (l) , we only need (1 -err (l) ) > 1/K, i.e., the accuracy of each weak classifier should be better than random guess (Hastie et al., 2009) . This can be met easily to guarantee the weights to be updated in the right direction. Then we adjust nodes' weights by increasing weights on incorrectly classified ones: w i ← w i • exp α (l) • I c i = f (l) θ (x i ) , i = 1, . . . , n After re-normalizing the weights, we then compute Âl+1 X = Â • ( Âl X) to sequentially extract knowledge from l+1-th hop of neighbors in the following base classifier f (l+1) θ . One crucial point of AdaGCN is that different from traditional AdaBoost, we only define one f θ , e.g. a two-layer fully connected neural network, which in practice is recursively optimized in each base classifier just similar to a recurrent neural network. This also indicates that the parameters from last base classifier are leveraged as the initialization of next base classifier, which coincides with our intuition that l + 1-th hop of neighbors are directly connected from l-th hop of neighbors. The efficacy of this kind of layer-wise training has been similarly verified in (Belilovsky et al., 2018) recently. Further, we combine the predictions from different orders of neighbors in an Adaboost way to obtain the final prediction C(A, X): C(A, X) = arg max k L l=0 α (l) f (l) θ ( Âl X) Finally, we obtain the concise form of AdaGCN in the following: Âl X = Â • ( Âl-1 X) Z (l) = f (l) θ ( Âl X) Z = AdaBoost(Z (l) ) (7) Note that f θ is non-linear, rather than linear in SGC (Wu et al., 2019) , to guarantee the representation power. As shown in Figure 1 , the architecture of AdaGCN is a variant of RNN with synchronous sequence input and output. Although the same classifier architecture is adopted for f (l) θ , their parameters are different, which is different from vanilla RNN. We provide a detailed description of the our algorithm in Section 3. Architectural Difference. As illustrated in Figure 1 and 2, there is an apparent difference among the architectures of GCN (Kipf & Welling, 2017) , SGC (Wu et al., 2019) , Jumping Knowledge (JK) (Xu et al., 2018b) and AdaGCN. Compared with these existing graph convolutional approaches that sequentially convey intermediate result Z (l) to compute final prediction, our AdaGCN transmits weights of nodes w i , aggregated features of different hops of neighbors Âl X. More importantly, in AdaGCN the embedding Z (l) is independent of the flow of computation in the network and the sparse adjacent matrix Â is also not directly involved in the computation of individual network because we compute Â(l+1) X in advance and then feed it instead of Â into the classifier f (l+1) θ , thus yielding significant computation reduction, which will be discussed further in Section 3.

GCN

Connection with PPNP and APPNP. We also established a strong connection between AdaGCN and previous state-of-the-art Personalized Propagation of Neural Predictions (PPNP) and Approximate PPNP (APPNP) (Klicpera et al., 2018) method that leverages personalized pagerank to reconstruct graph convolutions in order to use information from a large and adjustable neighborhood. The analysis can be summarized in the following Proposition 1. Proof can refer to Appendix A.3. Proposition 1. Suppose that γ is the teleport factor. Let matrix sequence {Z (l) } be from the output of each layer l in AdaGCN, then PPNP is equivalent to the Exponential Moving Average (EMA) with exponentially decreasing factor γ on {Z (l) } in a sharing parameters version, and its approximate version APPNP can be viewed as the approximated form of EMA with a limited number of terms. Proposition 1 illustrates that AdaGCN can be viewed as an adaptive form of APPNP, formulated as: Z = L l=0 α (l) f (l) θ ( Âl X) Specifically, the first discrepancy between AdaGCN and APPNP lies in the adaptive coefficient α (l) in AdaGCN determined by the error of l-th base classifier f (l) θ rather than fixed exponentially decreased weights in APPNP. In addition, AdaGCN employs classifier f (l) θ with different parameters to learn the embedding of different orders of neighbors, while APPNP shares these parameters in its form. We verified this benefit of our approach in our experiments shown in Section 4.2. Definition 1. General layer-wise Neighborhood Mixing: A graph convolution network has the ability to represent the layer-wise neighborhood mixing if for any b 0 , b 1 , ..., b L , there exists an injective mapping f with a setting of its parameters, such that the output of this graph convolution network can express the following formula: f L l=0 b l σ Âl X Proposition 2. AdaGCNs defined by our proposed approach (Eq. equation 7) are capable of representing general layer-wise neighborhood mixing, i.e., can meet the Definition 1. Albeit the similarity, AdaGCN distinguishes from MixHop in many aspects. Firstly, MixHop concatenates all outputs from each order of neighbors while we combines these predictions in an Adaboost way, which has theoretical generalization guarantee based on boosting theory Hastie et al. (2009) . Oono & Suzuki (2020) have recently derived the optimization and generalization guarantees of multi-scale GNNs, serving as the theoretical backbone of AdaGCN. Meantime, MixHop allows full linear mixing of different orders of neighboring features, while AdaGCN utilizes different nonlinear transformation f (l) θ among all layers, enjoying stronger expressive power.

3. ALGORITHM

In practice, we employ SAMME.R (Hastie et al., 2009) , the soft version of SAMME, in AdaGCN. SAMME.R (R for Real) algorithm (Hastie et al., 2009) leverages real-valued confidence-rated predictions, i.e., weighted probability estimates, rather than predicted hard labels in SAMME, in the prediction combination, which has demonstrated a better generalization and faster convergence than SAMME. We elaborate the final version of AdaGCN in Algorithm 1. We provide the analysis on the choice of model depth L in Appendix A.7, and then we elaborate the computational advantage of AdaGCN in the following. Analysis of Computational Advantage. Due to the similarity of graph convolution in Mix-Hop (Abu-El-Haija et al., 2019) , AdaGCN also requires no additional memory or computational complexity compared with previous GCN models. Meanwhile, our approach enjoys huge computational advantage compared with GCN-based models, e.g., PPNP and APPNP, stemming from excluding the additional computation involved in sparse tensors, such as the sparse tensor multiplication between Â and other dense tensors, in the forward and backward propagation of the neural network. Specifically, there are only L times sparse tensor operations for an AdaGCN model with L layers, i.e., Âl X = Â • ( Âl-1 X) for each layer l. This operation in each layer yields a dense tensor Obtain the weighted probability estimates p (l) ( X(l) ) for f (l) θ : p (l) k ( X(l) ) = Softmax(f (l) θ (c = k| X(l) )), k = 1, . . . , K 5: Compute the individual prediction h (l) k (x) for the current graph convolutional classifier f (l) θ : h (l) k ( X(l) ) ← (K -1) log p (l) k ( X(l) ) - 1 K k log p (l) k ( X(l) ) where k = 1, . . . , K.

6:

Adjust the node weights w i for each node x i with label y i on training set: w i ← w i • exp - K -1 K y i log p (l) (x i ) , i = 1, . . . , n 7: Re-normalize all weights w i .

8:

Update l+1-hop neighbor feature matrix X(l+1) : X(l+1) = Â X(l) 9: end for 10: Combine all predictions h (l) k ( X(l) ) for l = 0, ..., L. C(A, X) = arg max k L l=0 h (l) k ( X(l) ) 11: return Final combined prediction C(A, X). B l = Âl X for the l-th layer, which is then fed into the computation in a two-layer fully-connected network, i.e., f (l) θ (B l ) = ReLU(B l W (0) )W (1) . Due to the fact that dense tensor B l has been computed in advance, there is no other computation related to sparse tensors in the multiple forward and backward propagation procedures while training the neural network. By contrast, this multiple computation involved in sparse tensors in the GCN-based models, e.g., GCN: Â ReLU( ÂXW (0) )W (1) , is highly expensive. AdaGCN avoids these additional sparse tensor operations in the neural network and then attains huge computational efficiency. We demonstrate this viewpoint in the Section 4.3.

4. EXPERIMENTS

Experimental Setup. We select five commonly used graphs: CiteSeer, Cora-ML (Bojchevski & Günnemann, 2018; McCallum et al., 2000) , PubMed (Sen et al., 2008) , MS-Academic (Shchur et al., 2018) and Reddit. Dateset statistics are summarized in Table 1 . Recent graph neural networks suffer from overfitting to a single splitting of training, validation and test datasets (Klicpera et al., 2018) . To address this problem, inspired by (Klicpera et al., 2018) , we test all approaches on multiple random splits and initialization to conduct a rigorous study. Detailed dataset splittings are provided in Appendix A.6. (Buchnik & Cohen, 2018) and JK (jumping knowledge networks with concatenation) (Xu et al., 2018b) . In the computation part, we additionally compare AdaGCN with FastGCN (Chen et al., 2018) and GraphSAGE (Hamilton et al., 2017) . We refer to the result of baselines from (Klicpera et al., 2018) and the implementation of AdaGCN is adapted from APPNP. For AdaGCN, after the line search on hyper-parameters, we set h = 5000 hidden units for the first four datasets except Ms-academic with h = 3000, and 15, 12, 20 and 5 layers respectively due to the different graph structures. In addition, we set dropout rate to 0 for Citeseer and Cora-ML datasets and 0.2 for the other datasets and 5 × 10 -3 L 2 regularization on the first linear layer. We set weight decay as 1 × 10 -3 for Citeseer while 1 × 10 -4 for others. More detailed model parameters and analysis about our early stopping mechanism can be referred from Appendix A.6.

4.1. DESIGN OF DEEP GRAPH MODELS TO CIRCUMVENT OVERSMOOTHING EFFECT

It is well-known that GCN suffers from oversmoothing (Li et al., 2018) with the stacking of more graph convolutions. However, combination of knowledge from each layer to design deep graph models is a reasonable method to circumvent oversmoothing issue. In our experiment, we aim to explore the prediction performance of GCN, GCN with residual connection (Kipf & Welling, 2017), SGC and our AdaGCN with a growing number of layers. From Figure 3 , it can be easily observed that oversmoothing leads to the rapid decreasing of accuracy for GCN (blue line) as the layer increases. In contrast, the speed of smoothing (green line) of SGC is much slower than GCN due to the lack of ReLU analyzed in Section 2.1. Similarly, GCN with residual connection (yellow line) partially mitigates the oversmoothing effect of original GCN but fails to take advantage of information from different orders of neighbors to improve the prediction performance constantly. Remarkably, AdaGCN (red line) is able to consistently enhance the performance with the increasing of layers across the three datasets. This implies that AdaGCN can efficiently incorporate knowledge from different orders of neighbors and circumvent oversmoothing of original GCN in the process of constructing deep graph models. In addition, the fluctuation of performance for AdaGCN is much lower than GCN especially when the number of layer is large.

4.2. PREDICTION PERFORMANCE

We conduct a rigorous study of AdaGCN on four datasets under multiple splittings of dataset. The results from Table 2 suggest the state-of-the-art performance of our approach and the improvement compared with APPNP validates the benefit of adaptive form for our AdaGCN. More rigorously, p values under paired t test demonstrate the significance of improvement for our method. In the realistic setting, graphs usually have different labeled nodes and thus it is necessary to investigate the robust performance of methods on different number of labeled nodes. Here we utilize label rates to measure the different numbers of labeled nodes and then sample corresponding labeled nodes per class on graphs respectively. Table 3 presents the consistent state-of-the-art performance of AdaGCN under different label rates. An interesting manifestation from Table 3 is that AdaGCN yields more improvement on fewer label rates compared with APPNP, showing more efficiency on graphs with few labeled nodes. Inspired by the Layer Effect on graphs (Sun et al., 2019) , we argue that the increase of layers in AdaGCN can result in more benefits on the efficient propagation of label signals especially on graphs with limited labeled nodes. GCN on Cora-ML and MS-Academic datasets. Furthermore, with multiple power iterations involved in sparse tensors, APPNP unfortunately has relatively expensive computation cost. It should be noted that this computational advantage of AdaGCN is more significant when it comes to large datasets, e.g., Reddit. Table 4 demonstrates AdaGCN has the potential to perform much faster on larger datasets. Besides, we explore the computational cost of ReLU and sparse adjacency tensor with respect to the number of layers in the right part of Figure 4 . We focus on comparing AdaGCN with SGC and GCN as other GCN-based methods, such as GraphSAGE and APPNP, behave similarly with GCN. Particularly, we can easily observe that both SGC (green line) and GCN (red line) show a linear increasing tendency and GCN yields a larger slope arises from ReLU and more parameters. For SGC, stacking more layers directly is undesirable regarding the computation. Thus, a limited number of SGC layers is preferable with more advanced optimization techniques Wu et al. (2019) . It also shows that the computational cost involved sparse matrices in neural networks plays a dominant role in all the cost especially when the layer is large enough. In contrast, our AdaGCN (pink line) displays an almost constant trend as the layer increases simply because it excludes the extra computation involved in sparse tensors Â, such as  • • • Â ReLU( ÂXW (0) )W (1) θ with a fixed architecture in each layer while the layer-wise optimization, therefore displaying a nearly constant computation cost within each epoch although more epochs are normally needed in the entire layer-wise training. We leave the analysis of exact time and memory complexity of AdaGCN as future works, but boosting-based algorithms including AdaGCN is memory-efficient (Oono & Suzuki, 2020) .

5. DISCUSSIONS AND CONCLUSION

One potential concern is that AdaBoost (Hastie et al., 2009; Freund et al., 1999 ) is established on i.i.d. hypothesis while graphs have inherent data-dependent property. Fortunately, the statistical convergence and consistency of boosting (Lugosi & Vayatis, 2001; Mannor et al., 2003) can still be preserved when the samples are weakly dependent (Lozano et al., 2013) . More discussion can refer to Appendix A.5. In this paper, we propose a novel RNN-like deep graph neural network architecture called AdaGCNs. With the delicate architecture design, our approach AdaGCN can effectively explore and exploit knowledge from different orders of neighbors in an Adaboost way. Our work paves a way towards better combining different-order neighbors to design deep graph models rather than only stacking on specific type of graph convolution.

A APPENDIX A.1 RELATED WORKS ON DEEP GRAPH MODELS

A straightforward solution (Kipf & Welling, 2017; Xu et al., 2018b) inspired by ResNets (He et al., 2016) was by adding residual connections, but this practice was unsatisfactory both in prediction performance and computational efficiency towards building deep graph models, as shown in our experiments in Section 4.1 and 4.3. More recently, JK (Jumping Knowledge Networks (Xu et al., 2018b) ) introduced jumping connections into final aggregation mechanism in order to extract knowledge from different layers of graph convolutions. However, this straightforward change of GCN architecture exhibited inconsistent empirical performance for different aggregation operators, which cannot demonstrate the successful construction of deep layers. In addition, Graph powering-based method (Jin et al., 2019) implicitly leveraged more spatial information by extending classical spectral graph theory to robust graph theory, but they concentrated on defending adversarial attacks rather than model depth. LanczosNet (Liao et al., 2019) utilized Lanczos algorithm to construct low rank approximations of the graph Laplacian and then can exploit multi-scale information. Moreover, APPNP (Approximate Personalized Propagation of Neural Predictions, (Klicpera et al., 2018) ) leveraged the relationship between GCN and personalized PageRank to derive an improved global propagation scheme. Beyond these, DeepGCNs (Li et al., 2019) directly adapted residual, dense connection and dilated convolutions to GCN architecture, but it mainly focused on the task of point cloud semantic segmentation and has not demonstrated its effectiveness in typical graph tasks. Similar to our work, Deep Adaptive Graph Neural Network (DAGNN) (Liu et al., 2020 ) also focused on incorporating information from large receptive fields through the entanglement of representation transformation and propagation, while our work efficiently ensembles knowledge from large receptive fields in an Adaboost manner. Other related works based on global attention models (Puny et al., 2020) and sample-based methods (Zeng et al., 2019) are also helpful to construct deep graph models.

A.2 INSUFFICIENT REPRESENTATION POWER OF ADASGC

As illustrated in Figure 5 , with the increasing of layers, AdaSGC with only linear transformation has insufficient representation power both in extracting knowledge from high-order neighbors and combining information from different orders of neighbors while AdaGCN exhibits a consistent improvement of performance as the layer increases. 

A.3 PROOF OF PROPOSITION 1

Firstly, we further elaborate the Proposition 1 as follows, then we provide the proof. Suppose that γ is the teleport factor. Consider the output ZPPNP = γ(I -(1 -γ) Â) -1 f θ (X) in PPNP and ZAPPNP from its approxminated version APPNP. Let matrix sequence {Z (l) } be from the output of each layer l in AdaGCN, then PPNP is equivalent to the Exponential Moving Average (EMA) with exponentially decreasing factor γ, a first-order infinite impulse response filter, on {Z (l) } in a sharing parameters version, i.e., f (l) θ ≡ f θ . In addition, APPNP, which we reformulate in Eq. 10, can be viewed as the approximated form of EMA with a as a conference paper at ICLR 2021 limited number of terms.

ZAPPNP = (γ

L-1 l=0 (1 -γ) l Âl + (1 -γ) L ÂL )f θ (X) Proof. According to Neumann Theorem, ZPPNP can be expanded as a Neumann series: ZPPNP = γ(I -(1 -γ) Â) -1 f θ (X) = γ ∞ l=0 (1 -γ) l Âl f θ (X), where feature embedding matrix sequence {Z (l) } for each order of neighbors share the same parameters f θ . If we relax this sharing nature to the adaptive form with respect to the layer and put Âl into f θ , then the output Z can be approximately formulated as: ZPPNP ≈ γ ∞ l=0 (1 -γ) l f (l) θ ( Âl X) This relaxed version from PPNP is the Exponential Moving Average form of matrix sequence {Z (l) } with exponential decreasing factor γ. Moreover, if we approximate the EMA by truncating it after L -1 items, then the weight omitted by stopping after L -1 items is (1 -γ) L . Thus, the approximated EMA is exactly the APPNP form: ZAPPNP = (γ L-1 l=0 (1 -γ) l Âl + (1 -γ) L ÂL )f θ (X) A.4 PROOF OF PROPOSITION 2 Proof. We consider a two layers fully-connected neural network as f in Eq. 8, then the output of AdaGCN can be formulated as: Z = L l=0 α (l) σ( Âl XW (0) )W (1) Particularly, we set W (0) = b l sign(b l )α (l) I and W (1) = sign(b l )I where sign(b l ) is the signed incidence scalar w.r.t b l . Then the output of AdaGCN can be presented as: Z = L l=0 α (l) σ( Âl X b l sign(b l )α (l) I)sign(b l )I = L l=0 α (l) σ( Âl X) b l sign(b l )α (l) sign(b l ) = L l=0 b l σ Âl X The proof that GCNs-based methods are not capable of representing general layer-wise neighborhood mixing has been demonstrated in MixHop (Abu-El-Haija et al., 2019) . Proposition 2 proved. A.5 EXPLANATION ABOUT CONSISTENCY OF BOOSTING ON DEPENDENT DATA Definition 2. (β-mixing sequences.) Let σ j i = σ(W ) = σ(Wi, Wi+1, ..., Wj) be the σ-field generated by a strictly stationary sequence of random variables W = (Wi, Wi+1, ..., Wj). The β-mixing coefficient is defined by: βW (n) = sup k E sup P A|σ k 1 -P(A) : A ∈ σ ∞ k+n Then a sequence W is called β-mixing if limn→∞βW (n) = 0. Further, it is algebraically β-mixing if there is a positive constant r β such that βW (n) = O(n -r β ). Definition 3. (Consistency) A classification rule is consistent for a certain distribution P if E(L(hn)) = P {hn(X) = Y } → a as n → ∞ where a is a constant. It is strongly Bayes-risk consistent if limn→∞L(hn) = a almost surely. Under these definitions, the convergence and consistence of regularized boosting method on stationary βmixing sequences can be proved under mild assumptions. More details can be referred from (Lozano et al., 2013) .

A.6 EXPERIMENTAL DETAILS

Early Stopping on AdaGCN. We apply the same early stopping mechanism across all the methods as (Klicpera et al., 2018) for fair comparison. Furthermore, boosting theory also has the capacity to perfectly incorporate early stopping and it has been shown that for several boosting algorithms including AdaBoost, this regularization via early stopping can provide guarantees of consistency (Zhang et al., 2005; Jiang et al., 2004; Bühlmann & Yu, 2003) . Dataset Splitting. We choose a training set of a fixed nodes per class, an early stopping set of 500 nodes and test set of remained nodes. Each experiment is run with 5 random initialization on each data split, leading to a total of 100 runs per experiment. On a standard setting, we randomly select 20 nodes per class. For the two different label rates on each graph, we select 6, 11 nodes per class on citeseer, 8, 16 nodes per class on Cora-ML, 7, 14 nodes per class on Pubmed and 8, 15 nodes per class on MS-Academic dataset. Model parameters. For all GCN-based approaches, we use the same hyper-parameters in the original paper: learning rate of 0.01, 0.5 dropout rate, 5 × 10 -4 L2 regularization weight, and 16 hidden units. For FastGCN, we adopt the officially released code to conduct our experiments. PPNP and APPNP are adapted with best setting: K = 10 power iteration steps for APPNP, teleport probability γ = 0.1 on Cora-ML, Citeseer and Pubmed, γ = 0.2 on Ms-Academic. In addition, we use two layers with h = 64 hidden units and apply L2 regularization with λ = 5 × 10 -3 on the weights of the first layer and use dropout with dropout rate d = 0.5 on both layers and the adjacency matrix. The early stopping criterion uses a patience of p = 100 and an (unreachably high) maximum of n = 10000 epochs.The implementation of AdaGCN is adapted from PPNP and APPNP. Corresponding patience p = 300 and n = 500 in the early stopping of AdaGCN. Moreover, SGC is re-implemented in a straightforward way without incorporating advanced optimization for better illustration and comparison. Other baselines are adopted the same parameters described in PPNP and APPNP. Settings on Reddit dataset. By repeatedly tuning the parameters of these typical methods on Reddit, we finally choose weight decay rate as 10 -4 , hidden layer size 100 and epoch 20000 for AdaGCN. For APPNP, we opt weight decay rate as 10 -5 , dropout rate as 0 and epoch 500. V.GCN applies the same parameters in (Kipf & Welling, 2017) and we choose epoch as 500. All approaches have not deployed early stopping due to the expensive computational cost on the large Reddit dataset, which is also a fair comparison. A.7 CHOICE OF THE NUMBER OF LAYERS Different from the "forcible" behaviors in CNNs that directly stack many convolution layers, in our AdaGCN there is a theoretical guidance on the choice of model depth L, i.e., the number of base classifiers or layers, derived from boosting theory. Specifically, according to the boosting theory, the increasing of L can exponentially decreases the empirical loss, however, from the perspective of VC-dimension, an overly large L can yield overfitting of AdaGCN. It should be noted that the deeper graph convolution layers in AdaGCN are not always better, which indeed heavily depends on the the complexity of data. In practice, L can be determined via cross-validation. Specifically, we start a VC-dimension-based analysis to illustrate that too large L can yield overfitting of AdaGCN. For L layers of AdaGCN, its hypothesis set is FL = arg max k L l=1 α (l) f (l) θ : α (l) ∈ R, l ∈ [1, L] Then the VC-dimension of FT can be bounded as follows in terms of the VC-dimension d of the family of base hypothesis: VCdim (FL) ≤ 2(d + 1)(L + 1) log 2 ((L + 1)e), where e is a constant and the upper bounds grows as L increases. Combined with VC-dimension generalization bounds, these results imply that larger values of L can lead to overfitting of AdaBoost. This situation also happens in AdaGCN, which inspires us that there is no need to stack too many layers on AdaGCN in order to avoid overfitting. In practice, L is typically determined via cross-validation.



Figure 1: The RNN-like architecture of AdaGCN with each base classifier f (l)

Figure 2: Comparison of the graph model architectures. f a in JK network denotes one aggregation layer with aggregation function such as concatenation or max pooling.

AdaGCN based on SAMME.R AlgorithmInput: Features Matrix X, normalized adjacent matrix Â, a two-layer fully connected network f θ , number of layers L and number of classes K. Output: Final combined prediction C(A, X).1: Initialize the node weights w i = 1/n, i = 1, 2, ..., n on training set, neighbors feature matrix X(0) = X and classifier f

Figure 4: Left: Per-epoch training time of AdaGCN vs other methods under 5 runs on four datasets. Right: Per-epoch training time of AdaGCN compared with GCN and SGC with the increasing of layers and the digit after "k =" denotes the slope in a fitted linear regression.

Figure 5: AdaSGC vs AdaGCN.

Comparison with MixHop MixHop (Abu-El-Haija et al., 2019) applied the similar way of graph convolution by repeatedly mixing feature representations of neighbors at various distance. Proposition 2 proves that both AdaGCN and MixHop are able to represent feature differences among neighbors while previous GCNs-based methods cannot. Proof can refer to Appendix A.4. Recap the definition of general layer-wise Neighborhood Mixing (Abu-El-Haija et al., 2019) as follows:

Comparison of test accuracy of different models as the layer increases. We regard the l-th base classifier as the l-th layer in AdaGCN as both of them are leveraged to exploit the information from l-th order of neighbors for current nodes.

Average accuracy across different label rates with 20 splittings of datasets under 100 runs.



ACKNOWLEDGMENTS

Z. Lin is supported by NSF China (grant no.s 61625301 and 61731018), Major Scientific Research Project of Zhejiang Lab (grant no.s 2019KB0AC01 and 2019KB0AB02), Beijing Academy of Artificial Intelligence, and Qualcomm.

