Changes in the Hurst exponent of heartbeat intervals during physical activity

A fractional-order SIRD model with time-dependent memory indexes for encompassing the multi-fractional characteristics of the COVID-19

COVID-19 is a novel coronavirus affecting all the world since December last year. Up to date, the spread of the outbreak continues to complicate our lives, and therefore, several research efforts from many scientific areas are proposed. Among them, mathematical models are an excellent way to understand and predict the epidemic outbreaks evolution to some extent. Due to the COVID-19 may be modeled as a non-Markovian process that follows power-law scaling features, we present a fractional-order SIRD (Susceptible-Infected-Recovered-Dead) model based on the Caputo derivative for incorporating the memory effects (long and short) in the outbreak progress. Additionally, we analyze the experimental time series of 23 countries using fractal formalism. Like previous works, we identify that the COVID-19 evolution shows various power-law exponents (no just a single one) and share some universality among geographical regions. Hence, we incorporate numerous memory indexes in the proposed model, i.e., distinct fractional-orders defined by a time-dependent function that permits us to set specific memory contributions during the evolution. This allows controlling the memory effects of more early states, e.g., before and after a quarantine decree, which could be less relevant than the contribution of more recent ones on the current state of the SIRD system. We also prove our model with Italy's real data from the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University.

Entropy Analysis of RR-Time Series From Stress Tests

The RR-interval time series or tachograms obtained from electrocardiograms have been widely studied since they reflect the cardiac variability, and this is an indicative of the health status of a person. The tachogram can be seen as a highly non-linear and complex time series, and therefore, should be analyzed with non-linear techniques. In this work, several entropy measures, Sample Entropy (SampEn), Approximate Entropy (ApEn), and Fuzzy Entropy (FuzzyEn) are used as a measure of heart rate variability (HRV). Tachograms belonging to thirty-nine subjects were obtained from a cardiac stress test consisting of a rest period followed by a period of moderate physical activity. Subjects are grouped according to their physical activity using the IPAQ sedentary and active questionnaire, we work with youth and middle-aged adults. The entropy measures for each group show that for the sedentary subjects the values are high at rest and decrease appreciably with moderate physical activity, This happens for both young and middle-aged adults. These results are highly reproducible. In the case of the subjects that exercise regularly, an increase in entropy is observed or they tend to retain the entropy value that they had at rest. It seems that there is a possible correlation between the physical condition of a person with the increase or decrease in entropy during moderate physical activity with respect to the entropy at rest. It was also observed that entropy during longer physical activity tests tends to decrease as fatigue accumulates, but this decrease is small compared to the change that occurs when going from rest to physical activity.

Estimation of statistical properties of rough surface profiles from the Hurst exponent of speckle patterns

We applied the Hurst exponent technique to an experimental study of rough metallic surface profiles and the speckle patterns generated by them. Characterization of important statistical properties of the surface profile and speckle patterns were performed. We observed a clear correlation between the Hurst exponent of a surface profile and the one calculated from the associated speckle patterns. Therefore, in principle, information of the Hurst exponent of the profile can be obtained from the Hurst exponent of speckle patterns. Range and sampling analyses were performed in the Hurst exponent calculations showing the robustness of the method. As an additional application, we performed a basic simulation to show that the Hurst exponent is sensitive to surface waviness.

Nonlinear Methods Most Applied to Heart-Rate Time Series: A Review

The heart-rate dynamics are one of the most analyzed physiological interactions. Many mathematical methods were proposed to evaluate heart-rate variability. These methods have been successfully applied in research to expand knowledge concerning the cardiovascular dynamics in healthy as well as in pathological conditions. Notwithstanding, they are still far from clinical practice. In this paper, we aim to review the nonlinear methods most used to assess heart-rate dynamics. We focused on methods based on concepts of chaos, fractality, and complexity: Poincaré plot, recurrence plot analysis, fractal dimension (and the correlation dimension), detrended fluctuation analysis, Hurst exponent, Lyapunov exponent entropies (Shannon, conditional, approximate, sample entropy, and multiscale entropy), and symbolic dynamics. We present the description of the methods along with their most notable applications.

