Next Article in Journal
PV System Failures Diagnosis Based on Multiscale Dispersion Entropy
Previous Article in Journal
Work and Thermal Fluctuations in Crystal Indentation under Deterministic and Stochastic Thermostats: The Role of System–Bath Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Geometric Attractor Structure and Recurrence Analysis in Professional Dancers

by
Michalina Błażkiewicz
Faculty of Rehabilitation, Józef Piłsudski University of Physical Education in Warsaw, 00-809 Warszawa, Poland
Entropy 2022, 24(9), 1310; https://doi.org/10.3390/e24091310
Submission received: 8 July 2022 / Revised: 8 September 2022 / Accepted: 14 September 2022 / Published: 16 September 2022
(This article belongs to the Topic Recent Trends in Nonlinear, Chaotic and Complex Systems)

Abstract

:
Background: Human motor systems contain nonlinear features. The purpose of this study was to evaluate the geometric structure of attractors and analyze recurrence in two different pirouettes (jazz and classic) performed by 15 professional dancers. Methods: The kinematics of the body’s center of mass (CoM) and knee of the supporting leg (LKNE) during the pirouette were measured using the Vicon system. A time series of selected points were resampled, normalized, and randomly reordered. Then, every second time series was flipped to be combined with other time series and make a long time series out of the repetitions of a single task. The attractors were reconstructed, and the convex hull volumes (CHV) were counted for the CoM and LKNE for each pirouette in each direction. Recurrence quantification analysis (RQA) was used to extract additional information. Results: The CHVs calculated for the LKNE were significantly lower for the jazz pirouette. All RQA measures had the highest values for LKNE along the mediolateral axis for the jazz pirouette. This result underscores the high determinism, high motion recurrence, and complexity of this maneuver. Conclusions: The findings offer new insight into the evaluation of the approximation of homogeneity in motion control. A high determinism indicates a highly stable and predictive motion trajectory.

1. Introduction

According to Stergiou and Decker [1], the deviations observed during multiple repetitions of a given task are defined as movement variability. This variability is an integral feature of living organisms because the repetition of activities is associated with unique neuromotor patterns [2,3]. It implies that even an elite performer cannot perform the same task identically twice. According to Riley and Turvey [4] and van Emmerik et al. [5], human movement systems present stochastic behaviors potentially affected by deterministic and random processes. However, as van Mourik et al. [6] noted, most studies focused only on the repetitive, deterministic features of human motion. Nevertheless, research emphasizing the stochastic aspects of human movement has become more common in recent years. Despite this development, the deterministic and stochastic features of human movement are rarely evaluated together. Accordingly, human motion should be treated as a dynamic system in which a subject moves between states over time.
The dynamical systems approach focuses on how the system (1) maintains its current state (the stability problem), (2) changes or transitions between states (the task of variability and adaptability), and (3) regulates the complexity resulting from interactions between the central nervous system, body, and the surrounding environment. When dealing with dynamic systems, it is worth checking whether they exhibit nonlinear characteristics before applying nonlinear tools [2,7]. Nonlinear tools for evaluating the features of the system mentioned above include the Lyapunov exponent, Hurst exponent, fractal dimension, entropy families [8,9,10], the possibility of attractor reconstruction (phase space reconstruction) [11], and recurrence quantification analysis (RQA) [12]. Kędziorek and Błażkiewicz [13] have shown that sample entropy, fractal dimension, and the Lyapunov exponent are highly applicable to assessing postural control. On the other hand, RQA and the phase space concept in postural control assessment are less common in the literature.
One of the most powerful tools for identifying different behaviors of a dynamical system is the Lyapunov exponent (LyE), more precisely, its value (sign) [14]. It measures the exponential rate at which nearby trajectories converge or diverge. The Lyapunov exponent shows whether the observed dynamic system is stable, periodic, quasi-periodic, or chaotic [15]. Over the years, many algorithms have been developed for calculating Lyapunov exponents from time series [15,16,17,18]. The two most popular are the Wolf algorithm [15] and Rosenstein’s method [16]. According to Rispens et al. [19], the accuracy of Rosenstein’s algorithm is equal to or greater than Wolf’s algorithm. The latter appeared more sensitive to the number of observations and the amount of noise in the data. Thus, Wolf’s algorithm is more accurate for long-term data with low noise levels [16]. In chaos research, another parameter, known as the Hurst exponent, can be defined. The Hurst exponent measures the predictability of a time series and describes its level of statistical self-similarity and long-term memory [20]. The larger the value of H, the smoother the time series. H reaches values between 0 and 1. The properties of H can be summarized as follows: (1) H = 0.5 indicates random series (time series behaves like the Brownian process); (2) 0.5 < H ≤ 1 indicates persistent long-range power-law correlations; (3) 0 ≤ H < 0.5 indicates anti-persistent (short-term memory, anti-correlated) process [21].
The behavior of a dynamic system is typically represented in phase space. A phase space is a space in which all possible states of a system are represented, with each possible state corresponding to one unique point in the phase space. It seems appropriate at this point to mention the definition of an attractor. According to Devaney [22], an attractor is a set toward which most other points from the state space tend toward under iteration. In dimension two, such a set would consist of a collection of limit cycles, equilibrium points, and solutions connecting them. In higher dimensions, these attractors may be much stranger. The dynamics of the attractor itself may be chaotic. There are many possible definitions of chaos. One that is recognized comes from Devaney [22]: Let X be a metric space. A continuous map f:XX is said to be chaotic on X if (1) f is transitive; (2) the periodic points of f are dense in X; (3) f has sensitive dependence on the initial conditions. An attractor is strange if it displays sensitive dependence on the initial conditions (points initially close to each other on the attractor become exponentially separated with time). According to Takens [23], to reconstruct the state space, two variables are needed: (1) minimal embedding dimension (D) and (2) optimal time delay (lag) − τ (tau) (Table 1).
The embedding dimension is the minimum number of variables required to create a valid state space from a given time series. Vectors in the embedding space (a new space) are constructed from the time-delayed values of the scalar measurement [2]. Many methods (Table 1) have been developed for estimating embedding dimensions [25,26,27,28]. There is no rule for determining the minimum embedding dimension, D, and none of the published proposals are widely accepted. Of all the possible techniques, the correlation dimension and false nearest neighbors stand out the most [29]. The autocorrelation of a time series or mutual information is a common approach for time-lag selection [24]. It is worth noting that both the embedding dimension and the time delay are also necessary when using RQA. RQA quantifies the number and duration of recurrences of a dynamical system represented by its phase space trajectory [12]. RQA allows for a model of features, such as the determinism of a system, the dimensionality or complexity of its dynamics, and the amount of patterning it holds [30]. In this section, it is worth writing that attractor reconstruction and RQA have been successfully used to characterize human behaviors like heart rate variability, postural fluctuations, and the variability and stability of dynamic biological systems [31,32].
Som, et al. [33] proposed high-dimensional shape descriptors for the reconstructed attractors of the center-of-pressure (CoP) paths collected from subjects with Parkinson’s disease to assess their balance impairment. Gates and Dingwell [34] reconstructed a state space for the assessment of shoulder movement. Bradley and Stuart [35] used attractor geometry to describe the predefined motion sequences in dance. Peppoloni et al. [36] characterized the interference of neural control strategies for dynamic fingertip forces from attractor reconstruction. In addition, attractor reconstruction can determine changes in the shape and variability of periodic signals, such as arterial blood pressure signals [37] and electroencephalogram signals [38], providing a two-dimensional attractor with features like density and symmetry [39]. Recurrence quantification analysis was used to determine the structure of CoP data during a standing state [40]. Labini et al. [41] assessed the walking balance complexity of the head, trunk, and pelvis for the gait of normal and hypo vestibular subjects. Riley and Clark [42] examined how the availability of, and alterations in, sensory information influenced the amount, variability, and temporal structure of spontaneous postural sway in young, healthy adults. Moreover, RQA is often used to evaluate heart rate variability [43,44,45]. Until now, no paper has used attractor reconstruction for pirouettes and used RQA to evaluate their repeatability. The purpose of the study was to evaluate the geometric structure of the attractors and analyze the recurrence of two different pirouettes performed by professional dancers.

