An Empirical Evaluation of Time-Series Feature Sets

by   Trent Henderson, et al.
The University of Sydney

Solving time-series problems with features has been rising in popularity due to the availability of software for feature extraction. Feature-based time-series analysis can now be performed using many different feature sets, including hctsa (7730 features: Matlab), feasts (42 features: R), tsfeatures (63 features: R), Kats (40 features: Python), tsfresh (up to 1558 features: Python), TSFEL (390 features: Python), and the C-coded catch22 (22 features: Matlab, R, Python, and Julia). There is substantial overlap in the types of methods included in these sets (e.g., properties of the autocorrelation function and Fourier power spectrum), but they are yet to be systematically compared. Here we compare these seven sets on computational speed, assess the redundancy of features contained in each, and evaluate the overlap and redundancy between them. We take an empirical approach to feature similarity based on outputs across a diverse set of real-world and simulated time series. We find that feature sets vary across three orders of magnitude in their computation time per feature on a laptop for a 1000-sample series, from the fastest sets catch22 and TSFEL ( 0.1ms per feature) to tsfeatures ( 3s per feature). Using PCA to evaluate feature redundancy within each set, we find the highest within-set redundancy for TSFEL and tsfresh. For example, in TSFEL, 90 of the variance across 390 features can be captured with just four PCs. Finally, we introduce a metric for quantifying overlap between pairs of feature sets, which indicates substantial overlap. We found that the largest feature set, hctsa, is the most comprehensive, and that tsfresh is the most distinctive, due to its incorporation of many low-level Fourier coefficients. Our results provide empirical understanding of the differences between existing feature sets, information that can be used to better tailor feature sets to their applications.



There are no comments yet.


page 1

page 5

page 6


catch22: CAnonical Time-series CHaracteristics

Capturing the dynamical properties of time series concisely as interpret...

Highly comparative feature-based time-series classification

A highly comparative, feature-based approach to time series classificati...

The Wasserstein-Fourier Distance for Stationary Time Series

We introduce a novel framework for analysing stationary time series base...

MultiRocket: Effective summary statistics for convolutional outputs in time series classification

Rocket and MiniRocket, while two of the fastest methods for time series ...

Massive feature extraction for explaining and foretelling hydroclimatic time series forecastability at the global scale

Statistical analyses and descriptive characterizations are sometimes ass...

Hydroclimatic time series features at multiple time scales

A comprehensive understanding of the behaviours of the various geophysic...

Feature Extraction for Functional Time Series: Theory and Application to NIR Spectroscopy Data

We propose a novel method to extract global and local features of functi...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

Time series, or ordered measurements taken over time, are a ubiquitous form of data in the sciences and industry. Applications of time-series analysis are correspondingly diverse, from forecasting (the spread of disease through a population or future rainfall) to classification (from inferring cord pH from fetal heart-rate dynamics [1] to detecting abnormalities in silicon semiconductors [2]). As a result, there is a wide range of quantitative tools for analyzing time series.

In this work, we focus on methods for extracting informative summary statistics, or features, from univariate time series (sampled uniformly in time). Such methods are diverse, and capture properties of the time-series distribution of values, autocorrelation structure, entropy, model-fit statistics, nonlinear time-series analysis, stationarity, and many others [4, 3]. The first comprehensive effort to organize the interdisciplinary time-series analysis literature, using an approach termed highly comparative time-series analysis, encoded thousands of diverse analysis methods as features and compared their behavior on a wide range of time-series data [3]. This culminated in a software framework that includes implementations of over 7 000 time-series features, as the Matlab package hctsa [4, 5]. While hctsa is comprehensive in its coverage of time-series analysis methods, it involves a computational cost exceeding what is required of many real-world applications, for which a smaller number of simpler features is often sufficient. Additionally, hctsa requires access to the proprietary Matlab software, which limits broader use. To address this limitation, catch22 is a reduced set of 22 canonical time-series features, that were derived from hctsa via a procedure designed to both maximize classification performance across 93 tasks and reduce redundancy among the top-performing features [6]. These 22 features capture the main conceptual areas of time-series analysis covered in hctsa, including properties of distributional shape, linear and nonlinear autocorrelation, long-range scaling, symbolic transition rates, and others. The catch22 feature set is coded in C, with wrappers for Python and Matlab, and derived native packages in R (Rcatch22 [7]) and Julia (Catch22.jl [8]).

