summaries: This paper proposes InfoTS, a method for learning augmentations that improve contrastive learning of time series data. The core contribution is a learnable augmentation strategy that uses Concrete/Gumbel-Softmax distributions. The paper shows empirical results on several time-series forecasting and classification benchmarks.

summaries: This paper describes an information-aware approach to representation learning for time series. The formulation focuses on how to obtain effective data augmentations and addresses the underlying problem from information-theoretic viewpoints, leading to the two optimization criteria, namely, high fidelity and high variety. The experimental results on several time series datasets for forecasting and classification show improvements over the methods in comparison. Detailed comments are listed below.

summaries: This is a very good and solid paper that has both empirical elements and theoretical basis. As mentioned in the paper data augmentation in time series is a notoriously difficult problem for the simple reason that it can distort the time series completely - which in contrast to images - humans cannot verify this as it can happen with images. So the contribution of this paper is a new data augmentation approach based on information theory, a meta learning approach and an approach to select optimal data augmentation for contrastive learning.

summaries: The authors propose a contrastive learning framework for time series data, where data augmentations are adaptively being selected, given a fidelity and variety criterion. Additionally to these two criteria, a contrastive learning objective is applied both on local and global level. The model is tested on a time-series forecasting and classification task on multiple datasets.

abstractText: Various contrastive learning approaches have been proposed in recent years and achieve significant empirical success. While effective and prevalent, contrastive learning has been less explored for time series data. A key component of contrastive learning is to select appropriate augmentations imposing some priors to construct feasible positive samples, such that an encoder can be trained to learn robust and discriminative representations. Unlike image and language domains where “desired” augmented samples can be generated with the rule of thumb guided by prefabricated human priors, the ad-hoc manual selection of time series augmentations is hindered by their diverse and human-unrecognizable temporal structures. How to find the desired augmentations of time series data that are meaningful for given contrastive learning tasks and datasets remains an open question. In this work, we address the problem by encouraging both high fidelity and variety based upon information theory. A theoretical analysis leads to the criteria for selecting feasible data augmentations. On top of that, we employ the meta-learning mechanism and propose an information-aware approach, InfoTS, that adaptively selects optimal time series augmentations for contrastive representation learning. The meta-learner and the encoder are jointly optimized in an endto-end manner to avoid sub-optimal solutions. Experiments on various datasets show highly competitive performance with up to 11.4% reduction in MSE on forecasting task and up to 2.8% relative improvement in accuracy on classification task over the leading baselines.

authors:
Shaojie Bai
J Zico Kolter
Vladlen Koltun

title: An empirical evaluation of generic convolutional and recurrent networks for sequence modeling

venue: arXiv preprint arXiv:1803.01271,

title: Tristounet: triplet loss for speaker turn embedding

venue: In ICASSP,

authors:
Defu Cao
Yujing Wang
Juanyong Duan
Ce Zhang
Xia Zhu
Conguri Huang
Yunhai Tong
Bixiong Xu
Jing Bai
Jie Tong

title: Spectral temporal graph neural network for multivariate timeseries forecasting

venue: arXiv preprint arXiv:2103.07719,

authors:
Ting Chen
Simon Kornblith
Mohammad Norouzi
Geoffrey Hinton

title: A simple framework for contrastive learning of visual representations

venue: In ICML,

authors:
Ekin D Cubuk
Barret Zoph
Dandelion Mane
Vijay Vasudevan
Quoc V Le

title: Autoaugment: Learning augmentation strategies from data

venue: In CVPR,

authors:
Hoang Anh Dau
Anthony Bagnall
Kaveh Kamgar
Chin-Chia Michael Yeh
Yan Zhu
Shaghayegh Gharghabi
Chotirat Ann Ratanamahatana
Eamonn Keogh

title: The ucr time series archive

venue: IEEE/CAA Journal of Automatica Sinica,

authors:
Dheeru Dua
Casey Graff

title: UCI machine learning repository, 2017

venue: URL http://archive. ics.uci.edu/ml

authors:
Emadeldeen Eldele
Mohamed Ragab
Zhenghua Chen
Min Wu
Chee Keong Kwoh
Xiaoli Li
Cuntai Guan

title: Time-series representation learning via temporal and contextual contrasting

venue: arXiv preprint arXiv:2106.14112,

authors:
Haoyi Fan
Fengbin Zhang
Yue Gao

title: Self-supervised time series representation learning by inter-intra relational reasoning

venue: arXiv preprint arXiv:2011.13548,

authors:
Jean-Yves Franceschi
Aymeric Dieuleveut
Martin Jaggi

title: Unsupervised scalable representation learning for multivariate time series

venue: arXiv preprint arXiv:1901.10738,

authors:
Ryuichiro Hataya
Jan Zdenek
Kazuki Yoshizoe
Hideki Nakayama

title: Faster autoaugment: Learning augmentation strategies using backpropagation

venue: In ECCV,

authors:
Dan Hendrycks
Kevin Gimpel

title: Gaussian error linear units (gelus)

venue: arXiv preprint arXiv:1606.08415,

authors:
R Devon Hjelm
Alex Fedorov
Samuel Lavoie-Marchildon
Karan Grewal
Phil Bachman
Adam Trischler
Yoshua Bengio

title: Learning deep representations by mutual information estimation and maximization

venue: arXiv preprint arXiv:1808.06670,

authors:
Aapo Hyvarinen
Hiroshi Morioka

title: Unsupervised feature extraction by time-contrastive learning and nonlinear ica

venue: In NIPS,

authors:
Eric Jang
Shixiang Gu
Ben Poole

title: Categorical reparameterization with gumbel-softmax

venue: arXiv preprint arXiv:1611.01144,

authors:
Diederik P Kingma
Jimmy Ba

title: Adam: A method for stochastic optimization

venue: arXiv preprint arXiv:1412.6980,

authors:
Guokun Lai
Wei-Cheng Chang
Yiming Yang
Hanxiao Liu

title: Modeling long-and short-term temporal patterns with deep neural networks

venue: In SIGIR,

authors:
Arthur Le Guennec
Simon Malinowski
Romain Tavenard

title: Data Augmentation for Time Series Classification using Convolutional Neural Networks

venue: In ECML/PKDD Workshop on Advanced Analytics and Learning on Temporal Data,

authors:
Shiyang Li
Xiaoyong Jin
Yao Xuan
Xiyou Zhou
Wenhu Chen
Yu-Xiang Wang
Xifeng Yan

title: Enhancing the locality and breaking the memory bottleneck of transformer on time series forecasting

venue: In NeurIPS,

authors:
Yonggang Li
Guosheng Hu
Yongtao Wang
Timothy Hospedales
Neil M Robertson
Yongxin Yang

title: Dada: Differentiable automatic data augmentation

venue: arXiv preprint arXiv:2003.03780,

authors:
Dongsheng Luo
Wei Cheng
Dongkuan Xu
Wenchao Yu
Bo Zong
Haifeng Chen
Xiang Zhang

title: Parameterized explainer for graph neural network

venue: arXiv preprint arXiv:2011.04573,

authors:
Dongsheng Luo
Wei Cheng
Jingchao Ni
Wenchao Yu
Xuchao Zhang
Bo Zong
Yanchi Liu
Zhengzhang Chen
Dongjin Song
Haifeng Chen

title: Unsupervised document embedding via contrastive augmentation

venue: arXiv preprint arXiv:2103.14542,

authors:
Chris J Maddison
Andriy Mnih
Yee Whye Teh

title: The concrete distribution: A continuous relaxation of discrete random variables

venue: arXiv preprint arXiv:1611.00712,

authors:
Sebastian Nowozin
Botond Cseke
Ryota Tomioka

title: f-gan: Training generative neural samplers using variational divergence minimization, 2016

authors:
Aaron van den Oord
Yazhe Li
Oriol Vinyals