Local field potential dynamics in the primate cortex in relation to parkinsonism reveled by machine learning: A comparison between the primary motor cortex and the supplementary area

The present study compares the cortical local field potentials (LFPs) in the primary motor cortex (M1) and the supplementary motor area (SMA) of non-human primates rendered Parkinsonian with administration of dopaminergic neurotoxin, 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine. The dynamic of the LFPs was investigated under several mathematical frameworks and machine learning was used to discriminate the recordings based on these features between healthy, parkinsonian with off-medication and parkinsonian with on-medication states. The importance of each feature in the discrimination process was further investigated. The dynamic of the LFPs in M1 and SMA was affected regarding its variability (time domain analysis), oscillatory activities (frequency domain analysis) and complex patterns (non-linear domain analysis). Machine learning algorithms achieved accuracy near 0.90 for comparisons between conditions. The TreeBagger algorithm provided best accuracy. The relative importance of these features differed with the cortical location, condition and treatment. Overall, the most important features included beta oscillation, fractal dimension, gamma oscillation, entropy and asymmetry of amplitude fluctuation. The importance of features in discriminating between normal and pathological states, and on- or off-medication states depends on the pair-comparison and it is region-specific. These findings are discussed regarding the refinement of current models for movement disorders and the development of on-demand therapies.

Non-linear dynamics of heart rate variability during incremental cycling exercise

Within the last years complex models of cardiovascular regulation and exercise fatigue have implemented heart rate variability (HRV) as a measure of autonomic nervous system. Using detrended fluctuation analysis (DFA) to assess heart rate correlation properties, the present study examines the influence of exercise intensity on total variability and complexity in non-linear dynamics of HRV. Sixteen cyclists performed a graded exercise test on a bicycle ergometer. HRV time domain measures and fractal correlation properties were analyzed using short-term scaling exponent alpha1 of DFA. Amplitude and complexity of HRV parameters decreased significantly. DFA-alpha1 increased from rest to low exercise intensity and showed an almost linear decrease from higher intensities until exhaustion. These findings support a qualitative change in self-organized heart rate regulation from a complex autonomic control at rest and low intensities towards a breakdown of the interaction in control mechanisms with non-autonomic heart rate control dominating at high intensities.

Hidden Signals-The History and Methods of Heart Rate Variability

The understanding of heart rate variability (HRV) has increased parallel with the development of modern physiology. Discovered probably first in 1847 by Ludwig, clinical applications evolved in the second part of the twentieth century. Today HRV is mostly used in cardiology and research settings. In general, HRV can be measured over shorter (e.g., 5-10 min) or longer (12 or 24 h) periods. Since 1996, most measurements and calculations are made according to the standard of the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. As the first step, the series of times between successive R-peaks in the ECG are in milliseconds. It is crucial, however, to identify and remove extrasystoles and artifacts according to standard protocols. The series of QRS distances between successive heartbeats can be analyzed with simple or more sophisticated algorithms, beginning with standard deviation (SDNN) or by the square root of the mean of the sum of squares of differences between adjacent normal RR (rMSSD). Short-term HRV is frequently analyzed with the help of a non-parametric fast Fourier transformation quantifying the different frequency bands during the measurement period. In the last decades, various non-linear algorithms have been presented, such as different entropy and fractal measures or wavelet analysis. Although most of them have a strong theoretical foundation, their clinical relevance is still debated.

Aging effects on cardiac and respiratory dynamics in healthy subjects across sleep stages

STUDY OBJECTIVES

Respiratory and heart rate variability exhibit fractal scaling behavior on certain time scales. We studied the short-term and long-term correlation properties of heartbeat and breathing-interval data from disease-free subjects focusing on the age-dependent fractal organization. We also studied differences across sleep stages and night-time wake and investigated quasi-periodic variations associated with cardiac risk.

DESIGN

Full-night polysomnograms were recorded during 2 nights, including electrocardiogram and oronasal airflow.

SETTING

Data were collected in 7 laboratories in 5 European countries.

PARTICIPANTS

180 subjects without health complaints (85 males, 95 females) aged from 20 to 89 years.

INTERVENTIONS