Other independently developed packages for computing time-series features also exist in R, including feasts [9] and tsfeatures [10]. feasts is an updated and refined iteration of tsfeatures

, and both include sets of features tailored largely to time series commonly encountered in econometrics, such as seasonality and autoregressive processes. Examples of features included in these sets include those associated with generalized autoregressive conditional heteroscedasticity (GARCH) models, crossing points, seasonality, and Seasonal and Trend decomposition using Loess (STL).

In Python, three popular libraries for time-series feature extraction are tsfresh [11], TSFEL [12], and Kats [13]. tsfresh is a popular feature set that contains implementations of up to 1 558 features, including autocorrelation and Fourier decomposition methods, entropy, and distributional properties. TSFEL

includes features associated with autocorrelation properties and Fourier transforms, spectral quantities and wavelet decompositions, and distributional characteristics. The feature set included in

Kats, developed by Facebook’s Engineering team, was designed to be similar in composition to R’s existing tsfeatures package and includes features associated with crossing points, STL decomposition, sliding windows, autocorrelation properties, and Holt-Winters algorithms.

There are major barriers to leveraging these various feature-extraction packages, as they span programming languages and differ in required data formats and syntax. For example, Kats requires a DateTime

formatted time index and a corresponding numeric vector of values, whereas other packages represent uniformly sampled time-series data as an ordered vector of values. The R package,

theft: Tools for Handling Extraction of Features from Time series [15], addresses these difficulties, providing a standardized computational framework for time-series feature extraction, supporting the catch22, feasts, tsfeatures, tsfresh, TSFEL, and Kats feature sets.

While the time-series analysis community now has ready access to multiple time-series feature sets, how does one select an appropriate feature set to use for a given application? Does the behavior of features in different feature sets largely overlap, or are some distinctive? Different applications may require different trade-offs between feature comprehensiveness and computation speed—how do the feature-sets compare? While some sets contain many features, are these features largely capturing independent properties of the data, or are they highly correlated (and hence redundant), such that care will need to be taken to statistically account for the non-independence (or dimensionality reduction may be particularly useful). Addressing these important questions requires systematic comparison of the behavior of the existing tools for time-series feature extraction, which has not previously been done. In this work we aim to address this shortcoming by comparing seven time-series feature sets: hctsa (v1.06), catch22 (v0.1.12), feasts (v0.2.1), tsfeatures (v1.0.2), tsfresh (v0.18.0), TSFEL (v0.1.4), and Kats (initial release). We assess these feature sets in three different settings: (i) their computation time and its scaling with time-series length; (ii) the redundancy of features within a given feature set; and (iii) overlap and redundancy in feature behavior between feature sets. We take a novel empirical approach of assessing features using their behavior across a diverse set of real-world and model-generated time series. Note that tsfeatures includes multiple feature modules: we included all features here.

Ii Computation Time

In applications involving large datasets, where computational resources are limited, smaller feature sets that are fast to compute are advantageous. But for applications in which accuracy is prioritized, larger sets of complex time-series features may be slower to compute but can yield higher performance. For example, some feature sets, like hctsa, contain complex nonlinear time-series analysis methods, and tsfeatures includes GARCH modeling. These features are slower to compute than simple properties of the distribution (e.g., mean and variance) or linear correlation statistics (e.g., Fourier coefficients or linear autocorrelations). It is thus important to empirically benchmark the computation-time scaling of different time-series feature sets with time-series length.