title: Representation learning with contrastive predictive coding

venue: arXiv preprint arXiv:1807.03748,

authors:
Boris N Oreshkin
Dmitri Carpov
Nicolas Chapados
Yoshua Bengio

title: N-beats: Neural basis expansion analysis for interpretable time series forecasting

authors:
F. Pedregosa
G. Varoquaux
A. Gramfort
V. Michel
B. Thirion
O. Grisel
M. Blondel
P. Prettenhofer
R. Weiss
V. Dubourg
J. Vanderplas
A. Passos
D. Cournapeau
M. Brucher
M. Perrot
E. Duchesnay

title: Scikit-learn: Machine learning in Python

authors:
Ben Poole
Sherjil Ozair
Aaron Van Den Oord
Alex Alemi
George Tucker

title: On variational bounds of mutual information

venue: In International Conference on Machine Learning,

authors:
David Salinas
Valentin Flunkert
Jan Gasthaus
Tim Januschowski

title: Deepar: Probabilistic forecasting with autoregressive recurrent networks

venue: International Journal of Forecasting,

authors:
Yonglong Tian
Chen Sun
Ben Poole
Dilip Krishnan
Cordelia Schmid
Phillip Isola

title: What makes for good views for contrastive learning

venue: arXiv preprint arXiv:2005.10243,

authors:
Naftali Tishby
Fernando C Pereira
William Bialek

title: The information bottleneck method

venue: arXiv preprint physics/0004057,

authors:
Sana Tonekaboni
Danny Eytan
Anna Goldenberg

title: Unsupervised representation learning for time series with temporal neighborhood coding

venue: arXiv preprint arXiv:2106.00750,

authors:
Qingsong Wen
Liang Sun
Xiaomin Song
Jingkun Gao
Xue Wang
Huan Xu

title: Time series data augmentation for deep learning: A survey

venue: In AAAI,

authors:
Mark van der Wilk
Matthias Bauer
ST John
James Hensman

title: Learning invariances using the marginal likelihood

authors:
Qizhe Xie
Zihang Dai
Eduard Hovy
Minh-Thang Luong
Quoc V Le

title: Unsupervised data augmentation for consistency training

authors:
Qiang Yang
Xindong Wu

title: 10 challenging problems in data mining research

venue: International Journal of Information Technology & Decision Making,

authors:
Rex Ying
Dylan Bourgeois
Jiaxuan You
Marinka Zitnik
Jure Leskovec

title: Gnnexplainer: Generating explanations for graph neural networks

authors:
Yuning You
Tianlong Chen
Yongduo Sui
Ting Chen
Zhangyang Wang
Yang Shen

title: Graph contrastive learning with augmentations

authors:
Zhihan Yue
Yujing Wang
Juanyong Duan
Tianmeng Yang
Congrui Huang
Bixiong Xu

title: Ts2vec: Towards universal representation of time series

venue: arXiv preprint arXiv:2106.10466,

authors:
George Zerveas
Srideepika Jayaraman
Dhaval Patel
Anuradha Bhamidipaty
Carsten Eickhoff

title: A transformer-based framework for multivariate time series representation learning

venue: In SIGKDD,

authors:
Haoyi Zhou
Shanghang Zhang
Jieqi Peng
Shuai Zhang
Jianxin Li
Hui Xiong
Wancai Zhang

title: Informer: Beyond efficient transformer for long sequence time-series forecasting

heading: 1 INTRODUCTION

text:
Time series data in the real world is high dimensional, unstructured, and complex with unique properties, leading to challenges for data modeling (Yang & Wu, 2006). In addition, without human recognizable patterns, it is much harder to label time series data than images and languages in realworld applications. These labeling limitations hinder deep learning methods, which typically require a huge amount of labeled data for training, been applied on time series data (Eldele et al., 2021). Representation learning learns a fixed-dimension embedding from the original time series that keeps their inherent features. Comparing to the raw time series data, these representations are with better transferability and generalization capacity. To deal with labeling limitations, contrastive learning methods have been widely adopted in various domains for their soaring performance on representation learning, including vision, language, and graph-structured data (Chen et al., 2020; Xie et al., 2019; You et al., 2020). In a nutshell, contrastive learning methods typically train an encoder to map instances to an embedding space where dissimilar (negative) instances are easily distinguishable from similar (positive) ones and model predictions to be invariant to small noise applied to either input examples or hidden states.
Despite being effective and prevalent, contrastive learning has been less explored in the time series domain (Eldele et al., 2021; Franceschi et al., 2019; Fan et al., 2020; Tonekaboni et al., 2021). Existing contrastive learning approaches often involve a specific data augmentation strategy that creates novel and realistic-looking training data without changing its label to construct positive alternatives for any input sample. Their success relies on carefully designed rules of thumb guided by domain expertise. Routinely used data augmentations for contrastive learning are mainly designed for image and language data, such as color distortion, flip, word replacement, and back-translation (Chen et al., 2020; Luo et al., 2021). These augmentation techniques generally do not apply to time series data.
Recently, some researchers propose augmentations for time series to enhance the size and quality of the training data (Wen et al., 2021). For example, Eldele et al. (2021) and Fan et al. (2020) propose to adopt jittering, scaling, and permutation strategies to generate augmented instances. Franceschi et al. (2019) extracts subsequences for data augmentation. In spite of the current progress, existing methods have two main limitations. First, unlike images with human recognizable features, time series data are often associated with inexplicable underlying patterns. Strong augmentation such as permutation may ruin such patterns and consequently, the model will mistake the negative handcrafts for positive ones. While weak augmentation methods such as jittering may generate augmented instances that are too similar to the raw inputs to be informative enough for contrastive learning. On the other hand, time series datasets from different domains may have diverse nature. Adapting a universal data augmentation method, such as subsequence in Xie et al. (2019), for all datasets and tasks leads to sub-optimal performances. Other works follow empirical rules to select suitable augmentations from expensive trial-and-error. Akin to hand-crafting features, hand-picking choices of data augmentations are undesirable from the learning perspective. The diversity and heterogeneity of real-life time series data further hinder these methods away from wide applicability.
To address the challenges, we first introduce the criteria for selecting good data augmentations in contrastive learning. Data augmentation benefits generalizable, transferable, and robust representation learning by correctly extrapolating the input training space to a larger region (Wilk et al., 2018). The positive instances enclose a discriminative zone in which all the data points should be similar to the original instance. The desired data augmentations for contrastive representation learning should have both high fidelity and high variety. High fidelity encourages the augmented data to maintain the semantic identity that is invariant to transformations (Wilk et al., 2018). For example, if the downstream task is classification, then the generated augmentations of inputs should be class-preserving. Meanwhile, generating augmented samples with high variety benefits representation learning by increasing the generalization capacity (Chen et al., 2020). From the motivation, we theoretically analyze the information flows in data augmentations based upon information theory and derive the criteria for selecting desired time series augmentations. Due to the inexplicability in practical time series data, we assume that the semantic identity is presented by the target in the downstream task. Thus, high fidelity can be achieved by maximizing the mutual information between the downstream label and the augmented data. A one-hot pseudo label is assigned to each instance in the unsupervised setting when downstream labels are unavailable. These pseudo labels encourage augmentations of different instances to be distinguishable from each other. We show that data augmentations preserving these pseudo labels can add new information without decreasing the fidelity. Concurrently, we maximize the entropy of augmented data conditional on the original instances to increase the variety of data augmentations.
Based on the derived criteria, we propose an adaptive data augmentation method, InfoTS (as shown in Figure 1), by employing a meta-learning mechanism to avoid ad-hoc choices or painstakingly trial-and-error tuning. Specifically, we utilize a meta-network to learn the augmentation prior in tandem with contrastive learning. The meta-learner automatically selects optimal augmentations from candidate augmentations to generate feasible positive samples. Along with random sampled negative instances, augmented instances are then fed into a time series encoder to learn representations in a contrastive manner. With a reparameterization trick, the meta-network can be efficiently optimized with back-propagation based upon the proposed criteria. Therefore, the meta-network can automatically select data augmentations in a per dataset and per learning task manner without resorting to expert knowledge or tedious downstream validation. Our main contributions include:
• We propose criteria to guide the selection of data augmentations for contrastive time series representation learning without prefabricated knowledge.
• We propose a meta-learning based method to automatically select feasible data augmentations for different time series datasets, which can be efficiently optimized with backpropagation.
• We empirically verify the effectiveness of the proposed criteria to find optimal data augmentations. Extensive experiments demonstrate that InfoTS can achieve highly competitive performance with up to 11.4% reduction in MSE on forecasting task and up to 2.8% relative improvement in accuracy on classification task over the leading baselines.

