Recurrence quantification analysis

Recurrence quantification analysis (RQA) is a method of nonlinear data analysis (cf. chaos theory) for the investigation of dynamical systems. It quantifies the number and duration of recurrences of a dynamical system presented by its phase space trajectory.

Background

The recurrence quantification analysis was developed in order to quantify differently appearing recurrence plots (RPs) based on the small-scale structures therein. Recurrence plots are tools which visualise the recurrence behaviour of the phase space trajectory of dynamical systems. They mostly contain single dots and lines which are parallel to the mean diagonal ("line of identity", LOI) or which are vertical/horizontal. Lines parallel to the LOI are referred to as "diagonal lines" and the vertical structures as "vertical lines". Because an RP is usually symmetric, horizontal and vertical lines correspond to each other, and, hence, only vertical lines are considered. The lines correspond to a typical behaviour of the phase space trajectory: whereas the diagonal lines represent such segments of the phase space trajectory which run parallel for some time, the vertical lines represent segments which remain in the same phase space region for some time.

If only a time series is available, the phase space can be reconstructed by using a time delay embedding (see Takens' theorem):

:vec{x}(i) = (u(i), u(i+ au), ldots, u(i+ au(m-1)),

where u(i) is the time series, m the embedding dimension and au the time delay.

The RQA quantifies the small-scale structures of recurrence plots, which present the number and duration of the recurrences of a dynamical system. The measures introduced for the RQA were developed heuristically between 1992 and 2002 (Zbilut & Webber 1992; Webber & Zbilut 1994; Marwan et al. 2002). They are actually measures of complexity. The main advantage of the recurrence quantification analysis is that it can provide useful information even for short and non-stationary data, where other methods fail.

RQA can be applied to almost every kind of data. It is widely in physiology, but was also successfully applied on problems from engineering, chemistry, Earth sciences etc.

RQA measures

The simplest measure is the recurrence rate, which is the density of recurrence points in an RP:

:RR = frac{1}{N^2} sum_{i,j=1}^N mathbf{R}(i,j).

The recurrence rate corresponds with the probability that a specific state will recur. It is almost equal with the definition of the correlation sum, where the LOI is excluded from the computation.

The next measure is the percentage of recurrence points which diagonal lines in the RP of minimal length l_{min}:

:DET=frac{sum_{l=l_{min^N l, P(l)}{sum_{i,j=1}^N R(i,j)},

where P(l) is the frequency distribution of the lengths l of the diagonal lines. This measures is called determinism and is related with the predictability of the dynamical system, because white noise has an RP with almost only single dots and very few diagonal lines, whereas a deterministic process has an RP with very few single dots but many long diagonal lines.

The amount of recurrence points which form vertical lines can be quantified in the same way:

:LAM=frac{sum_{v=v_{min^{N}vP(v)}{sum_{v=1}^{N}vP(v)},

where P(v) is the frequency distribution of the lengths v of the vertical lines, which have at least a length of v_{min}. This measure is called laminarity and related with the amount of laminar phases in the system (intermittency).

The lengths of the diagonal and vertical lines can be measured as well. The averaged diagonal line length

:L = frac{sum_{l=l_{min^N l, P(l)}{sum_{l=l_{min^N P(l)}

is related with the "predictability time" of the dynamical systemand the trapping time, measuring the average lengthof the vertical lines,

:TT = frac{sum_{v=v_{min^{N} v P(v)} {sum_{v=v_{min^{N} P(v)}

is related with the "laminarity time" of the dynamical system, i.e. how long the system remains in a specific state.

Because the length of the diagonal lines is related on the time how long segments of the phase space trajectory run parallel, i.e. on the divergence behaviour of the trajectories, it was sometimes stated that the reciprocal of the maximal length of the diagonal lines (without LOI) would be an estimator for the positive maximal Lyapunov exponent of the dynamical system. Therefore, the maximal diagonal line length L_{max} or the divergence

:DIV = frac{1}{L_{max

are also measures of the RQA. However, the relationship between these measures with the positive maximal Lyapunov exponent is not as easy as stated, but even more complex (to calculate the Lyapunov exponent from an RP, the whole frequency distribution of the diagonal lines has to be considered). The divergence can have the trend of the positive maximal Lyapunov exponent, but not more. Moreover, also RPs of white noise processes can have a really long diagonal line, although very seldom, just by a finite probability. It is obvious that therefore the divergence cannot reflect the maximal Lyapunov exponent.

The probability p(l) that a diagonal line has exactly length l can be estimated from the frequency distribution P(l) with p(l) = frac{P(l)}{sum_{l=l_{min^N P(l)}. The Shannon entropy of this probability,

:ENTR = - sum_{l=l_{min^N p(l) ln p(l),

reflects the complexity of the deterministic structure in the system. However, this entropy depends sensitively on the bin number and, thus, may differ for different realisations of the same process, as well as for different data preparations.

The last measure of the RQA quantifies the paling of the RP. The trend is the number of the linear relationship between the density of recurrence points in a line parallel to the LOI and its distance to the LOI. More exact, let us consider the recurrence rate in a diagonal line parallel to LOI of distance k ("diagonal-wise recurrence rate"):

:RR_k = frac{1}{N-k} sum_{j-i=k}^{N-k} mathbf{R}(i,j),

then the trend is defined by

:TREND = frac{sum_{i=1}^ ilde{N} (i- ilde{N}/2)(RR_i - langle RR_i angle)}{sum_{i=1}^ ilde{N} (i- ilde{N}/2)^2},

with langle cdot angle as the average value and ilde{N} < N. This latter relation should ensure to avoid the edge effects of too low recurrence point densities in the edges of the RP. The measure "trend" provides information about the stationarity of the system.

Similar to the diagonal-wise defined recurrence rate, the other measures based on the diagonal lines (DET, L, ENTR) can be defined diagonal-wise. These definitions are useful to study interrelations or synchronisation between different systems (using recurrence plots or cross recurrence plots).

Time-dependent RQA

Instead of computing the RQA measures of the entire recurrence plot, they can be computed in small windows moving over the recurrence plot along the LOI. This provides time-dependent RQA measures which are able to detect, e.g., chaos-chaos transitions (Marwan et al. 2002). Note: the choice of the size of the window can strongly influence the measure "trend".

Example


References

*
*
*
*

See also

* Recurrence plot, a powerful visualisation tool of recurrences in dynamical (and other) systems.


Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Recurrence period density entropy — (RPDE) is a method, in the fields of dynamical systems, stochastic processes, and time series analysis, for determining the periodicity, or repetitiveness of a signal. Overview Recurrence period density entropy is useful for characterising the… …   Wikipedia

  • Recurrence plot — In descriptive statistics and chaos theory, a recurrence plot (RP) is a plot showing, for a given moment in time, the times at which a phase space trajectory visits roughly the same area in the phase space. In other words, it is a graph of… …   Wikipedia

  • RQA — Recurrence Quantification Analysis (Academic & Science » Mathematics) * Register of Qualified Aromatherapists (Medical » Physiology) * Rational QualityArchitect (Computing » Software) * Remote Query Accelerator (Computing » General) * Recursive… …   Abbreviations dictionary

  • Time series — Time series: random data plus trend, with best fit line and different smoothings In statistics, signal processing, econometrics and mathematical finance, a time series is a sequence of data points, measured typically at successive times spaced at …   Wikipedia

  • List of mathematics articles (R) — NOTOC R R. A. Fisher Lectureship Rabdology Rabin automaton Rabin signature algorithm Rabinovich Fabrikant equations Rabinowitsch trick Racah polynomials Racah W coefficient Racetrack (game) Racks and quandles Radar chart Rademacher complexity… …   Wikipedia

  • Correlation dimension — In chaos theory, the correlation dimension (denoted by ν) is a measure of the dimensionality of the space occupied by a set of random points, often referred to as a type of fractal dimension.[1][2][3] For example, if we have a set of random… …   Wikipedia

  • Lyapunov exponent — In mathematics the Lyapunov exponent or Lyapunov characteristic exponent of a dynamical system is a quantity that characterizes the rate of separation of infinitesimally close trajectories. Quantitatively, two trajectories in phase space with… …   Wikipedia

  • RQA — is a three letter acronym that may stand for:* The Raising Quality and Achievement programme, a UK scheme for coordination of standards of educational achievement and student performance in colleges * The Return of Qualified Afghans programme, an …   Wikipedia

  • Correlation sum — In chaos theory, the correlation sum is the estimator of the correlation integral, which reflects the mean probability that the states at two different times are close: where N is the number of considered states , ε is a threshold distance, a… …   Wikipedia

  • Correlation integral — In chaos theory, the correlation integral is the mean probability that the states at two different times are close: where N is the number of considered states , ε is a threshold distance, a norm (e.g. Euclidean norm) and …   Wikipedia


Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.