Title: DROPEDGE: TOWARDS DEEP GRAPH CONVOLU-TIONAL NETWORKS ON NODE CLASSIFICATION

Abstract: Over-fitting and over-smoothing are two main obstacles of developing deep Graph Convolutional Networks (GCNs) for node classification. In particular, over-fitting weakens the generalization ability on small dataset, while over-smoothing impedes model training by isolating output representations from the input features with the increase in network depth. This paper proposes DropEdge, a novel and flexible technique to alleviate both issues. At its core, DropEdge randomly removes a certain number of edges from the input graph at each training epoch, acting like a data augmenter and also a message passing reducer. Furthermore, we theoretically demonstrate that DropEdge either reduces the convergence speed of over-smoothing or relieves the information loss caused by it. More importantly, our DropEdge is a general skill that can be equipped with many other backbone models (e.g. GCN, ResGCN, GraphSAGE, and JKNet) for enhanced performance. Extensive experiments on several benchmarks verify that DropEdge consistently improves the performance on a variety of both shallow and deep GCNs. The effect of DropEdge on preventing over-smoothing is empirically visualized and validated as well. Codes are released on https://github.com/DropEdge/DropEdge. * Wenbing Huang is the corresponding author. 1 When counting the number of layers (or network depth) of GCN, this paper does not involve the input layer.

1. INTRODUCTION:

Graph Convolutional Networks (GCNs), which exploit message passing or equivalently certain neighborhood aggregation function to extract high-level features from a node as well as its neighborhoods, have boosted the state-of-the-arts for a variety of tasks on graphs, such as node classification (Bhagat et al., 2011; Zhang et al., 2018) , social recommendation (Freeman, 2000; Perozzi et al., 2014) , and link prediction (Liben-Nowell & Kleinberg, 2007) to name some. In other words, GCNs have been becoming one of the most crucial tools for graph representation learning. Yet, when we revisit typical GCNs on node classification (Kipf & Welling, 2017) , they are usually shallow (e.g. the number of the layers is 2 1 ). Inspired from the success of deep CNNs on image classification, several attempts have been proposed to explore how to build deep GCNs towards node classification (Kipf & Welling, 2017; Li et al., 2018a; Xu et al., 2018a; Li et al., 2019) ; nevertheless, none of them delivers sufficiently expressive architecture. The motivation of this paper is to analyze the very factors that impede deeper GCNs to perform promisingly, and develop method to address them.

We begin by investigating two factors: over-fitting and over-smoothing. Over-fitting comes from the case when we utilize an over-parametric model to fit a distribution with limited training data, where the model we learn fits the training data very well but generalizes poorly to the testing data. It does exist if we apply a deep GCN on small graphs (see 4-layer GCN on Cora in Figure 1 ). Over-smoothing, towards the other extreme, makes training a very deep GCN difficult. As first introduced by Li et al. (2018a) and further explained in Wu et al. (2019) ; Xu et al. (2018a) ; Klicpera et al. (2019) , graph convolutions essentially push representations of adjacent nodes mixed with each Figure 1 : Performance of Multi-layer GCNs on Cora. We implement 4-layer GCN w and w/o DropEdge (in orange), 8-layer GCN w and w/o DropEdge (in blue)2 . GCN-4 gets stuck in the over-fitting issue attaining low training error but high validation error; the training of GCN-8 fails to converge satisfactorily due to over-smoothing. By applying DropEdge, both GCN-4 and GCN-8 work well for both training and validation. other, such that, if extremely we go with an infinite number of layers, all nodes' representations will converge to a stationary point, making them unrelated to the input features and leading to vanishing gradients. We call this phenomenon as over-smoothing of node features. To illustrate its influence, we have conducted an example experiment with 8-layer GCN in Figure 1 , in which the training of such a deep GCN is observed to converge poorly. Both of the above two issues can be alleviated, using the proposed method, DropEdge. The term "DropEdge" refers to randomly dropping out certain rate of edges of the input graph for each training time. There are several benefits in applying DropEdge for the GCN training (see the experimental improvements by DropEdge in Figure 1 ). First, DropEdge can be considered as a data augmentation technique. By DropEdge, we are actually generating different random deformed copies of the original graph; as such, we augment the randomness and the diversity of the input data, thus better capable of preventing over-fitting. Second, DropEdge can also be treated as a message passing reducer. In GCNs, the message passing between adjacent nodes is conducted along edge paths. Removing certain edges is making node connections more sparse, and hence avoiding over-smoothing to some extent when GCN goes very deep. Indeed, as we will draw theoretically in this paper, DropEdge either reduces the convergence speed of over-smoothing or relieves the information loss caused by it.

We are also aware that the dense connections employed by JKNet (Xu et al., 2018a) are another kind of tools that can potentially prevent over-smoothing. In its formulation, JKNet densely connects each hidden layer to the top one, hence the feature mappings in lower layers that are hardly affected by over-smoothing are still maintained. Interestingly and promisingly, we find that the performance of JKNet can be promoted further if it is utilized along with our DropEdge. Actually, our DropEdge-as a flexible and general technique-is able to enhance the performance of various popular backbone networks on several benchmarks, including GCN (Kipf & Welling, 2017) , ResGCN (Li et al., 2019) , JKNet (Xu et al., 2018a) , and GraphSAGE (Hamilton et al., 2017) . We provide detailed evaluations in the experiments.

2. RELATED WORK:

GCNs Inspired by the huge success of CNNs in computer vision, a large number of methods come redefining the notion of convolution on graphs under the umbrella of GCNs. The first prominent research on GCNs is presented in Bruna et al. (2013) , which develops graph convolution based on spectral graph theory. Later, Kipf & Welling (2017); Defferrard et al. (2016) ; Henaff et al. (2015) ; Li et al. (2018b) ; Levie et al. (2017) apply improvements, extensions, and approximations on spectralbased GCNs. To address the scalability issue of spectral-based GCNs on large graphs, spatial-based GCNs have been rapidly developed (Hamilton et al., 2017; Monti et al., 2017; Niepert et al., 2016; Gao et al., 2018) . These methods directly perform convolution in the graph domain by aggregating the information from neighbor nodes. Recently, several sampling-based methods have been proposed for fast graph representation learning, including the node-wise sampling methods (Hamilton et al., 2017) , the layer-wise approach (Chen et al., 2018) and its layer-dependent variant (Huang et al., 2018) . Specifically, GAT (Velickovic et al., 2018) has discussed applying dropout on edge attentions. While it actually is a post-conducted version of DropEdge before attention computation, the relation to over-smoothing is never explored in Velickovic et al. (2018) . In our paper, however, we have formally presented the formulation of DropEdge and provided rigorous theoretical justification of its benefit in alleviating over-smoothing. We also carried out extensive experiments by imposing DropEdge on several popular backbones. One additional point is that we further conduct adjacency normalization after dropping edges, which, even simple, is able to make it much easier to converge during training and reduce gradient vanish as the number of layers grows.