heading: 2 METHODOLOGY

heading: 2.1 NOTATIONS AND PROBLEM DEFINITION

text: A time series instance x has dimension T × F , where T is the length of sequence and F is the dimension of features. Given a set of time series instances X, we aim to learn an encoder fθ(x) that maps each instance x to a fixed-length vector z ∈ RD, where θ is the learnable parameters of the encoder network and D is the dimension of representation vectors. In semi-supervised settings, each instance x in the labelled set XL ⊆ X is associated with a label y for the downstream task. Specially, XL = X holds in the fully supervised setting. In the work, we use the Sans-serif style lowercase letters, such as x, to denote random time series variables and italic lowercase letters, such as x, for sampled instances.

heading: 2.2 WHAT ARE GOOD AUGMENTATIONS FOR CONTRASTIVE LEARNING?

text:
The goal of data augmentation for contrastive learning is to create realistically rational instances that maintain semantics through different transformation approaches. Unlike instances in vision and language domains, the underlying semantics of time series data is not recognizable to human, making it hard, if not impossible, to include human knowledge to data augmentation for time series data. For example, rotating an image will not change its content or the label. While permuting a time series instance may ruin its signal patterns and generates a meaningless time series instance. In addition, the tremendous heterogeneity of real-life time series datasets further makes selections based on trial-and-errors impractical. Although multiple data augmentation methods have been proposed for time series data, there is less discussion on what is a good augmentation that is meaningful for a given learning task and dataset without prefabricated human priors. From our perspective, ideal data augmentations for contrastive representation should keep high fidelity, high variety, and adaptive to different datasets. The illustration and examples are shown in Figure 2.
High Fidelity. Augmentations with high fidelity maintain the semantic identity that is invariant to transformations. Considering the inexplicability in practical time series data, it is challenging to visually check the fidelity of augmentations. Thus, we assume that the semantic identity of a time series instance is presented by its label in the downstream task, which might be either available or unavailable during the training period. Here, we start our analysis from the supervised case and will extend it to the unsupervised case later. Based on the information theory (Tishby et al., 2000), we define the objective that keeps high fidelity as the large mutual information (MI) between augmentation v and the label y, i.e., MI(v; y).
We consider augmentation v as a probabilistic function of x and a random variable , that v = g(x; ). From the definition of mutual information, we have MI(v; y) = H(y)−H(y|v), where H(y) is the (Shannon) entropy of y and H(y|v) is the entropy of y conditioned on augmentation v. Since H(y) is irrelevant to data augmentations, the objective is equivalent to minimizing the conditional entropy H(y|v). Considering the efficient optimization, we follow Ying et al. (2019) and Luo et al. (2020) to approximate it with cross-entropy between y and ŷ, where ŷ is the prediction with augmentation v as the input and calculated via
v = g(x; ) z = fθ(v) ŷ = hw(z), (1)
where z is the representation and hw(·) is a prediction projector parameterized by w. The prediction projector is optimized by the classification objective. Then, the objective of high fidelity for
supervised or semi-supervised cases is to minimize
CE(y; ŷ) = − C∑ c=1 P (y = c) logP (ŷ = c), (2)
where C is the number of labels. In the unsupervised settings where y is unavailable, one-hot encoding ys ∈ R|X| is utilized as the pseudo label to replace y in Eq. (2). The motivation is that augmented instances are still distinguishable from other instances with the classifier. We theoretically show that augmentations that preserving pseudo labels have the following properties.
Property 1 (Preserving Fidelity). If augmentation v preserves the one-hot encoding pseudo label, the mutual information between v and the downstream task label y (although not visible to training) is equivalent to that between raw input x and y, i.e., MI(v; y) = MI(x; y).
Property 2 (Adding New Information). By preserving the one-hot encoding pseudo label, augmentation v contains new information comparing to the raw input x, i.e., H(v) ≥ H(x). Detailed proofs are shown in the Appendix A. These properties show that in the unsupervised setting, preserving the one-hot encoding pseudo label guarantees that the generated augmentations will not decrease the fidelity, regardless of the downstream tasks and variances inherent in the augmentations. Concurrently, it may introduce new information for contrastive learning.
Since the number of labels is equal to the number of instances in dataset X in an unsupervised case, direct optimization of Eq. (2) is inefficient and unscalable. Thus, we further relax it by approximating y with the batch-wise one-hot encoding yB , which decreases the number of labels C from the dataset size to the batch size.
High Variety. Sufficient variances in augmentations improve the generalization capacity of contrastive learning models. In the information theory, the uncertainty inherent in the random variable’s possible outcomes is described by its entropy. Considering that augmented instances are generated based on the raw input x, we maximize the entropy of v conditioned on x, H(v|x), to maintain a high variety of augmentations. From the definition of conditional entropy, we have
H(v|x) = H(v)−MI(v; x). (3)
We dismiss the first part since the unconstrained entropy of v can be dominated by meaningless noise. Considering the continuity of both v and x, we adopt a mutual information neural estimator, InfoNCE (Oord et al., 2018; Tian et al., 2020) to approximately compute the mutual information for its practical effectiveness. Other MI estimators, such as Jensen-Shannon (JSD) estimator (Nowozin et al., 2016), normalized temperature-scaled cross-entropy (NT-Xent) (Chen et al., 2020), and leaveone-out (Poole et al., 2019) can also conveniently be the plug-and-play component in our framework. Then, the objective to encourage high variety is to maximize the InfoNCE between v and x:
InfoNCE(v; x) = −Ex [ log
exp(sim(zx, zv))∑ x′∈X exp(sim(zx, zv′))
] , (4)
where v′ is an augmented instance of input instance x′. zx, zv , and zv′ are representations of instance x, v, and v′ respectively. sim(z1, z2) = zT1 z2 is the inner product of vectors z1 and z2.
Criteria. Combining the information aware definition of both high fidelity and variety, we propose the criteria for selecting good augmentations without prior knowledge,
min v −InfoNCE(v; x) + βCE(y;hw(fθ(v))), (5)
where β is a hyper-parameter to achieve the trade-off between fidelity and variety. Note that in the unsupervised settings, y is replaced by one-hot encoding pseudo label..
Relation to Information Bottleneck. Although the formation is similar to information bottleneck in data compression, minp(e|x) MI(x; e) − βMI(e; y), our criteria are different in the following aspects. First, e in the information bottleneck is a representation of input x, while v in Eq.(5) represents the augmented instances. Second, information bottleneck aims to keep minimal and sufficient information for data compression, while our criteria are designed for data augmentations in contrastive learning. Third, in information bottleneck, the compressed representation e is a deterministic function of input x with no variances. MI(e; y) and MI(e; x) are constraint by MI(x; y) and H(x) that MI(e; y) ≤ MI(x; y) and MI(e; x) = H(e), where H(e) is the entropy of e. In our criteria, v is a probabilistic function of input x. As a result, the variances of v makes the augmentation space much larger than the compression representation space in information bottleneck.
Relation to InfoMin. InfoMin is designed based on the information bottleneck that good views should keep minimal and sufficient information from the original input (Tian et al., 2020). Similar to the information bottleneck, InfoMin assumes that augmented views are functions of the input, which heavily constrains the variance of data augmentations. Besides, high fidelity property is dismissed in the unsupervised setting. It works for image datasets due to the availability of human knowledge. However, it may fail to generate reasonable augmentations for time series data.

