D EEP I NTERACTION P ROCESSES FOR T IME -E VOLVING G RAPHS
Anonymous authors
Paper under double-blind review
A BSTRACT
Time-evolving graphs are ubiquitous such as online transactions on an e-commerce platform and user interactions on social networks.
While neural approaches have been proposed for graph modeling, most of them focus on static graphs.
In this paper we present a principled deep neural approach that models continuous timeevolving graphs at multiple time resolutions based on a temporal point process framework.
To model the dependency between latent dynamic representations of each node, we deﬁne a mixture of temporal cascades in which a node’s neural representation depends on not only this node’s previous representations but also the previous representations of related nodes that have interacted with this node.
We generalize LSTM on this temporal cascade mixture and introduce novel time gates to model time intervals between interactions.
Furthermore, we introduce a selection mechanism that gives important nodes large inﬂuence in both k − depth subgraphs of nodes in an interaction.
To capture temporal dependency at multiple time-resolutions, we stack our neural representations in several layers and fuse them based on attention.
Meanwhile, our method can process unseen nodes and their interactions.
Experimental results on interaction prediction and classiﬁcation tasks – including a real-world ﬁnancial application – illustrate the effectiveness of the time gate, the selection and fusion mechanisms of our approach, as well as its superior performance over the alternative approaches.
I NTRODUCTION
Representation learning over graph data has become a core machine learning task with a wide range of applications including e-commerce, ﬁnance, social networks, and bioinformatics.
Various neural graph representations such as (Perozzi et al., 2014; Grover & Leskovec, 2016; Wang et al., 2016;
Kipf & Welling, 2017; Defferrard et al., 2016; Scarselli et al., 2009;
Ying et al., 2018; Hamilton et al., 2017b; Monti et al., 2017; Den Berg et al., 2017) have been proposed to learn from static graph data and successfully used for downstream tasks ( e.g. , classiﬁcation).
Graph data, however, are often dynamic in practice; nodes and interactions between them can grow and shrink.
A straightforward approach to handle dynamic graphs is to compress them into one or several static graphs.
The drawbacks of this approach are multifold; we not only blur temporal structural information but also miss time information that can be critical for real-world applications.
An illustrative example is given in ﬁgure 1.
To handle continuous time-evolving graph, we can approximate a by a sequence of snapshot graphs, each of which includes all interactions that occur during a user-speciﬁed discrete-time interval, as shown in (Goyal et al., 2018; Leskovec et al., 2007;
Zhou et al., 2018; Sankar et al., 2019).
This treatment reduces time resolution and it is tricky to specify the appropriate aggregation granularity.
To avoid these problems, Nguyen et al.
(2018) proposed continuous-time dynamic networks (CTDNE) that generalize deep walk methods to learn time-dependent network embedding.
As a transductive method, CTDNE cannot handle the growth of new nodes.
Dai et al.
(2016) applied temporal point processes to model time-evolving graphs and, as a nonparametric Bayesian approach, their approach can naturally cope with the growth of new nodes and interactions.
They used recurrent neural networks (RNNs) to deﬁne an intensity function in temporal point processes.
These RNN models are shallow and one-step unrolled, making it easy to compute but relatively limited in modeling power.
Trivedi et al.
(2019) extended this approach by modeling two-time scale and adopting temporal-attention mechanism.
[Figure]
The rest of the paper is organized as follows.
In Section 2 we give background on temporal point processes and in Section 3 we present the new DIP approach.
In Section 4 we discuss related works.
In Section 5 we report experimental results on multiple interaction prediction and classiﬁcation tasks including an important real-world anti-fraud ﬁnancial application, demonstrating superior performance of the new approach over the alternatives.
TEMPORAL POINT PROCESSES
We ﬁrst describe temporal point processes (a class of nonparametric Bayesian models) that our approach is based on.
Speciﬁcally, a temporal point process is a stochastic process that generates a sequence of discrete events localized at times { t i } Ni =1 in any given observed time window [0 , T ] , where N is the number of events.
An important way to characterize temporal point processes is via the conditional intensity function λ ( t | H t )
-the stochastic model for the next event time t given all historical events before time t , denoted as H t = { t i | t i < t } .
Formally, within a small time window [ t, t + dt ) , λ ( t | H t ) dt is the probability for the occurrence for a new event given the H t : λ ( t | H t ) dt
= P { event in [ t, t + dt ) | H t } .
From the survival analysis theory(Aalen et al., 2008), given the times of the past events { t 1 , t 2 , . . .
, t i } , the conditional density that an event occurs at t i +1 is given as follows: p (cid:0) t i +1 | H t i +1 (cid:1) =
λ (cid:0) t i +1 | H t i +1 (cid:1) exp (cid:110) − (cid:82) t
i +
t
i λ ( t | H t ) dt (cid:111) ,where the exponential part in the above equation means the conditional probability that no event happens during [ t i , t i +1 ) .
The functional forms of the conditional intensity function λ ( t | H t ) can represent certain forms of dependencies of the historical events.
For instance, for Poisson processes(Kingman, 2005) we set λ to be constant – making the assumption that the process is stationary and the temporal events in history are independent of each other.
For classical Hawkes processes(Hawkes, 1971), the intensity function λ is often set to be a sum of multiple exponential functions, assuming that the mutual excitation among events is positive, additive over the past events, and exponentially decaying with time.
Mei &
Eisner (2017a) removed these limiting assumptions using LSTM to learn λ from data.
D EEP I NTERACTION PROCESSES
In this section, we present a new neural non-parametric Bayesian approach over continuous-time evolving graphs.
First, we present a temporal dependency graph that is a mixture of the temporal cascades, to model interdependence between graph nodes (as well as latent node representations).
Then we present a novel deep model to learn dynamic node representations in the temporal dependency graph.
This model naturally generalizes a chain-structured LSTM to a temporal Graph-structured
LSTM equipped with time gates to handle interaction with irregular time intervals.
Furthermore, we propose the FUSION and SELECTION methods to enhance the dynamic interactive nodes representation.
Given the dynamic node representations, we deﬁne deep interaction processes that model potential interactions between any two nodes over time and we ﬁnally layout the maximum likelihood estimation method to optimize it.
A toy example is shown in Figure 2 to introduce the corresponding concepts and the whole procedure of computing enhanced dynamic representation is given in Figure 3.
3. T EMPORAL DEPENDENCY GRAPH
Consider a collection of people-movie interactions at different time points ( e.g. , David saw the movie
The Skull at t 6 .) as shown in ﬁgure 2(a).
The people and movies form a temporal dependency graph in which each person or movie is a node and interactions happen over time (in ﬁgure 2(b)).
After one interaction occurs, we update the neural representation of the two nodes linked to this interaction; e.g. , right after time t 6 , we update the representations for David and movie The Skull.
The new neural representation of David depends on both his current and previous interactions – as a result, depending on the representations of the two nodes associated with the previous interaction.
This naturally forms a dependency cascade.
Similarly we can obtain a dependency cascade for Lucy’s representations.
Because of the common movies David and Lucy saw, their dependency cascades overlap and form a cascade mixture.
Formally, we denote a dynamic interaction or link at time t by l u,v,t where u and v are two nodes associated with this interaction.
We denote the node u at time t by u ( t ) and the two nodes associated with u ’s precedent interaction at time t − as u 1 ( t ) and u 2 ( t ) .
For later usage, we denote the k − depth temporal dependency subgraph for u ( t ) as subgraph ( u ( t ) , k ) as shown in ﬁgure 2(b).
3. DIP NEURAL UNIT AND D YNAMIC R EPRESENTATION
Now we present a novel neural unit to update dynamic latent representations of nodes over the temporal dependency graph which is illustrated in ﬁgure 2(c) and ﬁgure 2(d).
First, let us denote node u ’s features or embedding (i.e., a static representation jointly learned from data) at time t by x u ( t ) and denote features of interaction l by x l .
The interaction feature can be empty if the interaction contains only the temporal information.
The concatenation of x u ( t ) and x l is denoted by ˆ x u ( t ) .
Let ∆ ( u,t ) =
t − t − be the time interval between two consecutive interactions involving u at time t and t − .
Our neural unit, i.e, DIP unit, generalizes a chain-structured LSTM unit to depict the temporal dependency on graph data.
The DIP unit has an update gate s , an input gate z , an output gate o and two forget gates f over ˆ x u ( t ) , dynamic representation of u i ( t ) , i.e, h u i ( t ) and cell states c u
i ( t ) ( i = 1 , 2) as shown in ﬁgure 2(c).
Additionally, we introduce time gates g to capture the impacts of irregular time interval ∆ ( u,t ) .
Speciﬁcally, h u ( t ) and c u ( t ) are updated as follows:
[Equation]
where σ, tanh and (cid:12) represent the sigmoid function, the hyperbolic tangent function, and
the
Hadamard product (pointwise multiplication), respectively, and parameters in the unit including the recurrent weights
R z i , R
o i ,
R s i , R f i and R g i , the projection matrices W z , W o , W s , W f i and W g
i , the bias vectors b z , b o , b s , b f i and b g i and the time weight matrix M g i are learned from data.
For convenience, we use DIP-UNIT ( · ) to summarize the above equations, then the dynamic representation of node u ( t ) is given as follows:
[Equation]
where Θ represent all the parameters.
Considering the computational cost in practice, when obtaining nodes’ dynamic representation given an happened interaction, we only utilize history information in the subgraph ( u ( t ) , k ) as illustrated in ﬁgure 2(d), which is similar to a chain-structured LSTM training unfolded with the max k steps.
The k is a hyper-parameter.
3. E NHANCED D YNAMIC R EPRESENTATION 3.3. S TACKING AND F USION
To model nonlinear dependency relationships at different temporal resolutions, we stack L layers of
DIP-UNIT together.
The output of the j -th layer is computed recursively as follows:
[Equation]
were h 0 u ( t )
=
ˆ x u ( t ) , j = 1 , . . .
, L .
To train deeper dynamic neural networks easily, we employ the residual connection as the following form: skip ( h j − 1 u ( t ) , h ju ( t ) )
=
W skip h j − 1 u ( t ) + h ju ( t ) where W skip is a weight matrix.
Motivated by ELMo (Peters et al., 2018), we fuse all internal dynamic representations from all the layers to achieve rich dynamic representations.
The fusion is a weighted summation of all layers deﬁned as follows: h u ( t ) = γ (cid:80)
Lj =
α j h ju ( t ) , where α j are softmax-normalized weights and γ is a scaling parameter.
They both are learned parameters.
3.3. SELECTION
Given an interaction l u,v,t , it is reasonable to assume that not all the historical interactive nodes have the equal importance for formalizing this interaction.
Thus we use a two-phase gating mechanism to select relevant nodes to learn dynamic representations and cell states of the current node.
Speciﬁcally, a co-attention mechanism is ﬁrst used to measure relevance of historical time-evolving patterns between subgraph ( u ( t ) , k ) and subgraph ( v ( t ) , k ) , [Equation]
[Equation]
) where H ju = (cid:104) h j 1 , . .
, h ja , . .
, h jm (cid:105) , a ∈ subgraph ( u ( t ) , k ) , H jv
= (cid:104)
h
j 1 , . .
, h je , . . .
, h jn (cid:105) , e ∈ subgraph ( v ( t ) , k ) , m and n are the numbers of nodes in the two corresponding subgraphs , W Q ∈
R d × d are the weight parameters.
The Q j is a co-attention afﬁnity matrix which captures the relevance information in subgraph ( u ( t ) , k ) and subgraph ( v ( t ) , k ).
The co-dependent global embedding p ju , p jv are obtained by the following equations.
[Equation]
where Max means max-pooling operation which is used to choose the most relevant information for the maximum inﬂuence (or afﬁnities) on nodes in the corresponding subgraph.
In addition, to adjust the importance of historical nodes, two adaptive gate functions are designed for nodes in subgraph ( u ( t ) , k ) and subgraph ( v
[Equation] [Equation]
respectively,
[Equation]
where the weights w p , and w h are shared by all the stacked layers.
Using these gates, we enhance the dynamic node representations as follows:
[Equation]
Similarly, we can compute ( h jv ( t ) , c jv ( t ) ) for v ( t ) based on the selection mechanism.
3. CONDITIONAL INTENSITY FUNCTION
We model the dynamic interactions as a multi-dimensional temporal point process.
Speciﬁcally, we deﬁne the conditional intensity function of the temporal point process at the dimension indexed by ( u, v ) , given its graph-structured history H u,v t where H u,v t = subgraph ( u 1 ( t ) , k ) ∪ subgraph ( u 2 ( t ) , k ) ∪ subgraph ( v 1 ( t ) , k ) ∪ subgraph ( v 2 ( t ) , k ) , as follows:
[Equation]
(8) where h u,vt =
[Equation]
[Equation]
=
[Equation]
∆ (
[Equation] t ) , ∆ ( [Equation] t ) (cid:3) (cid:62) , the scalar b λ can be viewed as a base intensity level for the occurrence of the next interaction, and the SoftPlus function is used to ensure the non-negativity of the intensity.
A key step for obtaining H u,vt is to get the k -depth subgraphs of u and v ’s direct dependants.
Please see Appendix.
A for more details about fast obtaining k -depth subgraphs .
3. M AXIMUM LIKELIHOOD PARAMETER ESTIMATION 3.5. INTERACTION PREDICTION
Given a set of interactions as I = { ( u i , v i , t i ) } i = N i =1 observed in a time window [0, T], we can learn the model by minimizing the negative joint log-likelihood of I as follows:
L 1 = − (cid:80) i log P u
i ,v i (cid:0) t i |
H u i ,v i t i (cid:1) where P u i ,v i ( t i |
H u i ,v i t i ) represents the probability of formalizing an interaction between u i and v i at time t i given the dependant history of non-chain structures H u
i ,v i t i .
Based on the intensity deﬁnition, we have L 1 = − (cid:80) i log λ
u i ,v i (cid:0) t i |
H u i ,v i t i (cid:1)
+ (cid:82) T 0 Λ( t ) dt , where Λ( t ) =
(cid:80) u,v λ u,v ( t ) .
Since the survival part does not have an analytic solution, we apply Monte
Carlo to do numerical integrations.
We follow the negative sampling approaches used by Dai et al.
(2016) and Trivedi et al.
(2019) to accelerate the survival term calculation.
3.5. INTERACTION CLASSIFICATION
An interaction sequence with markers is denoted as I (cid:48) =
{ ( u i , v i , t i , y i ) } i = N i =1 , where y i is a marker at time t i and usually is a discrete variable.
In practice, the markers have different meanings in distinct scenes.
A marker can be treated as a magnitude in modeling earthquakes and aftershocks.
For ﬁnancial transactions, a marker can be used to label whether a transaction is a fraudulent trading or not.
The joint conditional density of an interaction( u i , v i , t i ) with marker y i is given as
P u i ,v i (cid:16) t i , y i
| ˆ H u i ,v i [Equation] (cid:17) .
By applying the Bayesian rule , the joint conditional density can be written as:
P u i ,v i
[Equation]
[Equation]
i
[Equation]
v
[Equation]
,where P u i ,v i ( t i | ˆ H u
i ,
[Equation]
i [Equation] i ) has the same meaning as given in subsection 3.5.1, while P ( y i | t i , ˆ H u
i ,v i t i ) means the distribution of y i given the interaction happened at t i with interaction history
ˆ H u i ,v i t i .
It should be noted that the history ˆ H u i ,v
i t i contains the information of history markers and one can design a marker-speciﬁc intensity function like Mei & Eisner (2017b).
For simplicity, in our marked temporal point processes, we assume that
P u i ,v i (cid:16) t i , y i
| ˆ H u i ,v i t i (cid:17) is independent of historical markers.
We model P (cid:16)
y i | t i , ˆ H u i ,v i t i (cid:17) as
[Equation]
where h u i ,v i ,t
i is the concatenation of h u ( t i ) and h v ( t i ) which can be regarded as a dynamic repre- sentation for an interaction between u and v at t i , V y
i is the weight parameters for the y
i th class.
Then the overall cost function is L 2 =
L 1 + L cross − entropy , where L cross − entropy is a cross-entropy loss over marks:
[Equation]
R ELATED W ORK
Inspired by the Skip-gram (Mikolov et al., 2013) for word embedding, a series of node embedding methods based on the random walks on graphs have been proposed(Perozzi et al., 2014;
Tang et al., 2015; Grover & Leskovec, 2016; Wang et al., 2016; 2017).
GCN and its variants (Bruna et al., 2013; Hamilton et al., 2017a; Kipf & Welling, 2017) are a recent class of algorithms which extend convolutions from spatial domains to graph-structured domains.
Meanwhile they can efﬁciently generate node embeddings for previously unseen data.
All models above are designed for static graphs.
The intuitive and popular approaches for modeling dynamic graphs are based on a sequence for graph snapshots(Goyal et al., 2018;
Zhou et al., 2018; Seo et al., 2018; Yu et al., 2018), but it can be difﬁcult to specify the appropriate aggregation granularity.
Nguyen et al.
(2018) adds a temporal constraint on random walk sampling, but it can’t model the rich temporal information explicitly.
Temporal point processes (TPPs) are an another alternative to model dynamics(Daley & Vere-Jones,
Several dynamic graph modeling methods based on the TPPs (Dai et al., 2016; Trivedi et al., 2019) have been proposed.
Our method DIP differs from these TPP-based methods by the extension of the LSTM model over temporal dependency graphs, the multiple time resolution modeling via stacking, fusing and the selection mechanism.
A recently proposed method Kumar et al.
(2019) models interactions directly by predicting the next interaction embedding.
More detailed related work are included in Appendix.
C.
E XPERIMENTS
We evaluate the proposed DIP model for interaction prediction and interaction classiﬁcation on several real-world datasets.
5. B ASELINES
GraphSage (Hamilton et al., 2017a) is an inductive graph neural network framework consisting of three different aggregators which are GCN , Mean and LSTM aggregators respectively.
We report the best result among these three aggregators noted as Graphsage*.
What’s more, for comparing with
GAT (Veli ˇ ckovi ´ c et al., 2017)
we also implement a graph attention aggregator based on GraphSage.
CTDNE (Nguyen et al., 2018) is a newly-proposed temporal network embedding method and also a tranductive method like DeepWalk(Perozzi et al., 2014).
It incorporates temporal order constraint when sampling walks from time-continuous graphs.
DynGEM (Goyal et al., 2018) takes a sequence of static graph snapshots as inputs to learn node embeddings by a deep auto-encoder network.
DeepCoevolve (Dai et al., 2016) models dynamic interaction sequences with two co-evolution recurrent neural networks.
Hidden embeddings are learned for interactive nodes after each interaction.
DyREP (Trivedi et al., 2019) uses a two-time scale deep temporal point process model to capture dynamics of graphs.
JODIE (Kumar et al., 2019) models interaction processes in a novel way by predicting the next interaction embedding directly instead of modeling the intensity function.
5. E XPERIMENTAL S ETTING
We conduct all the experiments with a hyper-parameter grid search strategy.
For all methods, we search the dimension of embedding from {16, 32, 64, 128} and the learning rate from {0.01, 0.001, 0.0005, 0.0001, 0.00001}.
For our DIP model, we go through {1, 2, 3, 4} for k and {1, 2, 3} for
L .
For Graphsage, the maximum number of 1-hop and 2-hop neighbor nodes are set to 25 and 20 respectively.
All the models are trained for at most 50 epochs with an early-stop operation if the performance does not improve for 5 epochs.
For Graphsage, DynGEM ,DeepCoevolve and
JODIE, we use the open source codes provided by the authors.
We implement the CTNDE and GAT based on the Graphsage framework, and implement DyREP based on the pytorch implementation of DeepCoevolve.
After the best conﬁguration is found, we repeat the full experiments 5 times and report the mean results and standard deviation.
5. I NTERACTION P REDICTION 5.3. D ATASETS
CollegeMsg (Leskovec & Krevl, 2014) consists of sending message interactions on an online social network at the University of California, Irvine during 193 days.
Ubuntu (Leskovec & Krevl, 2014) is a temporal interaction dataset extracted from the stack exchange website.
An interaction between two users means one answered another’s questions or replied to his/her posts.
Amazon (McAuley et al., 2015) is composed of commodity rating data from amazon users.
We use the Clothing subset of this dataset.
MathOverﬂow (Leskovec & Krevl, 2014) is comprised of interactions of commenting an existing answer on the Math Overﬂow website.
Table 1 shows the detailed dataset statistics.
In this table, Sparsity indicates | Edge | / ( | u ( t ) | ∗ | v ( t ) | ) .
Repetition is the rate of repeated interaction in datasets.
For each dataset,we ﬁrst sort these interactions by occurrence time and split them to be training/validation/test sets.
The cold-start participants which only exist in validation set or test set are removed.
[Table] R ESULTS AND A NALYSIS
For each interaction l u,v,t in test set, we ﬁx the ﬁrst participant u , then replace the second participant v by all possible participant candidates.
The conditional density p u,v ( t ) = λ u,v ( t ) S u,v ( t ) for all candidates are ﬁrst computed and then sorted by a descending order.
The rank of the correct participant is ﬁnally stored and denoted as rank i for the i-th test interaction.
We report the Mean Rank deﬁned as Mean Rank =
[Table]
Test Interaction .
Figure 4 summarizes the Mean Rank performances of all the baselines and our DIP method.
On the whole, our DIP method consistently beats all baselines by 65.84%, 41.64%, 10.69% and 43.99% over the four diverse datasets.
The CTDNE method performs worst across all the datasets since the generated embedding is static and can’t be updated and evolved across validation and test datasets. Meanwhile we can see that there is no consistent winner among baselines and all the methods perform relatively better on the CollegeMsg , Ubuntu and
MathOverﬂow datasets than do on the Amazon-Clothing dataset. It is also noteworthy that the static methods GAT and GraphSage* perform competitive with dynamic baselines on these three datasets. This phenomenon could be explained that the CollegeMsg , Ubuntu and MathOverﬂow have lots of repetitive interactions (at least 75% for CollegeMsg and MathOverﬂow on test data, 68% for ubuntu on test data as shown in Table 1) and help predict future recurring interaction. The two static methods perform worse than the other dynamic methods(but JODIE ) on the Amazon-Clothing because no repetitive information can be reused and nodes representation can’t be updated across validation and test like CTDNE .
As for comparison with dynamic methods, our work performs better than the Dyrep , DeepCoevolve and JODIE which use the similar mutually-recursive RNNs to capture and update co-evolution states sequentially. We argue that the vallina RNNs can’t capture long-term dependency history well and have optimization problems which may lead to worse performances.
5. I NTERACTION CLASSIFICATION
We conduct this task on an industrial dataset: Huabei Trade Data .
For the performance measures, we employ KS (Kolmogorov, 1933) value which is a big concern of the loans provider in anti-fraud detection, as well as AUC(Area under the ROC Curve) score.
This dataset consists of about 150, transaction records processed by Huabei during August 2018.
Each transaction is initiated with three parties: the buyer, the seller and transaction details such as merchant category and transaction amount.
Around 15% of the transaction are fraudulent and is labeled by a complicated Ex-Post method.
In order to conducting an Ex-ante detection, we are required to ﬁnd out fraudulent transaction at the time of initiating using those basic transaction features.
For each interaction event, there are 11 context features that can be obtained when the trade request is created, including information about buyer types, seller types, purchased items’ categories and trading platform.
We use the ﬁrst 10 days data as training set, the following 10 days data as validation set, the rest as test set.
Note that, in this scenario, there are always users who only appear in validation/testing dataset.
Thus, traditional transductive methods are not applicable on this task.
Alternatively, we employ the XGBoost (Chen & Guestrin, 2016) as an additional baseline which is a popular model in the cash-out detection task(Hu et al., 2019).
Table 2 compares the results.
Obviously, our model outperforms all the baseline methods again.
The contributions of different modules of the proposed DIP models are given in ﬁgure 6. [Figure] [Table] [Figure]
[Table]
[Figure]
[Table]
[Figure]
BLATION S TUDY
As we described in Section 3, the DIP model consists of three important components: First, it uses a
Time Gate in the DIP neural unit to explicitly model the temporal information.
Second, the selection mechanism enables our model to select more important historical information for interactions.
Third, the Fusion of multi-layer DIP-UNIT’s hidden state vector helps to extract high level feature.
We investigate the contribution of each component by disabling each of them one by one, and compare the corresponding result to the full model.
ﬁgure 5 and ﬁgure 6 give the detailed ablation results.
FullModel in the two ﬁgures means all the three components are enabled.
• No Time Gate:
In this conﬁguration, the time gate in DIP-UNIT is disabled.
This leads to a signiﬁcant drop of the Mean Rank performance.
It provides a strong evidence for the effectiveness of the time gate.
• No Selection:
In this conﬁguration the selection mechanism is disabled.
Accordingly, all the historical node representations contribute equally, thus again leading to a performance drop.
• No Fusion: In this variant, we directly use the hidden state vector of the last layer.
Again, the performance degrades signiﬁcantly.
This demonstrates that a fusion of different layers’ representations gives richer information than the last layer only.
C ONCLUSIONS
In this paper, we have proposed a deep multidimensional point process approach, DIP, to learn dynamic graph representations.
We generalize LSTM over temporal dependency graphs and model [Figure]
[Figure]
A TEMPORAL DEPENDENCY GRAPH
The key step for efﬁcient training or inference is to fast obtain subgraph ( u ( t ) , k ) or subgraph ( v ( t ) , k ) for an interaction ( u, v, t ) .
Obviously, our deﬁnition of temporal dependency graph( TDG ) in Subsection.3 provides a recursive and incremental way for its construction.
As illustrated in ﬁgure 2, we can easily ﬁnd that subgraph ( u ( t ) , k ) is the union of subgraph ( u ( t ) , k ), subgraph ( w ( t − , k ) ) and the newly added dependency edges.
Based on this feature, to fast obtain k − depth subgraph , we design two algorithms, TDG- COLORING and C ONSTRUCTION OF K - DEPTH TDG orderly.
A. TDG- COLORING
The TDG- COLORING algorithm is shown in Algorithm 1, which is the pre-step of Algorithm 2.
This algorithm takes I
= { ( u i , v i , t i ) } i =
N i
=1 as input and generate a new sorted sequence of interactions ordered by color numbers in an ascending way.
The TDG- COLORING algorithm ensures that, interactions in the same color group are independent and interactions from groups with smaller color numbers are precedents of the larger ones.
Algorithm
Require:
I : A sequence of interaction with a chronological order.
Ensure:
ColorGpSeq :
A sequence of sorted interactions by color no.
1: Initialize ColorGroupArray [ x ] ← − x nodes
I initial color no - all nodes 2: Initialize LastNodeTime [ x ] ← − 1 (cid:46) Record the latest time when node x was involved in an interaction and initialize with3: for event in I do 4: cur _ u, cur _ v, cur _ time = event.u, event.v, event.t 5: if cur _ time >
LastNodeTime [ cur _ u ] then 6:
ColorGroupArray [ cur _ u ] + + 7:
LastNodeTime [ cur _ u ] =
cur _ time 8: end if 9: if cur _ time >
LastNodeTime [ cur _ v ] then 10:
ColorGroupArray [ cur _ v ] + + 11:
LastNodeTime [ cur _ v ] = cur _ time 12: end if 13: event.gn = max( ColorGroupArray [ cur _ u ] , ColorGroupArray [ cur _ v ])
(cid:46) gn is short for group no 14: end for 15: ColorGpSeq =
Sort ( I ) (cid:46) sort
I by an ascending order with assigned group color no.
16: return ColorGpSeq
A. C ONSTRUCTION OF K - HOP TDG
The details of C ONSTRUCTION OF K - DEPTH TDG are given in Algorithm 2.
In Algorithm 2, based on the deﬁnition of TDG and ColorGpSeq output by Algorithm 1, we can incrementally construct the k-depth subgraphs for any two nodes in a new interaction.
Speciﬁcally, we ﬁrst call
P REVIOUS function to get each node’s adjacent interaction in which it was involved before current interaction.
Then incrementally construct a graph( ugraph ) rooted at cur _ u (lines 6 to 10).
Here
H k_subgraph [ u _ preInteraction ] stores the subgraph of ugraph which was obtained before ugraph because u _ preInteraction ranks ahead of current event in ColorGpSeq .
Likely, we can get vgraph incrementally.
At last, we call the Breadth-First-Search algorithm with the traverse depth =
K to obtain the k-depth subgraphs for cur _ event and then store it into H k_subgraph [ cur _ event ] .
A. O BTAINING AND U PDATING K - DEPTH TDG
As for efﬁciency, fast obtaining k-depth subgraph of TDG for an interaction is both important for ofﬂine training and online inference.
Based on Algorithm 2, for an interaction l from the
Require:
ColorGpSeq : A seq of sorted interaction by group color no.
Ensure:
HashTable H k_subgraph : Map an interaction to its corresponding k_subgraph .
1: Initialize an empty HashTable H k_subgraph 2: for cur _ event in ColorGpSeq do 3: cur _ u, cur _ v = cur _ event.u, cur _ event.v 4: u _ preInteraction
=
Previous [ cur _ u ] (cid:46)
Previous Function returns last interaction in which cur _ u was involved 5: v _ preInteraction =
Previous [ cur _ v ] 6: if u _ preInteraction exits then 7: edge 1 = ( u _ preInteraction.u, cur _ u ) (cid:46)
The edge is deﬁned according to the deﬁnition of dependency graph 8: edge 2 = ( u _ preInteraction.v, cur _ u ) 9:
ugraph = edge 1 ∪ edge10: ugraph =
H k_subgraph [ u _ preInteraction ] ∪ ugraph (cid:46)
Incremental Update 11: else 12: ugraph = cur _ u 13: end if 14: if v _ preInteraction exits then 15: edge 3 = ( v _ preInteraction.u, cur _ v ) 16: edge 4 = ( v _ preInteraction.v, cur _ v ) 17: vgraph = edge 3 ∪ edge18: vgraph =
H k_subgraph [ v _ preInteraction ] ∪ vgraph (cid:46)
Incremental Update 19: else 20: vgraph = cur _ v 21: end if 22:
G K − depth =
BFS ( ugraph ∪ vgraph, depth = K ) (cid:46) Call Breadth-First-Search with max depth =
K 23:
H k_subgraph [ cur _ event ] =
G K − depth 24: end for 25: return H k_subgraph collected training interaction data, we can obtain its corresponding k-depth subgraphs directly using
H k_subgraph [ l ] with time complexity O
For a new incoming interaction l with two nodes cur _ u and cur _ v , we can easily reuse Algorithm 2’s code from line 3 to 22 to ﬁnd ugraph and vgraph , respectively.
Then we merge them to get k-depth subgraph for doing inference online.
At the same time, we can incrementally update TDG by H k_subgraph [ l ].
The time complexity here for updating
TDG mainly depends on union between ugraph and vgraph in line 22 of Algorithm 2 and is O ( n + e ) where n and e are the total number of nodes and edges for ugraph and vgraph .
B A DDITIONAL D ATA STATISTICS A ND E XPERIMENT RESULTS B.
[Table]
[Equation] [Table] [Figure]
[Figure] [Figure]
B. H IT P ERFORMANCE
Figure 7 and Figure 8 provide HIT@5 and Hit@1 results in addition to the Mean Rank results in
Section 5.3.2 of the main paper.
HIT@n is deﬁned as HIT @ n = (cid:80)
δ i
# of Test Interaction , where δ i =
if rank i < = n else
[Figure]
which measures the ability of top rank prediction
[Figure] [Figure]
[Figure]
B. E FFECT OF THE DEPTH K AND THE NUMBER OF LAYERS L
For a given node u ( t )
[Figure]
only the history information in k-depth temporal subgraph ,
i
[Figure]
e, subgraph ( u ( t ) , k ) , is considered to calculate the dynamic representation, so the number of depth k could affect the ﬁnal representation.
Meanwhile, our proposed fusion mechanism is based on deep model which can capture different information at different layers.
So, in this subsection
[Figure]
we investigate the effects of k and L for CollegeMsg dataset on the Interaction Prediction task and
Huabei Trade dataset on the Interaction Recognition, which is shown in Figure 9.
As we can see, stacking multiple layers does help improve the interaction prediction and interaction classiﬁcation performances
[Figure]
but it doesn’t mean the larger L will always give the better results.
And with the k increases which means we could utilize more history information, the trend of the performances go better for both tasks
[Figure]
B. E FFECT OF EMBEDDING SIZES
To investigate the effect of the dynamic embedding size on the interaction prediction and interaction classiﬁcation tasks, we vary the dynamic embedding dimension in {16, 32, 64, 128} and report the corresponding Mean Rank results for interaction prediction on the CollegeMsg dataset and KS results(AUC is similar) for interaction classiﬁcation on Huabei Trade Data in ﬁgure 10.
In ﬁgure(a), we see that the dynamic embedding size has little impact on DIP model, and it consistently outperforms all baselines.
The performances of several baselines drop more or less when the embedding size is set to 128.
The baseline model could overﬁt at this embedding size.
In ﬁgure 10 (b) we ﬁnd that when all the baselines are set with the same dynamic embedding sizes, our DIP model always outperforms the competitors.
C D ETAILED R ELATED W ORK
Graph representation learning, which is also known as graph embedding or network embedding, is a task aiming to learn low-dimensional dense vectors for vertex and edges that preserve the original graph structural information and network properties.
Before the era of deep learning, conventional graph representation learning methods usually adopts dimension reduction (Belkin & Niyogi, 2002;
Tenenbaum et al., 2000; Roweis & Saul, 2000) or matrix factorization(Ahmed et al., 2013) techniques.
While these methods usually suffer a heavy computational cost.
After the great success of deep learning methods in the ﬁeld of computer vision and natural language processing, especially Skipgram(Mikolov et al., 2013) in word representation learning, this problem has been widely investigated by neural network methods.
We review the development of this research direction from the following perspectives.
Static Graph Representation Learning :
Perozzi et al.
(2014) proposed the DeepWalk model utilizing random walks on graph to generate sequences of nodes, and feed them into the Skip-gram model to learn nodes’ low dimensional representation vector.
Then, Tang et al.
(2015) proposed the deﬁnition of ﬁrst-order proximity and second-order proximity in the model named LINE , that jointly models this two different level structure information.
Grover & Leskovec (2016) designed a biased random-walk sampling method to capture more ﬂexible neighborhood information which is helpful to learn richer representations.
And Wang et al.
(2016; 2017) extended these models by considering high-order proximity and community structure.
But these mentioned methods are somehow shallow.
Another type of method is the Graph Convolution Network (Bruna et al., 2013; Kipf & Welling, 2017).
They extended the convolution neural network to the graph spectral space, thus applied DEEP model on graph data.
While these methods are still transductive , and can not jointly models the network structure information and the attributed on nodes or edges of graph.
To overcome this problem, Hamilton et al. (2017a) proposed an inductive graph embedding framework, which divided the representation learning into two phrase: sampling and aggregation, such that it could incorporates node feature information and thus generate embedding from these features using the well-trained graph neural network for unseen nodes.
Veliˇckovi´c et al.
(2017) used self-attention to learn weights for neighborhood, then aggregated them by the self-adapted weight.
While, most of these approaches can only model static graph data.
Dynamic Graph Representation Learning: A popular approach for modeling dynamic graph data is considering the dynamics as a sequence of graph snapshot(Soundarajan et al., 2016).
Zhu et al.
(2016) uses an non-negative matrix factorization technique to embed social network in a temporal latent space.
Zhou et al.
(2018) models the speciﬁc triadics closure formation procedure over the snapshots.
Goyal et al.
(2018) adapts an graph autoencoder and learns a stable embedding over time.
Seo et al.
(2018) combines a CNN module and a RNN module to capture spatial characteristics and temporal characteristics.
NetWalk(Yu et al., 2018) also learns vertex representations from sequences of snapshots, and it is designed specially for anomaly detection.
In contrast, Nguyen et al.
(2018) adds temporal order constraint on random walk sampling to capture evolving neighborhood.
But, it can’t explicitly models the rich temporal information.
Temporal Point Process: Temporal Point Process is a powerful statistical tool for modeling sequences of events with unequal interval.
It has been wildly used in recent research with different intensity function: from parametric(Du et al., 2016) to recurrent neural networks(Dai et al., 2016; Mei & Eisner, 2017b), event reinforcement learning based(Xiao et al., 2017).
DeepCoevolve(Dai et al., 2016) designs a recurrent neural network to capture the co-evolution dynamics of interaction data.
But it just takes the interaction counterpart into the co-evolution module, which may loss temporal structural information.
DyREP(Trivedi et al., 2019) divides the interactions between nodes into two different types as communication and association and models them separately.
Although it uses the neighborhood of interaction counterpart to update node representation, the relevance between the different neighbors of counterpart and the node itself has not been modeled.