Here we estimated the computation-time scaling of each feature set across the length range of 100–1000 samples (indicative of time-series lengths in many applied problems). We computed the time taken to evaluate each of the seven feature sets on a set of noise time series generated as independent samples from a Gaussian distribution, across five lengths,

, and

. The median and interquartile range of computation time was returned after repeating the process for ten independent noise time series of each length. We confirmed that the results are not specific to our choice to use white noise time series for benchmarking; we found quantitatively similar results when repeating the analysis with noisy sinusoids (generated from:

, for i.i.d. and ). Calculations were run on a 2019 MacBook Pro with an Intel Core i7 2.6 GHz 6-Core CPU, and all feature sets were timed in their native programming language (hctsa was run in parallel mode, allowing access to all six CPU cores).

Feature sets varied considerably in the number of features they contain, from catch22 (22 features), feasts (43 features), tsfeatures (62 features), Kats (40 features), TSFEL (185 features for time series of length , 260 for length , 285 for length , and 390 for lengths ; feature numbers scale approximately linearly with time-series length until ), tsfresh (1 558 features), and hctsa (7 730 features). All features computed successfully for all lengths for catch22, feasts, tsfeatures, and TSFEL. For tsfresh, 392 features (%) did not compute successfully for , but no issues were noted for any other time-series lengths. These features all corresponded to fast Fourier transform (FFT) coefficients for real, imaginary, absolute, and angle components of 51–99. For Kats, 3 features (%) did not compute successfully across all time-series lengths. These three features related to level, trend, and seasonal parameters of a Holt–Winters model fit, which understandably would not fit successfully on a white noise time series. For hctsa, 706 features (%) did not compute successfully for , 469 features (%) for , 458 features (%) for , 470 features (%) for , and 443 features (%) for . Some features (such as stationarity on a certain fixed timescale) will not compute on short time series, and some TISEAN-based features [14], requiring compilation from C and Fortran code, were excluded from this analysis.

Median and interquartile ranges of computation times for each feature set are plotted as a function of time-series length in Fig. 1A. We found that feature-set computation times varied over approximately four orders of magnitude. For example, for a 1000-sample time series, feature-extraction times varied from  ms for the 22 C-coded features in catch22, to  s for the hctsa features. The other feature sets showed intermediate computation times, from TSFEL ( s) and Kats ( s), to feasts (0.47 s), tsfresh (2.53 s), and tsfeatures (6.18 s). To more fairly compare feature sets of different sizes (from the 22 features in catch22 to the features in hctsa), we next investigated the computation time per successfully computed feature, shown in Fig. 1B. We find that computation time per feature also varied over several orders of magnitude. For a 1000-sample time series, feature computation times were fastest for TSFEL ( ms per feature) and catch22 (0.3 ms per feature); intermediate for Kats (2 ms per feature), hctsa (2 ms per feature), tsfresh (2 ms per feature), and feasts (11 ms per feature); and slowest for tsfeatures (280 ms per feature).

Fig. 1: Computation times of different feature sets vary across orders of magnitude. A Median and interquartile range of computation time across ten Gaussian noise time series is plotted as a function of time-series length on logarithmic axes for lengths , and . B Median computation time per successfully computed feature is plotted for the same data as in A. Results are shown for catch22 (dark green), feasts (orange), hctsa (purple), Kats (pink), tsfeatures (light green), TSFEL (yellow), and tsfresh (brown).

Iii Evaluation Dataset

Having evaluated the differences in computation time of different feature sets, we next aimed to compare their behavior on real data, in order to assess the similarity of the constituent features, both within and between feature sets. To do this, we use the ‘Empirical 1000’ dataset (version 10) [16], which consists of a diverse set of 1000 real-world and model-generated time series [17]. Model-generated data spans a broad range of generative processes, including chaotic and non-chaotic deterministic dynamical systems (both discrete and continuous), and a range of stochastic processes. The real-world time series span a broad range of domains, including hydrology, meteorology, finance, astrophysics, medicine, and animal-sound audio. While the Empirical 1000 time series range in length from to samples, we analyzed the first 1000 samples for time series longer than 1000 points to minimize the computational expense. To demonstrate the range of temporal patterns contained in the Empirical 1000 dataset, six example time series (after truncation) are plotted in Fig. 2. This diversity of dynamical patterns reflects the corresponding diversity in the types of real-world and generative models studied in science and industry [17].

Given the nature of some of the time series in the Empirical 1000, some individual features failed to compute successfully. To minimize issues with further analysis, we employed a filtering process to retain the largest subset of features possible for each feature set while also removing consistently problematic time series from the Empirical 1000. First, all features were first normalized independently using a -score transformation. We then filtered each set to only the -score transformed features that had % missing values. This removed no features for catch22 or TSFEL, 2 features (%) for feasts (41 remaining), 3 features (%) for Kats (37 remaining), 34 features (%) for tsfresh (745 remaining—a maximum of 779 features were computed on the Empirical 1000 compared to the 1 558 features computed for the Gaussian noise computation time analysis), 3 features (%) for tsfeatures (59 remaining), and 22 features (%) for hctsa (7 708 remaining). After filtering for problematic features, we then filtered out any time series in the Empirical 1000 with missing values across any feature in any set. This removed 105 time series, leaving 895 time series remaining. Of the 105 removed problematic time series, 99 (%) were generated due to features in TSFEL, (%) were generated due to tsfresh, 1 (%) was generated due to feasts, and 1 (%) was generated due to tsfeatures. We use these 895 time series and the retained features for each set to characterize feature behaviour throughout the remainder of this paper.

Fig. 2: We assess time-series features based on their behavior on a dataset of 1000 diverse real-world and model-generated time series, a selection of which are plotted here [16, 17]. Six example time series from the Empirical 1000 dataset [16, 17] are plotted (after truncation to the first 1000 samples, as described in text): A audio of the Black-bellied Plover (AS_s2.1_f4_b8_l8280_135971), B the Trans Polar Index normalized pressure difference between Hobart and Stanley (CM_tpi_hob), C a simulation from the dissipative forced Brusselator ODE [18] (FL_fbruss_L300_N10000_IC_0.3_2_y), D a gait time series from the Parkinson’s disease gait dataset [19, 20] (MD_gaitpdb_GaCo05_01l_SNIP_12562955), E river flow (ME_winter_rf), and F a simulation from the dissipative predator-prey map [18] (MP_predprey_L5000_IC_0.5_0.55_x). Time-series names are provided in parentheses for reference to their corresponding data entries on [17].

Iv Within Feature-Set Redundancy

In constructing a feature set, choices need to be made in selecting the range of conceptually distinct types of time-series analysis algorithms to include, their parametric variations, and how to represent the method’s outputs as one or more real-valued features. For example, after computing the linear autocorrelation function, one could simply extract the raw values of the function at a chosen set of time-lags, (e.g., ), but could also capture different properties of the shape of the function (e.g., an exponential fit), or extract a timescale estimate, such as the function’s first or zero crossing. Similar choices need to be made when extracting features from the output of many other time-series analysis methods, including a Fourier or wavelet decomposition, the Lyapunov exponent spectrum, and many others. In this process of coding the outputs of time-series analysis methods as real-valued features, some feature sets that take a comprehensive approach, in which many features are included and thus use multiple features to capture similar properties, are expected to have higher feature–feature redundancy. Other feature sets may take a more restrained approach, aiming to include a smaller set of features that are less redundant with each other. When doing statistical analysis with a feature set, it is important to understand the effective number of independent features it measures, especially when there is a high degree of collinearity between the included features.

We therefore aimed to investigate feature-set redundancy by assessing the range of distinct behaviors exhibited by features in each of the seven feature sets. Within-set redundancy was measured by the ability of low-dimensional feature components to capture variance in the full feature set when applied to the Empirical 1000 dataset [16]