heading: 2.3 TIME SERIES META-CONTRASTIVE LEARNING

text: We aim to design a learnable augmentation selector that learns to select feasible augmentations in a data-driven manner. With such adaptive data augmentations, the contrastive loss is then used to train the encoder that learns representations from raw time series.

heading: 2.3.1 ARCHITECTURE

text:
The adopted encoder fθ(x) : RT×F → RD consists of two components, a fully connected layer, and a 10-layer dilated CNN module (Franceschi et al., 2019; Yue et al., 2021). To explore the inherent structure of time series, we include both global-wise (instance-level) and local-wise (subsequencelevel) losses in the contrastive learning framework to train the encoder.
Global-wise contrastive loss is designed to capture the instance level relations in a time series dataset. Formally, given a batch of time series instances XB ⊆ X, for each instance x ∈ XB , we generate an augmented instance v with an adaptively selected transformation, which will be introduced later. (x, v) is regarded as a positive pair and other (B−1) combinations {(x, v′)}, where v′ is an augmented instance of x′ and x′ 6= x, are considered as negative pairs. Following Chen et al. (2020); You et al. (2020), we design the global-wise contrastive loss based on InfoNCE (Hjelm et al., 2018). The batch-wise instance-level contrastive loss is
Lg = − 1 |XB | ∑ x∈XB log exp(sim(zx, zv))∑ x′∈XB exp(sim(zx, zv′)) . (6)
Local-wise contrastive loss is proposed to explore the intra-temporal relations in time series. For an augmented instance v of a time series instance x, we first split it into a set of subsequences S, each with length L. For each subsequence s ∈ S, we follow Tonekaboni et al. (2021) to generate a positive pair (s, p) by selecting another subsequence close to it. Non-neighboring samples, N̄s, are adopted to generate negative pairs. Detailed descriptions can be found in Appendix B. Then, the local-wise contrastive loss for an instance x is:
Lcx = − 1 |S| ∑ s∈S log exp(sim(zs, zp)) exp(sim(zs, zp)) + ∑ j∈N̄s exp(sim(zs, zj)) . (7)
Across all instances in a batch, we haveLc = 1|XB | ∑ x∈XB Lcx. Then, the final contrastive objective is: min θ Lg + αLc, (8)
where α is a hyper-parameter to achieve the trade-off between global and local contrastive losses.

heading: 2.3.2 META-LEARNER NETWORK

text:
Previous time series contrastive learning methods (Franceschi et al., 2019; Fan et al., 2020; Eldele et al., 2021; Tonekaboni et al., 2021) generate augmentations with either rule of thumb guided by prefabricated human priors or tedious trial-and-errors, which are designed for specific datasets and learning tasks. In this part, we discuss how to adaptively select the optimal augmentations with a meta-learner network based on the proposed information-aware criteria. We can regard its choice of optimal augmentation as a kind of prior selection. We first choose a set of candidate transformations T, such as jittering and time warping. Each candidate transformation ti ∈ T is associated with an important weight pi ∈ (0, 1), inferring the probability of selecting transformation ti. For an instance x, the augmented instance vi through transformation ti can be computed by:
ai ∼ Bernoulli(pi) vi = (1− ai)x+ aiti(x). (9)
Considering multiple transformations, we pad all vi to be with the same length. Then, the adaptive augmented instance can be achieved by combing candidate ones, v = 1|T| ∑ i vi.
To enable the efficient optimization with gradient-based methods, we approximate discrete Bernoulli processes with binary concrete distributions (Maddison et al., 2016). Specifically, we approximate ai in Eq. (9) with
∼ Uniform(0, 1) ai = σ((log − log(1− ) + log pi
1− pi )/τ), (10)
where σ(·) is the sigmoid function and τ is the temperature controlling the approximation. The rationality of such approximation is given in Appendix A. Moreover, with temperature τ > 0, the gradient ∂v∂pi is well-defined. Therefore, our meta-network is end-to-end differentiable. Detailed algorithm is shown in Appendix B.

heading: 3 RELATED WORK

heading: 3.1 CONTRASTIVE TIME SERIES REPRESENTATION LEARNING

text: Contrastive learning has been utilized widely in representation learning with superior performances in various domains (Chen et al., 2020; Xie et al., 2019; You et al., 2020). Recently, some efforts have been devoted to applying contrastive learning to the time series domain (Oord et al., 2018; Franceschi et al., 2019; Fan et al., 2020; Eldele et al., 2021; Tonekaboni et al., 2021; Yue et al., 2021). Time Contrastive Learning trains a feature extractor with a multinomial logistic regression classifier to discriminate all segments in a time series (Hyvarinen & Morioka, 2016). Franceschi et al. (2019) generates positive and negative pairs based on subsequences. TNC employs a debiased contrastive objective to ensure that in the representation space, signals in the local neighborhood are distinguishable from non-neighboring signals (Tonekaboni et al., 2021). SelfTime adopts multiple hand-crafted augmentations for unsupervised time series contrastive learning by exploring both inter-sample and intra-sample relations (Fan et al., 2020). TS2Vec learns a representation for each time stamp and conducts contrastive learning in a hierarchical way (Yue et al., 2021). However, data augmentations in these methods are either universal or selected by error-and-trail, hindering them away from been widely applied in complex real-life datasets.

heading: 3.2 TIME SERIES FORECASTING

text:
Forecasting is a critical task in time series analysis. Deep learning architectures used in the literature include Recurrent Neural Networks (RNNs) (Salinas et al., 2020; Oreshkin et al., 2019), Convolutional Neural Networks (CNNs) (Bai et al., 2018), Transformers (Li et al., 2019; Zhou et al., 2021), and Graph Neural Networks (GNNs) (Cao et al., 2021). N-BEATS deeply stacks fully-connected layers with backward and forward residual links for univariate times series forecasting (Oreshkin
et al., 2019). TCN utilizes a deep CNN architecture with dilated causal convolutions (Bai et al., 2018). Considering both long-term dependencies and short-term trends in multivariate time series, LSTnet combines both CNNs and RNNS in a unified model (Lai et al., 2018). LogTrans brings the Transformer model to time series forecasting with causal convolution in its attention mechanism Li et al. (2019). Informer further proposes a sparse self-attention mechanism to reduce the time complexity and memory usage (Zhou et al., 2021). StemGNN is a GNN based model that considers the intra-temporal and inter-series correlations simultaneously Cao et al. (2021). Unlike these works, we aim to learn general representations for time series data that can not only be used for forecasting but also other tasks, such as classification. Besides, the proposed framework is compatible with various architectures as encoders.

heading: 3.3 ADAPTIVE DATA AUGMENTATION

text: Data augmentation is an important component in contrastive learning. Existing researches reveal that the choices of optimal augmentation are dependent on downstream tasks and datasets (Chen et al., 2020; Fan et al., 2020). Some researchers have explored adaptive selections of optimal augmentations for contrastive learning in the vision field. AutoAugment automatically searches the combination of translation policies via a reinforcement learning method (Cubuk et al., 2019). Faster-AA improves the searching pipeline for data augmentation using a differentiable policy network (Hataya et al., 2020). DADA further introduces an unbiased gradient estimator for an efficient one-pass optimization strategy (Li et al., 2020). Within contrastive learning frameworks, Tian et al. (2020) applies the Information Bottleneck theory that optimal views should share minimal and sufficient information, to guide the selection of good views for contrastive learning in the vision domain. Considering the inexplicability of time series data, directly applying the InfoMin framework may keep insufficient information during augmentation. Different from Tian et al. (2020), we focus on the time series domain and propose an end-to-end differentiable method to automatically select the optimal augmentations for each dataset.