2. Materials and Methods

2.1. Participants and Data Collection

A total of 15 modern dance dancers, whose primary preparation was sports acrobatics, participated in this study (Table 2). Participants had no lower limb injuries or balance disorders that would affect the quality of their task performance. It is crucial to mention that the examined dancers were from one dance group and had been dancing with each other for more than four years. The study protocol was approved by the University Research Ethics Committee (SEK 01-09/2020).
The kinematic parameters of jazz and classic pirouettes were collected using a three-dimensional motion capture system (Vicon Motion Systems, Oxford, UK). The system consisted of 9 cameras operating at 100 Hz. The two force platforms (Kistler Holding AG, Winterthur, Switzerland) were synchronized with the Vicon system. The platforms operated at 1000 Hz. A total of 15 markers were placed on the subjects, according to the lower body Plug-In-Gait scheme (Figure 1). Each dancer, after warm-up, performed two types of single-turn pirouette en dehors (classic and jazz on demi-plié) (Figure 1). All participants were barefoot and performed one turn on the left foot. Dancers rotated all the time in the clockwise direction [46].

2.2. Phases of Pirouettes and Parameters

Each pirouette (classic and jazz) begins in the classical fourth position. From this position, the dancer performs plié to initiate a clockwise rotation on the support leg (the forward leg). The turn starts with pushing off the ground with the non-supporting leg and moving the foot up to the knee of the supporting leg. This position remains throughout the rotation. The dancer then pushes up onto the pointe or demi-pointe (the ball of the foot) with the supporting leg and closes both arms in front [46,47,48]. While in the rotation, the dancer maintains their body in the vertical rotational axis. The pirouette ends with a specific body pose, adopted by placing the non-supporting leg on the ground (Figure 1A).
The main difference between the two pirouettes is the range of knee flexion and plantar flexion of the ankle joint of the supporting leg and, thus, the CoM position of the body (Figure 1B). In the classic pirouette, the knee joint should have full extension. The ankle joint should have maximum plantar flexion. In contrast, in the jazz pirouette, the knee joint should remain flexed, and the plantar flexion should be less. Although entire pirouettes were recorded, the phases of turning were taken for analysis. The analyzed movements began with knee flexion, as the first movement of the rotation began and ended when the foot of the non-supporting leg first touched the ground. The three-dimensional CoM and LKNE time series were taken for further analysis (LKNE was the marker on the lateral epicondyle of the left knee).

2.3. Phase Space Reconstruction

In this paper, it was decided to assess the structure of the considered time series by checking the nonlinearity, non-stationarity, and signs of chaos, before reconstructing the phase space. All calculations described below were completed using MatLab R2021a software (MathWorks, Natick, MA, USA).

2.3.1. New Time Series Reconstruction

The collected data for both pirouettes, for a single individual, were not long enough to meet the basic requirements for the ability to reconstruct the phase space. A new time series has been reconstructed separately for LKNE and CoM for jazz and classic pirouettes to overcome this issue. The reconstruction steps are below.
(1)
The P i = x i 1 y i 1 z i 1 x i n y i n z i n , where P i —is the matrix of LKNE or CoM coordinates for the analyzed pirouettes of the i-th person (i = 1, …, 15); n—is the length of the recorded time series. The n values were in the range of 101 to 250 frames. The x-coordinate describes motion in the anterior-posterior direction, y—along a vertical axis (inferior-superior), and z—in the mediolateral direction.
(2)
The signals from each matrix Pi were resampled to obtain 300 samples. Its new length was longer than the maximum length of the recorded time series (250). The new signals were normalized by their maximal value.
P ˜ i = x ˜ i 1 y ˜ i 1 z ˜ i 1 x ˜ i 300 y ˜ i 300 z ˜ i 300 , where P ˜ i —matrix of LKNE or CoM after transformation, for the analyzed pirouettes of the i-th person (i =1, …, 15).
(3)
Next, the new time series for LKNE and CoM for the x, y, and z coordinates were created. Individuals in the new time series have been shuffled to avoid bias based on their order.
/ s h u f f l e   i n   t h e   m a t r i x / P ˜ 1 P ˜ 15 = / i n v e r s i o n   o f   i n d e x e s / x ˜ 1 1 y ˜ 1 1 z ˜ 1 1 x ˜ 1 300 y ˜ 1 300 z ˜ 1 300 x ˜ 2 300 y ˜ 2 300 z ˜ 2 300 x ˜ 2 1 y ˜ 2 1 z ˜ 2 1 x ˜ 15 300 y ˜ 15 300 z ˜ 15 300 x ˜ 15 1 y ˜ 15 1 z ˜ 15 1
Connecting the time series constituted adding additional time series from the individuals. The even indexes were flipped to create the loop for the signal. After applying the above steps, the new coordinates (x, y, z) for LKNE and CoM were 4500 points long. The quantiles’ function was then applied. This function returned the quantiles of the input data elements for probabilities in the range [0.1, 1]. This procedure reduced the impulsive peaks in the output trajectories, resulting in shorter time series (4048 points).

2.3.2. Hurst Exponent

One of the most common tools to measure autocorrelation (persistence and long memory) and the level of noise is the Hurst exponent. The Hurst exponent (H) is defined as follows:
R / S n = C · n H
where R / S —rescaled range; R—range of n cumulative deviations from the average; S—the standard deviation of n observations; n—the number of elements of the time series; C —positive constant, and H —Hurst exponent. To calculate the Hurst exponent, it is necessary to compute the mean value of R / S n for different n, and then, using ordinary linear regression, determine H from the equation:
l o g E R / S n = H l o g n + l o g c
where E R / S n —expected value of the rescaled range [20].
Due to the specific nature of the behavior that occurs at certain values of the exponent, it was possible to determine three ranges of its occurrence: (1) 0 ≤ H < 0.5—the time series is anti-persistent, characterized by high variability; (2) H = 0.5—represents a Brownian Motion. The examined time series does not have a dominant trend by which subsequent changes take on a random character; (3) 0.5 < H ≤ 1—the time series is persistent. This time series follows an orderly course, which has the effect of maintaining the current trend. The orderliness of this interval is greater when the value of the exponent is higher [49]. In the present paper, the Kugiumtzis and Tsimpiris [50] code published on https://www.mathworks.com/matlabcentral/fileexchange/27561-measures-of-analysis-of-time-series-toolkit-mats (accessed on 20 June 2022) was applied for LKNEx, LKNEy, LKNEz, and CoMx, CoMy, CoMz time series for jazz and classic pirouettes constructed as described in the previous section.

2.3.3. Test of Non-Stationarity