Deep GCNs Despite the fruitful progress, most previous works only focus on shallow GCNs while the deeper extension is seldom discussed. The attempt for building deep GCNs is dated back to the GCN paper (Kipf & Welling, 2017) , where the residual mechanism is applied; unexpectedly, as shown in their experiments, residual GCNs still perform worse when the depth is 3 and beyond. The authors in Li et al. (2018a) first point out the main difficulty in constructing deep networks lying in over-smoothing, but unfortunately, they never propose any method to address it. The follow-up study (Klicpera et al., 2019) solves over-smoothing by using personalized PageRank that additionally involves the rooted node into the message passing loop; however, the accuracy is still observed to decrease when the depth increases from 2. JKNet (Xu et al., 2018a) employs dense connections for multi-hop message passing which is compatible with DropEdge for formulating deep GCNs. Oono & Suzuki (2019) theoretically prove that the node features of deep GCNs will converge to a subspace and incur information loss. It generalizes the conclusion in Li et al. (2018a) by further considering the ReLu function and convolution filters. Our interpretations on why DropEdge can impede over-smoothing is based on the concepts proposed by Oono & Suzuki (2019) . A recent method (Li et al., 2019) has incorporated residual layers, dense connections and dilated convolutions into GCNs to facilitate the development of deep architectures. Nevertheless, this model is targeted on graph-level classification (i.e. point cloud segmentation), where the data points are graphs and naturally disconnected from each other. In our task for node classification, the samples are nodes and they all couple with each other, thus the over-smoothing issue is more demanded to be addressed. By leveraging DropEdge, we are able to relieve over-smoothing, and derive more enhanced deep GCNs on node classification.

3. NOTATIONS AND PRELIMINARIES:

Notations. Let G = (V, E) represent the input graph of size N with nodes v i ∈ V and edges (v i , v j ) ∈ E. The node features are denoted as

X = {x 1 , • • • , x N } ∈ R N ×C ,

and the adjacency matrix is defined as A ∈ R N ×N which associates each edge (v i , v j ) with its element A ij . The node degrees are given by d = {d 1 , • • • , d N } where d i computes the sum of edge weights connected to node i. We define D as the degree matrix whose diagonal elements are obtained from d.

GCN is originally developed by Kipf & Welling (2017) . The feed forward propagation in GCN is recursively conducted as

EQUATION

where

H (l+1) = {h (l+1) 1 , • • • , h (l+1) N

} are the hidden vectors of the l-th layer with h (l)

i as the hidden feature for node i; Â = D-1/2 (A + I) D-1/2 is the re-normalization of the adjacency matrix, and D is the corresponding degree matrix of A + I; σ(•) is a nonlinear function, i.e. the ReLu function; and W (l) ∈ R C l ×C l-1 is the filter matrix in the l-th layer with C l refers to the size of l-th hidden layer. We denote one-layer GCN computed by Equation 1 as Graph Convolutional Layer (GCL) in what follows.

4. OUR METHOD: DROPEDGE:

This section first introduces the methodology of the DropEdge technique as well as its layer-wise variant where the adjacency matrix for each GCN layer is perturbed individually. We also explain how the proposed DropEdge can prevent over-fitting and over-smoothing in generic GCNs. Particularly for over-smoothing, we provide its mathematical definition and theoretical derivations on showing the benefits of DropEdge.

4.1. METHODOLOGY:

At each training epoch, the DropEdge technique drops out a certain rate of edges of the input graph by random. Formally, it randomly enforces V p non-zero elements of the adjacency matrix A to be zeros, where V is the total number of edges and p is the dropping rate. If we denote the resulting adjacency matrix as A drop , then its relation with A becomes

EQUATION

where A is a sparse matrix expanded by a random subset of size V p from original edges E. Following the idea of Kipf & Welling (2017) , we also perform the re-normalization trick on A drop , leading to Âdrop . We replace Â with Âdrop in Equation 1 for propagation and training. When validation and testing, DropEdge is not utilized.

Preventing over-fitting. DropEdge produces varying perturbations of the graph connections. As a result, it generates different random deformations of the input data and can be regarded as a data augmentation skill for graphs. To explain why this is valid, we provide an intuitive understanding here. The key in GCNs is to aggregate neighbors' information for each node, which can be understood as a weighted sum of the neighbor features (the weights are associated with the edges). From the perspective of neighbor aggregation, DropEdge enables a random subset aggregation instead of the full aggregation during GNN training. Statistically, DropEdge only changes the expectation of the neighbor aggregation up to a multiplier p, if we drop edges with probability p. This multiplier will be actually removed after weights normalization, which is often the case in practice. Therefore, DropEdge does not change the expectation of neighbor aggregation and is an unbiased data augmentation technique for GNN training, similar to typical image augmentation skills (e.g. rotation, cropping and flapping) that are capable of hindering over-fitting in training CNNs. We will provide experimental validations in § 5.1.

Layer-Wise DropEdge. The above formulation of DropEdge is one-shot with all layers sharing the same perturbed adjacency matrix. Indeed, we can perform DropEdge for each individual layer. Specifically, we obtain Â(l) drop by independently computing Equation 2 for each l-th layer. Different layer could have different adjacency matrix Â(l) drop . Such layer-wise version brings in more randomness and deformations of the original data, and we will experimentally compare its performance with the original DropEdge in § 5.2.

Over-smoothing is another obstacle of training deep GCNs, and we will detail how DropEdge can address it to some extent in the next section. For simplicity, the following derivations assume all GCLs share the same perturbed adjacency matrix, and we will leave the discussion on layer-wise DropEdge for future exploration.

4.2. TOWARDS PREVENTING OVER-SMOOTHING:

By its original definition in Li et al. (2018a) , the over-smoothing phenomenon implies that the node features will converge to a fixed point as the network depth increases. This unwanted convergence restricts the output of deep GCNs to be only relevant to the graph topology but independent to the input node features, which as a matter of course incurs detriment of the expressive power of GCNs. Oono & Suzuki (2019) has generalized the idea in Li et al. (2018a) by taking both the non-linearity (i.e. the ReLu function) and the convolution filters into account; they explain over-smoothing as convergence to a subspace rather than convergence to a fixed point. This paper will use the concept of subspace by Oono & Suzuki (2019) for more generality.

We first provide several relevant definitions that facilitate our later presentations.

Definition 1 (subspace). Let M := {EC|C ∈ R M ×C } be an M -dimensional subspace in R N ×C , where E ∈ R N ×M is orthogonal, i.e. E T E = I M , and M ≤ N .

Definition 2 ( -smoothing). We call the -smoothing of node features happens for a GCN, if all its hidden vectors H (l) beyond a certain layer L have a distance no larger than ( > 0) with respect to a subspace M that is independent to the input features, namely,

EQUATION

where d M (•) computes the distance between the input matrix and the subspace M.3 Definition 3 (the -smoothing layer). Given the subspace M and , we call the minimal value of the layers that satisfy Equation 3 as the -smoothing layer, that is, l * (M, ) := min l {d M (H (l) ) < }.

Since conducting analysis exactly based on the -smoothing layer is difficult, we instead define the relaxed -smoothing layer which is proved to be an upper bound of l * .

Definition 4 (the relaxed -smoothing layer). Given the subspace M and , we call l(M, ) =

log( /d M (X)) log sλ

as the relaxed smoothing layer, where, • computes the ceil of the input, s is the supremum of the filters' singular values over all layers, and λ is the second largest eigenvalue of Â.

Besides, we have l ≥ l *4 .

According to the conclusions by the authors in Oono & Suzuki (2019) , a sufficiently deep GCN will certainly suffer from the -smoothing issue for any small value of under some mild conditions (the details are included in the supplementary material). Note that they only prove the existence of -smoothing in deep GCN without developing any method to address it.

Here, we will demonstrate that adopting DropEdge alleviates the -smoothing issue in two aspects: 1. By reducing node connections, DropEdge is proved to slow down the convergence of over-smoothing; in other words, the value of the relaxed -smoothing layer will only increase if using DropEdge. 2.

The gap between the dimensions of the original space and the converging subspace, i.e. N -M measures the amount of information loss; larger gap means more severe information loss. As shown by our derivations, DropEdge is able to increase the dimension of the converging subspace, thus capable of reducing information loss.

We summarize our conclusions as follows.

Theorem 1. We denote the original graph as G and the one after dropping certain edges out as G .

Given a small value of , we assume G and G will encounter the -smoothing issue with regard to subspaces M and M , respectively. Then, either of the following inequalities holds after sufficient edges removed.

• The relaxed smoothing layer only increases: l(M, ) ≤ l(M , );

• The information loss is decreased:

N -dim(M) > N -dim(M ).

The proof of Theorem 1 is based on the derivations in Oono & Suzuki (2019) as well as the concept of mixing time that has been studied in the random walk theory (Lovász et al., 1993) . We provide the full details in the supplementary material. Theorem 1 tells that DropEdge either reduces the convergence speed of over-smoothing or relieves the information loss caused by it. In this way, DropEdge enables us to train deep GCNs more effectively.

4.3. DISCUSSIONS:

This sections contrasts the difference between DropEdge and other related concepts including Dropout, DropNode, and Graph Sparsification.

DropEdge vs. Dropout The Dropout trick (Hinton et al., 2012) is trying to perturb the feature matrix by randomly setting feature dimensions to be zeros, which may reduce the effect of over-fitting but is of no help to preventing over-smoothing since it does not make any change of the adjacency matrix. As a reference, DropEdge can be regarded as a generation of Dropout from dropping feature dimensions to dropping edges, which mitigates both over-fitting and over-smoothing. In fact, the impacts of Dropout and DropEdge are complementary to each other, and their compatibility will be shown in the experiments.

DropEdge vs. DropNode Another related vein belongs to the kind of node sampling based methods, including GraphSAGE (Hamilton et al., 2017) , FastGCN (Chen et al., 2018) , and AS-GCN (Huang et al., 2018) . We name this category of approaches as DropNode. For its original motivation, DropNode samples sub-graphs for mini-batch training, and it can also be treated as a specific form of dropping edges since the edges connected to the dropping nodes are also removed. However, the effect of DropNode on dropping edges is node-oriented and indirect. By contrast, DropEdge is edge-oriented, and it is possible to preserve all node features for the training (if they can be fitted into the memory at once), exhibiting more flexibility. Further, to maintain desired performance, the sampling strategies in current DropNode methods are usually inefficient, for example, GraphSAGE suffering from the exponentially-growing layer size, and AS-GCN requiring the sampling to be conducted recursively layer by layer. Our DropEdge, however, neither increases the layer size as the depth grows nor demands the recursive progress because the sampling of all edges are parallel.

DropEdge vs. Graph-Sparsification Graph-Sparsification (Eppstein et al., 1997 ) is an old research topic in the graph domain. Its optimization goal is removing unnecessary edges for graph compressing while keeping almost all information of the input graph. This is clearly district to the purpose of DropEdge where no optimization objective is needed. Specifically, DropEdge will remove the edges of the input graph by random at each training time, whereas Graph-Sparsification resorts to a tedious optimization method to determine which edges to be deleted, and once those edges are discarded the output graph keeps unchanged.

5. EXPERIMENTS:

Datasets Joining the previous works' practice, we focus on four benchmark datasets varying in graph size and feature type: (1) classifying the research topic of papers in three citation datasets: Cora, Citeseer and Pubmed (Sen et al., 2008) ; (2) predicting which community different posts belong to in the Reddit social network (Hamilton et al., 2017) . Note that the tasks in Cora, Citeseer and Pubmed are transductive underlying all node features are accessible during training, while the task in Reddit is inductive meaning the testing nodes are unseen for training. We apply the full-supervised training fashion used in Huang et al. (2018) and Chen et al. (2018) on all datasets in our experiments.

The statics of all datasets are listed in the supplemental materials.

5.1. CAN DROPEDGE GENERALLY IMPROVE THE PERFORMANCE OF DEEP GCNS?:

In this section, we are interested in if applying DropEdge can promote the performance of current popular GCNs (especially their deep architectures) on node classification.