heading: 4 EXPERIMENTS

text: In this section, we compare InfoTS with SOTA baselines on time series forecasting and classification tasks. We also conduct case studies to show insights into the proposed criteria and meta-learning framework. Detailed experimental setups are shown in Appendix C. Full experimental results and extra experiments, such as parameter sensitivity studies, are presented in Appendix D.

heading: 4.1 TIME SERIES FORECASTING

text:
Time series forecasting aims to predict the future Ly time stamps, with the last Lx observations. We follow Yue et al. (2021) to train a linear model regularized with the L2 norm penalty to make predictions. The output has dimension Ly in the univariate case and Ly × F for the multivariate case, where F is the feature dimension.
Datasets and Baselines. Four benchmark datasets for time series forecasting are adopted, including ETTh1, ETTh2, ETTm1 (Zhou et al., 2021), and the Electricity dataset (Dua & Graff, 2017). These datasets are used in both univariate and multivariate settings. We compare unsupervised InfoTS to the SOTA baselines, including TS2Vec (Yue et al., 2021), Informer (Zhou et al., 2021), StemGNN (Cao et al., 2021), TCN (Bai et al., 2018), LogTrans (Li et al., 2019), LSTnet (Lai et al., 2018), and N-BEATS (Oreshkin et al., 2019). Among these methods, N-BEATS is merely designed for the univariate and StemGNN is for multivariate only. We refer to Yue et al. (2021) to set up baselines for a fair comparison. Standard metrics for a regression problem, Mean Squared Error (MSE), and Mean Absolute Error (MAE) are utilized for evaluation. Evaluation results of univariate time series forecasting are shown in Table 1, while multivariate forecasting results are reported in the Appendix D.4 due to the space limitation.
Performance. As shown in Tabel 1 and Tabel 4, comparison in both univariate and multivariate settings indicates that InfoTS consistently matches or outperforms the leading baselines. Some results of StemGNN are unavailable due to the out-of-memory issue (Yue et al., 2021). Specifically, we have the following observations. TS2Vec, another contrastive learning method with data augmentations, achieves the second-best performance in most cases. The consistent improvement of TS2Vec
over other baselines indicates the effectiveness of contrastive learning for time series representations learning. However, such universal data augmentations may not be the most informative ones to generate positive pairs. Comparing to TS2Vec, InfoTS decreases the average MSE by 11.4%, and the average MAE by 8.6% in the univariate setting. In the multivariate setting, the MSE and MAE decrease by 4.6% and 2.3%, respectively. The reason is that InfoTS can adaptively select the most suitable augmentations in a data-driven manner with high variety and high fidelity. Encoders trained with such informative augmentations learn representations with higher quality.

heading: 4.2 TIME SERIES CLASSIFICATION

text:
Following the previous setting, we evaluate the quality of representations on time series classification in a standard supervised manner (Franceschi et al., 2019; Yue et al., 2021). We train an SVM classifier with a radial basis function kernel on top of representations in the training split and then compare the prediction in the test set.
Datasets and Baselines. Two types of benchmark datasets are used for evaluations. The UCR archive (Dau et al., 2019) is a collection of 128 univariate time series datasets and the UEA archive (Bredin, 2017) consists of 30 multivariate datasets. We compare InfoTS with baselines including TS2Vec (Yue et al., 2021), T-Loss (Franceschi et al., 2019), TS-TCC (Eldele et al., 2021), TST (Zerveas et al., 2021), and DTW (Franceschi et al., 2019). For our methods, InfoTSs, training labels are only used to train the meta-network to select suitable augmentations, and InfoTS is with a purely unsupervised setting for representation learning.
Performance. The results on the UEA datasets are summarized in Table 2. Full results are provided in Appendix D.4. With the ground-truth label guiding the meta-network, InfoTSs substantially outperforms other baselines. On average, it improves the classification accuracy by 2.8% over the best baseline, TS2Vec, with an average rank value 2.133 on all 30 UEA datasets. Under the purely unsupervised setting, InfoTS preserves fidelity by adopting one-hot encoding as the pseudo labels. InfoTS achieves the second best average performance in Table 2, with an average rank value 2.633. Performances on the 128 UCR datasets are shown in Table 9 in Appendix D.4. These datasets are univariate with easily recognized patterns, where data augmentations have marginal or even negative effects (Yue et al., 2021). However, with adaptively selected augmentations for each dataset based on our criteria, InfoTSs and InfoTS still outperform the state-of-the-arts.

heading: 4.3 EVALUATION OF THE CRITERIA

text:
In Section 2.2, we propose criteria of data augmentations for time series based on information theory that good augmentations should have high variety and fidelity. With mutual information neural
estimation and cross-entropy as approximations, we get the criteria in Eq. (5). To empirically verify the effectiveness of the proposed criteria, we adopt the dataset CricketX from the UCR archive (Dau et al., 2019; Fan et al., 2020), and conduct augmentations with different configurations, such as jittering with noise from N (0, 0.3) and cutout 10% time stamps. For each configuration, we calculate the criteria score and the corresponding classification accuracy within the setting in Section 4.2. As shown in Figure 3, in general, accuracy performance is positively related to the proposed criteria in both supervised and unsupervised settings, verifying the correctness of using the criteria as the objective in the meta-network training.

heading: 4.4 ABLATION STUDIES

text: To present deep insights into the proposed method, we verify key components in InfoTS with multiple ablation studies on the Electricity dataset. Results are shown in Table 3. 1) We verify the effectiveness of the adopted contrastive learning objective, Eq. (8), by comparing InfoTS with its two variants. “w/o Local” removes the local contrastive loss between subsequences and “w/o Global” removes the instance-level contrastive loss. The comparison between these variants and InfoTS demonstrates that both instance-level and intra-temporal objectives are important in contrastive learning for time series data. 2) To demonstrate the advantage of adaptive selection of augmentations, we compare InfoTS with variants “Random” and “All”. “Random” randomly selects an augmentation from candidate transformation functions each time and “All” sequentially applies transformations to generate augmented instances. Performance decreases are observed in these variances, verifying the key role of adaptive selection in our method. 3) To show the effects of variety and fidelity objectives in meta-network training, we include two variants, “w/o Fidelity” and “w/o Variety”, which dismiss the fidelity or variety objective, respectively. The comparison between InfoTS and the variants empirically confirms both variety and fidelity are important for data augmentation in contrastive learning. Note that the adopted basic augmentations are manually tuned in a previous paper (Fan et al., 2020) with high qualities for contrastive learning. Including low quality basic augmentations, we conduct more ablation studies are shown in Appendix D.3.

heading: 5 CONCLUSIONS

text:
We propose criteria of data augmentations for time series data based on the information theory that good augmentations should preserve high variety and high fidelity. We approximate the criteria with mutual information neural estimation and cross-entropy estimation. Based on the approximated criteria, we present a meta-learning method to adaptively select optimal augmentations for contrastive representation learning. Comprehensive experiments show that representations produced by our method are high qualified and easy to use in various downstream tasks, such as time series forecasting and classification, with state-of-the-art performances.
STATEMENTS
Reproducibility Statement. The training algorithm of the proposed InfoTS is shown in Appendix B. Detailed experimental settings, including architectures and hyper-parameters settings, are given in Appendix C. Core source codes with an example of usages are attached in the supplementary material. Full version will be made available upon acceptance.

heading: A DETAILED PROOFS