A time series is non-stationary if its statistical structure changes over time. The code of Zhivomirov and Nedelchev [51] was applied to check the signal’s non-stationarity (https://www.mathworks.com/matlabcentral/fileexchange/75118-signal-stationarity-estimation-with-matlab) (accessed on 2 June 2022). The new CoM and LKNE time series for all x, y, and z coordinates, for both pirouettes, were non-stationary.

2.3.4. Nonlinearity of Time Series

The surrogate data test follows a null hypothesis that the signal is a realization of a linear Gaussian stochastic process (the data are a random series). Here, a surrogate data set is generated by transforming the analyzed data in a way so as to keep any characteristics recorded in the hypothesis but destroy any possible nonlinear features in the data. The iterative amplitude-adjusted Fourier transform (IAAFT) method was applied to generate the surrogate data [52]. New data matched the amplitude spectrum and signal distribution, according to Schreiber and Schmitz [53]. Then, one or more test statistics from the distribution of the surrogates and the original series were calculated. If the values of the test statistics differ significantly, the null hypothesis is rejected, suggesting that the series under consideration has a deterministic character. Here, the time series were nonlinear.

2.3.5. Detection of Chaos Based on Largest Lyapunov Exponent

Using the Largest Lyapunov exponent (LyE) to identify chaos in a system assumes that, if the average distance between two points increases exponentially, then the system is sensitive to a change in the initial conditions (the value of LyE is higher than zero). Thus, LyE can be described by the equation: d t = C e L y E t , where d(t) is the average divergence at time t, and C is a constant that normalizes the initial separation. Therefore, a positive LyE is considered a necessary and sufficient condition for the presence of chaos in the system. For stable boundary cycles, the Lyapunov exponent is zero, and for stable equilibrium points, it is negative [54].
In this paper, the largest Lyapunov exponent was computed for new time series using the lyapunovExponent function, which is available in MatLab. This algorithm is based on the paper by Rosenstein, Collins, and De Luca [16]. According to the MatLab website, the embedding dimension (D) and time delay (tau) values were input into the lyapunovExponent function, the calculation of which is described in the next section.

2.3.6. Embedding Dimension, Time Delay, and Phase Space Reconstruction

The embedding dimension (D) and time delay (tau) are key factors in phase space reconstruction. First, the time delay was estimated using the first minimum of the average mutual information function (AMIF), as discussed in [55]. For this purpose, the mdDelay function was used as follows:
tau = mdDelay(LKNE or CoM data, ‘maxLag’, 100, ‘plottype’, ‘all’).
In this code, the maximum time delay was set to 100 using the ‘maxLag’ parameter. The default value of 10 was not large enough, as shown by the output data in Figure 4 (in the Results section, the mutual information did not reach the minimum for those delays less than or equal to 10). The ‘plottype’ parameter was set to “all”, which means that the AMIF for each data dimension is presented in the plot. The mdDelay function used the threshold (1/e) criterion per default [55].
In the next step, the embedding dimension (D) was determined using the false nearest neighbors (FNN) function, which computes the percentage of false nearest neighbors for multidimensional input time series as a function of the embedding dimension [28]. To estimate the embedding dimension, the mdFnn function was used:
[fnnPerc, embTimes] = mdFnn(LKNE or CoM data, tau); where tau is the time lag calculated in the previous step for a given time series.
The functions presented above are well described in the paper by Wallot and Mønster [55] and are available on github.com/danm0nster/mdembedding (accessed on 20 June 2022). It is worth noting that the LKNE and CoM points are matrices constructed from x, y, and z coordinates describing the position of the knee joint and the center of mass for jazz and classic pirouettes. Therefore, from this point of view, they are dependent variables. The embedding dimension and time lag were calculated for LKNE and CoM for 3D space (named global), as for the Lorenz system [55,56], and also for each coordinate separately (named directional).
A total of 12 attractors were reconstructed after applying global and directional tau and D for the LKNE and CoM coordinates. Then, convex hulls were determined for each time series [57]. The convex hull of a sample of points is the minimum convex set enclosing them all [58]. Calculating the volume of the convex hull made it possible to compare the phase spaces [59]. The phase spaces were compared within the coordinates (x, y, z) of the LKNE and CoM points and between the pirouettes.

2.4. Recurrence Quantification Analysis (RQA)

Recurrence quantification analysis (RQA) was applied to quantify the duration and number of recurrences [5,60] presented within the generated state space for CoM and LKNE, separately for x, y, and z directions, for each pirouette. As was carried out previously, the analysis was performed for global and directional time lags and embedding dimensions. The recurrence plot is symmetrical with respect to the main diagonal. Thus, all quantitative features take place within the upper triangle. The main diagonal and the lower triangle (which provides only redundant information) are excluded from the analysis [61].
This paper uses the Toolbox of recurrence plot and recurrence quantification analysis (https://www.mathworks.com/matlabcentral/fileexchange/58246-tool-box-of-recurrence-plot-and-recurrence-quantification-analysis) (accessed on 2 June 2022) [62]. A good description of the introduced tool can be found in [63,64].
In this paper, six measures have been calculated as follows: % of recurrence (REC), % of determinism (DET), the length of the longest diagonal line segment in the plot (LMAX), Shannon information entropy (ENT), % of laminarity (LAM), trapping time (TT) [12]. The listed measures are the result of the density of the recurrence points and the structure of the diagonal and vertical lines of the recurrence plot.
The REC describes the percentage of recurrent points lying within a specified radius r(i). This variable can range from 0% (no recurrent points) to 100% (all points recurrent) [12,65].
R E C = 1 N 2 i , j = 1 N R i , j ;
where R i , j = Θ r i x i x j ; i, j = 1, 2, …, N; N—number of points on the phase space trajectory; r i —specified threshold, —the norm or metric (In this study, the Euclidean norm was calculated), Θ —the Heaviside-function, defined as:
Θ x = 1 ,   x 0 0 ,   x < 0
The 0 and 1 values of the Heaviside function are represented in white and black, respectively. According to Zbilut, et al. [66], the r(i) value was set at 1% of the maximum phase space diameter.
DET determines the proportion of recurrent points lying along diagonal line structures, except those within the main diagonal [61]. The diagonal line segments must have a certain minimum length to be not excluded. In this paper, according to Riley, et al. [67], the choice of 2 points was considered as the number of successive points defining a line segment.
D E T = l = l m i n N l P l i = 1 N l P l ;
where P l —histogram of the lengths l of the diagonal lines; N —number of points on the phase space trajectory, l m i n = 2 .
LMAX is the length of the longest diagonal line in the plot, excluding the main diagonal. The lower the LMAX value, the more chaotic the analyzed signal is [12,65].
L M A X = m a x l i ; i = 1 , N l ;
where l i —diagonal structures; N l —number of diagonal lines in the recurrence plot.
ENT is the Shannon entropy of the probability p l = P l N l to find a diagonal line of length l in the recurrence plot among the total number of diagonal lines N l [12,65]. ENT demonstrates the complexity of the recurrent trajectories versus the diagonal lines.
E N T = l = l m i n N p l l n p l .
LAM is similar to DET. LAM measures the percentage of repeated points containing vertical line structures, not diagonal lines like in the case of DET. TT is simply the average length of vertical line structures [12,60,65]. For LAM and TT, the value of the minimum line was 2 points, as for DET.
L A M = v = v m i n N v P v v = 1 N v P v ;   T T = v = v m i n N v P v v = v m i n N P v ;
where P v —is the histogram of the lengths of the vertical lines, N —number of points on the phase space trajectory.

3. Results

3.1. New Time Series Reconstruction

The reconstruction of the new time series proceeded for each of the x, y, and z coordinate of the LKNE marker and CoM for the jazz and classic pirouette, respectively. The recorded time series were short (101 to 250 points). Therefore, for each person, the series were resampled to 300 points (Figure 2A,C). Figure 2A,C show the LKNE and CoM curves for a single person for two pirouettes to help understand the creation of a new series and demonstrate their structure (Figure 2B,D). Importantly, all further calculations involved new time series.

3.2. Hurst Exponent Analysis and the Largest Lyapunov Exponent

The Hurst exponent was calculated to check if the new time series showed persistence and robustness and to find if there was any significant difference between the two pirouettes. The time series for the coordinates (x, y, z) of LKNE and CoM had H values of between 0.5 and 1, implying that the time series were persistent (Table 3). On average, for both pirouettes, the H values were 0.72.
The values of the largest Lyapunov exponents, which were positive and indicate the presence of chaos for the analyzed time series, are shown in Table 4. The average LyE is 1.64 and 1.78 for the classic and jazz pirouettes, respectively.

3.3. Determination of the Embedding Dimension and Time Delay

The global embedding dimension, calculated in 3D space (all D), was 2 for each time series (CoM and LKNE) except for the LKNE from the jazz pirouette, for which this value was equal to 3.
The directional embedding dimension calculated for the knee of the supported limb was 4 for the jazz and classic pirouettes in the x and y directions (Figure 3A,C). For the z-direction, the dimension was 1 greater for the jazz pirouette (Figure 3A,C). For CoM, the identical embedding dimension was found only in the y-direction (D = 3) (Figure 3B,D).
The global mutual information (all tau) was the average of the minimum values calculated separately for the LKNE and CoM (x, y, z) time series (Figure 4).
For the LKNE time series of the classic pirouette, the all tau value increased by 11% compared to the value calculated for the jazz pirouette (Figure 4A,C). For CoM, the situation was the opposite. The all tau value for the classic pirouette decreased by 16% compared to the value noted for jazz (Figure 4B,D).

3.4. Phases Space Reconstruction and Convex Hull Calculation

The reconstructed 3D phase space for LKNE and CoM in a group of 15 dancers was treated as a long time series (Figure 5). The reconstruction was carried out using the global embedding dimension (all D) and global time delay (all tau) (Figure 3 and Figure 4).
In addition, phase space reconstruction was also performed using embedding dimensions and time lags specific to the x, y, and z directions of the analyzed time series. Since the visualization of the phase space reconstruction had not changed, only the values of the convex hull volumes (calculated for the global tau and D and those directional ones) are included in Figure 5. It is worth highlighting that the two reconstruction methods yielded similar results.
The convex hull values calculated in both ways for the LKNE were always significantly higher for the classic pirouette than for the jazz pirouette (Figure 5A). The highest percentage differences were sequentially for LKNEz (global tau and D—854.55%, directional tau and D–864.51%), which are the mediolateral motion, LKNEx (anterior-posterior), with 132.87% (134.38%), and LKNEy (inferior-superior), with 7.69% (7.84%).
For the CoM, the situation was different (Figure 5B). The trend of the changes continued only for the z-direction (mediolateral), where the convex hull volume value was 225.3% (217.07%), which was higher for the classic pirouette. For the y-direction (inferior-superior), the volume values were equal. However, for the x-direction (anterior-posterior), the convex hull volume value was 4.41% (5.46%), which was higher for the jazz pirouette. Therefore, the classic pirouette appears to be the most compact for CoM x and y, as the convex hull volumes are the smallest in these directions. The jazz pirouette is more coherent for the knee joint. The knee joint for this pirouette has a slight change in the range of motion during rotation compared to the initial and final phases.

3.5. Recurrence Quantification Analysis

Recurrence is a measure of the reproducibility of a limit cycle. The recurrence values were calculated for global tau and D and for directional ones. The results are shown in Table 5.
It is worth noting that only the REC values calculated for the LKNE marker for the jazz pirouette in each direction were higher (by an average of 7.13%) regarding the global setting (Table 5). Overall, the mean percentage difference between recurrence measures calculated for the global, and directional settings were higher for the classic pirouette. DET and LAM are the parameters with the lowest changes, while the most significant differences were recorded for TT and REC (Table 6).
The shorter lengths of the diagonal lines within the recurrence plots indicate lower recurrence rates in the phase space. This is mainly evident for LKNEx, when compared to LKNEz and LKNEy in the classic pirouette (Figure 6). The highest recurrence rate (allREC = 90.41%; dREC = 82.46%) was for LKNE_z across the mediolateral axis in the jazz pirouette (Table 5). High values (allREC = 86.35%; dREC = 86.52%) were also recorded for CoM_z for the same pirouette.
The recurring points form easily visible diagonals that are parallel to the main diagonal (origin line) yet offset from it. This is typical of a periodic deterministic structure (Figure 6). The percentage determinism (DET), calculated for global and directional settings, was extremely high throughout the LKNE and CoM for all coordinates in both pirouettes (Table 5). The highest DET value was for the CoM and LKNE during the Jazz pirouette in the z-direction, similar to REC. The %LAM was similar to %DET. This measures the percentage of the repeated points containing vertical line structures rather than diagonal lines. Here, the LAM values (allLAM and dLAM) are high, showing that the vertical line structures are widespread in the analyzed pirouettes. The highest value was, as before, for CoMz and LKNEz for the jazz pirouette. TT is the average length of vertical line structures. The highest TT value was found for CoM and LKNE in the z-direction for the jazz pirouette, and the lowest was found for LKNE for the classic pirouette in the x-direction. As for entropy, its lowest values were for LKNEx, y, and z for the classic pirouette. The highest entropy value (allENT = 8.23 and dENT = 8.26) and, thus, the most complex signal was that of CoMz for the jazz pirouette. The entropy values calculated for LKNEz (allENT = 7.54 and dENT = 8.17) were also not much lower for the jazz pirouette.

4. Discussion

The pirouette is an essential skill for all gymnasts and dancers. It is one of the most popular movements used in training [68]. The purpose of this study was to evaluate the geometric structure of the related attractors and analyze “recurrence” in two different pirouettes (jazz and classic), as performed by professional dancers. According to Kantz and Schreiber [69], the phase space representation approach, rather than the time or frequency domain analysis, is a feature of nonlinear dynamical time series analysis. Moreover, recurrence quantification analysis (RQA) is an alternative nonlinear method for describing signal dynamics. It is remarkably robust for short, noisy, and non-stationary signals [70,71,72]. This approach was used to observe if the evaluated pirouettes in the study group were recurrent. In addition, it was possible to see which pirouette is more stable or complex.
Until now, no one has focused on attractor reconstruction for a group of dancers performing pirouettes. This study considered the most characteristic time series that differentiated the jazz and classic pirouettes. These included the center of body mass (CoM) and the position of the knee joint of the supporting limb (LKNE). The mentioned parameters are characteristic due to the manner in which the pirouettes are performed [46]. In the middle of the turn phase, the average minimum knee flexion value for the classic pirouette was 2.4°, while for the jazz pirouette, it was 22.8° [46]. In addition, Błażkiewicz [46] showed that pelvic alignment changes when performing these two pirouettes, so the positioning of the center of mass is different.
It is worth noting that typical pirouettes produce short time series, which contain about one to two revolutions. These time series are not long enough to use nonlinear measures and reconstruct the attractor, as most conditions will not be satisfied. The present study addressed this challenge by constructing a new CoM and LKNE time series, containing data from all the studied dancers. The dancers performed only single turns. It is crucial to mention that the dancers were from one dance group and had been dancing with each other for more than four years. Therefore, one may assume that the individual dancers were at the same level of motor development as far as the performance of the analyzed movements is concerned. The time series of the selected points were resampled, normalized, and randomly reordered. Then every second time series was flipped to be combined with other time series and make a long time series out of the repetitions of a single task. This approach made it possible to satisfy the conditions for the presence of chaos.
After lengthening the time series, it turned out that there are positive Lyapunov exponents in three directions for the CoM and LKNE for jazz and classic pirouettes, indicating the presence of chaos. In this study, the Rosenstein algorithm [16] was used to calculate the Lyapunov exponent. The values for LyE for CoM were within the ranges reported by Mehdizadeh [73] for the upper trunk during gait with different velocities. The ranges of LyE values reported by Mehdizadeh [73] were 0.15–3.30, 0.10–2.85, and 0.27–2.40 for anterior-posterior, mediolateral, and vertical directions, respectively. In this study, the average LyE calculated for CoM (x, y, z) was 1.62 for classic pirouettes and 1.84 for jazz pirouettes. In contrast, the average LyE calculated for LKNE (x, y, z) and CoM (x, y, z) was 1.64 and 1.78 for the classic and jazz pirouettes, respectively. Moreover, Hurst exponent analysis was applied to check the persistence of these new time series, which showed all out-coming time series to have a persistent behavior. On average, for both pirouettes, the H values were 0.72. The similar values for the Hurst exponent indicate that both pirouettes have a stable movement pattern for both the knee joint and the center of mass. It is worth highlighting that this is the first approach in the literature that presents such preparation of the time series. This preparation allows dancers to be treated as samples for the same dance moves, and thus, this time series can evaluate the global system’s dynamic behavior. Unfortunately, this approach does not allow for the individual assessment of the dancers. The global dynamic evaluation of classic and jazz pirouettes for the new time series was realized using RQA and phase space reconstruction. It is worth mentioning that the phase space reconstruction and RQA analysis were specific to embedding dimensions and time delay. Here, two types of these parameters were detected: global and directional. The directional parameters were calculated separately for the x, y, and z directions of the displacements of the analyzed points (LKNE and CoM) on the dancer’s body. However, assuming that the x, y, and z variables are dependent because they describe the movement of individual body parts, the so-called global parameters were then calculated. It is worth highlighting that both approaches yield similar results for RQA and phase space reconstruction, although the results for directional parameters were slightly higher.
The results of the RQA analysis showed that the percentage determinism (DET) calculated for both ways was exceptionally high (99.59–99.98%) for LKNE and CoM (for all coordinates) in both pirouettes. It indicates that the signals have deterministic features due to repeated (recurrent) time series at various frequencies (because they came from different individuals). A determinism of 100, meaning perfect cyclic motion, is near impossible and probably undesirable, even in professional athletes. This could be related to skill optimization in the learning process [74] and the functional variability required for flexible movement behavior [75]. The signals for the jazz pirouette had a very high recurrence. The recurrence was the highest for this pirouette along the mediolateral axis at the knee joint (REC = 90.41% (global), REC = 82.46% (directional)), and CoM (REC = 86.35% (global), REC = 86.52% (directional)). The lowest recurrence was for the x-direction of the knee in the classic pirouette (REC = 1.6% (global), REC = 1.62% (directional)). On the other hand, the entropy value for this movement was the lowest, which suggests little complexity in the distribution of periodic components [76]. The high values for the RQA measures, like trapping time (TT), laminarity (%LAM), and a maximum length of a diagonal line (LMAX), imply low complexity in the system’s dynamics. In the case of this paper, it is in the mediolateral direction for the knee joint for the jazz pirouette.
The attractor reconstruction process confirmed the results discussed for the recurrence quantification analysis. The convex hulls were calculated to compare the reconstructed attractors. According to Zhao [59], convex hulls are a common way to generalize the shape of an attractor and extract useful information, which can be used for comparison. The values for the volumes of the convex hulls were, on average, 51.41% (51.24%-calculated for directional tau and D), with this being smaller for the knee joint in the jazz pirouette compared to those calculated for the classic pirouette. This means that the jazz pirouette is more coherent for the knee joint of the supporting limb. On the other hand, for CoM, the average value of convex hulls volume was 4.06% (3.13%-calculated for directional tau and D), slightly smaller for the classic pirouette compared to Jazz. Therefore, the classic pirouette seems to be a little more compact for the CoM. The main direction causing these changes was the motion along the x-axis for both LKNE and CoM. The convex hull volume values were the highest for the x-direction (anterior-posterior) and lowest for the z-direction (mediolateral). Once again, referring to the paper of Błażkiewicz [46], it is worth emphasizing that the range of motion in all joints of the supporting leg in the sagittal plane, and the hip joint in the frontal plane for the non-supporting limb, was significantly higher for the classic pirouette compared to the jazz one. This makes the classic pirouette more challenging because it requires higher precision of movement, better neuromuscular control, and thus stability.
It is worth emphasizing that the discussed methodology is not popular. The author found only one paper [77] that dealt with similar issues. Burton et al. [77] evaluated long swings from the perspective of attractor dynamics. They proved that the senior gymnasts displayed increased recurrence characteristics in addition to a long, long swing duration and a lower angular velocity of their center of mass. Moreover, all groups of gymnasts had highly recurrent and predictable limit cycle characteristics. Burton et al. [77] showed that attractor reconstruction helps to understand the development of an athlete’s skills and physical preparation. However, the cited study had an adequate time series length for each participant, allowing for local and global evaluations of the system’s dynamics.

5. Conclusions

The approach shown in this paper made it possible to evaluate a group of dancers as a single, global organism, performing two types of pirouettes. The jazz pirouette, and more precisely, the movement of the knee joint of the supporting limb, which was assessed in the way presented in the paper, is more compact and reproducible (smaller convex hull volume values and large REC values) compared to the classic pirouette. The emergent dynamics of the dance movement, and the established similarities between the dancers, suggest that the global level of movement may somehow be disconnected from the lower-level neuromuscular dynamics that generate it. Thus, the results show that the global behavior of the pirouette, instead of the local, is individual to the behavior of the dancers. This approach has not been applied in the literature and might be controversial. Overall, the reconstruction of the attractors, the evaluation of the Lyapunov exponent and the Hurst exponent, and the application of RQA appear to provide a comprehensive view of the general motion presented by the dancers.

Funding

This research was funded by the Ministry of Science and Higher Education during 2020–2022 under Research Group No. 3 at the Józef Pilsudski University of Physical Education in Warsaw, “Motor system diagnostics in selected dysfunctions as a basis for planning the rehabilitation process”.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Józef Piłsudski University of Physical Education in Warsaw (protocol code SEK 01-09/2020, date of approval 01.09.2020).

Data Availability Statement

The measurement data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgments

I wanted to thank MohammadHesam Azadjou from University of Southern California, Valero Lab for the idea of calculating the convex hull for the attractor.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stergiou, N.; Decker, L.M. Human movement variability, nonlinear dynamics, and pathology: Is there a connection? Hum. Mov. Sci. 2011, 30, 869–888. [Google Scholar] [CrossRef] [PubMed]
  2. Stergiou, N. Nonlinear Analysis for Human Movement Variability; CRC Press: Boca Raton, FL, USA, 2016; 388p. [Google Scholar]
  3. Bernstein, N. A new method of mirror cyclographie and its application towards the study of labor movements during work on a workbench. Hyg. Saf. Pathol. Labor 1930, 5, 3–9. [Google Scholar]
  4. Riley, M.A.; Turvey, M.T. Variability and Determinism in Motor Behavior. J. Mot. Behav. 2002, 34, 99–125. [Google Scholar] [CrossRef] [PubMed]
  5. van Emmerik, R.E.A.; Ducharme, S.W.; Amado, A.C.; Hamill, J. Comparing dynamical systems concepts and techniques for biomechanical analysis. J. Sport Health Sci. 2016, 5, 3–13. [Google Scholar] [CrossRef]
  6. van Mourik, A.M.; Daffertshofer, A.; Beek, P.J. Deterministic and stochastic features of rhythmic human movement. Biol. Cybern. 2006, 94, 233–244. [Google Scholar] [CrossRef]
  7. Weber, I.; Oehrn, C.R. NoLiTiA: An Open-Source Toolbox for Non-linear Time Series Analysis. Front. Neuroinform. 2022, 16, 876012. [Google Scholar] [CrossRef]
  8. Donker, S.F.; Roerdink, M.; Greven, A.J.; Beek, P.J. Regularity of center-of-pressure trajectories depends on the amount of attention invested in postural control. Exp. Brain Res. 2007, 181, 1–11. [Google Scholar] [CrossRef]
  9. Roerdink, M.; Hlavackova, P.; Vuillerme, N. Center-of-pressure regularity as a marker for attentional investment in postural control: A comparison between sitting and standing postures. Hum. Mov. Sci. 2011, 30, 203–212. [Google Scholar] [CrossRef]
  10. Ramdani, S.; Tallon, G.; Bernard, P.L.; Blain, H. Recurrence quantification analysis of human postural fluctuations in older fallers and non-fallers. Ann. Biomed. Eng. 2013, 41, 1713–1725. [Google Scholar] [CrossRef]
  11. Bradley, E.; Kantz, H. Nonlinear time-series analysis revisited. Chaos 2015, 25, 097610. [Google Scholar] [CrossRef]
  12. Webber, C.; Zbilut, J. Recurrence quantification analysis of nonlinear dynamical systems. Tutor. Contemp. Nonlinear Methods Behav. Sci. 2005, 94, 26–94. [Google Scholar]
  13. Kędziorek, J.; Błażkiewicz, M. Nonlinear Measures to Evaluate Upright Postural Stability: A Systematic Review. Entropy 2020, 22, 1357. [Google Scholar] [CrossRef]
  14. Khayat, O.; Nowshiravan-Rahatabad, F. Complex feature analysis of center of pressure signal for age-related subject classification. Ann. Mil. Health Sci. Res. 2014, 12, e63518. [Google Scholar]
  15. Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov exponents from a time series. Physica 1985, 16, 285–317. [Google Scholar] [CrossRef]
  16. Rosenstein, M.T.; Collins, J.J.; De Luca, C.J. A practical method for calculating largest Lyapunov exponents from small data sets. Phys. D Nonlinear Phenom. 1993, 65, 117–134. [Google Scholar] [CrossRef]
  17. Liu, H.-F.; Dai, Z.-H.; Li, W.-F.; Gong, X.; Yu, Z.-H. Noise robust estimates of the largest Lyapunov exponent. Phys. Lett. A 2005, 341, 119–127. [Google Scholar] [CrossRef]
  18. Kantz, H. A robust method to estimate the maximal Lyapunov exponent of a time series. Phys. Lett. A 1994, 185, 77–87. [Google Scholar] [CrossRef]
  19. Rispens, S.M.; Pijnappels, M.; van Dieën, J.H.; van Schooten, K.S.; Beek, P.J.; Daffertshofer, A. A benchmark test of accuracy and precision in estimating dynamical systems characteristics from a time series. J. Biomech. 2014, 47, 470–475. [Google Scholar] [CrossRef]
  20. Hurst, H.E. Long-Term Storage Capacity of Reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–799. [Google Scholar] [CrossRef]
  21. Tarnopolski, M. On the relationship between the Hurst exponent, the ratio of the mean square successive difference to the variance, and the number of turning points. Phys. A Stat. Mech. Its Appl. 2016, 461, 662–673. [Google Scholar] [CrossRef]
  22. Devaney, R. An Introduction to Chaotic Dynamical Systems; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  23. Takens, F. Detecting strange attractors in turbulence. In Proceedings of the Dynamical Systems and Turbulence, Warwick 1980; Springer: Berlin/Heidelberg, Germany, 1981; pp. 366–381. [Google Scholar]
  24. Fraser, A.M.; Swinney, H.L. Independent coordinates for strange attractors from mutual information. Phys. Rev. A 1986, 33, 1134–1140. [Google Scholar] [CrossRef]
  25. Cao, L. Practical method for determining the minimum embedding dimension of a scalar time series. Phys. D Nonlinear Phenom. 1997, 110, 43–50. [Google Scholar] [CrossRef]
  26. Grassberger, P.; Procaccia, I. Measuring the strangeness of strange attractors. Phys. D Nonlinear Phenom. 1983, 9, 189–208. [Google Scholar] [CrossRef]
  27. Sauer, T.; Yorke, J.A.; Casdagli, M. Embedology. J. Stat. Phys. 1991, 65, 579–616. [Google Scholar] [CrossRef]
  28. Kennel, M.B.; Brown, R.; Abarbanel, H.D.I. Determining embedding dimension for phase-space reconstruction using a geometrical construction. Phys. Rev. A 1992, 45, 3403–3411. [Google Scholar] [CrossRef]
  29. de Pedro-Carracedo, J.; Fuentes-Jimenez, D.; Ugena, A.M.; Gonzalez-Marcos, A.P. Phase Space Reconstruction from a Biological Time Series: A Photoplethysmographic Signal Case Study. Appl. Sci. 2020, 10, 1430. [Google Scholar] [CrossRef]
  30. Brick, T.R.; Gray, A.L.; Staples, A.D. Recurrence Quantification for the Analysis of Coupled Processes in Aging. J. Gerontol. Ser. B 2018, 73, 134–147. [Google Scholar] [CrossRef]
  31. Harbourne, R.T.; Stergiou, N. Movement variability and the use of nonlinear tools: Principles to guide physical therapist practice. Phys. Ther. 2009, 89, 267–282. [Google Scholar] [CrossRef]
  32. Pellecchia, G.L.; Shockley, K. Application of recurrence quantification analysis: Influence of cognitive activity on postural fluctuations. In Tutorials in Contemporary Nonlinear Methods for the Behavioral Sciences; National Science Foundation: Alexandria, VA, USA, 2005; pp. 95–141. [Google Scholar]
  33. Som, A.; Krishnamurthi, N.; Venkataraman, V.; Turaga, P. Attractor-shape descriptors for balance impairment assessment in Parkinson’s disease. IEEE Eng. Med. Biol. Soc. Annu. Int. Conf. 2016, 2016, 3096–3100. [Google Scholar] [CrossRef]
  34. Gates, D.H.; Dingwell, J.B. Comparison of different state space definitions for local dynamic stability analyses. J. Biomech. 2009, 42, 1345–1349. [Google Scholar] [CrossRef]
  35. Bradley, E.; Stuart, J. Using chaos to generate variations on movement sequences. Chaos 1998, 8, 800–807. [Google Scholar] [CrossRef] [PubMed]
  36. Peppoloni, L.; Lawrence, E.L.; Ruffaldi, E.; Valero-Cuevas, F.J. Characterization of the disruption of neural control strategies for dynamic fingertip forces from attractor reconstruction. PLoS ONE 2017, 12, e0172025. [Google Scholar] [CrossRef] [PubMed]
  37. Aston, P.J.; Nandi, M.; Christie, M.I.; Huang, Y.H. Comparison of Attractor Reconstruction and HRV Methods for Analysing Blood Pressure Data. Comput. Cardiol. 2014, 41, 437–440. [Google Scholar] [CrossRef] [Green Version]
  38. Lyle, J.; Charlton, P.; Bonet-Luz, E.; Chaffey, G.; Christie, M.; Nandi, M.; Aston, P. Beyond HRV: Analysis of ECG Signals Using Attractor Reconstruction. In Proceedings of the 2017 Computing in Cardiology (CinC), Rennes, France, 24–27 September 2017. [Google Scholar]
  39. Thanaj, M.; Chipperfield, A.J.; Clough, G.F. Attractor Reconstruction Analysis for Blood Flow Signals. IEEE Eng. Med. Biol. Soc. Annu. Int. Conf. 2019, 2019, 2281–2284. [Google Scholar] [CrossRef]
  40. van den Hoorn, W.; Hodges, P.W.; van Dieën, J.H.; Kerr, G.K. Reliability of recurrence quantification analysis of postural sway data. A comparison of two methods to determine recurrence thresholds. J. Biomech. 2020, 107, 109793. [Google Scholar] [CrossRef] [PubMed]
  41. Labini, F.S.; Meli, A.; Ivanenko, Y.P.; Tufarelli, D. Recurrence quantification analysis of gait in normal and hypovestibular subjects. Gait Posture 2012, 35, 48–55. [Google Scholar] [CrossRef]
  42. Riley, M.A.; Clark, S. Recurrence analysis of human postural sway during the sensory organization test. Neurosci. Lett. 2003, 342, 45–48. [Google Scholar] [CrossRef]
  43. Zimatore, G.; Gallotta, M.C.; Innocenti, L.; Bonavolontà, V.; Ciasca, G.; De Spirito, M.; Guidetti, L.; Baldari, C. Recurrence quantification analysis of heart rate variability during continuous incremental exercise test in obese subjects. Chaos 2020, 30, 033135. [Google Scholar] [CrossRef]
  44. Dimitriev, D.; Saperova, E.V.; Dimitriev, A.; Karpenko, Y. Recurrence Quantification Analysis of Heart Rate During Mental Arithmetic Stress in Young Females. Front. Physiol. 2020, 11, 40. [Google Scholar] [CrossRef]
  45. Billeci, L.; Marino, D.; Insana, L.; Vatti, G.; Varanini, M. Patient-specific seizure prediction based on heart rate variability and recurrence quantification analysis. PLoS ONE 2018, 13, e0204339. [Google Scholar] [CrossRef]
  46. Błażkiewicz, M. Joint loads and muscle force distribution during classical and jazz pirouettes. Acta Bioeng. Biomech. 2021, 23, 3–13. [Google Scholar] [CrossRef]
  47. Bläsing, B.; Schack, T. Mental Representation of Spatial Movement Parameters in Dance. Spat. Cogn. Comput. 2012, 12, 111–132. [Google Scholar] [CrossRef]
  48. Vaganova, A. Grundlagen des Klassischen Tanzes [Basic Principles of Classical Ballet]; Henschel: Berlin, Germany, 2002. [Google Scholar]
  49. Kleinow, T. Testing Continuous Time Models in Financial Markets. Ph.D Thesis, Humboldt-Universität zu Berlin, Wirtschaftswissenschaftliche Fakultät, Berlin, Germany, 2002. [Google Scholar]
  50. Kugiumtzis, D.; Tsimpiris, A. Measures of Analysis of Time Series (MATS): A MATLAB Toolkit for Computation of Multiple Measures on Time Series Data Bases. J. Stat. Softw. 2010, 33, 1–30. [Google Scholar] [CrossRef] [Green Version]
  51. Zhivomirov, H.; Nedelchev, I. A Method for Signal Stationarity Estimation. Rom. J. Acoust. Vib. 2020, 17, 149–155. [Google Scholar]
  52. Mallika, M.C.; Suriya Prabhaa, S.; Asokan, K.; Anil Kumar, K.S.; Ramamohan, T.R.; Kumar, K.S. Testing for nonlinearity in nonstationary time series: A network-based surrogate data test. Phys. Rev. E 2021, 104, 054217. [Google Scholar] [CrossRef]
  53. Schreiber, T.; Schmitz, A. Improved Surrogate Data for Nonlinearity Tests. Phys. Rev. Lett. 1996, 77, 635–638. [Google Scholar] [CrossRef]
  54. Coloma, W.A.; Costa, A.L.; Pereira, C.; da Silva, C.A.M. Time Series Analysis for BWR Stability Studies. Nucl. Technol. 2020, 206, 554–564. [Google Scholar] [CrossRef]
  55. Wallot, S.; Mønster, D. Calculation of Average Mutual Information (AMI) and False-Nearest Neighbors (FNN) for the Estimation of Embedding Parameters of Multidimensional Time Series in Matlab. Front. Psychol. 2018, 9, 1679. [Google Scholar] [CrossRef]
  56. Lorenz, E.N. Deterministic nonperiodic flow. J. Atmos. Sci. 1963, 20, 130–141. [Google Scholar] [CrossRef]
  57. Barber, C.B.; Dobkin, D.P.; Huhdanpaa, H. The quickhull algorithm for convex hulls. ACM Trans. Math. Softw. 1996, 22, 469–483. [Google Scholar] [CrossRef]
  58. Baíllo, A.; Chacón, J.E. Chapter 1—Statistical outline of animal home ranges: An application of set estimation. In Handbook of Statistics; Srinivasa Rao, A.S.R., Rao, C.R., Eds.; Elsevier: Amsterdam, The Netherlands, 2021; Volume 44, pp. 3–37. [Google Scholar]
  59. Zhao, F. Extracting and representing qualitative behaviors of complex systems in phase space. Artif. Intell. 1994, 69, 51–92. [Google Scholar] [CrossRef]
  60. Marwan, N.; Carmen Romano, M.; Thiel, M.; Kurths, J. Recurrence plots for the analysis of complex systems. Phys. Rep. 2007, 438, 237–329. [Google Scholar] [CrossRef]
  61. Jackson, E.S.; Tiede, M.; Riley, M.A.; Whalen, D.H. Recurrence Quantification Analysis of Sentence-Level Speech Kinematics. J. Speech Lang. Hear. Res. JSLHR 2016, 59, 1315–1326. [Google Scholar] [CrossRef]
  62. Yang, H. Tool Box of Recurrence Plot and Recurrence Quantification Analysis. Available online: https://www.mathworks.com/matlabcentral/fileexchange/58246-tool-box-of-recurrence-plot-and-recurrence-quantification-analysis (accessed on 2 June 2022).
  63. Yang, H. Multiscale recurrence quantification analysis of spatial cardiac vectorcardiogram signals. IEEE Trans. Biomed. Eng. 2011, 58, 339–347. [Google Scholar] [CrossRef]
  64. Chen, Y.; Yang, H. Multiscale recurrence analysis of long-term nonlinear and nonstationary time series. Chaos Solitons Fractals 2012, 45, 978–987. [Google Scholar] [CrossRef]
  65. Nayak, S.K.; Bit, A.; Dey, A.; Mohapatra, B.; Pal, K. A Review on the Nonlinear Dynamical System Analysis of Electrocardiogram Signal. J. Healthc. Eng. 2018, 2018, 6920420. [Google Scholar] [CrossRef]
  66. Zbilut, J.P.; Zaldivar-Comenges, J.-M.; Strozzi, F. Recurrence quantification based Liapunov exponents for monitoring divergence in experimental data. Phys. Lett. A 2002, 297, 173–181. [Google Scholar] [CrossRef]
  67. Riley, M.A.; Balasubramaniam, R.; Turvey, M.T. Recurrence quantification analysis of postural fluctuations. Gait Posture 1999, 9, 65–78. [Google Scholar] [CrossRef]
  68. Biringen, E.K. Analysis of pirouette execution for improved performance. Med. Probl. Perform. Artist. 2010, 25, 136–137. [Google Scholar] [CrossRef]
  69. Kantz, H.; Schreiber, T. Nonlinear Time Series Analysis; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  70. Zbilut, J.P.; Giuliani, A.; Webber, C.L. Detecting deterministic signals in exceptionally noisy environments using cross-recurrence quantification. Phys. Lett. A 1998, 246, 122–128. [Google Scholar] [CrossRef]
  71. Zbilut, J.P.; Giuliani, A.; Webber, C.L. Recurrence quantification analysis and principal components in the detection of short complex signals. Phys. Lett. A 1998, 237, 131–135. [Google Scholar] [CrossRef]
  72. Webber, C.L.; Marwan, N.; Facchini, A.; Giuliani, A. Simpler methods do it better: Success of Recurrence Quantification Analysis as a general purpose data analysis tool. Phys. Lett. A 2009, 373, 3753–3756. [Google Scholar] [CrossRef]
  73. Mehdizadeh, S. The largest Lyapunov exponent of gait in young and elderly individuals: A systematic review. Gait Posture 2018, 60, 241–250. [Google Scholar] [CrossRef] [PubMed]
  74. Newell, K.; Vaillancourt, D. Dimensional change in motor learning. Hum. Mov. Sci. 2001, 20, 695–715. [Google Scholar] [CrossRef]
  75. Stergiou, N.; Harbourne, R.T.; Cavanaugh, J.T. Optimal Movement Variability: A New Theoretical Perspective for Neurologic Physical Therapy. J. Neurol. Phys. Ther. 2006, 30, 120–129. [Google Scholar] [CrossRef] [PubMed]
  76. Curtin, P.; Curtin, A.; Austin, C.; Gennings, C.; Tammimies, K.; Bölte, S.; Arora, M. Recurrence quantification analysis to characterize cyclical components of environmental elemental exposures during fetal and postnatal development. PLoS ONE 2017, 12, e0187049. [Google Scholar] [CrossRef] [PubMed]
  77. Burton, S.; Vicinanza, D.; Exell, T.; Newell, K.M.; Irwin, G.; Williams, G.K.R. Attractor dynamics of elite performance: The high bar longswing. Sports Biomech. 2021, 1–14. [Google Scholar] [CrossRef]
Figure 1. Participant model during a pirouette. (A) Motion sequences; (B) the location of the center of mass (CoM) and position of the marker on the knee’s lateral epicondyle of the supporting leg (LKNE) for the classic and jazz pirouette. The coordinate system, where: x—horizontal axis (anterior-posterior), y—vertical axis (inferior-superior), and z—horizontal axis (mediolateral).
Figure 1. Participant model during a pirouette. (A) Motion sequences; (B) the location of the center of mass (CoM) and position of the marker on the knee’s lateral epicondyle of the supporting leg (LKNE) for the classic and jazz pirouette. The coordinate system, where: x—horizontal axis (anterior-posterior), y—vertical axis (inferior-superior), and z—horizontal axis (mediolateral).
Entropy 24 01310 g001
Figure 2. Time series for LKNE marker and CoM for x, y, and z coordinates for jazz and classic pirouettes (A,C) for a single individual after the resampling procedure; (B,D) reconstruction of the time series before applying the quantiles function.
Figure 2. Time series for LKNE marker and CoM for x, y, and z coordinates for jazz and classic pirouettes (A,C) for a single individual after the resampling procedure; (B,D) reconstruction of the time series before applying the quantiles function.
Entropy 24 01310 g002
Figure 3. The global (all D) and directional (xD, yD, zD) embedding dimensions as an output of the mdFnn function for the (A) LKNE marker for jazz pirouettes; (B) CoM for jazz pirouettes; (C) LKNE marker for classic pirouettes, and (D) CoM for classic pirouettes.
Figure 3. The global (all D) and directional (xD, yD, zD) embedding dimensions as an output of the mdFnn function for the (A) LKNE marker for jazz pirouettes; (B) CoM for jazz pirouettes; (C) LKNE marker for classic pirouettes, and (D) CoM for classic pirouettes.
Entropy 24 01310 g003
Figure 4. The global (all tau) and directional (xtau, ytau, ztau) time delays as the output of the mdDelay function for (A) LKNE marker for jazz pirouette; (B) CoM for jazz pirouette; (C) LKNE marker for classic pirouette; (D) CoM for classic pirouette. The threshold was set at 1/e.
Figure 4. The global (all tau) and directional (xtau, ytau, ztau) time delays as the output of the mdDelay function for (A) LKNE marker for jazz pirouette; (B) CoM for jazz pirouette; (C) LKNE marker for classic pirouette; (D) CoM for classic pirouette. The threshold was set at 1/e.
Entropy 24 01310 g004
Figure 5. The reconstructed phase space (odd lines) with a global time lag (τ = all tau), convex hulls, and their volumes (even lines) for the (A) LKNE marker and (B) CoM x, y, and z coordinates for jazz and classic pirouettes (x-anterior-posterior, y-inferior-superior, z-mediolateral). The tables show convex hull volumes for global and directional parameters (D and tau).
Figure 5. The reconstructed phase space (odd lines) with a global time lag (τ = all tau), convex hulls, and their volumes (even lines) for the (A) LKNE marker and (B) CoM x, y, and z coordinates for jazz and classic pirouettes (x-anterior-posterior, y-inferior-superior, z-mediolateral). The tables show convex hull volumes for global and directional parameters (D and tau).
Entropy 24 01310 g005
Figure 6. The illustration of the multiscale recurrence analysis for the (A) LKNE and (B) CoM time series separately for the x-anterior-posterior, y-inferior-superior, z-mediolateral directions for jazz and classic pirouette. Graphs obtained for the global time delay and embedding dimensions.
Figure 6. The illustration of the multiscale recurrence analysis for the (A) LKNE and (B) CoM time series separately for the x-anterior-posterior, y-inferior-superior, z-mediolateral directions for jazz and classic pirouette. Graphs obtained for the global time delay and embedding dimensions.
Entropy 24 01310 g006
Table 1. Methods for calculating the time delay τ (tau) and the embedding dimension (D) required for phase space reconstruction.
Table 1. Methods for calculating the time delay τ (tau) and the embedding dimension (D) required for phase space reconstruction.
Phase Space Reconstruction State Vector:
y(t) = (x(t), x(t + τ), x(t + 2τ), …, x(t + (D − 1)τ)) from Original Time Series Data x(t) by Using Time Delay (τ) and Dimension of the Attractor D
Possible ways of τ selection [24]:
(1) Autocorrelation function
(2) Mutual information
Possible ways of D selection:
(1) Principal Component Analysis (PCA) [25]
(2) Correlation dimension [26]
(3) Box-counting [27]
(4) False Nearest Neighbor (FNN) [28]
Table 2. Mean and standard deviation of anthropometric parameters describing the dancers.
Table 2. Mean and standard deviation of anthropometric parameters describing the dancers.
GroupAge [Years]Body Mass [kg]Body Height [m]Training Period [Years]
N = 1522.13 ± 2.7357.56 ± 6.761.68 ± 0.6212.19 ± 3.04
Table 3. The Hurst exponent values calculated for the x, y, and z coordinates of the reconstructed LKNE and CoM time series for classic and jazz pirouettes.
Table 3. The Hurst exponent values calculated for the x, y, and z coordinates of the reconstructed LKNE and CoM time series for classic and jazz pirouettes.
LKNExLKNEyLKNEzCoMxCoMyCoMz
Classic0.750.710.710.750.720.70
Jazz0.760.660.700.730.690.75
Table 4. The Lyapunov exponent values calculated for the x, y, and z coordinates of the reconstructed LKNE and CoM time series for classic and jazz pirouettes (for the directional dimensions and lags from Figure 3 and Figure 4).
Table 4. The Lyapunov exponent values calculated for the x, y, and z coordinates of the reconstructed LKNE and CoM time series for classic and jazz pirouettes (for the directional dimensions and lags from Figure 3 and Figure 4).
LKNExLKNEyLKNEzCoMxCoMyCoMz
Classic1.531.432.052.111.441.3
Jazz1.431.132.572.041.791.69
Table 5. Recurrence variables calculated for global and directional time lag and the embedding dimension for the LKNE marker and CoM motion in x-, y-, and z-directions for J—jazz and C—classic pirouettes. REC—% of recurrence, DET—%determinism, LMAX—the length of the longest diagonal line segment in the plot, ENT—Shannon information entropy, LAM—%laminarity, TT—trapping time. The prefix ‘all’ means global, and d denotes directional.
Table 5. Recurrence variables calculated for global and directional time lag and the embedding dimension for the LKNE marker and CoM motion in x-, y-, and z-directions for J—jazz and C—classic pirouettes. REC—% of recurrence, DET—%determinism, LMAX—the length of the longest diagonal line segment in the plot, ENT—Shannon information entropy, LAM—%laminarity, TT—trapping time. The prefix ‘all’ means global, and d denotes directional.
allRECdRECallDETdDETallLMAXdLMAXallENTdENTallLAMdLAMallTTdTT
LKNEx_J8.618.0899.5999.60101810195.755.8399.7499.7716.5419.95
LKNEy_J43.3441.299.899.85221222436.97.199.999.9354.0865.31
LKNEz_J90.4182.4699.9399.96393539727.548.1799.9799.98163.39223.12
LKNEx_C1.61.6299.699.61124613055.695.7598.3399.659.715.36
LKNEy_C16.320.199.7299.86210922106.467.0599.8799.9443.4577.38
LKNEz_C10.518.5199.8199.92242525146.447.2399.999.9736.3667.22
CoMx_J22.8799.7399.769159456.296.3099.7499.8516.220.82
CoMy_J11.9718.1699.9199.93393939677.237.3799.9599.9737.346.92
CoMz_J86.3586.5299.9899.98393939638.238.2699.9999.99261.59262
CoMx_C1.121.6399.6599.74153616006.256.3299.6499.8512.3620.96
CoMy_C9.8412.7499.8799.93392739876.837.1499.9599.9838.0455.04
CoMz_C32.5942.1899.9399.94392739547.157.3299.9799.9876.289.77
Table 6. The mean percentage difference between recurrence measures calculated for the global and directional settings.
Table 6. The mean percentage difference between recurrence measures calculated for the global and directional settings.
RECDETLMAXENTLAMTT
Jazz19.470.021.172.520.0322.07
Classic34.220.073.275.080.2958.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Błażkiewicz, M. Evaluation of Geometric Attractor Structure and Recurrence Analysis in Professional Dancers. Entropy 2022, 24, 1310. https://doi.org/10.3390/e24091310

AMA Style

Błażkiewicz M. Evaluation of Geometric Attractor Structure and Recurrence Analysis in Professional Dancers. Entropy. 2022; 24(9):1310. https://doi.org/10.3390/e24091310

Chicago/Turabian Style

Błażkiewicz, Michalina. 2022. "Evaluation of Geometric Attractor Structure and Recurrence Analysis in Professional Dancers" Entropy 24, no. 9: 1310. https://doi.org/10.3390/e24091310

APA Style

Błażkiewicz, M. (2022). Evaluation of Geometric Attractor Structure and Recurrence Analysis in Professional Dancers. Entropy, 24(9), 1310. https://doi.org/10.3390/e24091310

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop