L INK PREDICTION IN H YPERGRAPHS USING G RAPH CONVOLUTIONAL NETWORKS
Anonymous authors
Paper under double-blind review
A BSTRACT
Link prediction in simple graphs is a fundamental problem in which new links between nodes are predicted based on the observed structure of the graph.
However, in many real-world applications, there is a need to model relationships among nodes which go beyond pairwise associations.
For example, in a chemical reaction, relationship among the reactants and products is inherently higher-order.
Additionally, there is need to represent the direction from reactants to products.
Hypergraphs provide a natural way to represent such complex higher-order relationships.
Even though Graph Convolutional Networks (GCN) have recently emerged as a powerful deep learning-based approach for link prediction over simple graphs, their suitability for link prediction in hypergraphs is unexplored – we ﬁll this gap in this paper and propose Neural Hyperlink Predictor (NHP).
NHP adapts GCNs for link prediction in hypergraphs.
We propose two variants of NHP – NHP-U and NHP-D – for link prediction over undirected and directed hypergraphs, respectively.
To the best of our knowledge, NHP-D is the ﬁrst method for link prediction over directed hypergraphs.
Through extensive experiments on multiple real-world datasets, we show NHP’s effectiveness.
I NTRODUCTION
The problem of link prediction in graphs has numerous applications in the ﬁelds of social network analysis (Liben-Nowell & Kleinberg, 2003), knowledge bases (Nickel et al., 2016), bioinformatics (Lü & Zhou, 2011) to name a few.
However, in many real-world problems relationships go beyond pairwise associations.
For example, in chemical reactions data the relationship representing a group of chemical compounds that can react is inherently higher-order and similarly, the co-authorship relationship in a citation network is higher-order etc.
Hypergraphs provide a natural way to model such higher-order complex relations.
Hyperlink prediction is the problem of predicting such missing higher-order relationships in a hypergraph.
Besides the higher-order relationships, modeling the direction information between these relationships is also useful in many practical applications.
For example, in the chemical reactions data, in addition to predicting groups of chemical compounds which form reactants and/or products, it is also important to predict the direction between reactants and products, i.e., a group of reactants react to give a group of products.
Directed hypergraphs (Gallo et al., 1993) provide a way to model the direction information in hypergraphs.
Similar to the undirected hypergraphs, predicting the missing hyperlinks in a directed hypergraph is also useful in practical settings.
Figure 1 illustrates the difference between modeling the chemical reactions data using undirected and directed hypergraphs.
Most of the previous work on hyperlink prediction (Zhou et al., 2006; Zhang et al., 2018) focus only on undirected hypergraphs.
In this work we focus both on undirected and directed hypergraphs.
Recently, Graph Convolutional Networks (GCNs) (Kipf & Welling, 2017) have emerged as a powerful tool for representation learning on graphs.
GCNs have also been successfully applied for link prediction on normal graphs (Zhang & Chen, 2018; van den Berg et al., 2018; Schlichtkrull et al., 2018; Kipf & Welling, 2016).
Inspired by the success of GCNs for link prediction in graphs and deep learning in general Wang et al. (2017), we propose a GCN-based framework for hyperlink prediction which works for both undirected and directed hypergraphs.
We make the following contributions:
• We propose a Graph Convolutional Networks (GCN)-based framework called Neural Hyperlink Predictor (NHP) for the problem of hyperlink prediction.
To the best of our knowledge, this is the ﬁrst ever deep learning based approach for this problem.
• We extend the proposed NHP for the problem of hyperlink prediction in directed hypergraphs.
To the best of our knowledge, this is the ﬁrst ever attempt at the problem of link prediction in directed hypergraphs.
• Through extensive experiments on multiple real-world datasets, we show the effectiveness of proposed NHP for link prediction in both undirected and directed hypergraphs.
We have released NHP’s source code at this anonymous location: https://anonymous.4open.
science/repository/7d86231e-f6ba-4795-ae51-ac28d89f1521/ .
R ELATED WORK
In this section, we brieﬂy review related work in deep learning on graphs and link prediction on hypergraphs.
Learning representations on graphs : The key advancements in learning low-dimensional node representations in graphs include matrix factorisation-based methods, random-walk based algorithms, and deep learning on graphs (Hamilton et al., 2017).
Our work is based on deep learning on graphs.
Geometric deep learning (Bronstein et al., 2017) is an umbrella phrase for emerging techniques attempting to generalise (structured) deep neural network models to non-Euclidean domains such as graphs and manifolds.
The earliest attempts to generalise neural networks to graphs embed each node in an Euclidean space with a recurrent neural network (RNN) and use those embeddings as features for classiﬁcation or regression of nodes or graphs (Gori et al., 2005; Scarselli et al., 2009).
A CNN-like deep neural neural network on graphs was later formulated in the spectral domain in a pioneering work (Bruna et al., 2014) by a mathematically sound deﬁnition of convolution on graph employing the analogy between the classical Fourier transforms and projections onto the eigen basis of the graph Laplacian operator (Hammond et al., 2011).
Initial works proposed to learn smooth spectral multipliers of the graph Laplacian, although at high computational cost (Bruna et al., 2014;
Henaff et al., 2015).
To resolve the computational bottleneck and avoid the expensive computation of eigenvectors, the ChebNet framework (Defferrard et al., 2016) learns Chebyshev polynomials of the graph Laplacian
(hence the name ChebNet).
The graph convolutional network (GCN) (Kipf &
Welling, 2017) is a simpliﬁed ChebNet framework that uses simple ﬁlters operating on 1-hop local neighborhoods of the graph.
A second formulation of convolution on graph is in the spatial domain (or equivalently in the vertex domain) where the localisation property is provided by construction.
One of the ﬁrst formulations of a spatial CNN-like neural network on graph generalised standard molecular feature extraction methods based on circular ﬁngerprints (Duvenaud et al., 2015).
Subsequently, all of the above types
(RNN, spectral CNN, spatial CNN on graph) were uniﬁed into a single message passing neural network (MPNN) framework (Gilmer et al., 2017) and a variant of MPNN has been shown to achieve state-of-the-art results on an important molecular property prediction benchmark.
The reader is referred to a comprehensive literature review (Bronstein et al., 2017) and a survey (Hamilton et al., 2017) on the topic of deep learning on graphs and learning representation on graphs respectively.
Below, we give an overview of related research in link prediction on hypergraphs where relationships go beyond pairwise.
Link Prediction on hypergraphs : Machine learning on hypergraphs was introduced in a seminal work (Zhou et al., 2006) that generalised the powerful methodology of spectral clustering to hypergraphs and further inspired algorithms for hypergraph embedding and semi-supervised classiﬁcation of hypernodes.
Link prediction on hypergraph (hyperlink prediction) has been especially popular for social networks to predict higher-order links such as a user releases a tweet containing a hashtag (Li et al., 2013) and to predict metadata information such as tags, groups, labels, users for entities (images from Flickr)
(Arya & Worring, 2018).
Techniques for hyperlink prediction on social networks include ranking for link proximity information (Li et al., 2013) and matrix completion on the (incomplete) incidence matrix of the hypergraph (Arya & Worring, 2018;
Monti et al., 2017).
Hyperlink prediction has also been helpful to predict multi-actor collaborations (Sharma et al., 2014).
In other works, a dual hypergraph has been constructed from the intial (primal) hypergraph to cast the hyperlink prediction as an instance of vertex classiﬁcation problem on the dual hypergraph (Lugo-Martinez & Radivojac, 2017).
Coordinated matrix maximisation (CMM) predicts hyperlinks in the adjacency space with non-negative matrix factorisation and least square matching performed alternately in the vertex adjacency space (Zhang et al., 2018).
CMM uses expectation maximisation algorithm for optimisation for hyperlink prediction tasks such as predicting missing reactions of organisms’ metabolic networks.
P RELIMINARIES
Undirected hypergraph is an ordered pair H =
( V, E ) where V = { v 1 , · · · , v n } is a set of n hypernodes and E = { e 1 , · · · , e m } ⊆ 2 V is a set of m hyperlinks.
The problem of link prediction in an incomplete undirected hypergraph H involves predicting missing hyperlinks from ¯ E = 2 V − E based on the current set of observed hyperlinks E .
The number of hypernodes in any given hyperlink e ∈ E can be any integer between 1 and 2 n .
This variable cardinality of a hyperlink makes traditional graph-based link prediction methods infeasible because they are based on exactly two input features (those of the two nodes potentially forming a link).
The variable cardinality problem also results in an exponentially large inference space because the total number of potential hyperlinks is O (2 n ) .
However, in practical cases, there is no need to consider all the hyperlinks in ¯ E as most of them can be easily ﬁltered out (Zhang et al., 2018).
For example, for the task of ﬁnding missing metabolic reactions, we can restrict hyperlink prediction to all feasible reactions because the infeasible reactions seldom have biological meanings.
In other cases such as predicting multi-author collaborations of academic/technical papers, hyperlinks have cardinalities less than a small number, as papers seldom have more than 6 authors.
The number of restricted hyperlinks in such practical cases is not exponential and hence hyperlink prediction on the restricted set of hyperlinks becomes a feasible problem.
Formally, a hyperlink prediction problem (Zhang et al., 2018) is a tuple ( H, E ) , where H =
( V, E ) is a given incomplete hypergraph and E is a set of (restricted) candidate hyperlinks with E ⊆ E .
The problem is to ﬁnd the most likely hyperlinks missing in H from the set of hyperlinks E − E .
Directed hypergraph (Gallo et al., 1993) is an ordered pair H
=
( V, E ) where V = { v 1 , · · · , v n } is a set of n hypernodes and E = { ( t 1 , h 1 ) , · · · , ( t m , h m ) } ⊆
2 V is a set of m directed hyperlinks.
Each e ∈ E is denoted by ( t, h ) where t ⊆
V is the tail and
h ⊆ V is the head with t (cid:54) =
Φ , h (cid:54) = Φ .
As shown in ﬁgure 1, chemical reactions can be modeled by directed hyperlinks with chemical substances forming the set V .
Observe that this model captures and is general enough to subsume previous graph models:
• an undirected hyperlink is the special case when t = h
• a directed simple link (edge) is the special case when | t | = | h | =
Similar to the undirected case, the directed hyperlink prediction problem is a tuple ( H, E ) , where
H =
( V, E ) is a given incomplete directed hypergraph and E is a set of candidate hyperlinks with
E ⊆ E .
The problem is to ﬁnd the most likely hyperlinks missing in H from the set of hyperlinks
E − E .
P ROPOSED F RAMEWORK NHP
In this section we discuss the proposed approach for link prediction in hypergraphs.
Our proposed
NHP can predict both undirected and directed hyperlinks in a hypergraph.
We start with how NHP can be used to predict undirected hyperlinks and then explain how NHP can be extended to the directed case.
4. NHP-U: L INK P REDICTION IN U NDIRECTED H YPERGRAPHS
Given an undirected hypergraph H =
( V, E ) , as a ﬁrst step NHP constructs a dual hypergraph
H ∗ =
( V ∗ , E ∗ ) of H deﬁned as follows.
Hypergraph Duality (Scheinerman & Ullman, 2011)
The dual hypergraph of a hypergraph H =
( V, E ) with V
= { v 1 , . .
. , v n } and E = { e 1 , . . .
, e m } , denoted by H ∗ =
( V ∗ , E ∗ ) is obtained by taking V ∗ =
E as the set of hypernodes and E ∗ =
{ e ∗ 1 , . . .
, e ∗ n } such that e ∗ i = { e ∈ E : v i ∈
e } with e ∗ i corresponding to v i for i = 1 , . . . , n .
The vertex-hyperedge incidence matrices of H and H ∗ are transposes of each other.
The problem of link prediction in H can be posed as a binary node classiﬁcation problem in H ∗ (Lugo-Martinez & Radivojac, 2017).
A label +1 on a node in H ∗ indicates the presence of a hyperlink in H and a label of − 1 indicates the absence.
For the problem of semi-supervised node classiﬁcation on the dual hypergraph H ∗ , we use Graph Convolutional Networks (GCN) on the graph obtained from the clique expansion of H ∗ .
Clique expansion (Zhou et al., 2006) is a standard and a simple way of approximating a hypergraph into a planar graph by replacing every hyperlink of size s with an s -clique
(Feng et al., 2018).
Graph Convolutional Network (Kipf & Welling, 2017) Let G =
( V , E ) , with N = |V| , be a simple undirected graph with the adjacency matrix A ∈ R N × N , and let the data matrix be X ∈ R N × p .
The data matrix has p -dimensional real-valued vector representations for each node in the graph.
The forward model for a simple two-layer GCN takes the following simple form:
[Equation]
where ¯ A = ˜ D − 12 ˜
A ˜ D − 12 , ˜
A =
A + I, and ˜ D ii = (cid:80) Nj =
˜
A ij .
Θ (0) ∈
R p
×
h is an input-to-hidden weight matrix for a hidden layer with h hidden units and Θ (1)
∈ R h × r is a hidden-to-output weight matrix.
The softmax activation function is deﬁned as softmax ( x i ) =
exp ( x i ) (cid:80) j exp ( x j ) and applied row-wise.
For semi-supervised multi-class classiﬁcation with q classes, we minimise the cross-entropy error over the set of labeled examples, V L .
[Equation]
The weights of the graph convolutional network, viz.
Θ (0) and Θ (1) , are trained using gradient descent.
We note that Lugo-Martinez & Radivojac (2017) also follow a similar approach of constructing dual graph followed by node classiﬁcation for hyperlink prediction.
However, ours is a deep learning based approach and (Lugo-Martinez & Radivojac, 2017) do not perform link prediction in the directed hypergraphs. the more details section
The cross-entropy objective of GCN in equation 2 inherently assumes the presence of labeled data from at least two different classes and hence cannot be directly used for the positive unlabeled setting.
A challenging variant of semi-supervised learning is positive-unlabelled learning (Elkan & Noto, 2008) which arises in many real-world applications such as text classiﬁcation, recommender systems, etc.
In this setting, only a limited amount of positive examples are available and the rest are unlabelled examples.
In positive unlabelled learning framework, we construct a plausible negative sampling set based on data similarity (Han & Shen, 2016).
We calculate the average similarity of each unlabelled point u ∈ E − E , to the positive labelled examples in E :
[Equation]
where S represents a similarity function and f represents a map that maps an example to a ddimensional emebedding space f : E → R d .
We then rank the unlabelled training examples in E − E in the ascending order of their average similarity scores.
We then select the top ones in the ascending order (i.e. with the lowest similarity values) to construct the set of plausible negative examples
F ⊆ E − E .
The intuition here is that the set of plausible negative examples contain examples most dissimilar to the positive examples.
The GCN on the hypergraph can subsequently be run by minimising the objective 2 over the positive examples in E and the plausible negative examples in F i.e. V L = E ∪ F .
4. NHP-D: L INK PREDICTION IN DIRECTED H YPERGRAPHS
As explained in Section 3, hyperlink prediction in directed hypergraphs is the problem of predicting directing links which are tail-head set pairs.
However, in practice the collection of the tail and head sets is also incomplete.
Therefore, the problem of link prediction in directed hypergraphs also requires to predict missing tail and head sets of nodes besides predicting directed links among them.
The tail and head sets of nodes can be thought of as undirected hyperlinks and the directed hyperlink is between these undirected hyperlinks.
A straight forward approach for this problem would be to predict the undirected hyperlinks ﬁrst and then followed by predicting the directed links between pairs.
However, this sequential approach may not produce desirable results as the error during the training for directed part would not have any impact on the training for the undirected part.
Therefore, we propose the following joint learning scheme,
[Equation]
[Table] [Table]
We denote W + L = { ( t, h ) ∈ E } to be a set of positively labelled directed pairs in loss L d .
In other words d ij 1 = 1 and d ij 0 = 0 for all ( i, j ) ∈ W + L .
The tail hyperlink and the corresponfding head hyperlink are separate hypernodes in the dual and form directed hyperlinks in the primal (with the direction from t to h ).
The set W + L consists of directed hyperlinks that currently exist in the given directed hypergraph.
Note that, for loss L u , the set of positively labelled hypernodes will be,
[Equation]
We sample |W + L | = | E | hypernodes (in the dual) from the unlabelled data using the positive unlabelled approach of 3 to get the set of W − L pairs.
We label these pairs negative i.e. d ij 1 =and d ij 0 = 1 for all ( i, j ) ∈ W − L .
We construct V − L =
(cid:83) ( t,h ) ∈ W − L t ∪ h similarly as in 5.
The sets
V L = V + L ∪ V − L and W L = W + L ∪ W − L are used to minimise the objective 6.
To explain how D is computed, we rewrite equation 1 as:
[Equation]
We use D = g ( x 1 , x 2 ) with g being a function that takes the dual hypernode representations x 1 ∈ X and x 2 ∈ X and is parameterised by for example a simple neural network.
In the experiments, we used a simple 2 -layer multilayer perceptron on the concatenated embeddings x 1 || x 2 i.e.
[Equation]
We train Θ (0) , Θ (1) , and Θ end-to-end using backpropagation.
E XPERIMENTS FOR UNDIRECTED HYPERLINK PREDICTION
In this section, we evaluate NHP on hyperlink prediction in undirected hypergraphs.
We performed two different sets of experiments whose motivations and setups are as follows.
Predicting reactions of metabolic networks: Reconstructed metabolic networks are important tools for understanding the metabolic basis of human diseases, increasing the yield of biologically engineered systems, and discovering novel drug targets.
We used four datasets of (Zhang et al., 2018) and we show the statistics of the datasets used in table 1.
For each dataset, we randomly generated fake reactions according to the substance distribution of the existing reactions.
So, the candidate reactions contain already existing ones and the randomly generated fake ones.
The number of fake reactions generated is equal to the number of already existing ones.
Given a small number of reactions the task is to predict the other reactions.
[Table] [Table]
hyperlinks.
NHP achieves consistently superior performance over its baselines for all the datasets.
Refer to section 5 for more details. Predicting multi-author collaborations in coauthorship networks: Research collaborations in scientiﬁc community have been extensively studied to understand team dynamics in social networks (Newman, 2001; Barabási et al., 2002; Bai et al., 2018; Chen et al., 2015).
Coauthorship data provide a means to analyse research collaborations.
We used cora and dblp for coauthorship data.
The statistics are shown in table 2 and the construction of the datasets is pushed to the supplementary.
A coauthorship hypergraph (primal) contains each author as a hypernode and each paper represents a hyperlink connecting all the authors of the paper.
The corresponding dual hypergraph considers each author as a hyperlink connecting all the papers (hypernodes) coauthored by the author.
The hyperlink prediction problem is given a small number of collaborations, to essentially predict other collaborations among a given set of candidate collaborations.
5. E XPERIMENTAL SETUP
For each dataset we randomly sampled 10% of the hypernodes to get E and then we sampled an equal number of negatively-labelled hypernodes from E − E in the positive-unlabelled setting of 3.
To get the feature matrix X , we used random 32 − dimensional Gaussian features ( p = 32 in equation 1) for the metabolic networks and bag-of-word features shown in table 2 for the coauthorship datasets.
For fake papers we generated random Gaussian bag-of-word features.
We used node2vec (Grover & Leskovec, 2016) to learn low dimensional embedding mapping, f : E → R d with d = 128 .
We used the clique expansion of the dual hypergraph as input graph to node2vec and cosine similarity to compute similarity between two embeddings.
We compared NHP against the following state-of-the-art baselines for the same E as constructed above.
• Spectral hypergraph Clustering (SHC) (Zhou et al., 2006) :
SHC outputs classiﬁcation scores by f = ( I − ξ Θ) − 1 y .
We used SHC on the dual hypergraph.
• node2vec (Grover & Leskovec, 2016) :
One of the most popular node embedding approaches.
We note that (Grover & Leskovec, 2016) have shown node2vec to be superior to DeepWalk (Perozzi et al., 2014) and LINE (Tang et al., 2015) and hence we compared against only node2vec.
We used node2vec to embed the nodes of the clique expansion (of the dual).
We then used an MLP on the embeddings with the semi-supervised objective of equation 2 in the positive unlabelled setting of equation 3.
• Co-ordinated Matrix Maximisation (CMM) (Zhang et al., 2018) :
The matrix factorisation-based CMM technique uses the EM algorithm to determine the presence or absence of candidate hyperlinks.
• GCN on star expansion (Ying et al., 2018) :
PinSage (Ying et al., 2018) is a GCN-based method designed to work on the (web-scale) bipartite graph of Pinterest.
The Pinterest graph can be seen as the star expansion of a hypergraph (Agarwal et al., 2006) with pins
(hypernodes) on one side of the partition and boards (hyperlinks) on the other side. We have compared NHP against star expansion. This is essentially approximating the hypergraph with its star exapansion and then running a GCN over it (instead of the clique expansion of NHP).
5. E XPERIMENTAL RESULTS
Similar to (Zhang et al., 2018) [Table] we report mean AUC over 10 trials in table 3 and the mean number of hyperlinks recovered over 10 trials in the top ranked | ∆ E | ones in table 4.
Note that ∆ E ⊂ E is the set of missing hyperlinks with | ∆ E | = 0 .
9 ∗ | E | .
As we can observe, we consistently outperform the baselines in both the metrics.
We believe this is because of the powerful non-linear feature extraction capability of GCNs [Table]
E XPERIMENTS FOR DIRECTED HYPERLINK PREDICTION
[Table]
[Table] [Table]
We used the same four metabolic networks to construct directed hyperlinks.
The metabolic reactions are encoded by stoichiometric matrices.
The negative entries in a stoichiometric matrix indicate reactants and positive entries indicate products.
We extracted only those reactions which have at least two substances in each of the reactant side and the product side and the statistics are shown in table 5.
We labelled randomly sampled 10% of the hyperlinks in the data and use the remaining 90% unlabelled data for testing.
Tables 6 and 7 show the results on the datasets.
NHP-D (joint) is the model proposed in 4.2.
On the other hand, NHP-D (sequential) is the model which treats undirected hyperlink prediction and direction prediction separately.
NHP-D (sequential) ﬁrst runs a GCN on the clique expansion of the undirected hypergraph to get the node embeddings (without softmax) and then runs a multi-layer perceptron on the concatenated emebeddings to predict the directions [Table] We compared against the following baselines with a multi-layer perceptron (MLP): • node2vec
(Grover & Leskovec, 2016) + MLP : We used node2vec for the undirected hyperlink prediction part (as explained in section 5 and a 2 -layer perceptron to predict the direction between hyperlinks with the joint objective of equation 4. • Co-ordinated Matrix Maximisation (CMM) (Zhang et al., 2018) + MLP :
The matrix factorisation-based CMM technique uses the EM algorithm to determine the presence or absence of candidate hyperlinks for the following optimisation problem:
[Equation]
subject to λ i ∈ { 0 , 1 } , i = 1 ,
[Equation]
, m W ≥We used a 2 − layer perceptron in a sequential manner.
To get a representation for a hyperlink e , we took the sum of the representations of all the hypernodes v ∈ e from the matrix W above.
Note CMM works on the primal hypergraph.
• GCN on star expansion (Ying et al., 2018)
+ MLP : We used the star expansion to approximate the input hypergraph and then run a GCN on it (instead of the clique expansion of NHP).
A 2 -layer perceptron is used to predict the direction between hyperlinks with the joint objective of equation 4.
As we see in the tables, both NHP-D (joint) and NHP-D (sequential) perform similarly.
This can be attributed to the fact that training data to predict directions between hyperlinks is sparse and hence the learned hypernode representations of both the models are similar.
Please note that existing apporaches for link prediction on directed simple graphs cannot be trivially adopted for this problem because of the sparsity in the training data.
Comparison to baselines NHP-D outperforms the baselines on 3 out of 4 datasets.
The dataset iHN637 seems to be a very challenging dataset on which each model recovers less than half the number of missing hyperlinks.
C ONCLUSION AND FUTURE WORK
We have introduced NHP, a novel neural approach for hyperlink prediction in both undirected and directed hypergraphs.
To the best of our knowledge, this is the ﬁrst neural method for hyperlink prediction in undirected hypergraphs.
NHP is also the ﬁrst method for hyperlink prediction in directed hypergraphs.
Through extensive experiments on multiple real-world datasets, we have demonstrated NHP’s effectiveness over state-of-the art baselines.
Approaches that augment GCNs with attention (Veliˇckovi´c et al., 2018), self-training and co-training with random walks (Li et al., 2018), edge-feature learning in a dual-primal setup (Monti et al., 2018) have been recently proposed on graph-based semi-supervised learning tasks.
Our NHP framework provides the ﬂexibility to incorporate these approaches for more improved performance.
An interesting future direction is predicting hyperlinks in partial-order hypergraphs (Feng et al., 2018).
We leave extending NHP framework to inductive settings as part of future work.
A PPENDIX 7. H YPERPARAMETER DETAILS
We used the same hyperparameters as Kipf & Welling (2017) for the 2 -layer GCN model for all the datasets in all the experiments.
Additionally for the 2 -layer multi-layer perceptron used for the directed hyperlink prediction experiments, we used 16 hidden units with a dropout rate of 0 .
hyperparameter value number of hidden units number of hidden layers dropout rate 0 .
L2 regularisation 5 × 10 −learning rate 0 .non-linearity
ReLU
C OAUTHORSHIP DATASETS CONSTRUCTION DETAILS
Cora: We used the author data 1 to get the co-authorship hypergraph for cora 2 .
• DBLP:
We used the DBLP database v4 3 .
We ﬁltered out papers without abstracts, and processed each abstract by tokenizing it and removing stop-words removal.
Further, we ﬁltered out papers with one author only.
This left 540532 papers.
In order to ensure that the hypergraph formed would be sufﬁciently dense
we found the number of papers authored by each author and took the top 1000 authors as ‘selected authors’.
Then we ﬁltered out the papers that were not authored by at least three of the selected authors.
Finally, we were left with 1590 papers by 685 of the original 1000 selected authors.
To extract word features from each of these abstracts, we took all words appearing in
these abstracts
with a frequency
greater than 50. Each abstract was thus represented by a 602-dimensional bag-of-
words representation. For both datasets, we randomly sample | E | fake papers according to the author distribution of the existing non-fake papers ( 2708 and 1590 for CORA and DBLP respectively).
We randomly generated Gaussian p dimensional features for these fake papers ( 1433 and 602 for CORA and DBLP respectively).