Implementations We consider five backbones: GCN (Kipf & Welling, 2017), ResGCN (He et al., 2016; Li et al., 2019) , JKNet (Xu et al., 2018a) , IncepGCN5 and GraphSAGE (Hamilton et al., 2017) with varying depth from 2 to 64. 6 Since different structure exhibits different training dynamics on different dataset, to enable more robust comparisons, we perform random hyper-parameter search for each model, and report the case giving the best accuracy on validation set of each benchmark. The searching space of hyper-parameters and more details are provided in Table 4 in the supplementary material. Regarding the same architecture w or w/o DropEdge, we apply the same set of hyperparameters except the drop rate p for fair evaluation.

Overall Results Table 1 summaries the results on all datasets. We only report the performance of the model with 2/8/32 layers here due to the space limit, and provide the accuracy under other different depths in the supplementary material. It's observed that DropEdge consistently improves the testing accuracy for all cases. The improvement is more clearly depicted in Figure 2a , where we have computed the average absolute improvement over all backbones by DropEdge on each dataset under different numbers of layers. On Citeseer, for example, DropEdge yields further improvement for deeper architecture; it gains 0.9% average improvement for the model with 2 layers while achieving a remarkable 13.5% increase for the model with 64 layers. In addition, the validation losses of all 4-layer models on Cora are shown in Figure 2b . The curves along the training epoch are dramatically pulled down after applying DropEdge, which also explains the effect of DropEdge on alleviating over-fitting. Another valuable observation in Table 1 is that the 32-layer IncepGCN without DropEdge incurs the Out-Of-Memory (OOM) issue while the model with DropEdge survives, showing the advantage of DropEdge to save memory consuming by making the adjacency matrix sparse.

Comparison with SOTAs We select the best performance for each backbone with DropEdge, and contrast them with existing State of the Arts (SOTA), including GCN, FastGCN, AS-GCN and GraphSAGE in Table 2 ; for the SOTA methods, we reuse the results reported in Huang et al. (2018) .

We have these findings: (1) Clearly, our DropEdge obtains significant enhancement against SOTAs; particularly on Reddit, the best accuracy by our method is 97.02%, and it is better than the previous best by AS-GCN (96.27%), which is regarded as a remarkable boost considering the challenge on this benchmark. (2) For most models with DropEdge, the best accuracy is obtained under the depth beyond 2, which again verifies the impact of DropEdge on formulating deep networks.

(3) As mentioned in § 4.3, FastGCN, AS-GCN and GraphSAGE are considered as the DropNode extensions of GCNs. The DropEdge based approaches outperform the DropNode based variants as shown in Table 2 , which somehow confirms the effectiveness of DropEdge. Actually, employing DropEdge upon the DropNode methods further delivers promising enhancement, which can be checked by revisiting the increase by DropEdge for GraphSAGE in Table 1 .

5.2. HOW DOES DROPEDGE HELP?:

This section continues a more in-depth analysis on DropEdge and attempts to figure out why it works. Due to the space limit, we only provide the results on Cora, and defer the evaluations on other datasets to the supplementary material.

Note that this section mainly focuses on analyzing DropEdge and its variants, without the concern with pushing state-of-the-art results. So, we do not perform delicate hyper-parameter selection. We employ GCN as the backbone in this section. Here, GCN-n denotes GCN of depth n. random seed is fixed. We train all models with 200 epochs. Unless otherwise mentioned, we do not utilize the "withloop" and "withbn" operation (see their definitions in Table 4 in the appendix).

5.2.1. ON PREVENTING OVER-SMOOTHING:

As discussed in § 4.2, the over-smoothing issue exists when the top-layer outputs of GCN converge to a subspace and become unrelated to the input features with the increase in depth. Since we are unable to derive the converging subspace explicitly, we measure the degree of over-smoothing by instead computing the difference between the output of the current layer and that of the previous one. We adopt the Euclidean distance for the difference computation. Lower distance means more serious over-smoothing. Experiments are conducted on GCN-8.

Figure 3 (a) shows the distances of different intermediate layers (from 2 to 6) under different edge dropping rates (0 and 0.8). Clearly, over-smoothing becomes more serious in GCN as the layer grows, which is consistent with our conjecture. Conversely, the model with DropEdge (p = 0.8) reveals higher distance and slower convergent speed than that without DropEdge (p = 0), implying the importance of DropEdge to alleviating over-smoothing. We are also interested in how the oversmoothing will act after training. For this purpose, we display the results after 150-epoch training in Figure 3 (b). For GCN without DropEdge, the difference between outputs of the 5-th and 6-th layers is equal to 0, indicating that the hidden features have converged to a certain stationary point. On the contrary, GCN with DropEdge performs promisingly, as the distance does not vanish to zero when Here, we provide the experimental evaluation on assessing its effect. As observed from Figure 4b , the LW DropEdge achieves lower training loss than the original version, whereas the validation value between two models is comparable.

It implies that LW DropEdge can facilitate the training further than original DropEdge. However, we prefer to use DropEdge other than the LW variant so as to not only avoid the risk over-fitting but also reduces computational complexity since LW DropEdge demands to sample each layer and spends more time.

6. CONCLUSION:

We A APPENDIX: PROOF OF THEOREM 1

To prove Theorem 1, we need to borrow the following definitions and corollaries from Oono & Suzuki (2019) . First, we denote the maximum singular value of W l by s l and set s := sup l∈N+ s l . We assume that W l of all layers are initialized so that s ≤ 1. Second, we denote the distance that induced as the Frobenius norm from X to M by d M (X) := inf Y ∈M ||X -Y || F . Then, we recall Corollary 3 and Proposition 1 in Oono & Suzuki (2019) as Corollary 1 below.

Corollary 1. Let λ 1 ≤ • • • ≤ λ N be the eigenvalues of Â, sorted in ascending order. Suppose the multiplicity of the largest eigenvalue λ N is M (≤ N ), i.e., λ N -M < λ N -M +1 = • • • = λ N and the second largest eigenvalue is defined as

EQUATION

Let E to be the eigenspace associated with λ N -M +1 , • • • , λ N . Then we have λ < λ N = 1, and

EQUATION

where M := {EC|C ∈ R M ×C }. Besides, s l λ < 1, implying that the output of the l-th layer of GCN on G exponentially approaches M.

We also need to adopt some concepts from Lovász et al. (1993) in proving Theorem 1. Consider the graph G as an electrical network, where each edge represents an unit resistance. Then the effective resistance, R st from node s to node t is defined as the total resistance between node s and t.