. A small number of low-dimensional components will explain a large amount of variance in a feature set that contains lots of redundant features (that produce highly correlated outputs across a diverse time-series dataset). But feature sets containing linearly independent features that display distinctive behavior on time-series data will require many more components to explain a high percentage of variance across the full set of features. Here we used the linear dimensionality-reduction method, principal components analysis (PCA), on the 895 time series

filtered feature matrix obtained from performing feature extraction for a given feature set on the 895 time series in the Empirical 1000 dataset. Extracted principal components were ordered by their percentage of explained variance.

To put feature sets of different sizes on a comparable scale, we expressed the number of principal components (PCs) as a proportion of the total number of PCs computed on the dataset. This ranged from 22 for catch22, up to a maximum of 895 (the number of time series in the dataset). We then summarized redundancy in a feature set using a simple summary statistic: the proportion of PCs required to explain 90% of the variance in the feature space. Results for this with full cumulative variance curves are plotted in Fig. 3A and the overall summary statistic is plotted in Fig. 3B. The results reveal a substantial variation in feature redundancy across feature sets. We observed the highest level of within-set redundancy for TSFEL and tsfresh. For TSFEL, 90% of the variance across all 390 features can be captured with just 4 PCs (1.0%), and for tsfresh, 90% of the variance of all 779 features was captured with 55 PCs (7.4%). To investigate this result, we inspected the features that make up these feature sets. We found that a large proportion of features in both feature sets were derived from fast Fourier transform (FFT) coefficients: 256/390 (or 66%) of TSFEL features, and 400/779 (or 51%) of tsfresh features. This large proportion of features derived from the FFT, in different forms of the raw coefficients, contributes to the high redundancy of features in TSFEL and tsfresh.

The most diverse feature set for its size was catch22, for which 11 PCs (50%) were required to capture 90% of the variance across its 22 features. This is consistent with the construction of catch22 as a reduction from 1000s of initial features that included a clustering procedure to reduce redundancy [6]. An intermediate level of redundancy is seen in feasts, tsfeatures, and Kats, which have overlaps between their constituent features (the Kats feature set was based on tsfeatures, and feasts includes five of the same features as tsfeatures). Despite its goal of being as extensive as possible in surveying time-series analysis algorithms, hctsa exhibited an intermediate level of redundancy, with 24.9% of the PCs required to capture 90.0% of the variance across all features.

Fig. 3: Feature sets vary substantially in their within-set redundancy. We assess within-set redundancy in terms of the ability of low-dimensional principal components to capture large amounts of variance in the feature space comprising a diverse set of 895 real-world and model-generated time series [16]. A Cumulative variance explained is plotted as a function of the percentage of the number of principal components included for the seven labeled feature sets. The dashed horizontal line indicates the 90% cumulative variance threshold used in B. B Percentage of the number of principal components required to explain 90% of the variance in the feature space. Results are plotted for the seven feature sets, ordered from least to most redundancy.

V Between-Set Overlap and Redundancy

Feature sets vary both in their coverage of time-series analysis methods and in their style of distilling these methods in the form of features. Do some sets contain many unique features (that are not contained in other sets), or do some feature sets reproduce the behavior of features contained in other feature sets? In this section, we aimed to address these questions by investigating the similarity between features contained in different feature sets. We used the same set of well-behaved features and 895 time series that were used above (see Sec. III): catch22 (22 features), TSFEL (390), feasts (41), Kats (37), tsfresh (745), tsfeatures (59), and hctsa (7708).

We developed a method to quantify the redundancy of a feature set, given another feature set. To achieve this for a given pair of feature sets, we first labeled one as a benchmark, . We then aimed to capture how similar features in the second set, labeled test, , were to the features in the benchmark set. If all individual features in the test feature set, , behave similarly to a feature that is contained in the benchmark set, , then we conclude that is redundant given . To quantify the similarity between two features, we used the absolute value of the Spearman correlation coefficient, , between their outputs on 895 diverse time series from the Empirical 1000 dataset [16]. In this way, features with highly correlated outputs on these diverse data are judged as having similar behavior. We defined an ordered similarity metric, , to judge the redundancy of a feature set, , given a benchmark set, , based on the similarities of their constituent features as


