Interpolation and sampling effects on recurrence quantification measures.

The methods of recurrence plots (RPs) and recurrence quantification analysis (RQA) have been used to investigate the tribosystem. The morphology of RPs and RQA measures are strongly dependent on the embedding parameters of the recursive matrix and the segment sizes of the time-series. To improve the calculation accuracy of recursive characteristics analysis, the influences of the embedding parameters and segment sizes on the morphology of RPs and RQA measures have been studied in this letter. Three kinds of theoretical chaotic time-series and measured coefficient of friction (COF) signals during the running-in process were chosen as research objects, and the morphology of RPs and RQA measures were obtained using CRP toolbox afterward. The results indicate that no embedding was actually needed if the data sets are to be qualitatively analyzed using RPs and RQA. Additionally, the morphology of RPs and RQA measures are sensitive to the segment sizes for theoretical chaotic time-series, while the RQA measures of COF signal in the steady-state period are rather stable due to its self-similarity. Finally, according to the guidelines of the parameter settings, the dynamical evolution of measured COF signals during the running-in process have been investigated. It is indicated that recursive characteristics of COF signals could reveal the tribological behaviors’ evolution and conduct the running-in status identification. The results in this paper are significant for improving the calculation accuracy and saving computational time when using the method of recursive characteristics analysis on the tribological behaviors.

Recurrence quantification analysis (RQA) was developed in order to quantify differently appearing recurrence plots (RPs) based on their small-scale structures, which generally indicate the number and duration of recurrences in a dynamical system. Although RQA measures are traditionally employed in analyzing complex systems and identifying transitions, recent work has shown that they can also be used for pairwise dissimilarity comparisons of time series. We explain why RQA is not only a modern method for nonlinear data analysis but also is a very promising technique for various time series mining tasks.

Purpose: To investigate if there exists a difference in breathing patterns between patients treated with SBRT and IMRT using RQA and DC (using K2: correlation entropy and D2: correlation dimension) measures. Methods: 9 patients treated with SBRT and 8 treated with IMRT were scanned with 4D CT and the breathing patterns were acquired. One of the SBRT patients was scanned with and without meditation. Each breathing signal consisted of a scalar time series and recurrence quantification analysis (RQA) was utilized to determine the following measures: Periodicity of the system as percentage of recurrence points (%RR), determinism (DET), maximum diagonal line length (L_max) whose inverse the divergence (DIV) is measure for how fast trajectories diverge from each other, the average diagonal line length (L) that can be interpreted as the mean prediction time of the signal, and the entropy (ENTR) a measure of information complexity. In addition the invariant measures of K2 and D2 were also estimated. A locally nonlinear forecast was applied to predict future breathing signals of N time step ahead for the patient with and without meditation. Results: Our results showed %RR has significant correlation with L_max and has inverse correlation with DIV. DET has significant correlation with Lmax, L, ENTR and DIV. Independent t test suggests there is no difference between the SBRT and IMRT groups in terms of the RQA measures and K2. Patient that had undergone meditation showed improvement in %RR, L_max, DIV, K2 and had an estimated correlation dimension of 1.7. Prediction showed similar results for one and three time step ahead but meditation one had better prediction horizon when time step was higher. Conclusion: RQA is a powerful tool that allows one to analyze the dynamic nature of breathing pattern. No significant difference was found in the dynamical complexity of SBRT and IMRT patients.

In the paper, application of the Recurrence Plots (RP) method to damage detection and localization in the beam structure is proposed. The RP is a well-established method of nonlinear data analysis dedicated to investigation into changes in dynamic behaviour of mechanical systems during their normal work. Interpretation of information provided by the recurrence plots can be supported by the estimation of the measures of Recurrence Quantification Analysis (RQA).

This paper presents efficient low-power adaptive filter architecture for electroencephalogram (EEG) signal and epileptic seizure detection using recurrence quantification analysis (RQA). The preprocessing of EEG is done using notch, wavelet and adaptive filter. The comparison of signal-to-noise ratio of the filter outputs proves that the adaptive filter provides better performance and a parallel interleaved sample of direct form adaptive FIR filter architecture is implemented and low-power issues are addressed. An innovative compressor-based addition technique is utilized in the filter implementation to reduce the area and power consumption. The design is developed using Verilog HDL and mapped to 65-nm technological node. The power results are compared with conventional architecture of adaptive filter. The major advantage of choosing RQA is that it provides better information even for short non-stationary and nonlinear signals where other methods fail to provide good results. And it requires no conventions about data set size or dispersal of the data. The algorithm is applied on epileptic EEG signal from CHB-MIT database. The RQA measures are determined from the recurrence plot, its performance is measured in terms of sensitivity and specificity as 97.4 and 93.5 %, respectively, and leakage power is reduced to 10 %.