According to Corollary 3.3 and Theorem 4.1 (i) in Lovász et al. (1993) , we can build the connection between λ and R st for each connected component via commute time as the following inequality.

EQUATION

Prior to proving Theorem 1, we first derive the lemma below. Lemma 2. The -smoothing happens whenever the layer number satisfies

EQUATION

where • computes the ceil of the input. It means l ≥ l * .

Proof. We start our proof from Inequality 5, leading to

d M (H (l) ) ≤ s l λd M (H (l-1) ) ≤ ( l i=1 s i )λ l d M (X) ≤ s l λ l d M (X)

When it reaches -smoothing, the following inequality should be satisfied as

EQUATION

)

Since 0 ≤ sλ < 1, then log sλ < 0. Therefore, the Inequality 8 becomes

EQUATION

Clearly, we have l ≥ l * since l * is defined as the minimal layer that satisfies -smoothing. The proof is concluded. Now, we prove Theorem 1.

Proof. Our proof relies basically on the connection between λ and R st in Equation ( 6). We recall Corollary 4.3 in Lovász et al. (1993) that removing any edge from G can only increase any R st , then according to (6), the lower bound of λ only increases if the removing edge is not connected to either s or t (i.e. the degree d s and d t keep unchanged). Since there must exist a node pair satisfying R st = ∞ after sufficient edges (except self-loops) are removed from one connected component of G, we have the infinite case λ = 1 given in Equation ( 6) that both 1/d s and 1/d t are consistently bounded by a finite number,i.e. 1. It implies λ does increase before it reaches λ = 1. As l is positively related to λ (see the right side of Equation ( 7) where log(sλ) <0), we have proved the first part of Theorem 1, i.e., l(M, ) ≤ l(M , ) after removing sufficient edges.

When there happens R st = ∞, the connected component is disconnected into two parts, which leads to the increment of the dimension of M by 1 and proves the second part of Theorem 1. i.e., the information loss is decreased:

N -dim(M) > N -dim(M ).

B APPENDIX: MORE DETAILS IN EXPERIMENTS B.1 DATASETS STATISTICS:

Datasets The statistics of all datasets are summarized in Table 3 . Backbones Other than the multi-layer GCN, we replace the CNN layer with graph convolution layer to implement three popular backbones recasted from image classification. They are residual network (ResGCN) (He et al., 2016; Li et al., 2019) , inception network (IncepGCN) (Szegedy et al., 2016) and dense network (JKNet) (Huang et al., 2017; Xu et al., 2018b) . Figure 5 shows the detailed architectures of four backbones. Furthermore, we employ one input GCL and one output GCL on these four backbones. Therefore, the layers in ResGCN, JKNet and InceptGCN are at least 3 layers. All backbones are implemented in Pytorch (Paszke et al., 2017) . For GraphSAGE, we utilize the Pytorch version implemented by DGL (Wang et al., 2019) . Self Feature Modeling We also implement a variant of graph convolution layer with self feature modeling (Fout et al., 2017) :

EQUATION

where

W (l) self ∈ R C l ×C l-1 .

Hyper-parameter Optimization We adopt the Adam optimizer for model training. To ensure the re-productivity of the results, the seeds of the random numbers of all experiments are set to the same. We fix the number of training epoch to 400 for all datasets. All experiments are conducted on a NVIDIA Tesla P40 GPU with 24GB memory.

Given a model with n ∈ {2, 4, 8, 16, 32, 64} layers, the hidden dimension is 128 and we conduct a random search strategy to optimize the other hyper-parameter for each backbone in § 5.1. The decryptions of hyper-parameters are summarized in Table 4 . Table 5 depicts the types of the normalized adjacency matrix that are selectable in the "normalization" hyper-parameter. For GraphSAGE, the aggregation type like GCN, MAX, MEAN, or LSTM is a hyper-parameter as well.

For each model, we try 200 different hyper-parameter combinations via random search and select the best test accuracy as the result. Table 6 summaries the hyper-parameters of each backbone with the best accuracy on different datasets and their best accuracy are reported in Table 2 . 

:

To check the efficacy of DropEdge more clearly, here we have removed bias in all GCN layers, while for the experiments in § 5, the bias are kept.

The definition of dM(•) is provided in the supplementary material.

All detailed definitions and proofs are provided in the appendix.

The formulation is given in the appendix.

For Reddit, the maximum depth is 32 considering the memory bottleneck.

Figures and Tables:

Figure3: Analysis on over-smoothing. Smaller distance means more serious over-smoothing.

Figure 4 the number of layers grows; it probably has successfully learned meaningful node representations after training, which could also be validated by the training loss in Figure 3 (c).

Figure 5: The illustration of four backbones. GCL indicates graph convolutional layer.

Figure 6: The validation loss on different backbones w and w/o DropEdge. GCN-n denotes PlainGCN of depth n; similar denotation follows for other backbones.

Figure 7

Testing accuracy (%) comparisons on different backbones w and w/o DropEdge.