None.

MEASUREMENTS AND RESULTS

Short-term correlations in heartbeat intervals measured by the detrended fluctuation analysis (DFA) exponent alpha1 show characteristic age dependence with a maximum around 50-60 years disregarding the dependence on sleep and wake states. Long-term correlations measured by alpha2 differ in NREM sleep when compared with REM sleep and wake, besides weak age dependence. Results for respiratory intervals are similar to those for alpha2 of heartbeat intervals. Deceleration capacity (DC) decreases with age; it is lower during REM and deep sleep (compared with light sleep and wake).

CONCLUSION

The age dependence of alpha1 should be considered when using this value for diagnostic purposes in post-infarction patients. Pronounced long-term correlations (larger alpha2) for heartbeat and respiration during REM sleep and wake indicate an enhanced control of higher brain regions, which is absent during NREM sleep. Reduced DC possibly indicates an increased cardiovascular risk with aging and during REM and deep sleep.

Stratification pattern of static and scale-invariant dynamic measures of heartbeat fluctuations across sleep stages in young and elderly

Cardiac dynamics exhibit complex variability characterized by scale-invariant and nonlinear temporal organization related to the mechanism of neuroautonomic control, which changes with physiologic states and pathologic conditions. Changes in sleep regulation during sleep stages are also related to fluctuations in autonomic nervous activity. However, the interaction between sleep regulation and cardiac autonomic control remains not well understood. Even less is known how this interaction changes with age, as aspects of both cardiac dynamics and sleep regulation differ in healthy elderly compared to young subjects. We hypothesize that because of the neuroautonomic responsiveness in young subjects, fractal and nonlinear features of cardiac dynamics exhibit a pronounced stratification pattern across sleep stages, while in elderly these features will remain unchanged due to age-related loss of cardiac variability and decline of neuroautonomic responsiveness. We analyze the variability and the temporal fractal organization of heartbeat fluctuations across sleep stages in both young and elderly. We find that independent linear and nonlinear measures of cardiac control consistently exhibit the same ordering in their values across sleep stages, forming a robust stratification pattern. Despite changes in sleep architecture and reduced heart rate variability in elderly subjects, this stratification surprisingly does not break down with advanced age. Moreover, the difference between sleep stages for some linear, fractal, and nonlinear measures exceeds the difference between young and elderly, suggesting that the effect of sleep regulation on cardiac dynamics is significantly stronger than the effect of healthy aging. Quantifying changes in this stratification pattern may provide insights into how alterations in sleep regulation contribute to increased cardiac risk.

Endogenous circadian rhythm in human motor activity uncoupled from circadian influences on cardiac dynamics

The endogenous circadian pacemaker influences key physiologic functions, such as body temperature and heart rate, and is normally synchronized with the sleep/wake cycle. Epidemiological studies demonstrate a 24-h pattern in adverse cardiovascular events with a peak at approximately 10 a.m. It is unknown whether this pattern in cardiac risk is caused by a day/night pattern of behaviors, including activity level and/or influences from the internal circadian pacemaker. We recently found that a scaling index of cardiac vulnerability has an endogenous circadian peak at the circadian phase corresponding to approximately 10 a.m., which conceivably could contribute to the morning peak in cardiac risk. Here, we test whether this endogenous circadian influence on cardiac dynamics is caused by circadian-mediated changes in motor activity or whether activity and heart rate dynamics are decoupled across the circadian cycle. We analyze high-frequency recordings of motion from young healthy subjects during two complementary protocols that decouple the sleep/wake cycle from the circadian cycle while controlling scheduled behaviors. We find that static activity properties (mean and standard deviation) exhibit significant circadian rhythms with a peak at the circadian phase corresponding to 5-9 p.m. ( approximately 9 h later than the peak in the scale-invariant index of heartbeat fluctuations). In contrast, dynamic characteristics of the temporal scale-invariant organization of activity fluctuations (long-range correlations) do not exhibit a circadian rhythm. These findings suggest that endogenous circadian-mediated activity variations are not responsible for the endogenous circadian rhythm in the scale-invariant structure of heartbeat fluctuations and likely do not contribute to the increase in cardiac risk at approximately 10 a.m.