text:
Properties of data augmentations that preserve pseudo labels.
We assume that vi and vj are two augmented instances of inputs xi and xj , respectively. Preserving pseudo labels defined in one-hot encoding requires that the map between variable x and v is one-tomany. Formally, xi 6= xj → vi 6= vj . This can be proved by contradiction. If we have a pair of (i, j) that i 6= j and vi = vj , then we have fθ(vi) = fθ(vj), showing that augmentations cannot preserve pseudo labels.
Property 1 (Preserving Fidelity). If augmentation v preserves the one-hot encoding pseudo label, the mutual information between v and downstream task label y (although not visible to training) is equivalent to that between raw input x and y, i.e., MI(v; y) = MI(x; y).
Proof. From the definition of mutual information, we have MI(v; y) = H(y)−H(y|v) = H(y) + ∑ v,y p(v, y) log p(v, y) p(v)
= H(y) + ∑ x,y ∑ v∈V(x) p(v, y) log p(v, y) p(v) ,
where V(x) is the set of augmented instances of a time series instance x. In the unsupervised setting where the ground-truth label y is unknown, we assume that the augmentation v is a (probabilistic) function of x only. The only qualifier means p(v|x, y) = p(v|x). Since the mapping from x to v is one-to-many. For each v ∈ V(x) we have p(v, y) = p(v, x, y) and p(v) = p(v|x)p(x). Thus, we have
p(v, y)
p(v) =
p(v, x, y) p(v|x)p(x) = p(v|x, y)p(x, y) p(v|x)p(x) = p(v|x)p(x, y) p(v|x)p(x) = p(x, y) p(x)
MI(v; y) = H(y) + ∑ x,y ∑ v∈V(x) p(v, y) log p(x, y) p(x)
= H(y) + ∑ x,y [ ∑ v∈V(x) p(v, y)] log p(x, y) p(x)
= H(y) + ∑ x,y p(x, y) log p(x, y) p(x)
= MI(x; y).
Property 2 (Adding New Information). By preserving the one-hot encoding pseudo label, augmentation v contains new information comparing to the raw input x, i.e., H(v) ≥ H(x). In information theory, entropy describes the amount of information of a random variable. For simplicity, we assume a finite number of augmented instances for each input, and each augmented instance is generated independently. Then, we have p(x) = ∑ v∈V(x) p(v). Then we have that the entropy of variable x is no larger than the entropy of v.
H(x) = − ∑ x p(x) log p(x) = − ∑ x [ ∑ v∈V(x) p(v)] log[ ∑ v∈V(x) p(v)]
= − ∑ x ∑ v∈V(x) p(v) log[ ∑ v∈V(x) p(v)]
≤ − ∑ x ∑ v∈V(x) p(v) log p(v)
= − ∑ v p(v) log p(v) = H(v)
Rationality of approximation of Bernoulli distribution with binary concrete distribution in Eq. (10).
In the binary concrete distribution, parameter τ controls the temperature that achieves the trade-off between binary output and continuous optimization. When τ → 0, we have limτ→0 P (ai = 1) = pi, which is equivalent to the Bernoulli distribution.
Proof. lim τ→0 P (ai = 1) = lim τ→0 P (σ((log − log(1− ) + log pi 1− pi )/τ) = 1)
= P (log − log(1− ) + log pi 1− pi > 0)
= P (log 1− − log 1− pi pi > 0)
= P ( 1− > 1− pi pi )
Since , and pi are both in (0, 1), and function x1−x are monotonically increasing in this region. Thus, we have
lim τ→0
P (ai = 1) = P ( > 1− pi) = pi

heading: B ALGORITHMS

text:
B.1 TRAINING ALGORITHM
The training algorithm of InfoTS under both supervised and unsupervised settings is described in Algorithm 1. We first randomly initiate parameters in the encoder, meta-network, and classifier (line 2). Given a batch of training instances XB ⊆ X, for each candidate transformation ti, we utilize binary concrete distribution to get parameters ai (lines 6-7), which indicates whether the transformation should be applied (line 8). V(i)B denotes the batch of augmented instances generated from transformation function ti. The final augmented instances are generated by adaptively considering all candidates transformations (line 10). Parameters θ in the encoder are updated by minimizing the contrastive objective (lines 11-13). Meta-network is then optimized with information-aware criteria (lines 14-16). Then, classifier hw is optimized with the classification objective (line 17).
Algorithm 1: Algorithm for InfoTS in both supervised and unsupervised settings 1: Input: time series dataset X, the label set Y (supervised setting), a set of candidate transformations T,
hyper-parameters α and β, 2: Initialize the encoder fθ , parameters in meta-network {qi}|T|i=1, and the classifier hw. 3: for each epoch do 4: for each training batch XB ⊆ X do 5: for each transformation ti ∈ T do 6: pi ← σ(qi) 7: ai ← binaryConcrete(pi) 8: V(i)B ← transform each x ∈ XB with Eq. (9) 9: end for
10: Get VB by averaging {V(i)B } |T| i=1. 11: Compute global contrastive loss Lg with Eq. (6) 12: Compute local contrastive loss Lc with Eq. (7) 13: Update parameters θ in the encoder with Eq. (8) 14: Compute fidelity loss with Eq. (2) 15: Compute variety loss with Eq. (4) 16: Update parameters {qi}|T|i=1 in the meta-network with Eq. (5) 17: Update parametersw in the classifier hw with the cross-entropy loss 18: end for 19: end for
B.2 IMPLEMENTATION OF LOCAL-WISE CONTRASTIVE
Local-wise contrastive loss aims to capture the intra-temporal relations in each time series instance. For an augmented instance v, we first split it into multiple subsequences, as shown in Figure 4.
Each subsequence has length L. For each subsequence s, the neighboring subsequences within window size 1 are considered as positive samples. If s locates at the end of v, then we choose the subsequence in front of s as the positive pair p. Otherwise, we choose the subsequence following s instead. Subsequences out of window size 1 are considered as negative samples.

heading: C EXPERIMENTAL SETTINGS

text:
C.1 DATA AUGMENTATIONS
We follow Fan et al. (2020) to set up candidate data augmentations, including jittering, scaling, cutout, time warping, window slicing, window warping and subsequence augmentation (Franceschi et al., 2019). Detailed descriptions are listed as follows.
• Jittering augmentation adds the random noise sampled from a Gaussian distribution N (0, 0.3) to the input time series.
• Scaling augmentation multiplies the input time series by a scaling factor sampled from a Gaussian distribution N (0, 0.5).
• Cutout operation replaces features of 10% randomly sampled time stamps of the input with zeros.
• Time warping random changes the speed of the timeline1. The number of speed changes is 100 and the maximal ratio of max/min speed is 10. If necessary, over-sampling or sampling methods are adopted to ensure the length of the augmented instance is the same as the original one.
• Window slicing randomly crops half the input time series and then linearly interpolates it back to the original length (Le Guennec et al., 2016).
• Window warping first randomly selects 30% of the input time series along the timeline and then warps the time dimension by 0.5 or 2. Finally, we adopt linear interpolation to transform it back to the original length (Le Guennec et al., 2016).
• Subsequence operation random selects a subsequence from the input time series (Yue et al., 2021).
With the first 100 time stamps in the univariate Electricity dataset as an example, we visualize the original time series and the augmented ones in Figure 5.
C.2 HARDWARE AND IMPLEMENTATIONS
All experiments are conducted on a Linux machine with 4 NVIDIA GeForce RTX 2080 Ti GPUs, each with 11GB memory. CUDA version is 10.1 and Driver Version is 418.56. Our method InfoTS is implemented with Python 3.7.7 and Pytorch 1.7.1.
1https://tsaug.readthedocs.io/
C.3 HYPERPARAMETERS
We train and evaluate our methods with the following hyperparameters and configurations.
• Optimizer: Adam optimizer (Kingma & Ba, 2014) with learning rate and decay rates setting to 0.001 and (0.9,0.999), respectively.
• SVM: scikit-learn implementation (Pedregosa et al., 2011) with penalty C ∈ {10i|i ∈ [−4, 4] ∪∞} Franceschi et al. (2019).
• Encoder architecture: We follow Yue et al. (2021) to design the encoder. Specifically, the output dimension of the linear projection layer is set to 64, the same for the number of channels in the following dilated CNN module. In the CNN module, GELU (Hendrycks & Gimpel, 2016) is adopted as the activation function, and the kernel size is set to 3. The dilation is set to 2i in the i-the block.
• Classifier architecture: a fully connected layer that maps the representations to the label is adopted.
• Trade-off hyperparameters: β in Eq. (5) and α in Eq. (8) are searched in [0.1, 0.5, 1.0, 5, 10]. Parameter sensitivity studies are shown in Section D.1.
• Temperature in binary concrete distribution: we follow the practice in (Jang et al., 2016) to adopt the strategy by starting the training with a high temperature and anneal to a small value with a guided schedule.