<p>Grassland ecosystems are extremely complex and set up intricate structures, whose characteristics and dynamic properties are greatly influenced by climate and meteorological patterns. Climate change and global warming are factors that could impact negatively in the quality and productivity of these ecosystems.</p><p>Remote sensing techniques have been demonstrated as a powerful tool for monitoring extensive areas. In this study, two semi-arid grassland plots were selected in the centre of Spain. This region is characterized by low precipitation and moderate productivity per unit. Through scientific research, spectral vegetation indices (VIs) have been developed to characterize vegetation cover. The most common VI is the Normalized Difference Vegetation Index (NDVI). However, in vegetation scarcity conditions, bare soil reflectance is increased, and the feasibility of NDVI is reduced. This study aims to perform a method to compare soil and agro-climatic variables effect on vegetation time-series indices.</p><p>The construction of the time series was based on multispectral images of MODIS TERRA (MOD09A1.006) product acquired from 2002 till 2018. Three pixels with a temporal resolution of 8 days and a spatial resolution of 500 x 500 m were chosen in each area. To estimate and analyse VIs series, Red (620-670 nm) and Near Infrared (841-876 nm) channels were extracted and filtered by the quality of pixel. All spectral bands showed statistically significant differences confirming that both areas presented different soil properties. Moreover, average annual precipitation was different in each area of study.</p><p>NDVI calculation is only based on NIR and RED bands. To improve the estimation of vegetation in semi-arid areas, several indices have been developed to minimize the soil effect. Each one of them incorporates soil influence in a different way, i.e., Soil Adjusted Vegetation Index (SAVI) adds a constant soil adjustment factor (L), whereas, MSAVI, incorporate an L variable and dependant on soil characteristics.</p><p>Recurrence plots (RP) and recurrence quantification analysis (RQA) were computed to characterize the influence of agro-climatic variables in vegetation index dynamics. Characterization was based on various RQA measures, such as Determinism (DET), average diagonal length (LT) or entropy (ENT).</p><p>Our results showed different RPs depending on the area, VI utilized and precipitation. MSAVI patterns were further distinct, meanwhile, NDVI showed a noisy pattern. LT values in MSAVI were higher than in SAVI implying that MSAVI recurrent events are much longer than SAVI. Simultaneously, LT and DET values in ZSO, with a higher rain, were above ZEA values in MSAVI.</p><p>This indicates that incorporating more detailed information of soil and precipitation reinforce vegetation index estimation and allow to obtain a more distinct pattern of the time series. Therefore, in arid-semiarid grasslands, they should be considered.</p><p><strong>ACKNOWLEDGEMENTS</strong></p><p>The authors acknowledge support from Project No. PGC2018-093854-B-I00 of the Spanish <em>Ministerio de Ciencia Innovación y Universidades</em> of Spain and the funding from the Comunidad de Madrid (Spain) and Structural Funds 2014-2020 512 (ERDF and ESF), through project AGRISOST-CM S2018/BAA-4330, are highly appreciated.</p>

In this paper the dynamical properties of ARP data for anomaly of internet traffic detection are presented. The local network traffic dynamics for single workstation has been observed experimentally on networks with 42 workstations. In the paper RQA (Recurrence Quantification Analysis) measures have been analyzed. Three of them have been explored using RP windowing method: recurrence rate (RR), divergence (DIV) and the longest vertical line (VMAX). The RP windowing algorithm samples the explored ARP time series, as a result each RP window characterized the disintegration of recurrence plots. The RR and VMAX defined the low activity of ARP dataflow. Inversely, the DIV identified the high activity of ARP time series. The authors have proved that both, the recurrence plot technique and the RP windowing method as the tools in the analysis of network traffic anomaly detection are effective and accurate.

The purpose of this study was to apply recurrence plots and recurrence quantification analysis (RQA) on event related potentials (ERPs) recorded during memory recognition tests. Data recorded during memory retrieval in four scalp region was used. Tow most important ERP's components corresponding to memory retrieval, FN400 and LPC, were detected in recurrence plots computed for single-trial EEGs. In addition, the RQA was used to quantify changes in signal dynamic structure during memory retrieval, and measures of complexity as RQA variables were computed. Given the stimulus, amplitude of the RQA variables increases around 400ms, corresponding to dimension reduction of system. Furthermore, after 800ms these amplitudes decreased which can be as a consequence of an increase in system dimension and back to its basic state. The mean amplitude of Old items was more than New one. Using this method, we found its ability to detect memory components of EEG signals and do a distinction between Old/ New items. In contrast with linear techniques recurrence plots and RQA do not need large number of recorded trials, and they can indicate changes in even single-trial EEGs. RQA can also show differences between old and new events in a memory process.