Power Law Analysis Estimates of Analyte Concentration and Particle Size in Highly Scattering Granular Samples from Photon Time-of-Flight Measurements

Optical measurements of particle size and composition in granular samples are difficult to make due to complex light scattering from particles. These multiple scattering events bias absorption estimates and complicate the calculation of scattering and absorption coefficients used to estimate sample properties. Time series data, such as chromatograms and photon time-of-flight (TOF) profiles, contain self-repeating (fractal) characteristics. Power law analysis of photon TOF profiles allows the determination of absorption coefficients and particle sizes in a single experiment. A correlation dimension algorithm was used on photon TOF data from scattering samples. MLR models were then obtained from correlation dimension plots for the estimation of sample properties. Estimates of particle sizes and absorption coefficients were shown to agree well with theoretical values when compared using independent validation sets. Results show close to a 3-fold and up to a 5-fold decrease in the errors of estimation of dye concentration and particle size, respectively, as compared to steady-state measurements. The power law approach provides a useful means of determining sample properties in highly scattering media.

Fractal scale-invariant and nonlinear properties of cardiac dynamics remain stable with advanced age: a new mechanistic picture of cardiac control in healthy elderly

Heart beat fluctuations exhibit temporal structure with robust long-range correlations, fractal and nonlinear features, which have been found to break down with pathologic conditions, reflecting changes in the mechanism of neuroautonomic control. It has been hypothesized that these features change and even break down also with advanced age, suggesting fundamental alterations in cardiac control with aging. Here we test this hypothesis. We analyze heart beat interval recordings from the following two independent databases: 1) 19 healthy young (average age 25.7 yr) and 16 healthy elderly subjects (average age 73.8 yr) during 2 h under resting conditions from the Fantasia database; and 2) 29 healthy elderly subjects (average age 75.9 yr) during approximately 8 h of sleep from the sleep heart health study (SHHS) database, and the same subjects recorded 5 yr later. We quantify: 1) the average heart rate (<R-R>); 2) the SD sigma(R-R) and sigma(DeltaR-R) of the heart beat intervals R-R and their increments DeltaR-R; 3) the long-range correlations in R-R as measured by the scaling exponent alpha(R-R) using the Detrended Fluctuation Analysis; 4) fractal linear and nonlinear properties as represented by the scaling exponents alpha(sgn) and alpha(mag) for the time series of the sign and magnitude of DeltaR-R; and 5) the nonlinear fractal dimension D(k) of R-R using the fractal dimension analysis. We find: 1) No significant difference in (P > 0.05); 2) a significant difference in sigma(R-R) and sigma(DeltaR-R) for the Fantasia groups (P < 10(-4)) but no significant change with age between the elderly SHHS groups (P > 0.5); and 3) no significant change in the fractal measures alpha(R-R) (P > 0.15), alpha(sgn) (P > 0.2), alpha(mag) (P > 0.3), and D(k) with age. Our findings do not support the hypothesis that fractal linear and nonlinear characteristics of heart beat dynamics break down with advanced age in healthy subjects. Although our results indeed show a reduced SD of heart beat fluctuations with advanced age, the inherent temporal fractal and nonlinear organization of these fluctuations remains stable. This indicates that the coupled cascade of nonlinear feedback loops, which are believed to underlie cardiac neuroautonomic regulation, remains intact with advanced age.

The "Chaos Theory" and nonlinear dynamics in heart rate variability analysis: does it work in short-time series in patients with coronary heart disease?

BACKGROUND

Dynamic analysis techniques may quantify abnormalities in heart rate variability (HRV) based on nonlinear and fractal analysis (chaos theory). The article emphasizes clinical and prognostic significance of dynamic changes in short-time series applied on patients with coronary heart disease (CHD) during the exercise electrocardiograph (ECG) test.

METHODS

The subjects were included in the series after complete cardiovascular diagnostic data. Series of R-R and ST-T intervals were obtained from exercise ECG data after sampling digitally. The range rescaled analysis method determined the fractal dimension of the intervals. To quantify fractal long-range correlation's properties of heart rate variability, the detrended fluctuation analysis technique was used. Approximate entropy (ApEn) was applied to quantify the regularity and complexity of time series, as well as unpredictability of fluctuations in time series.