heading: D MORE EXPERIMENTAL RESULTS

text:
D.1 PARAMETER SENSITIVITY STUDIES
In this part, we adopt the electricity dataset to analyze the effects of two important hyper-parameters in our method InfoTS. The hyperparameter α in Eq. (8) controls the trade-off between local and global contrastive losses when training the encoder. β in Eq. (5) achieves the balance between high variety and high fidelity when training the meta-network. We tune these parameters in range [0.1, 0.5, 1.0, 5, 10] and show the results in Figure 6. These figures show that our method achieves high performance with a wide range of selections, demonstrating the robustness of the proposed method. In general, setting trade-off parameters to 0.5 or 1 achieves good performance.
D.2 CASE STUDY ON SIGNAL DETECTION
In this part, we show the potential usage of InfoTS to detect the informative signals in the time series. We adopt the CricketX dataset as an example for the case study. Subsequence augmentations
on 0-100, 100-200, and 200-300 periods are adopted as candidate transformations. We observe that the one operated on 100-200 period has high fidelity and variety, leading to better accuracy performance, which is consistent with the visualization results in Figure 7.
D.3 MORE ABLATION STUDIES
D.3.1 UPDATING PROCESS OF INFOTS
To show that our InfoTS can adaptively detect the most effective augmentation based on the data distribution, we follow the setting in Section 4.4 and conduct more ablation studies to investigate comprehensively into the proposed model. We compare performances of variants that each applies a single transformation to generate augmented instances in Table 5. From the table, we know that augmentation with subsequence benefits the most for the Electricity dataset. We visualize the weight
updating process of InfoTS in Figure 8, with each line representing the normalized importance score of the corresponding transformation. The weight for subsequence increase with the epoch, showing that InfoTS tends to adopt subsequence as the optimal transformation. Consistency between accuracy performance and weight updating process demonstrates the effectiveness of InfoTS to adaptively select feasible transformations. Besides, as shown in Table 5, InfoTS outperforms the variant that uses subsequence only. This comparison shows that the meta-network learns to consider the combinations, which is better than any (single) candidate augmentation.
D.3.2 INFOTS WITH CONTROLLABLE AUGMENTATIONS.
We observe that in Table 3, the improvement of InfoTS over its two variants without Variety and Fidelity objectives is not significant. The reason is that basic augmentations are adopted from a previous paper (Fan et al., 2020), which have been manually tuned to generate relatively good performances for contrastive learning. To provide deep insights into the proposed criteria and model, in this part, we consider more augmentations with controllable hyper-parameters for the Electricity dataset.
Two groups of augmentations, subsequence augmentations with different lengths and jitter augmentations with different standard deviations, are considered in this part. Subsequence augmentations work on the temporal dimension, and jitter augmentations work on the feature dimension. For the subsequence augmentations, we range the ratio of subsequences r in the range [0.01, 0.99]. The subsequence augmentation with ratio r is denoted by Subr, such as Sub0.01. For the jitter augmentations, the standard deviations are chosen from the range [0.01, 3.0]. The jitter augmentation with standard deviation std is denoted by Jitterstd, such as Jitter0.01.
Intuitively, with r increasing, Subr generates augmented instances with lower variety and higher fidelity. For example, with r = 0.01, Subr generates subsequences that only keep 1% time stamps from the original input, leading to high variety but extremely low fidelity. Similarly, for jitter augmentations, with std increasing, Jitterstd generates augmented instances with higher variety but lower fidelity. We show their forecasting performances in Table 6 and Table 7. We observed that augmentations that only keeps high variety, such as Sub0.01, Jitter3.0 achieve worse performances in terms of MSE and MAE compared to the ones considering both criteria, such as Sub0.3 and Jitter0.3. On the other hand, augmentations that only keep high fidelity, such as Sub0.99, Jitter0.01, and Jitter0.1 cannot lead to satisfactory performances either. Thus, we empirically verify the importance of both variety and fidelity for the selection of optimal data augmentations for contrastive learning. In Fig-

heading: 24 0.262 0.292 0.249 0.278 0.278 0.277 0.252 0.279 0.253 0.279

heading: 48 0.309 0.321 0.300 0.310 0.298 0.306 0.302 0.309 0.305 0.311

heading: 168 0.417 0.395 0.310 0.384 0.409 0.378 0.418 0.384 0.414 0.381

heading: 336 0.550 0.475 0.550 0.467 0.553 0.460 0.562 469 0.558 0.463

heading: 24 0.252 0.277 0.253 0.278 0.250 0.276 0.249 0.250 0.281 0.277

heading: 48 0.303 0.314 0.303 0.310 0.295 0.304 0.300 0.302 0.314 0.306

heading: 168 0.416 0.385 0.413 0.383 0.401 0.373 0.418 0.418 0.391 0.378

heading: 336 0.559 0.473 0.554 0.472 0.545 0.458 0.559 0.570 0.476 0.460