Table content: <table><tr><td/><td/><td>2 layers</td><td/><td>8 layers</td><td/><td>32 layers</td><td/></tr><tr><td>Dataset</td><td>Backbone</td><td colspan="6">Orignal DropEdge Orignal DropEdge Orignal DropEdge</td></tr><tr><td/><td>GCN</td><td>86.10</td><td>86.50</td><td>78.70</td><td>85.80</td><td>71.60</td><td>74.60</td></tr><tr><td/><td>ResGCN</td><td>-</td><td>-</td><td>85.40</td><td>86.90</td><td>85.10</td><td>86.80</td></tr><tr><td>Cora</td><td>JKNet IncepGCN</td><td>--</td><td>--</td><td>86.70 86.70</td><td>87.80 88.20</td><td>87.10 87.40</td><td>87.60 87.70</td></tr><tr><td/><td>GraphSAGE</td><td>87.80</td><td>88.10</td><td>84.30</td><td>87.10</td><td>31.90</td><td>32.20</td></tr><tr><td/><td>GCN</td><td>75.90</td><td>78.70</td><td>74.60</td><td>77.20</td><td>59.20</td><td>61.40</td></tr><tr><td/><td>ResGCN</td><td>-</td><td>-</td><td>77.80</td><td>78.80</td><td>74.40</td><td>77.90</td></tr><tr><td>Citeseer</td><td>JKNet IncepGCN</td><td>--</td><td>--</td><td>79.20 79.60</td><td>80.20 80.50</td><td>71.70 72.60</td><td>80.00 80.30</td></tr><tr><td/><td>GraphSAGE</td><td>78.40</td><td>80.00</td><td>74.10</td><td>77.10</td><td>37.00</td><td>53.60</td></tr><tr><td/><td>GCN</td><td>90.20</td><td>91.20</td><td>90.10</td><td>90.90</td><td>84.60</td><td>86.20</td></tr><tr><td/><td>ResGCN</td><td>-</td><td>-</td><td>89.60</td><td>90.50</td><td>90.20</td><td>91.10</td></tr><tr><td>Pubmed</td><td>JKNet IncepGCN</td><td>--</td><td>--</td><td>90.60 90.20</td><td>91.20 91.50</td><td>89.20 OOM</td><td>91.30 90.50</td></tr><tr><td/><td>GraphSAGE</td><td>90.10</td><td>90.70</td><td>90.20</td><td>91.70</td><td>41.30</td><td>47.90</td></tr><tr><td/><td>GCN</td><td>96.11</td><td>96.13</td><td>96.17</td><td>96.48</td><td>45.55</td><td>50.51</td></tr><tr><td/><td>ResGCN</td><td>-</td><td>-</td><td>96.37</td><td>96.46</td><td>93.93</td><td>94.27</td></tr><tr><td>Reddit</td><td>JKNet IncepGCN</td><td>--</td><td>--</td><td>96.82 96.43</td><td>97.02 96.87</td><td>OOM OOM</td><td>OOM OOM</td></tr><tr><td/><td>GraphSAGE</td><td>96.22</td><td>96.28</td><td>96.38</td><td>96.42</td><td>96.43</td><td>96.47</td></tr></table>

Accuracy (%) comparisons with SOTAs. The number in parenthesis denotes the network depth for the models with DropEdge.

Table content: <table><tr><td colspan="5">Published as a conference paper at ICLR 2020</td><td/><td/><td/></tr><tr><td>0% 2% 4% 6% 8% 10% 12% 14% 16%</td><td>Cora</td><td>Citeseer</td><td colspan="2">Pubmed 2 layers 8 layers 32 layers</td><td colspan="2">Reddit 4 layers 16 layers 64 layers</td><td>Validation Loss</td><td>0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2</td><td>Epoch 0 25 50 75 100 125 150 175 200 Cora ResGCN-4 GCN-4 InceptGCN-4 JKNet-4 ResGCN-4+DropEdge GCN-4+DropEdge InceptGCN-4+DropEdge JKNet-4+DropEdge</td></tr><tr><td/><td colspan="6">(a) The average absolute improvement by DropEdge.</td><td colspan="2">(b) The validation loss on different</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td colspan="2">backbones w and w/o DropEdge.</td></tr><tr><td/><td/><td/><td/><td colspan="3">Figure 2</td><td/></tr><tr><td/><td/><td/><td/><td/><td/><td>Transductive</td><td/><td>Inductive</td></tr><tr><td/><td/><td/><td/><td>Cora</td><td/><td>Citeseer</td><td colspan="2">Pubmed</td><td>Reddit</td></tr><tr><td/><td>GCN</td><td/><td/><td colspan="2">86.64</td><td>79.34</td><td/><td>90.22</td><td>95.68</td></tr><tr><td/><td>FastGCN</td><td/><td/><td colspan="2">85.00</td><td>77.60</td><td/><td>88.00</td><td>93.70</td></tr><tr><td/><td>ASGCN</td><td/><td/><td colspan="2">87.44</td><td>79.66</td><td/><td>90.60</td><td>96.27</td></tr><tr><td/><td colspan="2">GraphSAGE</td><td/><td colspan="2">82.20</td><td>71.40</td><td/><td>87.10</td><td>94.32</td></tr><tr><td/><td colspan="2">GCN+DropEdge</td><td/><td colspan="2">87.60(4)</td><td>79.20(4)</td><td colspan="2">91.30(4)</td><td>96.71(4)</td></tr><tr><td/><td colspan="2">ResGCN+DropEdge</td><td/><td colspan="5">87.00(4) 79.40(16) 91.10(32) 96.48(16)</td></tr><tr><td/><td colspan="2">JKNet+DropEdge</td><td/><td colspan="2">88.00(16)</td><td colspan="3">80.20(8) 91.60(64)</td><td>97.02(8)</td></tr><tr><td/><td colspan="2">IncepGCN+DropEdge</td><td/><td colspan="2">88.20(8)</td><td>80.50(8)</td><td colspan="2">91.60(4)</td><td>96.87(8)</td></tr><tr><td/><td colspan="3">GraphSAGE+DropEdge</td><td colspan="2">88.10(4)</td><td>80.00(2)</td><td colspan="2">91.70(8)</td><td>96.54(4)</td></tr></table>

have presented DropEdge, a novel and efficient technique to facilitate the development of deep Graph Convolutional Networks (GCNs). By dropping out a certain rate of edges by random, DropEdge includes more diversity into the input data to prevent over-fitting, and reduces message passing in graph convolution to alleviate over-smoothing. Considerable experiments on Cora, Citeseer, Pubmed and Reddit have verified that DropEdge can generally and consistently promote the performance of current popular GCNs, such as GCN, ResGCN, JKNet, IncepGCN, and GraphSAGE. It is expected that our research will open up a new venue on a more in-depth exploration of deep GCNs for broader potential applications.

Table content: <table/>

Dataset Statistics

Table content: <table><tr><td colspan="2">Datasets Nodes</td><td>Edges</td><td colspan="4">Classes Features Traing/Validation/Testing Type</td></tr><tr><td>Cora</td><td>2,708</td><td>5,429</td><td>7</td><td>1,433</td><td>1,208/500/1,000</td><td>Transductive</td></tr><tr><td>Citeseer</td><td>3,327</td><td>4,732</td><td>6</td><td>3,703</td><td>1,812/500/1,000</td><td>Transductive</td></tr><tr><td>Pubmed</td><td>19,717</td><td>44,338</td><td>3</td><td>500</td><td>18,217/500/1,000</td><td>Transductive</td></tr><tr><td>Reddit</td><td colspan="2">232,965 11,606,919</td><td>41</td><td>602</td><td>152,410/23,699/55,334</td><td>Inductive</td></tr><tr><td colspan="3">B.2 MODELS AND BACKBONES</td><td/><td/><td/><td/></tr></table>

Hyper-parameter Description