Recurrence quantification analysis (RQA) is a time-series analysis method that uses autocorrelation properties of typing data to detect regularities within the writing process. The following paper first gives a detailed introduction to RQA and its application to time series data. We then apply RQA to keystroke logging data of first and foreign language writing to illustrate how outcome measures of RQA can be understood as skill-driven constraints on keyboard typing performance. Forty native German students performed two prompted writing assignments, one in German and one in English, a standardized copy task, and a standardized English placement test. We assumed more fluent and skilled writing to reveal more structured typing time series patterns. Accordingly, we expected writing in a well-mastered first language to coincide with higher values in relevant RQA measures as compared to writing in a foreign language. Results of mixed model ANOVAs confirmed our hypothesis. We further observed that RQA measures tend to be higher, thus indicating more structured data, whenever parameters of pause, burst, and revision analyses indicate more fluent writing. Multiple regression analyses revealed that, in addition to typing skills, language proficiency significantly predicts outcomes of RQA. Thus, the present data emphasize RQA being a valuable resource for studying time series data that yields meaningful information about the effort a writer must exert during text production.

Detecting intrinsic dynamics of traffic flow with recurrence analysis and empirical mode decomposition

Our contribution focuses on the characterization of sleep apnea from a cardiac rate point of view, using Recurrence Quantification Analysis (RQA), based on a Heart Rate Variability (HRV) feature selection process. Three parameters are crucial in RQA: those related to the embedding process (dimension and delay) and the threshold distance. There are no overall accepted parameters for the study of HRV using RQA in sleep apnea. We focus on finding an overall acceptable combination, sweeping a range of values for each of them simultaneously. Together with the commonly used RQA measures, we include features related to recurrence times, and features originating in the complex network theory. To the best of our knowledge, no author has used them all for sleep apnea previously. The best performing feature subset is entered into a Linear Discriminant classifier. The best results in the “Apnea-ECG Physionet database” and the “HuGCDN2014 database” are, according to the area under the receiver operating characteristic curve, 0.93 (Accuracy: 86.33%) and 0.86 (Accuracy: 84.18%), respectively. Our system outperforms, using a relatively small set of features, previously existing studies in the context of sleep apnea. We conclude that working with dimensions around 7–8 and delays about 4–5, and using for the threshold distance the Fixed Amount of Nearest Neighbours (FAN) method with 5% of neighbours, yield the best results. Therefore, we would recommend these reference values for future work when applying RQA to the analysis of HRV in sleep apnea. We also conclude that, together with the commonly used vertical and diagonal RQA measures, there are newly used features that contribute valuable information for apnea minutes discrimination. Therefore, they are especially interesting for characterization purposes. Using two different databases supports that the conclusions reached are potentially generalizable, and are not limited by database variability.

Grasslands in the Iberian Peninsula are some of the most valuable ecosystems in Europe and are vulnerable because of their location in arid-semiarid regions. Remote sensing techniques have the potential for monitoring grasslands using vegetation indices (VIs), which can reveal bare soil and non-photosynthetic vegetation reflectance in these regions. The temporal variability of the VI time-series is commonly measured as the standard deviation of the records, insufficient to study the system dynamics. Recurrence plots (RP) and recurrence quantification analysis (RQA) allow us to visualize and quantify system dynamics based on topology. These advanced analyses calculate stochasticity through determinism (DET) and predictability degree based on the average length of diagonal structures (LT). This study aims to evaluate RPs, Cross Recurrence Plots (CRP), and RQA to visualize and quantify VIs and climatic series and their anomalies responses dynamics.The approach was tested at two sites with different agro-climatic characteristics in Central Spain: Guadalix de la Sierra (ZGU) and Soto del Real (ZSO). The Normalized Difference Vegetation Index (NDVI) and the Modified Soil Adjusted Vegetation Index (MSAVI), which includes a soil factor, were calculated from MODIS images. Temperature (TEMP) and precipitation (PCP) series were acquired from two meteorological stations close to each site. All the anomalies series were obtained seasonally adjusted to the original ones.RPs showed the different VIs responses in both areas. The MSAVI patterns showed a clear structure, while NDVI presented a noisy pattern showing a higher random behaviour. To quantify these visual differences with RQA, four parameters were calculated, including DET and LT. The MSAVI series and the corresponding anomaly series presented higher DET and LT than the NDVI in both sites. CRPs of VIs with TMIN shows a higher DET in ZGU, suggesting that time synchronization is higher than in ZSO. When CRPs were applied over anomalies series, we could observe that precipitation is more synchronized with VIs anomalies, most likely because of the influence of soil moisture. LT of CRPs in all the cases was low and close to 2, indicating low predictability. Overall, our results suggest that VI time series show a more evident pattern and differentiate both sites better when a variable soil adjustment factor is included. Therefore, in arid-semiarid grasslands, a soil factor should be considered in VI calculations.

