T IME
D
EPENDENT R EPRESENTATION FOR N EURAL
E VENT S EQUENCE P REDICTION
Anonymous authors
Paper under double-blind review
A BSTRACT
Existing sequence prediction methods are mostly concerned with time-independent sequences, in which the actual time span between events is irrelevant and the distance between events is simply the difference between their order positions in the sequence.
While this time-independent view of sequences is applicable for data such as natural languages, e.g., dealing with words in a sentence, it is inappropriate and inefﬁcient for many real world events that are observed and collected at unequally spaced points of time as they naturally arise, e.g., when a person goes to a grocery store or makes a phone call.
The time span between events can carry important information about the sequence dependence of human behaviors.
In this work, we propose a set of methods for using time in sequence prediction.
Because neural sequence models such as RNN are more amenable for handling token-like input, we propose two methods for time-dependent event representation, based on the intuition on how time is tokenized in everyday life and previous work on embedding contextualization.
We also introduce two methods for using next event duration as regularization for training a sequence prediction model.
We discuss these methods based on recurrent neural nets.
We evaluate these methods as well as baseline models on ﬁve datasets that resemble a variety of sequence prediction tasks.
The experiments revealed that the proposed methods offer accuracy gain over baseline models in a range of settings.
I NTRODUCTION
Event sequence prediction is a task to predict the next event 1 based on a sequence of previously occurred events.
Event sequence prediction has a broad range of applications, e.g., next word prediction in language modeling (Józefowicz et al., 2016), next place prediction based on the previously visited places, or next app to launch given the usage history.
Depending on how the temporal information is modeled, event sequence prediction often decomposes into the following two categories: discrete-time event sequence prediction and continuous-time event sequence prediction.
Discrete-time event sequence prediction primarily deals with sequences that consist of a series of tokens (events) where each token can be indexed by its order position in the sequence.
Thus such a sequence evolves synchronously in natural unit-time steps.
These sequences are either inherently time-independent, e.g, each word in a sentence, or resulted from sampling a sequential behavior at an equally-spaced point in time, e.g., busy or not busy for an hourly trafﬁc update.
In a discrete-time event sequence, the distance between events is measured as the difference of their order positions.
As a consequence, for discrete-time event sequence modeling, the primary goal is to predict what event will happen next.
Continuous-time event sequence prediction mainly attends to the sequences where the events occur asynchronously.
For example, the time interval between consecutive clinical visits of a patient may potentially vary largely.
The duration between consecutive log-in events into an online service can change from time to time.
Therefore, one primary goal of continuous-time event sequence prediction is to predict when the next event will happen in the near future.
Although these two tasks focus on different aspects of a future event, how to learn a proper representation for the temporal information in the past is crucial to both of them.
More speciﬁcally, even though for a few discrete-time event sequence prediction tasks (e.g., neural machine translation), they do not involve an explicit temporal information for each event (token), a proper representation of the position in the sequence is still of great importance, not to mention the more general cases where each event is particularly associated with a timestamp.
For example, the next destination people want to go to often depends on what other places they have gone to and how long they have stayed in each place in the past.
When the next clinical visit (Choi et al., 2016a) will occur for a patient depends on the time of the most recent visits and the respective duration between them.
Therefore, the temporal information of events and the interval between them are crucial to the event sequence prediction in general.
However, how to effectively use and represent time in sequence prediction still largely remains under explored.
A natural and straightforward solution is to bring time as an additional input into an existing sequence model (e.g., recurrent neural networks).
However, it is notoriously challenging for recurrent neural networks to directly handle continuous input that has a wide value range, as what is shown in our experiments.
Alternatively, we are inspired by the fact that humans are very good at characterizing time span as high-level concepts.
For example, we would say "watching TV for a little while" instead of using the exact minutes and seconds to describe the duration.
We also notice that these high-level descriptions about time are event dependent.
For example, watching movies for 30 minutes might feel much shorter than waiting in the line for the same amount of time.
Thus, it is desirable to learn and incorporate these time-dependent event representations in general.
Our paper offers the following contributions:
• We propose two methods for time-dependent event representation in a neural sequence prediction model: time masking of event embedding and event-time joint embedding.
We use the time span associated with an event to better characterize the event by manipulating its embedding to give a recurrent model additional resolving power for sequence prediction.
• We propose to use next event duration as a regularizer for training a recurrent sequence prediction model.
Speciﬁcally, we deﬁne two ﬂavors of duration-based regularization: one is based on the negative log likelihood of duration prediction error and the other measures the cross entropy loss of duration prediction in a projected categorical space.
• We evaluated these proposed methods as well as several baseline methods on ﬁve datasets (four are public).
These datasets span a diverse range of sequence behaviors, including mobile app usage, song listening pattern, and medical history.
The baseline methods include vanilla RNN models and those found in the recent literature.
These experiments offer valuable ﬁndings about how these methods improve prediction accuracy in a variety of settings.
B ACKGROUND
In recent years, recurrent neural networks (RNN) especially with Long-Short Term Memory (LSTM) (Hochreiter & Schmidhuber, 1997) have become popular in solving a variety of discretetime event sequence prediction problems, including neural machine translation (Bahdanau et al., 2014), image captioning (Xu et al., 2015) and speech recognition (Soltau et al., 2016).
In a nutshell, given the sequence of previously occurred events, { e 1 , e 2 , ..., e t } , the conditional probability
P ( e t +1 |{ e 1 , e 2 , ..., e t } ) =
P ( e t +1 | h t , θ ) of the next event e t +1 is estimated by using a recurrent neural network with parameters θ and the hidden state vector h t =
f
( h t − 1 , e t , θ ) which is assumed to encode the information of the past events.
To feed an event into a recurrent neural network, the event, often described as a categorical variable, needs to be represented in a continuous vector space.
A common way to achieve this is to use embedding (Bengio et al., 2003)
x t = 1( e t )
E x where 1( e t ) is a one-hot vector.
For the j th event in the vocabulary V , e j , its one-hot vector has 0 s for all the entries except the j th entry being 1 .
E x ∈ R | V |× E is the embedding matrix, where | V | is the number of unique events (the vocabulary size) and E is the embedding dimension.
The use of embedding provides a dense representation for an event that improves learning (Turian et al., 2010).
Through training, the embedding vector of an event encodes its meaning relative to other events.
Events that are similar tend to have embedding vectors closer to each other in the embedding space than those that are not.
On the other hand, temporal point processes are a mathematical abstraction for the continuous-time event sequence prediction task by explicitly modeling the inter-event interval as a continuous random variable.
Since the occurrence of an event may be triggered by what happened in the past, we can essentially specify different models for the timing of the next event given what we have already known so far.
Very recently, (Du et al., 2016; Mei & Eisner, 2017; Xiao et al., 2017a;b) focus on expanding the ﬂexibility of temporal point processes using recurrent neural networks where the prediction of the next event time is based on the current hidden state h t of RNN.
However, all of these work use the direct concatenation between the inter-event interval and the respective event embedding as the input to the recurrent layer where the representation of the temporal information is limited.
In contrast, in this work, we focus on a more principled way to learn a proper representation for encoding temporal information in a sequence, which can be useful for general event sequence prediction.
We also remark that although we focus on improving discrete-time event prediction in the following experiments, our temporal embedding can be generalized to the continuous-time setting straightforwardly.
T IME -D EPENDENT E VENT R EPRESENTATION
There are two notions about time spans in a sequential behavior: duration and intervals.
Duration is how long an event lasts, e.g., listening to music for an half hour, and an interval is the time span between two adjacent events.
To unify both types of time spans, we treat the idle period when no event is occurring (e.g., the person is not using any app for an app usage sequence) as a special event.
Thus, duration becomes an inherent property of an event–the interval between two events is the duration of an idle event (see Figure 1).
With this, h t = f ( h t − 1 , e t , d t ; θ ) where d t is the duration of event e t .
We here propose two methods to bring continuous time, d t , into a neural sequence prediction model.
Both achieve timedependent event representation by manipulating event embedding vectors using time.
Our methods are schematically illustrated in Figure 2.
3. "C ONTEXTUALIZING " E VENT E MBEDDING WITH T IME M ASK
Recent work by (Choi et al., 2016b) revealed that in neural machine translation the embedding vector of a word encodes multiple meanings of the word.
As a result, it requires a recurrent layer to sacriﬁce its capacity to disambiguate a word based on its context, instead of focusing on its main task for learning the higher-level compositional structure of a sentence.
To address this problem, they used a mask computed based on all the words in a sentence to contextualize the embedding of a target word.
Based on this recent work, we propose a method to learn a time mask to "contextualize" event embedding, by which we hope a time-dependent embedding would give the recurrent layer additional resolving power.
Similar to the word mask proposed by Choi et al.
(Choi et al., 2016b), we ﬁrst compute a time context vector for duration, c d .
[Equation]
φ is a nonlinear transformation of d t and is implemented as a feedforward neural network parameterized by θ .
d t is log transformed before it is fed to φ to effectively cover the wide numerical range of duration values, e.g., it can range from seconds to hours for app usage events.
[Figure]
We compute a time mask by linearly transforming c d with weights W d ∈ R C × E and bias b d ∈ R E , which is followed by a sigmoid nonlinear activation, σ , to generate a mask m d ∈ R E and R E → [0 , 1] .
C is the size of the time context vector, and E is the event embedding dimension.
[Equation]
We then apply the mask to an event embedding by performing an element-wise multiplication, (cid:12) , between the embedding vector and the mask.
Finally, the product is fed to the recurrent layer.
[Equation]
3. E VENT -T IME J OINT E MBEDDING
Humans developed many ways to tokenize continuous time in everyday life.
For example, we would say "talk to someone brieﬂy" instead of using exact minutes and seconds to characterize the length of the conversation.
Such a kind of tokenization is extensively used in natural languages.
In addition, our perception about the duration also depends on the speciﬁc event that we are experiencing.
Based on these intuitions, we propose a method to ﬁrst encode the duration of an event using soft one-hot encoding and then use the encoding to form the joint embedding with the event.
To do so, we ﬁrst project the scalar duration value onto a vector space, where W d ∈ R 1 × P is the weight matrix, b d ∈ R P is the bias vector, and P is the projection size.
[Equation]
We then compute the soft one-hot encoding, s d , of a duration value by applying a softmax function to the projection vector, p d .
Softmax has been typically used in the output layer (Graves, 2012) and in the attention mechanisms (Bahdanau et al., 2014; Xu et al., 2015) for selecting one out of many.
The i th entry of the encoding vector is calculated as the following and p di is the i th entry in p d .
[Equation]
All the entries in the soft one-hot encoding are positive.
Similar to a regular one-hot encoding, (cid:80)
Pi =1 s di = 1 .
We then project the soft one-hot encoding onto a time embedding space, g d .
It has the same dimension as the event embedding.
E s ∈ R P
× E is the embedding matrix.
[Equation]
Embedding for a regular one-hot encoding essentially takes a single row of the embedding matrix that is corresponding to the non-zero entry as the embedding vector.
In contrast, embedding for a soft one-hot encoding computes a weighted sum over all the rows in the embedding matrix.
Finally, we form the joint embedding of an event and its duration by taking the mean of their embedding vectors, which is then fed to the recurrent layer.
[Equation]
N EXT E VENT D URATION AS A R EGULARIZER
While our goal here is to predict next event, it can help learning by introducing an additional loss component based on the prediction of the next event duration (see Figure 2).
The duration prediction of the next event at step t , d (cid:48) t +1 , is computed from a linear transformation of the recurrent layer.
A loss deﬁned on the prediction error of d (cid:48)
t +1 provides additional information during back propagation, acting like a regularizer.
Optionally, one can use the concatenation of the recurrent layer output and a hidden layer on the path for event prediction to regularize more layers.
We discuss two alternatives for the loss function over d (cid:48) t +1 .
4. N EGATIVE L OG L IKELIHOOD OF T IME P REDICTION E RROR
A common way for the loss over a continuous value is to use the squared error.
Here, it is ( d (cid:48)
t +1 − d t +1 ) 2 where d t +1 is the observed duration of the next event.
However, such a loss needs to be at the same scale as that of of event prediction, which is typically a log likelihood of some form.
Hinton and Van Camp (Hinton & van Camp, 1993) have shown that minimizing the squared error can be in fact formulated as maximizing the probability density of a zero-mean Gaussian distribution.
Note that this does not require duration to obey a Gaussian distribution but rather the prediction error.
We deﬁne our regularizer, R Nt , as the negative log likelihood of duration prediction error at step t .
[Equation]
The variance, σ i , is seeded with an initial value (e.g., the variance of duration values in the training data) and updated iteratively during training based on the duration prediction error distribution of the learned model at each update i .
4. C ROSS E NTROPY L OSS ON T IME P ROJECTION
In Section 3.2, we proposed to use softmax to project a continuous duration value onto a categorical space.
Using the same technique, by projecting both d (cid:48) t +1 and d t +1 onto a categorical space, we can then compute a cross entropy loss based on the two projections as another regularizer R Xt .
[Equation]
[Equation]
Proj is the softmax projection process we deﬁned in Equation 4 and 5, Proj k is the k th entry in the projection vector.
When event-time joint embedding and R Xt are both used, the embedding and the regularizer can use the same projection function, i.e., sharing the same projection weights (Equation 4).
E XPERIMENTS
In this section, we evaluate the effectiveness of our proposed approaches on the following ﬁve real-world datasets across a diverse range of domains.
• Electrical Medical Records.
MIMIC II medical dataset is a collection of de-identiﬁed clinical visit records of Intensive Care Unit patients for seven years.
The ﬁltered dataset released by (Du et al., 2016) include 650 patients and 204 diseases.
The goal is to predict which major disease will happen to a given patient.
• Stack Overﬂow Dataset.
The Stack Overﬂow dataset includes two years of user awards on a question-answering website.
The awarded badges are treated as the events.
(Du et al., 2016) collected 6,000 users with a total of 480,000 events.
The goal is to predict the next badge a user will receive.
• Financial Transaction Dataset.
(Du et al., 2016) collected a long stream of high frequency transactions for a single stock from NYSE where the events correspond to the "buy" and "sell" actions.
The task is to predict the next action a user might take.
• App Usage Dataset.
Mobile users often use a large number of apps, ranging from tens to hundreds.
It is time consuming to ﬁnd a target app on mobile devices.
One promising way to address this problem is to predict the next app a user will use based on their app usage history.
Being able to predict next apps also allows the mobile platform to preload an app in memory to speed up its startup.
We have collected 5,891 app usage sequences comprising of 2.8 million app usage events.
The task is to predict the next app that will be used for a given user.
• Music Recommendation.
The music dataset represents the longitudinal listening habits of 992 users (Last. FM, 2009; Celma, 2010) involving millions of listening events. The goal is to predict the next ﬁve unique songs that the user has not listened given the user’s listen history.
5. D ATA P REPARATION
For the MIMIC II, Stack Overﬂow, and Financial data, we follow (Du et al., 2016) to pre-process the data and seek to predict every single held-out event from the history.
We evaluate the prediction accuracy with the binary 0-1 loss.
For the app usage data, to avoid users who participated in the data collection only brieﬂy, we exclude sequences that have fewer than 50 app launches or if the time span of the sequence is shorter than a week.
This resulted in 5,891 app usage sequences, one from each unique user.
These sequences include 2,863,095 app usage events and the longest sequence spanned 551 days.
We split the dataset on users into the training (80%), validation (10%) and test (10%) such that each user is only in one of these partitions.
Hence there is no intersection of users between training, validation and test sets.
For an event that has fewer than 5 occurrences in the training dataset, we assign it the OOV id for out of vocabulary.
In total, there are 7,327 events in the vocabulary, including 7,325 unique apps, the idle event and the OOV (out of vocabulary).
In practice, predicting the next 5 apps is often desired so we use Precision@K to evaluate the performance.
For the music recommendation, each listen event in the original dataset has a timestamp.
We removed sequences that are shorter than 50 and songs that have fewer than 50 listens.
Based on these listen habits, we generate a collection of examples where each example consists of a listen history and a set of 5 unique songs to recommend.
To do so, we split each original listen sequence into segments (examples).
To extract an example, we ﬁrst take the 40 events out in order from the beginning of the sequence as the listen history, and then take more events out from the beginning of the sequence until we ﬁnd 5 unique songs that have not occurred in the listen history.
We do so repeatedly to extract each example until we exhaust all the original sequences.
This data processing resulted in 221,  sequence examples with 71,619 unique songs (the vocabulary size).
We then allocate these sequence examples for the training (80%), validation (10%) and test (10%).
Because the original dataset does not have the duration information for each listen event, we did not inject the additional idle event in the sequence to differentiate duration versus intervals.
Because in practice, the ranking order of the recommended music often matters, we further use MAP@K and Prevision@K to evaluate the performance.
5. M ODEL C ONFIGURATIONS
We compare with the following ﬁve models: NoTime in which simple LSTM sequence model is used;
TimeConcat in which we feed time (log transformed) directly into the recurrent layer along the event embedding; TimeMask (Section 3.1) and TimeJoint (Section 3.2) for generating time-dependent event embedding as input to the recurrent layer; and RMTPP for the model introduced previously by (Du et al., 2016).
Moreover, we also include four regularized models based on R Xt and R Nt deﬁned earlier.
For TimeMask , the size of the time context vector is C = 32 , and we use ReLu for the activation function in φ in Equation 2.
For TimeJoint , we chose the projection size, P = 30 (Equation 4).
For all these models, we use a two-layer hierarchical softmax (Morin & Bengio, 2005) for the output layer.
For the MIMIC II, Stack Overﬂow, and Financial data, we follow (Du et al., 2016) for RMTPP ’s model parameters.
For the app usage data, we determined the parameters of each model based on the training and the validation datasets on a distributed parallel tuning infrastructure.
We used
LSTM units (Hochreiter & Schmidhuber, 1997) for the recurrent layer, and Rectiﬁed Linear Units (ReLu) (Nair & Hinton, 2010) for the activation function in the nonlinear projection layer.
The event embedding dimension, the number of LSTM units, and the nonlinear projection layer size are all set to 128.
For the music recommendation data, we use a setting similar to the app prediction experiment where we chose the embedding size as 128 and LSTM size as 256.
We did not use the nonlinear projection layer after the LSTM layer for this task because it does not seem to help.
We implemented all the models in TensorFlow (TensorFlow, 2017).
5. T RAINING AND T ESTING
For the experiments based on MIMIC II, Stack Overﬂow and Financial Transaction datasets, we use the same training and testing strategy of (Du et al., 2016).
For App Usage and Music Recommendation tasks, we selected the model architecture and hyper parameters with early stopping based on the validation dataset of each task, and report the performance of each model based on the test dataset.
For the App Usage experiment, we used truncated back-propagation through time with the number of unroll to be 30.
We used an adaptive gradient descent optimizer (Zeiler, 2012), using a learning rate of 0.024 with a threshold for gradient clipping of 1.0, and a batch size of 32.
We decided not to use dropout as it did not seem to improve accuracy on this task.
For the Music Recommendation experiment, we used the full sequence back-propagation through time with 2% dropout ratio on the recurrent layer for better generalization.
We used the Adam optimizer by (Kingma & Ba, 2014) for adaptive learning with a learning rate of 0.00005 and a gradient clipping threshold at 1.0.
The mini-batch size is 256.
We trained the models by minimizing the cross-entropy loss, plus the regularization loss if the duration regularizer is used, over all the sequences in the training dataset.
The training for App Usage and
Music Recommendation was conducted on a distributed learning infrastructure (Dean et al., 2012) with 50 GPU cores where updates are applied asynchronously across multiple replicas.
5. E XPERIMENTAL R ESULTS
Effectiveness of Temporal Representation .
Figure 3 presents the comparisons between all baselines.
First, we can observe a consistent performance gain with using the proposed methods for timedependent event embedding compared to the NoTime baseline and the simple TimeConcat approach.
Second, we should remark that our primary goal here is to improve the predictive performance of the next event by utilizing the temporal information.
In comparison, the RMTPP seeks to predict both the next event and its occurrence time, which essentially learns a balance between these two goals.
Methods
A cc u r a cy %
Methods
NoTimeTimeConcatRMTPPTimeMaskTimeJoint 46.
Methods
A cc u r a cy %
Methods
NoTimeTimeConcatRMTPPTimeMaskTimeJoint 62.
Methods
A cc u r a cy %
Methods
NoTimeTimeConcatRMTPPTimeMaskTimeJoint
[Figure]
[Table] [Table]
TimeMask 11.
TimeMask + R Xt 11.
TimeMask + R Nt 11.
TimeJoint 11.
TimeJoint + R Xt 12.
TimeJoint + R Nt 11.
Therefore, if our major task is the next event prediction, we can also observe a performance gain over
RMTPP as well.
Effectiveness of Event Duration Regularization .
Finally, we demonstrate the performance boosting gained from our proposed temporal regularization in Table 1 and 2, respectively.
We can observe that our proposed cross-entropy regularizer is able to give consistent performance gain with the temporal embedding approaches.
C ONCLUSIONS
We proposed a set of methods for leveraging the temporal information for event sequence prediction.
Based on our intuition about how humans tokenize time spans as well as previous work on contextual representation of words, we proposed two methods for time-dependent event representation.
They transform a regular event embedding with learned time masking and form time-event joint embedding based on learned soft one-hot encoding.
We also introduced two methods for using next duration as a way of regularization for training a sequence prediction model.
Experiments on a diverse range of real data demonstrate consistent performance gain by blending time into the event representation before it is fed to a recurrent neural network.
A PPENDIX