Table content: <table><tr><td colspan="2">Hyper-parameter Description</td></tr><tr><td>lr</td><td>learning rate</td></tr><tr><td>weight-decay</td><td>L2 regulation weight</td></tr><tr><td colspan="2">sampling-percent edge preserving percent (1 -p)</td></tr><tr><td>dropout</td><td>dropout rate</td></tr><tr><td>normalization</td><td>the propagation models (Kipf &amp; Welling, 2017)</td></tr><tr><td>withloop</td><td>using self feature modeling</td></tr><tr><td>withbn</td><td>using batch normalization</td></tr></table>

The normalization / propagation models

Table content: <table><tr><td>Description</td><td>Notation</td><td>A</td></tr><tr><td>First-order GCN</td><td>FirstOrderGCN</td><td>I + D -1/2 AD -1/2</td></tr><tr><td>Augmented Normalized Adjacency</td><td>AugNormAdj</td><td>(D + I) -1/2 (A + I)(D + I) -1/2</td></tr><tr><td>Augmented Normalized Adjacency with Self-loop</td><td>BingGeNormAdj</td><td>I + (D + I) -1/2 (A + I)(D + I) -1/2</td></tr><tr><td>Augmented Random Walk</td><td>AugRWalk</td><td>(D + I) -1 (A + I)</td></tr></table>

The hyper-parameters of best accuracy for each backbone on all datasets.

Table content: <table><tr><td>Dataset</td><td>Backbone</td><td colspan="2">nlayers Acc.</td><td>Hyper-parameters</td></tr><tr><td/><td>GCN</td><td>4</td><td colspan="2">0.876 lr:0.010, weight-decay:5e-3, sampling-percent:0.7, dropout:0.8, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:FirstOrderGCN</td></tr><tr><td/><td>ResGCN</td><td>4</td><td colspan="2">0.87 lr:0.001, weight-decay:1e-5, sampling-percent:0.1, dropout:0.5, nor-</td></tr><tr><td>Cora</td><td>JKNet</td><td>16</td><td colspan="2">malization:FirstOrderGCN 0.88 lr:0.008, weight-decay:5e-4, sampling-percent:0.2, dropout:0.8, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:AugNormAdj</td></tr><tr><td/><td>IncepGCN</td><td>8</td><td colspan="2">0.882 lr:0.010, weight-decay:1e-3, sampling-percent:0.05, dropout:0.5,</td></tr><tr><td/><td/><td/><td/><td>normalization:AugNormAdj</td></tr><tr><td/><td>GraphSage</td><td>4</td><td colspan="2">0.881 lr:0.010, weight-decay:5e-4, sampling-percent:0.4, dropout:0.5, ag-</td></tr><tr><td/><td/><td/><td/><td>gregator:mean</td></tr><tr><td/><td>GCN</td><td>4</td><td colspan="2">0.792 lr:0.009, weight-decay:1e-3, sampling-percent:0.05, dropout:0.8,</td></tr><tr><td/><td/><td/><td/><td>normalization:BingGeNormAdj, withloop, withbn</td></tr><tr><td/><td>ResGCN</td><td>16</td><td colspan="2">0.794 lr:0.001, weight-decay:5e-3, sampling-percent:0.5, dropout:0.3, nor-</td></tr><tr><td>Citeseer</td><td>JKNet</td><td>8</td><td colspan="2">malization:BingGeNormAdj, withloop 0.802 lr:0.004, weight-decay:5e-5, sampling-percent:0.6, dropout:0.3, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:AugNormAdj, withloop</td></tr><tr><td/><td>IncepGCN</td><td>8</td><td colspan="2">0.805 lr:0.002, weight-decay:5e-3, sampling-percent:0.2, dropout:0.5, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:BingGeNormAdj, withloop</td></tr><tr><td/><td>GraphSage</td><td>2</td><td colspan="2">0.8 lr:0.001, weight-decay:1e-4, sampling-percent:0.1, dropout:0.5, ag-</td></tr><tr><td/><td/><td/><td/><td>gregator:mean</td></tr><tr><td/><td>GCN</td><td>4</td><td colspan="2">0.913 lr:0.010, weight-decay:1e-3, sampling-percent:0.3, dropout:0.5, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:BingGeNormAdj, withloop, withbn</td></tr><tr><td/><td>ResGCN</td><td>32</td><td colspan="2">0.911 lr:0.003, weight-decay:5e-5, sampling-percent:0.7, dropout:0.8, nor-</td></tr><tr><td>Pubmed</td><td>JKNet</td><td>64</td><td colspan="2">malization:AugNormAdj, withloop, withbn 0.916 lr:0.005, weight-decay:1e-4, sampling-percent:0.5, dropout:0.8, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:AugNormAdj, withloop,withbn</td></tr><tr><td/><td>IncepGCN</td><td>4</td><td colspan="2">0.916 lr:0.002, weight-decay:1e-5, sampling-percent:0.5, dropout:0.8, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:BingGeNormAdj, withloop, withbn</td></tr><tr><td/><td>GraphSage</td><td>8</td><td colspan="2">0.917 lr:0.007, weight-decay:1e-4, sampling-percent:0.8, dropout:0.3, ag-</td></tr><tr><td/><td/><td/><td/><td>gregator:mean</td></tr><tr><td/><td>GCN</td><td colspan="3">4 0.9671 lr:0.005, weight-decay:1e-4, sampling-percent:0.6, dropout:0.5, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:AugRWalk, withloop</td></tr><tr><td/><td>ResGCN</td><td colspan="3">16 0.9648 lr:0.009, weight-decay:1e-5, sampling-percent:0.2, dropout:0.5, nor-</td></tr><tr><td>Reddit</td><td>JKNet</td><td colspan="3">malization:BingGeNormAdj, withbn 8 0.9702 lr:0.010, weight-decay:5e-5, sampling-percent:0.6, dropout:0.5, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:BingGeNormAdj, withloop,withbn</td></tr><tr><td/><td>IncepGCN</td><td colspan="3">8 0.9687 lr:0.008, weight-decay:1e-4, sampling-percent:0.4, dropout:0.5, nor-</td></tr><tr><td/><td/><td/><td/><td>malization:FirstOrderGCN, withbn</td></tr><tr><td/><td>GraphSAGE</td><td colspan="3">4 0.9654 lr:0.005, weight-decay:5e-5, sampling-percent:0.2, dropout:0.3, ag-</td></tr><tr><td/><td/><td/><td/><td>gregator:mean</td></tr></table>

Accuracy (%) comparisons on different backbones with and without DropEdge