text:
ure 9, we show the relationship between forecasting performance and our proposed criteria. The results are consistent with the conclusion drawn from classification results that forecasting performance, evaluated with both MSE and MAE, is positively related to the proposed criteria.
With the proposed criteria, we empirically show the advantage of the developed meta-network on learning optimal augmentations. The overall observations are consistent with ones from Section 4.4. Here, we provide more details. 1) “Random” is a variant that randomly selects a candidate augmentation method to generate augmented instances. “All” combines all candidate augmentations evenly to generate augmented instances. These two variances are routinely utilized in previous works (Chen et al., 2020; You et al., 2020). They work well when candidate augmentations are manually verified. Otherwise, their performances are heavily affected by the low-quality candidate augmentations, as shown in Tabel 8. On the other hand, by adaptively selecting optimal augmentations, InfoTS significantly improve the forecasting performances even with low-quality candidate augmentations. 2) By comparing InfoTS to its variant w/o Fidelity and w/o Variety, we further empirically demonstrate that both variety and fidelity make contributions to the superior performances of InfoTS.
D.4 FULL RESULTS OF TIME SERIES PREDICTION AND CLASSIFICATION
The full results of multivariate time series forecasting are shown in Tabel 4. Results of StemGNN with Ly = 720 are not available due to the out-of-memory error (Yue et al., 2021). Full results of univariate time series classification on 128 UCR datasets are shown in Tabel 9. Results of T-Loss, TS-TCC, and TNC are not reported on several datasets because they are not able to deal with missing observations in time series data. These unavailable accuracy scores are dismissed when computing average accuracy and considered as 0 when calculating the average rank. Results of multivariate
classification on 30 UEA datasets are listed in Tabel 10. Computations of average accuracy scores and ranks follow the ones in 128 UCR datasets.
InfoTSs InfoTS TS2Vec T-Loss TNC TS-TCC TST DTW
Adiac 0.795 0.788 0.775 0.675 0.726 0.767 0.550 0.604 ArrowHead 0.874 0.874 0.857 0.766 0.703 0.737 0.771 0.703 Beef 0.900 0.833 0.767 0.667 0.733 0.600 0.500 0.633 BeetleFly 0.950 0.950 0.900 0.800 0.850 0.800 1.000 0.700 BirdChicken 0.850 0.900 0.800 0.850 0.750 0.650 0.650 0.750 Car 0.900 0.883 0.883 0.833 0.683 0.583 0.550 0.733 CBF 1.000 0.999 1.000 0.983 0.983 0.998 0.898 0.997 ChlorineConcentration 0.825 0.822 0.832 0.749 0.760 0.753 0.562 0.648 CinCECGTorso 0.896 0.928 0.827 0.713 0.669 0.671 0.508 0.651 Coffee 1.000 1.000 1.000 1.000 1.000 1.000 0.821 1.000 Computers 0.720 0.748 0.660 0.664 0.684 0.704 0.696 0.700 CricketX 0.780 0.774 0.805 0.713 0.623 0.731 0.385 0.754 CricketY 0.774 0.774 0.769 0.728 0.597 0.718 0.467 0.744 CricketZ 0.792 0.787 0.792 0.708 0.682 0.713 0.403 0.754 DiatomSizeReduction 0.997 0.997 0.987 0.984 0.993 0.977 0.961 0.967 DistalPhalanxOutlineCorrect 0.808 0.801 0.775 0.775 0.754 0.754 0.728 0.717 DistalPhalanxOutlineAgeGroup 0.763 0.763 0.727 0.727 0.741 0.755 0.741 0.770 DistalPhalanxTW 0.720 0.727 0.698 0.676 0.669 0.676 0.568 0.590 Earthquakes 0.821 0.821 0.748 0.748 0.748 0.748 0.748 0.719 ECG200 0.950 0.930 0.920 0.940 0.830 0.880 0.830 0.770 ECG5000 0.945 0.945 0.935 0.933 0.937 0.941 0.928 0.924 ECGFiveDays 1.000 1.000 1.000 1.000 0.999 0.878 0.763 0.768 ElectricDevices 0.691 0.702 0.721 0.707 0.700 0.686 0.676 0.602 FaceAll 0.929 0.929 0.805 0.786 0.766 0.813 0.504 0.808 FaceFour 0.864 0.818 0.932 0.920 0.659 0.773 0.511 0.830 FacesUCR 0.917 0.913 0.930 0.884 0.789 0.863 0.543 0.905 FiftyWords 0.809 0.793 0.774 0.732 0.653 0.653 0.525 0.690 Fish 0.949 0.937 0.937 0.891 0.817 0.817 0.720 0.920 FordA 0.925 0.915 0.948 0.928 0.902 0.930 0.568 0.555 FordB 0.795 0.785 0.807 0.793 0.733 0.815 0.507 0.620 GunPoint 1.000 1.000 0.987 0.980 0.967 0.993 0.827 0.907 Ham 0.848 0.838 0.724 0.724 0.752 0.743 0.524 0.467 HandOutlines 0.946 0.946 0.930 0.922 0.930 0.724 0.735 0.881 Haptics 0.545 0.546 0.536 0.490 0.474 0.396 0.357 0.377 Herring 0.703 0.656 0.641 0.594 0.594 0.594 0.594 0.531 InlineSkate 0.420 0.424 0.415 0.371 0.378 0.347 0.287 0.384 InsectWingbeatSound 0.664 0.639 0.630 0.597 0.549 0.415 0.266 0.355 ItalyPowerDemand 0.971 0.966 0.961 0.954 0.928 0.955 0.845 0.950 LargeKitchenAppliances 0.851 0.853 0.875 0.789 0.776 0.848 0.595 0.795 Lightning2 0.934 0.934 0.869 0.869 0.869 0.836 0.705 0.869 Lightning7 0.863 0.877 0.863 0.795 0.767 0.685 0.411 0.726 Mallat 0.967 0.974 0.915 0.951 0.871 0.922 0.713 0.934 Meat 0.967 0.967 0.967 0.950 0.917 0.883 0.900 0.933 MedicalImages 0.920 0.820 0.793 0.750 0.754 0.747 0.632 0.737 MiddlePhalanxOutlineCorrect 0.859 0.859 0.838 0.825 0.818 0.818 0.753 0.698 MiddlePhalanxOutlineAgeGroup 0.662 0.662 0.636 0.656 0.643 0.630 0.617 0.500 MiddlePhalanxTW 0.636 0.617 0.591 0.591 0.571 0.610 0.506 0.506 MoteStrain 0.854 0.871 0.863 0.851 0.825 0.843 0.768 0.835 NonInvasiveFetalECGThorax1 0.941 0.941 0.930 0.878 0.898 0.898 0.471 0.790 NonInvasiveFetalECGThorax2 0.943 0.944 0.940 0.919 0.912 0.913 0.832 0.865 OliveOil 0.933 0.933 0.900 0.867 0.833 0.800 0.800 0.833 OSULeaf 0.760 0.760 0.876 0.760 0.723 0.723 0.545 0.591 PhalangesOutlinesCorrect 0.826 0.826 0.823 0.784 0.787 0.804 0.773 0.728 Phoneme 0.272 0.281 0.312 0.276 0.180 0.242 0.139 0.228 Plane 1.000 1.000 1.000 0.990 1.000 1.000 0.933 1.000 ProximalPhalanxOutlineCorrect 0.924 0.927 0.900 0.859 0.866 0.873 0.770 0.784 ProximalPhalanxOutlineAgeGroup 0.883 0.883 0.844 0.844 0.854 0.839 0.854 0.805 ProximalPhalanxTW 0.849 0.844 0.824 0.771 0.810 0.800 0.780 0.761 RefrigerationDevices 0.624 0.624 0.589 0.515 0.565 0.563 0.483 0.464 ScreenType 0.510 0.493 0.411 0.416 0.509 0.419 0.419 0.397 ShapeletSim 0.856 0.856 1.000 0.672 0.589 0.683 0.489 0.650 ShapesAll 0.855 0.852 0.905 0.848 0.788 0.773 0.733 0.768 SmallKitchenAppliances 0.773 0.773 0.733 0.677 0.725 0.691 0.592 0.643 SonyAIBORobotSurface1 0.921 0.927 0.903 0.902 0.804 0.899 0.724 0.725 SonyAIBORobotSurface2 0.953 0.953 0.890 0.889 0.834 0.907 0.745 0.831 StarLightCurves 0.973 0.973 0.971 0.964 0.968 0.967 0.949 0.907 Strawberry 0.978 0.978 0.965 0.954 0.951 0.965 0.916 0.941 SwedishLeaf 0.954 0.950 0.942 0.914 0.880 0.923 0.738 0.792 Symbols 0.979 0.979 0.976 0.963 0.885 0.916 0.786 0.950 SyntheticControl 1.000 1.000 0.997 0.987 1.000 0.990 0.490 0.993 ToeSegmentation1 0.930 0.934 0.947 0.939 0.864 0.930 0.807 0.772 ToeSegmentation2 0.923 0.915 0.915 0.900 0.831 0.877 0.615 0.838 Trace 1.000 1.000 1.000 0.990 1.000 1.000 1.000 1.000 TwoLeadECG 0.999 0.998 0.987 0.999 0.993 0.976 0.871 0.905 TwoPatterns 1.000 1.000 1.000 0.999 1.000 0.999 0.466 1.000 UWaveGestureLibraryX 0.820 0.819 0.810 0.785 0.781 0.733 0.569 0.728 UWaveGestureLibraryY 0.745 0.736 0.729 0.710 0.697 0.641 0.348 0.634 UWaveGestureLibraryZ 0.768 0.768 0.770 0.757 0.721 0.690 0.655 0.658 UWaveGestureLibraryAll 0.966 0.967 0.934 0.896 0.903 0.692 0.475 0.892