RESULTS

It was found that the short-term fractal scaling exponent (alpha(1)) is significantly lower in patients with CHD (0.93 +/- 0.07 vs 1.09 +/- 0.04; P < 0.001). The patients with CHD had higher fractal dimension in each exercise test program separately, as well as in exercise program at all. ApEn was significant lower in CHD group in both RR and ST-T ECG intervals (P < 0.001).

CONCLUSIONS

The nonlinear dynamic methods could have clinical and prognostic applicability also in short-time ECG series. Dynamic analysis based on chaos theory during the exercise ECG test point out the multifractal time series in CHD patients who loss normal fractal characteristics and regularity in HRV. Nonlinear analysis technique may complement traditional ECG analysis.

Nonlinear analysis of discharge patterns in monkey basal ganglia

Spontaneous discharge of basal ganglia neurons is often analyzed with time- or frequency-domain methods. However, it has been shown that sequences of inter-spike interval series are not fully described by such linear procedures. We therefore carried out a characterization of the nonlinear features of spontaneous discharge of neurons in the primate basal ganglia. We studied the spontaneous activity of neurons in the subthalamic nucleus (22 cells), as well as neurons in the external and internal pallidal segments (53 and 39 cells, respectively), recorded with standard extracellular recording methods in two awake Rhesus monkeys. As a measure of the statistical irregularity of neuronal discharge, we compared the approximate entropy of inter-spike interval sequences with that of shuffled representations of the same data. In all three basal ganglia structures, approximately 95% of the original data showed lower approximate entropy values than the shuffled data, suggesting a temporal organization in the original sequence. Fano factor analysis confirmed the presence of a temporal organization of inter-spike interval sequences, and indicated the presence of self-similarity in the great majority of them. In addition, Hurst exponent analysis showed that the inter-spike interval series are persistent. Hurst exponents often differ between short and long scaling ranges. Subsequent principal component analyses allowed us to identify three distinct patterns of the temporal evolution of inter-spike interval sequences in the phase space. These types were found in varying distributions in all three nuclei. Our analyses demonstrate that the discharge of most neurons in the basal ganglia of awake monkeys has nonlinear features that may be important for information coding in the basal ganglia.

Cross-correlation of instantaneous phase increments in pressure-flow fluctuations: applications to cerebral autoregulation

We investigate the relationship between the blood flow velocities (BFV) in the middle cerebral arteries and beat-to-beat blood pressure (BP) recorded from a finger in healthy and post-stroke subjects during the quasisteady state after perturbation for four different physiologic conditions: supine rest, head-up tilt, hyperventilation, and CO2 rebreathing in upright position. To evaluate whether instantaneous BP changes in the steady state are coupled with instantaneous changes in the BFV, we compare dynamical patterns in the instantaneous phases of these signals, obtained from the Hilbert transform, as a function of time. We find that in post-stroke subjects the instantaneous phase increments of BP and BFV exhibit well-pronounced patterns that remain stable in time for all four physiologic conditions, while in healthy subjects these patterns are different, less pronounced, and more variable. We propose an approach based on the cross-correlation of the instantaneous phase increments to quantify the coupling between BP and BFV signals. We find that the maximum correlation strength is different for the two groups and for the different conditions. For healthy subjects the amplitude of the cross-correlation between the instantaneous phase increments of BP and BFV is small and attenuates within 3-5 heartbeats. In contrast, for post-stroke subjects, this amplitude is significantly larger and cross-correlations persist up to 20 heartbeats. Further, we show that the instantaneous phase increments of BP and BFV are cross-correlated even within a single heartbeat cycle. We compare the results of our approach with three complementary methods: direct BP-BFV cross-correlation, transfer function analysis, and phase synchronization analysis. Our findings provide insight into the mechanism of cerebral vascular control in healthy subjects, suggesting that this control mechanism may involve rapid adjustments (within a heartbeat) of the cerebral vessels, so that BFV remains steady in response to changes in peripheral BP.