Recurrence quantification analysis (RQA) is a widely used tool for studying complex dynamical systems, but its standard implementation requires computationally expensive calculations of recurrence plots (RPs) and line length histograms. This study introduces strategies to compute RQA measures directly from time series or phase space vectors, avoiding the need to construct RPs. The calculations can be further accelerated and optimised by applying a random sampling procedure, in which only a subset of line structures is evaluated. These modifications result in shorter run times, less memory use and access, and lower overall energy consumption during analysis whilst maintaining accuracy. This makes them especially appealing for large-scale data analysis and machine learning applications. The ideas are not limited to diagonal line measures, but can likewise be applied to vertical line-based measures and to recurrence network measures. By lowering computational costs, the proposed strategies contribute to energy saving and sustainable data analysis, and broaden the applicability of recurrence-based methods in modern research contexts.

An on‐line method is developed for monitoring of mean particle size in fluidized beds using pressure fluctuations (PFs) and acoustic emissions (AE) signal by recurrence plot (RP) and recurrence quantification analysis (RQA). PFs and AE signals of a lab‐scale fluidized bed were measured simultaneously at various superficial gas velocities and mean particle sizes. Although the AE signals are often very complicated due to many different acoustic sources in the bed, applying RP analyses showed that small changes in mean particle size can be detected by visual comparison of AE‐RP structures, while this cannot be distinguished by graphical RP analysis of PFs. Moreover, the hydrodynamics of the bed was inspected through RQA analysis of both signals. For this purpose, recurrence rate, determinism, laminarity, average length of diagonal and vertical lines were extracted from RPs showing the effect of an increase in the mean particle size. © 2016 American Institute of Chemical Engineers AIChE J, 62: 3547–3561, 2016

Several previous attempts have been made to divide the human amygdala into smaller subregions based on the unique functional properties of the subregions. Although these attempts have provided valuable insight into the functional heterogeneity in this structure, the possibility that spatial patterns of functional characteristics can quickly change over time has rarely been considered in previous studies. In the present study, we explicitly account for the dynamic nature of amygdala activity. Our goal was not only to develop another parcellation method but also to augment existing methods with novel information about amygdala subdivisions. We performed state-specific amygdala parcellation using resting-state fMRI (rsfMRI) data and recurrence quantification analysis (RQA). RsfMRI data from 102 subjects were acquired with a 3T Trio Siemens scanner. We analyzed values of several RQA measures across all voxels in the amygdala and found two amygdala subdivisions, the ventrolateral (VL) and dorsomedial (DM) subdivisions, that differ with respect to one of the RQA measures, Shannon's entropy of diagonal lines. Compared to the DM subdivision, the VL subdivision can be characterized by a higher value of entropy. The results suggest that VL activity is determined and influenced by more brain structures than is DM activity. To assess the biological validity of the obtained subdivisions, we compared them with histological atlases and currently available parcellations based on structural connectivity patterns (Anatomy Probability Maps) and cytoarchitectonic features (SPM Anatomy toolbox). Moreover, we examined their cortical and subcortical functional connectivity. The obtained results are similar to those previously reported on parcellation performed on the basis of structural connectivity patterns. Functional connectivity analysis revealed that the VL subdivision has strong connections to several cortical areas, whereas the DM subdivision is mainly connected to subcortical regions. This finding suggests that the VL subdivision corresponds to the basolateral subdivision of the amygdala (BLA), while the DM subdivision has some characteristics typical of the centromedial amygdala (CMA). The similarity in functional connectivity patterns between the VL subdivision and BLA, as well as between the DM subdivision and CMA, confirm the utility of our parcellation method. Overall, the study shows that parcellation based on BOLD signal dynamics is a powerful tool for identifying distinct functional systems within the amygdala. This tool might be useful for future research on functional brain organization.