Table content: <table><tr><td>64</td><td>DropEdge</td><td>53.20</td><td>84.80</td><td>87.90</td><td>88.20</td><td>31.90</td><td>45.60</td><td>75.30</td><td>80.00</td><td>79.90</td><td>25.10</td><td>79.00</td><td>90.20</td><td>91.60</td><td>90.00</td><td>62.30</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>32</td><td>DropEdge Orignal</td><td>74.60 52.00</td><td>86.80 79.80</td><td>87.60 86.30</td><td>87.70 85.30</td><td>32.20 31.90</td><td>61.40 44.60</td><td>77.90 21.20</td><td>80.00 76.70</td><td>80.30 79.00</td><td>53.60 16.90</td><td>86.20 79.70</td><td>91.10 87.90</td><td>91.30 90.60</td><td>90.50 OOM</td><td>47.90 40.70</td><td>50.51 -</td><td>94.27 -</td><td>OOM -</td><td>OOM -</td><td>96.47 -</td></tr><tr><td>16</td><td>DropEdge Orignal</td><td>84.30 71.60</td><td>86.90 85.10</td><td>88.00 87.10</td><td>87.70 87.40</td><td>84.50 31.90</td><td>76.80 59.20</td><td>79.40 74.40</td><td>80.10 71.70</td><td>80.20 72.60</td><td>74.50 37.00</td><td>90.30 84.60</td><td>91.00 90.20</td><td>91.50 89.20</td><td>91.30 OOM</td><td>87.80 41.30</td><td>90.54 45.55</td><td>96.48 93.93</td><td>96.78 OOM</td><td>OOM OOM</td><td>96.18 96.43</td></tr><tr><td>8</td><td>DropEdge Orignal</td><td>85.80 82.10</td><td>86.90 85.30</td><td>87.80 86.20</td><td>88.20 87.10</td><td>87.10 84.10</td><td>77.20 65.20</td><td>78.80 78.20</td><td>80.20 78.80</td><td>80.50 78.50</td><td>77.10 72.90</td><td>90.90 88.10</td><td>90.50 89.60</td><td>91.20 89.90</td><td>91.50 90.80</td><td>91.70 83.50</td><td>96.48 67.11</td><td>96.46 96.34</td><td>97.02 OOM</td><td>96.87 OOM</td><td>96.42 96.15</td></tr><tr><td>4</td><td>DropEdge Orignal</td><td>87.60 78.70</td><td>87.00 85.40</td><td>87.70 86.70</td><td>87.90 86.70</td><td>88.10 84.30</td><td>79.20 74.60</td><td>78.80 77.80</td><td>80.20 79.20</td><td>79.90 79.60</td><td>79.20 74.10</td><td>91.30 90.10</td><td>90.70 89.60</td><td>91.30 90.60</td><td>91.60 90.20</td><td>91.20 90.20</td><td>96.71 96.17</td><td>96.33 96.37</td><td>96.75 96.82</td><td>96.77 96.43</td><td>96.54 96.38</td></tr><tr><td>2</td><td>DropEdge Orignal</td><td>86.50 85.50</td><td>-86.00</td><td>-86.90</td><td>-85.60</td><td>88.10 87.10</td><td>78.70 76.70</td><td>-78.90</td><td>-79.10</td><td>-79.50</td><td>80.00 77.30</td><td>91.20 88.70</td><td>-90.70</td><td>-90.50</td><td>-89.90</td><td>90.70 89.40</td><td>96.13 96.62</td><td>-96.13</td><td>-96.54</td><td>-96.48</td><td>96.28 96.45</td></tr><tr><td/><td>Orignal</td><td>86.10</td><td>-</td><td>-</td><td>-</td><td>87.80</td><td>75.90</td><td>-</td><td>-</td><td>-</td><td>78.40</td><td>90.20</td><td>-</td><td>-</td><td>-</td><td>90.10</td><td>96.11</td><td>-</td><td>-</td><td>-</td><td>96.22</td></tr><tr><td/><td>Backbone</td><td>GCN</td><td>ResGCN</td><td>JKNet</td><td>IncepGCN</td><td>GraphSAGE</td><td>GCN</td><td>ResGCN</td><td>JKNet</td><td>IncepGCN</td><td>GraphSAGE</td><td>GCN</td><td>ResGCN</td><td>JKNet</td><td>IncepGCN</td><td>GraphSAGE</td><td>GCN</td><td>ResGCN</td><td>JKNet</td><td>IncepGCN</td><td>GraphSAGE</td></tr><tr><td/><td>Dataset</td><td/><td/><td colspan="2">Cora</td><td/><td/><td/><td colspan="2">Citeseer</td><td/><td/><td/><td colspan="2">Pubmed</td><td/><td/><td/><td colspan="2">Reddit</td><td/></tr></table>

Peer review (iclr2020):
This paper studied the problem of *deep* GCNs where the goal is to develop training methods that can make GCN becomes deeper while maintaining good test accuracy. The authors proposed a new method called *DropEdge*, where they randomly drop out the edges of the input graphs and demonstrate in experiments that this technique can indeed boost up the testing accuracy of deep GCN compared to other baselines. This paper is clearly well-written and the authors conducted a comprehensive study on deep GCNs. I also like the discussion in sec 4.3 where the authors explicitly clarify what are the difference between DropEdge, Dropout and DropNode, as the other two are the methods that will pop up during reading this paper. The extensive experiment results also show that for deeper GCNs, DropEdge always win over other baselines (see Tab 1) despite most of them are marginal except the backbone being GraphSAGE on Citeseer. Can you explain why this is the case? Why other backbones seem to have similar performance even with DropEdge (i.e. most of the accuracy increase are less than 3 %). Question: 1. When looking at Tab 1, it looks like most of the time, 2-layers networks are already the best (or close to the best) and are clearly better than 32 layers. Therefore, this makes me wonder: why do we need deeper networks at all if the shallow networks can already achieve a good (almost best) performance and it is also much similar and efficient in training? Can you please clarify why do we care to train a deeper network at all under this scenario? Are there any reasons that we would like to use deeper network as opposed to shallower networks? 2. It is less clear to me regarding this sentence: *DropEdge either retards the convergence speed of over-smoothing or relieves the information loss caused by it* Overall, I think this paper presents an interesting study on making deeper GCNs comparable to shallow network performance, but since the the boosted performance doesn*t really outperform most of the 2-layer networks, I would like to hear the justification of why we need the deeper networks for this node classification task.