π

π

π

π

catch22 Features

Searchβ¦

Nonlinear Autocorrelation

These features capture nonlinear autocorrelation properties of a time series

This statistic computes computes the average across the time series of the cube of successive time-series differences. It will be close to zero for time series for which the distribution of successive decreases in the time series matches the distribution of successive increases, but will be positive if increases tend to be larger in magnitude and negative if decreases tend to be larger in magnitude.

The statistic is computed as:

$\langle (x_{i+1} - x_i)^3\rangle_t\,,$

for time-series values, *x*, averaged across all possible time points*, t *(i.e., from index 1 to N-1, for a time series of length *N*).

It is based on a statistic used in nonlinear time-series analysis, cf. Surrogate time series, Schreiber and Schmitz, *Physica D*, **142: **346 (2000).

The original *hctsa* implementation is the code *z*-scored time series,

`CO_trev(x_z,1)`

, for a `x_z`

.For example, this Chua map time series has increments that are approximately symmetric (blue histogram), as are the cubic increments (orange histogram), and so this statistic has a value (the mean of the orange distribution) of approximately zero.

This simple flow has a multimodal distribution of successive increments (blue), with an asymmetry towards larger increases compared to smaller decreases. This is accentuated after taking a cube (orange), yielding a positive value of this statistic:

This Lozi map is an example of the opposite behavior: time-series decreases can be sudden and large in magnitude, whereas increases are relatively gradual. This leads to an asymmetry towards large negative values, leading to a negative value for this statistic.

(NB: renamed in current *hctsa *library, as the feature named

`CO_HistogramAMI_even_5bin_ami2`

)This feature is a nonlinear version of the autocorrelation function: using a nonlinear correlation metric (mutual information) instead of a conventional linear correlation metric, evaluated using a histogram with 5 bins and at a time delay *Ο* = 2 (from the *hctsa* code

`CO_HistogramAMI(x_z,2,'even',5)`

).This feature gives high values to time series like this Chaotic Web map:

, which has clear dependence structure of the time-series value at the current point, *, *yielding a high value for this feature of 1.25:

$x_t$

, and the value two time points ahead, $x_{t+2}$

This Tent map also has a clear (nonlinear) dependence of present time-series value with that two-steps into the future (and a moderate value of this feature, 0.473):

as seen in its embedding:

And a moderate value (0.25) for this heart rhythm:

ββ

But a low value (0.02) for this earthquake time series:

ββ

β

β

β

Last modified 6mo ago

Copy link

On this page

CO_trev_1_num

CO_HistogramAMI_even_2_5

IN_AutoMutualInfoStats_40_gaussian_fmmi