where is the maximum Spearman rank correlation between each feature, , in , and all features, , in . The average, , is then taken over all features, , in , yielding , the mean absolute maximum correlation between features in , given those that exist in the benchmark, .

Importantly, is asymmetric, allowing us to understand the relative overlaps and redundancies between feature sets. For example, consider a comprehensive feature set, , containing many diverse and useful time-series features, and another simple set, , that only contains autocorrelation-related features. Then will be high (indicating that there is, on average, a feature in that is similar to each feature in ), but will be low (indicating that the simple features in cannot capture the wider range of behavior of the features in ). In the case that and are sets containing very similar features, then both and will be high. We computed for all pairs of feature sets, treating each feature set as both a benchmark and a test. The result was represented as a pairwise matrix, with each benchmark feature set, , shown as rows, and each test feature set, , shown as columns, as shown in Fig. 4.

The matrix reveals a range of interesting overlaps and redundancies between feature sets. Overall, there are many high values, , indicating substantial overlap between the types of features included in the feature sets. Looking across a row, feature sets with many high values, , are ‘comprehensive’ in the sense that they contain features that exhibit similar behavior as those in the other feature sets. As expected, hctsa, the package with the most numerous and comprehensive coverage of time-series analysis algorithms, was the strongest benchmark set: on average, there exists an feature in hctsa that is highly correlated to features in catch22 (), feasts (), Kats (), tsfeatures (), and TSFEL (). The exception is tsfresh, for which

, indicating that, on average, for a feature selected from

hctsa, the most correlated tsfresh feature will have . Indeed, tsfresh is characterized by its distinctiveness overall: when considered as a test set, (inspecting the tsfresh column in Fig. 4), values are even lower () for all other feature sets treated as benchmarks, . This indicates that tsfresh contains many unique time-series features that are not present in the other feature sets.

To investigate this uniqueness, we examined the individual features in tsfresh that were the least correlated to features in hctsa; in particular, those with (the maximum taken across all hctsa features). We found 128 tsfresh features with to features in hctsa. Of these, 118 (or 92%) were Fourier coefficients (real and imaginary components, and angles). While hctsa contains summary statistics derived from the Fourier power spectrum, they do not output the set of all raw coefficients from the fast Fourier transform (FFT). This appears to drive the unique behavior of tsfresh, and also contributes to its relatively high within-set redundancy, cf. Fig. 3. Although TSFEL computes the absolute value of a large number of FFT coefficients, it does not also return the real and imaginary components and angles as features, and thus has a low .

As expected, given the history of their construction, tsfeatures and feasts were strongly overlapping in their behavior (in both directions, ), with high similarity scores also with the similarly constructed Kats ().

Fig. 4: Time-series feature sets show a wide variation of overlap and redundancy. We assessed the redundancy of a given feature set by treating it as a ‘test’ feature set, , given another feature set, which is treated as a ‘benchmark’, . The resulting directed similarity metric, , Eq. (1) (a mean of the maximum absolute Spearman correlation coefficients between feature outputs across a diverse time-series dataset), gives high values of when features in reproduce the behavior of features that exist in , and low values when features in are on average uncorrelated (distinctive) from the features in . Here we plot as a heatmap, where all feature sets have been treated as both a benchmark feature set (rows) and a test feature set (columns).

Vi Discussion

