R EPRESENTATION L EARNING FOR S EQUENCE D ATA WITH D EEP A UTOENCODING P REDICTIVE C OMPONENTS
Anonymous authors
Paper under double-blind review
A BSTRACT
We propose Deep Autoencoding Predictive Components (DAPC) – a selfsupervised representation learning method for sequence data, based on the intuition that useful representations of sequence data should exhibit a simple structure in the latent space.
We encourage this latent structure by maximizing an estimate of predictive information of latent feature sequences, which is the mutual information between past and future windows at each time step.
In contrast to the mutual information lower bound commonly used by contrastive learning, the estimate of predictive information we adopt is exact under a Gaussian assumption.
Additionally, it can be computed without negative sampling.
To reduce the degeneracy of the latent space extracted by powerful encoders and keep useful information from the inputs, we regularize predictive information learning with a challenging masked reconstruction loss.
We demonstrate that our method recovers the latent space of noisy dynamical systems, extracts predictive features for forecasting tasks, and improves automatic speech recognition when used to pretrain the encoder on large amounts of unlabeled data.
I NTRODUCTION
Self-supervised representation learning methods aim at learning useful and general representations from large amounts of unlabeled data, which can reduce sample complexity for downstream supervised learning.
These methods have been widely applied to various domains such as computer vision (Oord et al., 2018;
Hjelm et al., 2018; Chen et al., 2020; Grill et al., 2020), natural language processing (Peters et al., 2018;
Devlin et al., 2019; Brown et al., 2020), and speech processing (Schneider et al., 2019; Pascual et al., 2019b; Chung & Glass, 2020; Wang et al., 2020; Baevski et al., 2020).
In the case of sequence data, representation learning may force the model to recover the underlying dynamics from the raw data, so that the learnt representations remove irrelevant variability in the inputs, embed rich context information and become predictive of future states.
The effectiveness of the representations depends on the self-supervised task which injects inductive bias into learning, and the design of self-supervision has been an active research area.
One notable approach for self-supervised learning is based on maximizing mutual information between learned representation and inputs.
The most commonly used estimate of mutual information is based on contrastive learning.
A prominant example of this approach is CPC (Oord et al., 2018), where the representation of each time step is trained to distinguish between positive samples which are inputs from the near future, and negative samples which are inputs from distant future or other sequences.
The performance of contrastive learning heavily relies on the nontrivial selection of positive and negative samples, which lacks universal principles across different scenarios (He et al., 2020; Chen et al., 2020; Misra & Maaten, 2020).
Recent works have suspected that the mutual information lower bound estimate used by contrastive learning might be loose and may not be the sole reason for its success (Ozair et al., 2019; Tschannen et al., 2019).
In this paper, we leverage an estimate of information speciﬁc to sequence data, known as predictive information (PI, Bialek et al., 2001), which measures the mutual information between the past and future windows in the latent space.
The estimate is exact if the past and future windows have a joint Gaussian distribution, and is shown by prior work to be a good proxy for the true predictive information in practice (Clark et al., 2019).
We can thus compute the estimate with sample windows of the latent sequence (without sampling negative examples), and obtain a well-deﬁned objective for learning the encoder for latent representations.
However, simply using the mutual information as the learning objective may lead to degenerate representations, as PI emphasizes simple structures in the latent space and a powerful encoder could achieve this at the cost of ignoring information between latent representations and input features.
To this end, we adopt a masked reconstruction task to enforce the latent representations to be informative of the observations as well.
Similar to Wang et al.
, we mask input dimensions as well as time segments of the inputs, and use a decoder to reconstruct the masked portion from the learnt representations; we also propose variants of this approach to achieve superior performance.
Our method, Deep Autoencoding Predictive Components (DAPC), is designed to capture the above intuitions.
From a variational inference perspective, DAPC also has a natural probabilistic interpretation.
We demonstrate DAPC on both synthetic and real datasets of different sizes from various domains.
Experimental results show that DAPC can recover meaningful low dimensional dynamics from high dimensional noisy and nonlinear systems, extracts predictive features for forecasting tasks, and obtains state-of-the-art accuracies for Automatic Speech Recognition (ASR) with much lower cost, by pretraining encoders that are later ﬁnetuned with limited amounts of supervised data.
[Equation]
[Figure]
The main intuition behind Deep Autoencoding Predictive Components is to maximize the predictive information of latent representation sequence.
To ensure the learning process is tractable and nondegenerate, we make Gaussian assumption and regularize the learning with masked reconstruction.
In the following, we elaborate on how we estimate the predictive information and how we design the masked reconstruction task.
A probabilistic interpretation of DAPC is also provided to show the connections to deep generative models.
2. P REDICTIVE I NFORMATION
Given a sequence of observations X
= { x 1 , x 2 , ... }
where x i ∈ R n , we extract the corresponding latent sequence Z = { z 1 , z 2 , ... }
where z i ∈ R d with an encoder function e ( X ) , e.g., recurrent neural nets or transformers (Vaswani et al., 2017).
1 Let T > 0 be a ﬁxed window size, and denote Z pastt = { z t − T +1 , ..., z t } , Z futuret = { z t +1 , ..., z t + T } for any time step t .
The predictive information (PI) is deﬁned as the mutual information (MI) between Z pastt and Z futuret :
[Equation]
where H is the entropy function.
Intuitively, PI measures how much knowing Z pastt reduces the uncertainty about Z futuret (and vice versa).
PI reaches its minimum value
if Z pastt and Z futuret are independent, and it is maximized if Z futuret is a deterministic function of Z pastt .
Different from the MI estimate used by contrastive learning, which measures the MI between representation at each single time step and its future inputs, predictive information measures the MI between two windows of T time steps collectively.
The window size T used in PI estimation reﬂects the time resolution for which the time series is more or less stationary.
PI was designed as a general measure of the complexity of underlying dynamics which persists for a relatively long period of time (Li & Vit´anyi, 2008).
Furthermore, PI is aware of temporal structures: different dynamics could lead PI to converge or diverge even if they look similar.
These virtues of
PI contribute to the versatility of this measure.
The use of PI beyond a static complexity measure (and as a learning objective) is done only recently in machine learning by Clark et al.
(2019), which proposes to learn a linear dimensionality reduction method named Dynamical Component Analysis (DCA) to maximize the PI of the projected latent sequence.
One approach for estimating PI is through estimating the joint density P ( Z pastt , Z futuret ) , which can be done by density estimation methods such as k -NN and binning (Dayan & Abbott, 2001; Kraskov et al., 2004).
However, such estimates heavily rely on hyperparameters, and it is more challenging to come up with differentiable objectives based on them that are compatible with deep learning frameworks.
Our approach for estimating PI is the same as that of DCA.
Assume that every 2 T consecutive time steps { z t − T +1 , ..., z t , ..., z t + T } in the latent space form a stationary, multivariate
Gaussian distribution.
Σ 2 T ( Z ) is used to denote the covariance of the distribution, and similarly Σ T ( Z ) the covariance of T consecutive latent steps.
Under the stationarity assumption, H ( Z pastt ) remains the same for any t so we can omit the subscript t , and H ( Z past ) is equal to H ( Z future ) as well.
Using the fact that H ( Z past ) = 12 ln(2 πe ) dT | Σ T ( Z ) | , PI for the time series
z reduces to
[Equation]
Detailed derivations can be found in Appendix A.
It is then straightforward to collect samples of the consecutive 2 T -length windows and compute the sample covariance matrix for estimating Σ 2 T ( Z ) , and an empirical estimate of Σ T ( Z ) corresponds to the upper left sub-matrix of Σ 2 T ( Z ) .
Recall that, under the Gaussian assumption, the conditional distribution P ( Z future | Z past ) is again Gaussian, whose mean is a linear transformation of Z past .
Maximizing I T has the effect of minimizing the entropy of this conditional Gaussian, and thus reducing the uncertainty of future given past.
Though our estimation formula for PI is exact only under the Gaussian assumption, it was observed by Clark et al.
(2019) that the Gaussian-based estimate is positively correlated with a computationally intensive estimate based on non-parametric density estimate, and thus a good proxy for the full estimate.
We make the same weak assumption, so that optimizing the Gaussian-based estimate improves the true PI.
Our empirical results show that representations learnt with the Gaussian PI have strong predictive power of future (see Sec 4.2).
Furthermore, we ﬁnd that a probabilistic version of
DAPC (described in Sec 2.3) which models ( Z past , Z future ) with a Gaussian distribution achieves similar performance as this deterministic version (with the Gaussian assumption).
We now describe two additional useful techniques that we develop for PI-based learning.
Multi-scale PI One convenient byproduct of this formulation and estimation for I T is to reuse Σ 2 T ( Z ) for estimating I T/ 2 , I T/ 4 and so on, as long as the window size is greater than 1.
Since the upper left sub-matrix of Σ 2 T ( Z ) approximates Σ T ( Z ) , we can extract Σ T ( Z ) from Σ 2 T ( Z ) without extra computation, and similarly for Σ T/ 2 ( Z ) .
We will show that multi-scale PI, which linearly combines PI at different time scales, boosts the representation quality in ASR pretraining.
Orthogonality penalty Observe that the PI estimate in (2) is invariant to invertible linear transformations in the latent space.
To remove this degree of freedom, we add the penalty to encourage latent representations to have identity covariance, so that each of the d latent dimensions will have unit scale and different dimensions are linearly independent and thus individually useful.
This penalty is similar to the constraint enforced by deep canonical correlation analysis (Andrew et al., 2013), which was found to be useful in representation learning (Wang et al., 2015).
2. M ASKED R ECONSTRUCTION AND I TS S HIFTED V ARIATION
The PI objective alone can potentially lead to a degenerate latent space, when the mapping from input sequence to latent sequence is very powerful, as the latent representations can be organized in a way that increases our PI estimate at the cost of losing useful structure from the input.
This is also observed empirically in our experiments (see Sec 4.1).
To regularize PI-based learning, one simple idea is to force the learnt latent representations to be informative of the corresponding input observations.
For this purpose, we augment PI-based learning with a masked reconstruction task.
Masked reconstruction was ﬁrst proposed in BERT (Devlin et al., 2019), where the input text is fed to a model with a portion of tokens masked, and the task is to reconstruct the masked portion.
Wang et al.
(2020) extended the idea to continuous vector sequence data (spectrograms).
The authors found that randomly masking input dimensions throughout the sequence yields further performance gain, compared to masking only consecutive time steps.
We adopt their formulation in DAPC to handle continuous time series.
Given an input sequence X of length L and dimensionality n , we randomly generate a binary mask
M ∈ R n × L , where M i,j = 0 indicates X i,j is masked with value 0 and M i,j = 1 indicates X i,j is kept the same.
We feed the masked inputs to the encoder e ( · ) to extract representations (in R d ) for each time step, and use a feed-forward network g ( · ) to reconstruct the masked input observations.
e ( · ) and g ( · ) are trained jointly.
The masked reconstruction objective can be deﬁned as
[Equation]
Figure 1 gives an illustration for the masked spectrogram data.
For time axis, we randomly generate n T time masks each with width up to w T , and similarly n F frequency masks each with width up to w F .
In our experiments, we observe that input dimension masking makes the reconstruction task more challenging and yields higher representation quality.
Therefore, this strategy is useful for general time series data beyond audio.
We introduce one more improvement to masked reconstruction.
Standard masked reconstruction recovers the masked inputs for the same time step.
Inspired by the success of Autoregressive Predictive Coding (Chung & Glass, 2020), we propose a shifted variation of masked reconstruction, in which the latent state z
i is decoded to reconstruct a future frame x i + s (than x i ).
Formally, the shifted masked reconstruction loss R s is deﬁned as
[Equation]
where → s indicates right-shifting s time frames while the input dimensions remain unchanged.
When s = 0 , R s reduces to the standard masked reconstruction objective, and in the ASR experiments we ﬁnd that a nonzero s value helps.
We ensure no information leakage by enforcing that the portion to be reconstructed is never presented in the inputs.
As indicated by Chung et al. (2019b), predicting a future frame encourages more global structure and avoids the simple inference from local smoothness in domains like speech, and therefore helps the representation learning.
To sum up, our overall loss function is deﬁned as the combination of the losses described above:
[Equation]
where α, β, γ are tradeoff weights and R ortho =
|| Σ 1 − I d
|| 2 fro is the orthonormality penalty discussed in Sec. 2.1, with Σ 1 ∈ R d × d corresponding to the top left sub-matrix of Σ 2 T estimated from the latent sequence Z =
e ( X (cid:12) M ) .
The whole framework of DAPC is illustrated in Figure 1.
2. A P ROBABILISTIC I NTERPRETATION OF DAPC
We now discuss a probabilistic interpretation of DAPC in the Variational AutoEncoder (VAE) framework (Kingma & Welling, 2014).
Let X =
( X p , X f ) and Z =
( Z p , Z f ) , where the subscript p and f denote past and future respectively.
Consider a generative model, where the prior distribution is p ( Z p , Z f ) ∼ N ( 0 , Σ) .
One can write down explicitly p ( Z f | Z p ) , which is a Gaussian with
[Equation]
The linear dynamics in latent space are completely deﬁned by the covariance matrix Σ .
Large predictive information implies low conditional entropy H ( Z f | Z p ) .
Let ( Z p , Z f ) generate ( X p , X f ) with a stochastic decoder g
( X | Z ) .
We only observe X and would like to infer the latent Z by maximizing the marginal likelihood of X .
Taking a VAE approach, we parameterize a stochastic encoder e
( Z | X ) for the approximate posterior, and derive a lower bound for the maximum likelihood objective.
Different from standard VAE, here we would not want to parameterize the prior to be a simple Gaussian, in which case the Z p and Z f are independent and have zero mutual information.
Instead we encourage the additional structure of high predictive information for the prior.
This gives us an over all objective as follows:
[Equation]
where ˆ p ( X ) is the empirical distribution over training data, the ﬁrst term corresponds to the reconstruction loss, the second term measures the KL divergence between approximate posterior and the prior, and the last term is the PI deﬁned in (2).
The challenge is how to parameterize the covariance Σ .
We ﬁnd that simply parameterizing it as a positive deﬁnite matrix, e.g., Σ = AA T , does not work well in our experiments, presumably because there is too much ﬂexibility with such a formulation.
What we ﬁnd to work better is the pseudo-input technique discussed in VampPrior (Tomczak & Welling, 2018): given a set of pseudo-sequences X ∗ which are learnable parameters (initialized with real training sequences), we compute the sample covariance from e ( Z | X ∗ ) as Σ .
This approach yields an overall objective very similar to (5), with the beneﬁt of a well-deﬁned generative model (and the Gaussian assumption being perfectly satisﬁed), which allows us to borrow learning/inference techniques developed in the VAE framework.
For example, masking the input for the encoder can be seen as amortized inference regularization (Shu et al., 2018).
We show experimental results on this probabilistic DAPC in Appendix B and C.
In general, probabilistic
DAPC performs similarly to the deterministic counterpart, though the training process is more time and memory intensive.
On the other hand, these empirical results show that deviating from the
Gaussian assumption, as is the case for deterministic DAPC, does not cause signiﬁcant issues for representation learning in practice if proper regularization is applied.
Related to this interpretations are VAE-base sequential models (Chung et al., 2015; Hsu et al., 2017;
Li & Mandt, 2018; Anonymous, 2021) that also use reconstruction and enforce different structures/dynamics in the latent space.
Most of them are designed for the purpose of generating high quality sequence data, while the quality of their latent representation are mostly not shown for downstream tasks.
R ELATED W ORK
Mutual information (MI) maximization is a principled approach for representation learning (Bell &
Sejnowski, 1995), where the objective is to maximize the MI estimate between learnt representations and inputs.
The currently dominant approach for estimating MI is based on contrastive learning.
For sequence data, CPC (Oord et al., 2018) uses representations at current time as a classiﬁer to discriminate inputs of nearby frames (positive samples) from inputs of far-away steps or inputs from other sequences (negative samples) with a cross-entropy loss; this leads to the noise-contrastive estimation (NCE, Gutmann & Hyv ¨ arinen, 2010).
Deep InfoMax (DIM, (Hjelm et al., 2018)) generalizes the NCE estimator with a few variants, and proposes to maximize MI between global summary features and local features from intermediate layers (rather than the inputs as in CPC).
SimCLR (Chen et al., 2020) extends the contrastive loss to use a nonlinear transformation of the representation (than the representation itself) as classiﬁer for measuring MI.
Contrastive Multiview Coding (Tian et al., 2019) generalizes the contrastive learning frame to multiple views.
Momentum Contrast (He et al., 2020) saves memory with a dynamic dictionary and momentum encoder.
Meanwhile, there have been concerns about the contrastive learning framework.
One concern is that postive and negative sample selection is sometimes time and memory consuming.
To address this issue, BYOL (Grill et al., 2020) proposes to get rid of negative samples by learning a target network in an online fashion and gradually bootstrapping the latent space.
Another concern is regarding the
MI estimation.
Though contrastive learning has an MI backbone, Tschannen et al.
(2019) suggests that the inductive bias of the feature extractor and parametrization of estimators might contribute more than the MI estimate itself.
Ozair et al.
; McAllester & Stratos (2020) raise the concern that the MI lower bound used by contrastive learning might be too loose, and propose to use an estimate based on Wasserstein distance.
Unlike prior work, our principle for sequence representation learning is to maximize the MI between past and future latent representations, rather than the MI between representations and inputs (or shallow features of inputs).
Partially motivated by the above concerns, our mutual information estimate requires no sampling and is exact for Gaussian random variables.
To keep useful information from input, we use a masked reconstruction loss which has been effective for sequence data (text and speech), with an intuition resembling that of denoising autoencoders (Vincent et al., 2010). [Figure]
Note that by the data processing inequality, methods that maximize mutual information between current representation and future inputs also implicitly maximizes an upper bound of mutual information between high level representations, since MI ( Z past , Z future ) ≤ MI ( Z past , X future ) .
Our method explicitly maximizes the mutual information between high level representations itself, while having another regularization term (masked reconstruction) that maximizes information between current input and current representations.
Our results indicate that explicitly modeling the trade-off between the two can be advantageous.
In the audio domain where we will demonstrate the applicability of our method, there has been signiﬁcant interest in representation learning for reducing the need for supervised data.
Both contrastive learning based (Schneider et al., 2019; Baevski et al., 2019; Jiang et al., 2019) and reconstructionbased (Chorowski et al., 2019; Chung et al., 2019a; Song et al., 2019; Wang et al., 2020; Chung &
Glass, 2020; Ling et al., 2020; Liu et al., 2020) methods have been studied, as well as methods that combine multiple tasks (Pascual et al., 2019a; Ravanelli et al., 2020).
Our work promotes the use of a different MI estimate and combines different intuitions synergistically.
E XPERIMENTS 4. N OISY L ORENZ A TTRACTOR
Lorenz attractor (“butterﬂy effect”, see Appendix B) is a 3D time series depicting a chaotic system (Pchelintsev, 2014), as visualized in Figure 2.
We design a challenging dimension reduction task for recovering the Lorenz attractor from high dimensional noisy measurements.
We ﬁrst lift the 3D clean signals to 30D with a neural network of 2 hidden layers, each with 128 elu units (Clevert et al., 2015).
This lifting network has weights and biases drawn randomly from N (0 , 0 . 2) .
In addition, we corrupt the 30D lifted signals with white noise to obtain three different signal-to-noise ratio (SNR) levels, 0.3, 1.0, and 5.0, and use the noisy 30D measurements (see bottom left of Figure 2) as input for representation learning methods to recover the true 3D dynamics.
We generate a long Lorenz attractor trajectory using the governing differential equations, and chunk the trajectory into segments of 500 time steps.
We use 250, 25 and 25 segments for training, validation, test splits respectively.
After model selection on the validation split based on the R 2 regression score which measures the similarity between recovered and ground truth trajectories, our model is applied to the test split.
The optimal model uses T = 4 , s = 0 , α = 0 , γ = 0 . 1 , and β = 0 . 1 is set to balance the importance of PI and the reconstruction error.
We compare DAPC to representative unsupervised methods including DCA (Clark et al., 2019),
CPC (Oord et al., 2018), pure PI learning which corresponds to DAPC with β = 0 , and masked reconstruction (MR, Wang et al., 2020) which corresponds to DAPC without the PI term.
Except for DCA which corresponds to maximizing PI with a linear feedforward network, the other methods use bidirectional GRUs (Chung et al., 2014) for mapping the inputs into feature space (although uniGRU performs similarly well).
A feedforward DNN is used for reconstruction in MR and DAPC.
We show the latent representations of different methods in Figure 2 (right panel).
DCA fails completely since its encoder network has limited capacity to invert the nonlinear lifting process.
CPC is able to recover the 2 lobes, but the recovered trajectory is chaotic.
Maximizing PI alone largely ignores the global structure of the data.
MR is able to produce smoother dynamics for high SNRs, but its performance degrades quickly in the more noisy scenarios.
DAPC recovers a latent representation which has overall similar shapes to the ground truth 3D Lorenz attractor, and exhibits smooth dynamics enforced by the PI term.
In Appendix B, we provide the R 2 scores for different methods.
These results quantitatively demonstrate the advantage of DAPC across different noise levels.
4. F ORECASTING WITH L INEAR R EGRESSION
We then demonstrate the predictive power of learnt representations in downstream forecasting tasks on 3 real-world datasets used by Clark et al.
(2019), involving multi-city temperature time series data (Beniaguev, 2017), dorsal hippocampus study (Glaser et al., 2020), and motor cortex (O’Doherty et al., 2018).
For each model, unsupervised representation learning is performed on the training set with a uni-directional GRU, which prevents information leakage from the future.
After that, we freeze the model and use it as a feature extractor.
The representations at each time step are used as inputs for predicting the target at a future time step.
As an example, we can extract a representation for today’s weather based on past weather only (as the encoder is uni-directional), and use it to predict future temperature which is lag days away (a larger lag generally leads to a more difﬁcult forecasting task).
Following Clark et al.
(2019), the predictor from the extracted feature space to the target is a linear mapping, trained on samples of paired current feature and future target, using a least squares loss.
We use the same feature dimensionality as in their work, for each dataset.
These forecasting tasks are evaluated by the R 2 regression score, which measures the linear predictability.
More details can be found in Appendix C.
Besides DCA, CPC and MR, we further include PCA and SFA (Wiskott & Sejnowski, 2002) (similar to DCA with T = 1 for PI estimation), which are commonly used linear dimension reduction methods in these ﬁelds.
PCA is served as the baseline and we report R 2 score improvements from other methods.
Figure 3 gives the performances of different methods for the three datasets, with three different lag s (the number of time steps between current and future for the forecasting task): 5, 10, and 15.
DAPC consistently outperforms the other methods.
4. P RETRAINING FOR A UTOMATIC S PEECH R ECOGNITION (ASR)
A prominent usage of representation learning in speech processing is to pretrain the acoustic model with an unsupervised objective, so that the resulting network parameters serve as a good initialization for the supervised training phase using labeled examples (Hinton et al., 2012).
As supervised ASR techniques have been improved signiﬁcant over the years, recent works start to focus on pre-training with large amounts of unlabeled audio data, followed by ﬁnetuning on much smaller amounts of supervised data, so as to reduce the cost of human annotation.
We demonstrate different methods on two commonly used speech corpora for this setup: Wall Street
Journal (Paul & Baker, 1992) and LibriSpeech (Panayotov et al., 2015).
For WSJ, we pretrain on si284 partition (81 hours), and ﬁnetune on si84 partition (15 hours) or the si284 partition itself.
For
Librispeech, we pretrain on the train 960 partition (960 hours) and ﬁnetune on the train cleanpartition (100 hours).
Standard dev and test splits for each corpus are used for validation and testing. [Table]
[Table] (Liu & Huang, 2020) [Table] In the experiments, we largely adopt the transformers-based recipe from ESPnet (Watanabe et al., 2018), as detailed in Karita et al.
(2019), for supervised ﬁnetuning.
We provide the details regarding model architecture and data augmentation in Appendix D. Note that we have spent effort in building strong ASR systems, so that our baseline (without pretraining) already achieves low WERs and improving over it is non-trivial.
This can be seen from the results table where our baseline is often stronger than the best performance from other works.
For pretraining, we pretrain an encoder oftransformer layers, which will be used to initialize the ﬁrst 14 layers of ASR model.
For masked reconstruction, we use 2 frequency masks as in ﬁnetuning, but found more time masks can improve pretraining performance.
We set the number of time masks to 4 for WSJ, and 8 for LibriSpeech which has longer utterances on average.
The hyperparameters we tune include T, s, α, β and γ from our learning objective.
We select hyperparameters which give the best dev set WER, and report the corresponding test set WER.
In the end, we use T = 4 for estimating the PI term, γ = 0 .
05 , β = 0 . 005 and set s = 2 for WSJ and s = 1 for LibriSpeech if we use shifted reconstruction.
Since the pretraining objective is a proxy for extracting the structure of data and not fully aligned with supervised learning, we also tune the number of pretraining epochs, which is set to 5 for WSJ and 1 for LibriSpeech.
Other parameters will be shown in the ablation studies presented in Appendix D.
We perform an ablation study for the effect of different variants of DAPC (MR+PI) on the WSJ dataset, and give dev/test WERs in Table 1.
We tune the hyperparameters for multi-scale PI and shifted reconstruction for the 15 hours ﬁnetuning setup, and observe that each techniques leads to improvement over DAPC, while combining them delivers the best performance.
The same hyperparameters are used for the 81 hours ﬁnetuning setup, and we ﬁnd that with more supervised training data, the baseline without pretraining obtains much lower WERs, and pure masked reconstruction only slightly improves over the baseline, while the strongest DAPC variant still achieves 7.9% and 11.7% relative improvements on dev93 and eval92 respectively.
In Table 2, we provide a more thorough comparison with other representation learning methods on the LibriSpeech dataset.
We compare with CPC-type mutual information learning methods including wav2vec (Schneider et al., 2019), and vq-wav2vec which performs CPC-type learning with discrete tokens followed by BERT-stype learning (Baevski et al., 2019).
We also compare with two reconstruction-type learning approaches DeCoAR (Ling et al., 2020) and TERA (Liu et al., 2020); comparisons with more methods are given in Appendix D. Observe that DAPC and its variant achieve lower WER than MR: though our baseline is strong, DAPC still reduces WER by 8%, while MR only improves by 1.76%.
This shows the beneﬁt of PI-based learning in addition to masked reconstruction.
While our method does not outperform discrete BERT+vq-wav2vec, the size of our DAPC-based models (30M parameters) is only 1/5 of that of discrete BERT+vq-wav2vec (150M parameters).
In Appendix D, we provide additional experimental results where the acoustic model targets are subwords.
DAPC achieves 15% relative improvement over the baseline (without pretraining), showing that our method is generally effective for different types of ASR systems.
C ONCLUSIONS
In this work, we proposed a novel representation learning method, DAPC, for sequence data.
Our learnt latent features capture the essential dynamics of the underlying data, contain rich information of both input observations and context states, and are shown to be useful in a variety of tasks.
As future work, we may investigate other predictive information estimators that further alleviate the Gaussian assumption, and are accurate and efﬁcient to compute.
On the other hand, advanced variational inference techniques may be applied to the probabilistic version of DAPC to boost the performance.
DAPC provides a general alternative for mutual information-based learning for sequence data and we may investigate its potential usage in other domains such as NLP, biology, etc.
A D ERIVATION OF P REDICTIVE I NFORMATION
In this section, we give a self-contained derivation of the predictive information.
A multivariate Gaussian random variable X ∈ R N has the following PDF:
[Equation]
where µ is the mean and Σ is the covariance matrix with Σ
= E
[( X − E [ X ])( X − E
[ X ])] .
From the deﬁnition of entropy
[Equation]
we can derive the entropy formula for multivariate Gaussian
[Equation]
Consider the joint Gaussian distribution
[Equation]
where X ∈ R p , and Y ∈ R q .
We can plug in the entropy in (6) and obtained
[Equation]
For a latent sequence Z = { z 1 , z 2 , ... } where z i ∈ R
d , we deﬁne Z pastt = { z t − T +1 , ..., z t } , and
Z futuret = { z t +1 , ..., z t + T } .
Based on our stationarity assumption, all the length- 2 T windows of states within the range are drawn from the same Gaussian distribution with covariance Σ 2 T ( Z ) , and similarly for all the length- T windows.
As a result and under the stationary assumption, H ( Z pastt ) =
H ( Z futuret ) = 12 ln((2 πe )
Td | Σ T ( Z ) | ) , H ( Z pastt , Z futuret ) = 12 ln((2 πe ) 2 Td | Σ 2 T ( Z ) | ) , and
I T
=
MI
( Z past t , Z future t
[Equation]
futuret
[Equation]
futuret
[Equation]
The predictive information I T only depend on T but not speciﬁc time index t .
B L ORENZ A TTRACTOR
The Lorenz attractor system (also called “butterﬂy effect”) is generated by the following differential equations: (Strogatz, 2018; Clark et al., 2019):
[Equation]
where ( x, y, z ) are the 3D coordinates, and we use σ = 10 , β = 8 / 3 , ρ = 28 .
The integration step for solving the system of equations is 5 × 10 − 3 .
We lift the 3D trajectory into 30D using a nonlinear neural network with 2 hidden layers, each with 128 neurons and the Exponential Linear Unit (Clevert et al., 2015).
We further perturb the 30D trajectory with additive Gaussian noise to obtain datasets of three different Signal-to-Noise Ratios (SNRs, ratio between the power of signal and the power of noise): 0.3, 1.0, and 5.0.
The smaller the
SNR is, the more challenging it is to recover the clean 3D trajectory (see the left panel of Figurefor the comparison between noisy 30D trajectory and the clean 30D trajectory).
[Table]
[Figure]
[Table]
[Figure]
We compare both deterministic and probabilistic DAPC against representative unsupervised methods including DCA (Clark et al., 2019), CPC (Oord et al., 2018), pure PI learning which corresponds to DAPC with β = 0 , and masked reconstruction (MR) (Wang et al., 2020).
Except for DCA which correspond to maximizing PI with a linear orthogonal feedforward net, the other methods use bidirectional GRU (Chung et al., 2014) for mapping the inputs into feature space (although uni-GRU performs similarly well).
A feedforward DNN is used for reconstruction in MR and DAPC.
More speciﬁcally, CPC, MR, DAPC all use the bidirectional GRU where the learning rate is 0.001, and dropout rate is 0.7.
Our GRU has 4 encoding layers with hidden size 256. The batch size is (20, 500, 30). For CPC, the temporal lag k=4. For DAPC, β = 0 . 1 , T = 4 , s = 0 , α = 0 , γ = 0 . 1 .
For masked reconstruction, we use at most 2 masks on the frequency axis with width up to 5, and at most 2 mask on the time axis with width up to 40.
The DNN decoder has 3 hidden layers with size 512.
DCA’s setup is completely adopted from Clark et al.
The same architectures are used in the forecasting tasks in Appendix C.
Figure 4 provides qualitatively results for the recovered 3D trajectories by different methods (Figure 2 in the main text contains a subset of the results shown here).
Observe that DCA fails in this scenario since its feature extraction network has limited capacity to invert the nonlinear lifting process.
CPC is able to recover the 2 lobes, but the recovered signals are chaotic.
Masked reconstruction is able to produce smoother dynamics for high SNRs, but its performance degrades quickly in the more noisy scenarios.
Both deterministic and probabilistic DAPC recover a latent representation which has overall similar shapes to the ground truth 3D Lorenz attractor, and exhibits smooth dynamics enforced by the PI term.
We quantitatively measure the recovery performance with the R 2 score, which is deﬁned as coefﬁcient of determination.
R 2 score normally ranges from 0 to 1 where 1 means the perfect ﬁt.
Negative scores indicate that the model ﬁts the data worse than a horizontal hyperplane.
The Rresults are given in Table 3.
Our results quantitatively demonstrate the clear advantage of DAPC across different noise levels. [Table] [Table] [Figure]
[Table] We also give an ablation study on several components of DAPC in Table 4, where we attempt full reconstruction without masking, masked reconstruction with unidirectional encoder uni-GRU [Table] and the regular setup (masked reconstruction + bi-GRU).
Using full reconstruction yields worse results than using masked reconstruction at all noise levels, while uni-GRU degrades the performance less.
We show in Table 5 and Figure 5 how PI can improve both full reconstruction and masked reconstruction.
In Table 5, when SNR=0.3, PI can greatly boost the performance of full reconstruction.
We also tuned the temporal lag parameter k w.r.t. both quantitative and qualitative results (Tableand Figure 6). CPC performance starts to deteriorate after k=8, while k=4, 6, 8 have similar results. Based on the R 2 scores, we select k=4 as our ﬁnal temporal lag.
Similar tuning is also performed for CPC on the downstream forecasting experiments in Appendix C.
[Figure]
C D OWNSTREAM F ORECASTING WITH L INEAR R EGRESSION
For each of the three datasets used in downstream forecasting tasks, we divide the original dataset into training/validation/testing splits.
Unsupervised representation learning is performed on the training split, validated on the validation split and learn an encoder e ( · ) .
Denote the test sequence X = { x 1 , x 2 , ..., x L } .
The learnt e ( · ) transforms X into a sequence
Z = { z 1 , z 2 , ..., x L } of the same length.
Note that we use uni-directional GRU for representation learning so that and no future information is leaked.
For the forecasting tasks, z i will be used to predict a target y i which corresponds to an event of interest.
For the multi-city temperature dataset (Beniaguev, 2017), y i represents future multi-city temperatures, i.e., y i
=
x i + lag with lag > 0 .
For the hippocampus study (Glaser et al., 2020), x i is the multi-neuronal spiking activity of 55 single units recorded in rat hippocampal CA1 and y
i is a future location of the rat.
In the motor cortex dataset (O’Doherty et al., 2018), x i ’s are collected from multi-neuronal spiking activity ofsingle units recorded in monkey primary motor cortex (M1), and y i ’s are future behavior variables such as cursor kinematics.
The problems tend to be more challenging with larger lag .
[Equation]
8: On temperature dataset, we analyze the performances of different latent dimensions (dimension of z i ): 5, 10, 15, 20.
∆ R 2 corresponds to R 2 score improvement over PCA.
The performance for forecasting tasks is measured by the linear predictability from z i to y i .
Specifically, we solve the linear regression problem with inputs being the ( z i , y i ) pairs, and measure the R 2 score between the prediction and ground-truth target.
We use the R 2 score from the PCA projection as a baseline, and provide the improvements over PCA obtained by different representation learning methods in Figure 7 (so PCA’s ∆ R 2 is 0 ).
Both deterministic and probabilistic
DAPC consistently outperform other methods across all three datasets, with the deterministic version slightly outperforming the probabilistic counterpart.
Additionally, we provide sensitivity of the latent dimensionality for all methods in Figure 8, and DAPC outperforms others consistently across different dimensionalities.
Table 7 demonstrates the power of adding PI to full or masked reconstruction.
PI can improve both full reconstruction (4.09% in average) and masked reconstruction (17.22% in average), which shows
PI’s ability of generalization to different reconstruction methods.
[Table]
D A DDITIONAL RESULTS FOR A UTOMATIC S PEECH R ECOGNITION (ASR)
The acoustic model is trained with a multi-task objective (Watanabe et al., 2017) which combines attention Chorowski et al.
; Chan et al.
(2016) and CTC (Graves et al., 2006) losses, for predicting the output character sequence.
We extract 80D fbank features plus 3D pitch features from audio, with a frame size of 25ms and hop size of 10ms.
Every 3 consecutive frames are stacked to obtain the input sequence to acoustic model.
During ASR ﬁnetuning, the encoder shared by both attention and CTC consists of 14 transformer layers for WSJ and 16 layers for LibriSpeech, while the decoder consists of 6 transformer layers.
All attention operations use 4 heads of 64 dimensions each, and the output of multi-head attention goes through a one-hidden-layer position-wise feedforward network of 2048 ReLU units, before it is fed to the next layer.
During ﬁnetuning, we apply
SpecAugment (Park et al., 2019) to reduce overﬁtting, with max time warp set to 5 (frames), two frequency masks of widths up to 30 frequency bins, and two time masks of widths up to 40 frames.
We use the Adam optimizer with a warmup schedule for the learning rate.
Weight parameters of last 10 ﬁnetuning epochs is averaged to obtain the ﬁnal model.
For word-level decoding, we use a word RNNLM trained on the language model training data of each corpus, with a vocabulary size of 65K for WSJ, and 200K for LibriSpeech.
We use the lookahead scores derived from word RNNLM during beam search, for selecting promising character tokens at each step, as done by Hori et al.
A beam size of 20 is used for decoding.
In Table 8, we provide a thorough comparison with other representation learning methods on the
LibriSpeech dataset.
We compare with CPC-type mutual information learning methods including wav2vec (Schneider et al., 2019), and vq-wav2vec which performs CPC-type learning with discrete tokens followed by BERT-stype learning (Baevski et al., 2019).
We also compare with two reconstruction-type learning approaches DeCoAR (Ling et al., 2020) and TERA (Liu et al., 2020); comparisons with more methods are given in Appendix D. Note that TERA is quite similar to MR
(Wang et al., 2020) in performing masked reconstruction, although it adds recurrent layers to transformer-learned representations for ﬁnetuning, while our implementation of MR uses a pure transformer-based architecture throughout; we think the advantage of MR over TERA mainly comes from the acoustic model (attention vs. CTC) and a stronger language model (RNNLM vs. n-gram). [Table] [Table]
[Table] Observe that DAPC and its variant achieve lower WER than MR: though our baseline is strong,
DAPC still reduces WER by 8%, while MR only improves by 1.76%.
This shows the beneﬁt of
PI-based learning in addition to masked reconstruction.
While our method does not outperform discrete BERT+vq-wav2vec, the size of our DAPC-based models (30M parameters) is only 1/5 of that of discrete BERT+vq-wav2vec (150M parameters).
Compared to Table 2, here we include more comparison with recent related works.
Futhermore, to show our improvement is robust to details of ASR recipe, we include the comparison between baseline (without pretraining), MR, and DAPC when the ASR recipe uses 5000 unigrams as token set and decodes with token-level RNNLM.
These results are denoted with “(sub-word)” in Table 8.
The relative merits between methods are consistent with those observed for the character recipe.
DAPC obtains a relative improvement of 15% over the baseline on test clean ( 6 . 81% → 5 . 79% ).
Furthermore, we compare DAPC with other state-of-the-art methods on WSJ in Table 9.
Comparisons among different methods based on their key features are shown in Table 10.