The feature-based approach to time-series analysis, in which a range of algorithms underpinned by theory can be used to distill a large time-series dataset down to a set of interpretable summary statistics is a promising way to address important time-series analysis problems in science and industry. Here we characterized the behavior of, and overlaps between, the most prominent existing software libraries for time-series feature extraction: hctsa, catch22, tsfeatures, feasts, Kats, tsfresh, and TSFEL. We took an empirical approach of characterizing features in terms of their pattern of outputs across a set of over 800 diverse real-world and model-generated time series, allowing us to distinguish ‘similar’ features as those with highly correlated behavior on data. We found that TSFEL and catch22 were the fastest to evaluate for time series in the range 100–1000 samples ( s per feature). However, our low-dimensional analysis highlighted a large amount of redundancy in particularly TSFEL and tsfresh, with a large proportion of highly correlated features (including many raw Fourier coefficients) in both feature sets. By contrast, catch22, a feature set designed to reduce within-set redundancy, displayed the least within-set redundancy. We also analyzed the overlap between feature sets by introducing a directed pairwise metric of feature-set redundancy, that treats one set as a ‘benchmark’ and measures how correlated features in a second ‘test’ feature set are to those in the benchmark set. As expected, hctsa was distinguished as the most comprehensive feature set, with features in catch22, feasts, Kats, tsfeatures, and TSFEL all having a highly correlated match () to a feature contained in hctsa. The exception was tsfresh, which was marked as the most unique feature set, largely driven by a large number of features that capture the raw outputs of a fast Fourier transform as real and imaginary components, magnitudes, and phase angles (whereas many other packages extract summary statistics from the Fourier power spectrum).

Our results also highlight different choices made in feature composition and behavior across the seven time-series feature extraction libraries. It is important to understand the characteristics of different feature sets, as different applications of time-series feature extraction call for different choices, such as prioritizing speed or accuracy. Different feature sets take different approaches to feature-set construction, from the large hctsa feature set, which was developed with the aim of representing a comprehensive survey of the time-series analysis literature, to its distillation as the fast and less correlated features in catch22. We also found a large amount of redundancy in many of the feature sets, such as TSFEL and tsfresh, where the variance across a large number of included features can be well captured using a small proportion (

%) of reduced components. This redundancy was driven by large numbers of features quantifying similar types of time-series structure (such as Fourier coefficients or quantiles of the distribution of values). Feature comprehensiveness comes at the cost of redundancy, and could be beneficial for some applications, e.g., involving large datasets and for which raw Fourier coefficients (including phase information of specific components) is relevant. When applying statistical models to data, it is also important to be aware of and properly account for high collinearity of features within a feature set, which is substantial in all feature sets analyzed here.

While this work focused on the similarity of feature behavior on diverse time-series data (from the Empirical 1000 dataset [16, 17]

), it will be important for future work to properly evaluate the relative performance of the different feature sets on different types of problems involving time series. For example, the most common application of feature-based time-series analysis has been to time-series classification tasks—comparing the performance of different feature sets for different types of classification (and other) tasks will allow for clearer understanding of the relative strengths and weaknesses of these feature sets. Other tasks which may yield informative differences between feature sets include large-scale anomaly detection

[21] and feature-based forecasting [22]. Our results may also depend on our choice to quantify feature similarity using a wide variety of data from the Empirical 1000 dataset [16]. While features that are highly correlated across such diverse data can be considered to be measuring similar time-series properties in similar ways, feature–feature relationships may differ on more specific datasets, or on time series of lengths outside of the length range () considered here. We note that the R package, theft [15], provides easy access to the six Python and R feature sets analyzed here (excluding the Matlab-based hctsa), allowing researchers to leverage combinations of features from different feature sets. In the future, we hope that comparisons of time-series analysis methods across disciplines and feature sets will highlight fruitful directions for leveraging the combined strengths of diverse features, towards optimising an analysis approach for a given task.

Vii Available Code

Code to reproduce all analysis included in this paper (including downloading and processing the Empirical 1000 dataset [16]) is available at


  • [1] B. D. Fulcher, A. E. Georgieva, C. W. G. Redman, and N. S. Jones, “Highly comparative fetal heart rate analysis,” 34th Annual International Conference of the IEEE EMBS, 2012.
  • [2] B. D. Fulcher and N. S. Jones, “Highly comparative feature-Based time-series classification,” IEEE Transactions on Knowledge and Data Engineering, vol. 26, no. 12, pp. 3026–3037, 1 Dec. 2014.
  • [3]

    B. D. Fulcher, “Feature-based time-series analysis,” Feature engineering for machine learning and data analytics, pp. 87–-116, CRC Press, 2018.

  • [4] B. D. Fulcher, M. A. Little, and N. S. Jones, “Highly comparative time-series analysis: the empirical structure of time series and their methods,” J. Roy. Soc. Interface 10, 20130048, 2013.
  • [5] B. D. Fulcher and N. S. Jones. “hctsa: A computational framework for automated time-series phenotyping using massive feature extraction,” Cell Systems, 5, 527, 2017.
  • [6] C. H. Lubba, S. S. Sethi, P. Knaute, S. R. Schultz, B. D. Fulcher, and N. S. Jones, “catch22: CAnonical Time-series CHaracteristics.,” Data Min Knowl Disc, 33, pp. 1821-–1852, 2019.
  • [7] T. Henderson, “Rcatch22: Calculation of 22 CAnonical Time-Series CHaracteristics,” R package version 0.1.12, Zenodo,, 2021.
  • [8] B. Harris. “brendanjohnharris/Catch22.jl: v0.2.1,” Zenodo,, 2021.
  • [9] M. O’Hara-Wild, R. J. Hyndman, and E. Wang, “feasts: Feature Extraction and Statistics for Time Series,” R package version 0.2.1, 2021.
  • [10] R. J. Hyndman, Y. Kang, P. Montero-Manso, T. Talagala, E. Wang, Y. Yang, and M. O’Hara-Wild, “tsfeatures: Time series feature extraction,” R package version 1.0.2, 2020.
  • [11] M. Christ, N. Braun, J. Neuffer, and A. W. Kempa-Liehr, “Time Series FeatuRe Extraction on basis of Scalable Hypothesis tests (tsfresh – A Python package),” Neurocomputing, 307, pp.72-77, 2018.
  • [12] M. Barandas, D. Folgado, Fernandes, L, et al., “TSFEL: Time Series Feature Extraction Library,” SoftwareX 11, 2020.
  • [13]

    Facebook Infrastructure Data Science. “Kats,”
  • [14] R. Hegger, H. Kantz, and T. Schreiber, “Practical implementation of nonlinear time series methods: The TISEAN package,” CHAOS, 9, 413, 1999.
  • [15] T. Henderson, “theft: Tools for Handling Extraction of Features from Time series,” R package version 0.1.32, Zenodo,, 2021.
  • [16] B. D. Fulcher, “1000 Empirical Time series,” figshare, Dataset,, 2017.
  • [17] B. D. Fulcher, C. H. Lubba, S. S. Sethi, and N. S. Jones, “A self-organizing, living library of time-series data,” Scientific Data, 7, 213, 2020.
  • [18] J. C. Sprott, “Chaos and Time-Series Analysis,” Oxford University Press, Inc., USA, 2003.
  • [19] A. Goldberger, L. Amaral, L. Glass, J. Hausdorff, P. C. Ivanov, R. G. Mark, J. E. Mietus, G. B. Moody, C. K. Peng, and H. E. Stanley, “PhysioBank, PhysioToolkit, and PhysioNet: Components of a new research resource for complex physiologic signals,” Circulation [Online], 101 (23), pp. e215-–e220, 2000.
  • [20] J. Hausdorff, “Gait in Parkinson’s Disease v1.0.0,”, 2008.
  • [21] R. J. Hyndman, E. Wang, and N. Laptev, “Large-Scale Unusual Time Series Detection,” 2015 IEEE International Conference on Data Mining Workshop (ICDMW), pp. 1616-1619, 2015.
  • [22] P. Montero-Manso, G. Athanasopolous, R. J. Hyndman, and T. S. Talagala, “FFORMA: Feature-based forecast model averaging,” International Journal of Forecasting, 36, pp. 86-92, 2020.