Recurrence Quantification Analysis Based Methodology in Automatic Aerobic Threshold Detection: Applicability and Accuracy across Age Groups, Exercise Protocols and Health Conditions

Abstract

A new method based on the Recurrence Quantification Analysis (RQA) of the heart rate (HR) offers an objective, efficient alternative to traditional methods for Aerobic Threshold (AerT) identification that have practical limitations due to the complexity of equipment and interpretation. This study aims to validate the RQA-based method’s applicability across varied demographics, exercise protocols, and health status. Data from 123 cardiopulmonary exercise tests were analyzed, and participants were categorized into four groups: athletes, young athletes, obese individuals, and cardiac patients. Each participant’s AerT was assessed using both traditional ventilatory equivalent methods and the automatic RQA-based method. Ordinary Least Products (OLP) regression analysis revealed strong correlations (r > 0.77) between the RQA-based and traditional methods in both oxygen consumption (VO₂) and HR at the AerT. Mean percentage differences in HR were below 2.5%, and the Technical Error for HR at AerT was under 8%. The study validates the RQA-based method, directly applied to HR time series, as a reliable tool for the automatic detection of the AerT, demonstrating its accuracy across diverse age groups and fitness levels. These findings suggest a versatile, cost-effective, non-invasive, and objective tool for personalized exercise prescription and health risk stratification, thereby fulfilling the study’s goal of broadening the method’s applicability.

1. Introduction

Identifying the aerobic threshold (AerT) is vital for crafting individualized exercise programs and optimizing both health and athletic performance [1,2,3]. The AerT is defined as the exercise intensity at which there is a shift from predominantly aerobic to anaerobic energy production [2]. Traditionally, the AerT has been determined through visually detected and manually uncomfortable methods during physical activities, including physiological variables like gas exchange, heart rate, blood lactate concentration, and rate of perceived exertion, respectively [4,5]. However, these approaches are often time-consuming, and their result is subject to operator-dependent variability during interpretation.

To address these limitations, we recently introduced an objective, efficient, nonlinear method based on Recurrence Quantification Analysis (RQA) that employs heart rate variability (HRV) for automatic AerT detection [6,7,8]. The autonomic nervous system, which controls the body’s physiological reactions to exercise, is reflected in HRV. Consequently, a more thorough and dynamic examination of the physiological reactions to different exercise intensities is provided by utilizing HRV to measure AerT. This connection is particularly significant in clinical settings, where monitoring and management of training, particularly in patients with heart problems or on medications that affect HRV, depend heavily on a thorough understanding of physiological responses to workload. Although HRV-based techniques might not always be more effective, they do have some clear benefits in terms of simplicity of usage and reduced bulky. We hope to simplify the evaluation procedure while preserving the association with recognized physiological measures by employing HRV. This approach allows for a more practical and user-friendly method of AerT evaluation, which could potentially enhance its applicability in various settings. Our previous study tested the HRV-based algorithm on a homogenous sample in terms of age (15.26 ± 1.87 years) and health status (healthy athletes) undergoing the cardiovascular test while performing a single type of activity (cycling) [6].

Indeed, the method’s initial focus was primarily on optimizing athletic performance, which underscores the significance of accurately identifying the AerT not only in sports but also in managing health conditions. Our first study demonstrated a strong correlation (r > 0.77) between the RQA-based method and traditional techniques when applied to a homogenous group of athletes, with a technical error of less than 5% for HR detection at AerT [6]. The precise detection of AerT is critical for assessing the risk of cardiovascular and metabolic disorders in obese individuals [9,10]. For cardiac patients, particularly those on beta-blockers, accurately determining AerT is vital because it is significantly influenced by medications that affect HRV, a crucial indicator of autonomic nervous system function [11,12,13]. Given the foundational importance in both athletic and health contexts, this study aims to expand the utility of our automatic AerT detection method. We are broadening its application to include various physical activities and a more diverse demographic encompassing different ages, genders, fitness levels, and, most importantly, specific health conditions such as obesity and cardiovascular issues treated with beta-blockers. Specifically, we focus on four case studies: athletes, young athletes ages 10–17, obese individuals, and cardiac patients. These groups were selected because they are paradigmatic for evaluating metabolic thresholds.

One of the main advantages of this algorithm is its ability to be integrated into wearable systems for real-time detection, which is becoming increasingly feasible [14]. Flexible bioelectronics, like wearable and epidermal electronics that can be attached to the skin, clothing, or even the body itself, would be ideal for an objective evaluation of physical activity [15]. Furthermore, the rapidly evolving framework of future smart gadgets for tracking physical activity may be built based on small multifunctional sensors [16]. With the help of these new technologies, this automatic algorithm could be integrated into a wearable device, advancing real-time detection. These improvements would make fitness levels easily adjustable and accessible, being beneficial for diseases like obesity [7] and heart problems [8].

2. Materials and Methods

2.1. Participants

This work expands and extends the one presented in [6] by encompassing groups with a broader range of ages, varied training states, and different health conditions and employing a graduated test protocol. We then analyzed 123 cardiopulmonary tests (CPET) from four distinct groups: 27 athletes (Group A [6]), 20 obese individuals (Group O, [7]), 19 cardiac patients (Group C, [8]), and 57 young athletes (Group Y, [17]). The participants (25.35 ± 16.62 years) were between 11 and 70 years old and were divided based on age, fitness level, and health conditions. The athletes and the young athletes’ participants were free from serious health conditions like neuropathy, autonomic dysfunction, and cardiovascular disorders, while the obese individuals and cardiac patients met specific health criteria. The cardiac patients underwent testing before and after rehabilitation sessions. The study was conducted in accordance with the Declaration of Helsinki and received approval from various ethical committees: the CAR-IRB—University of Rome ‘Foro Italico’ Committee (Approval No. CAR 37/2020), the Sapienza University of Rome Ethical Committee (Approval No. 70/11, 2011) for groups A and O, the Ethical Committee of Tours (France) for group C, and the Institutional Review Board at Appalachian State University (IRB no. 18-0147) for group Y.

Participants provided informed consent, and clinical examinations were conducted to ensure safety and eligibility for physical activity. The combined analysis offers a comprehensive view of the physiological responses across a wide age range and varying health and fitness levels, employing time series data collected during incremental exercise tests. The methodologies and protocols were carefully designed to ensure data reliability and ethical compliance, providing valuable insights into cardiopulmonary responses under different physical conditions.

Groups underwent different cardiopulmonary tests (CPET), which are described in Appendix A. During these tests, time-series data related to gas exchange (VO₂, VCO₂, oxygen saturation SpO₂, etc.) and heart rate (HR) were collected. In this study, we have focused exclusively on analyzing the HR and VO₂ time series, comparing them against the gold standard for validation purposes.

2.2. Exclusion Criteria

HR time series were excluded from the analysis under the following criteria:

(a)

Insufficient consecutive HR data points (fewer than 150, corresponding to a test duration of 5 min)

(b)

Excessive difference (>10 beats per minute) between consecutive HR data points

(c)

HR values lower than the average of the preceding 10 data points within the time series

(d)

The presence of outliers in the HR time series is defined as values exceeding 4 standard deviations from the mean of that series.

2.3. Preprocessing

The first and the last steps of the graded exercise were discarded. In the context of a graded exercise test, whether to discard the first and last steps of the exercise for the analysis of the aerobic threshold depends on various factors, including the specific protocol of the test, the purpose of the analysis, and the characteristics of the individual being tested. Graded exercise typically involves a progressive increase in exercise intensity, and it is used to assess cardiovascular fitness and determine aerobic and anaerobic thresholds. Here are some considerations regarding the first and last steps.

The initial stage of the graded exercise is often a warm-up phase where the individual is just beginning to exercise, and the body is transitioning from a state of rest to one of physical activity. Moreover, physiological responses (like HR and oxygen consumption) might not have stabilized yet. This could lead to less reliable data for determining the aerobic threshold.

The final stage of the graded exercise usually represents the point of maximal or near-maximal effort. Here, the individual is often in a state of anaerobic metabolism. For aerobic threshold analysis, this phase may not be as relevant because the focus is on the point where the body transitions from aerobic to anaerobic metabolism, which typically occurs at a lower intensity level. Data from the last step can be more relevant for determining the anaerobic threshold or maximal aerobic capacity, which is one of our study’s possible developments. In this study, which aimed to generalize the algorithm in [6] between heterogeneous groups and protocols, we decided to discard the first and last steps for every group and protocol (see Figure 1 in [6]).

2.4. Data Cleaning

To cleanse the time series data of Oxygen Consumption (VO₂) and Heart Rate from outliers resulting from malfunctions in the CPET (Cardiopulmonary Exercise Testing) equipment and the Garmin heart rate monitor while minimally altering the raw data, we employed a smoothing technique based on a moving average, which is particularly useful in smoothing out short-term fluctuations and highlighting longer-term trends or cycles (see in details Appendix B). An example of a moving average applied to an HR time series is reported in Figure 1.

2.5. Gas Exchange Method (Gex Method) for AerT Detection

The gold standard method for the identification of critical physiological thresholds, such as the AerT, is the Gas Exchange Method (Gex method), which requires a Cardiopulmonary Exercise Test (CPET), an advanced and comprehensive test that measures pulmonary and cardiovascular response during exercise [18,19]. During a CPET, individuals typically perform graded exercises on devices such as those manufactured by Cosmed^®. These devices are engineered to meticulously monitor and record various physiological parameters. The experimental setup diagram and a detailed description of the test system are reported in Figure 2.

The Gex method specifically examines the relationship between the ventilatory equivalent for oxygen (VE/VO₂) and oxygen consumption (VO₂). The ventilatory equivalent for oxygen is a dimensionless number that represents the ratio of the volume of air ventilated to the amount of oxygen consumed. In practice, the method involves graphing VE/ VO₂ against VO₂ values to identify the AerT. The Aerobic Threshold is discerned (by visual inspection) at the point where the VE/VO₂ ratio reaches its minimum value before it begins to rise again during a progressively strenuous exercise protocol.

The Gex procedure for determining AerT was validated using a second method known as V-slope [2]. This method involves plotting a graph with VCO₂ on the y-axis and VO₂ on the x-axis. Two regression lines are then fitted to represent the upper and lower segments of this relationship. The intersection of these lines corresponds to the Aerobic Threshold determined by the Gex procedure.

To accurately determine the HR and VO₂ at the Aerobic Threshold, a precise method is employed. The average values of HR and VO₂ are calculated over the final 30 s of the exercise stage just before the AerT is reached. This approach is taken because these last 30 s are seen as reflective of the individual’s physiological state at the threshold. Capturing the average values during this critical period is thought to provide a stable and reliable set of data points. These data points are crucial for subsequent statistical analysis, offering a solid foundation for comparing various conditions or interventions that might affect the Aerobic Threshold.

2.6. RQA-Based Method for AerT Automatic Detection

RQA is a powerful method used to analyze nonlinear dynamical systems. It involves reconstructing time series data into a phase space using delay embedding techniques, where vectors are created based on an embedding dimension m and time delay τ. A matrix R is then formed by comparing points in this phase space using a predefined threshold distance ϵ. The resulting matrix is graphically represented as a Recurrence Plot (RP), showing patterns of recurrence, such as diagonal lines indicating repeating behaviors. Key parameters like ϵ, τ, and m are chosen carefully to suit the characteristics of the data and study goals. Measures derived from the recurrence plot, such as Determinism (DET) and Laminarity (LAM), quantify the predictability and structure of the system based on the lengths of diagonal and vertical (or horizontal) lines in the plot, respectively. For a more detailed overview of this nonlinear analysis, see Appendix C.

RQA has been used to provide a straightforward framework for describing the coupling of cardio-respiratory physiology under increasing physical exertion in terms of both complexity and dynamical transitions. This RQA-based methodology has been validated using DET and LAM to identify ventilatory thresholds, as outlined in references [6,7].

In particular, the assessment of the AerT is based on identifying a minimum in DET (DET_min), which is corroborated by a corresponding minimum in LAM (LAM_min).

Automatic Detection Algorithm

The RQA-based method has its limitations. Indeed, these thresholds have been identified so far by visually inspecting where DET/LAM underwent a sudden change or reached saturation. The need now is to determine these thresholds without the requirement of external user intervention but rather in an automated manner. To achieve this result, we developed a novel algorithm written in Python (version 5.3.3) that can automatically detect the most convex minima of a time series. This algorithm has been already presented in [6]. For more details, see Appendix D.

Here we report the main steps:

• Data cleaning and preprocessing on HR and VO₂ time series.
• The HR time series is divided into partially overlapping epochs of width = 100 points, shifted by 1 point, and for every epoch, RQA is performed, and DET values are calculated.
• The new DET time series is then analyzed to identify the most convex minima. After calculating the first and second derivatives, a point is considered a minimum if its second derivative is greater than or equal to the mean of the second derivative plus two times its standard deviation:

$$f″\left(t\right)>mean\left(f″\left(t\right)\right)+2\times SD\left(f″\left(t\right)\right)$$

(1)

The value $mean\left(f″\left(t\right)\right)+2\times SD\left(f″\left(t\right)\right)$ is named cut-off or radius.

This condition selects points where the curvature is significantly convex. Finally, the function returns a list of indices (minima) corresponding to these minima points.

4.

The AerT is identified by locating the most convex minimum of DET using the formula:

$${\mathrm{DET}}_{\min }=\max \left({\displaystyle \frac{d^2\mathrm{DET}\left(t\right)}{d{t}^2}}\right)$$

(2)

This selection is among all relative minima of DET. The point at which AerT occurs is termed DETmin, identified at the time point tVT1.

2.7. Statistics

In this context, variables are named following a specific convention to denote the method used for their calculation. Variables appended with “_RQA” are those calculated using the RQA-based method. This suffix indicates that the analysis and values of these variables are derived from RQA techniques. On the other hand, variables with the suffix “_Gex” are calculated using the Gas Exchange (Gex) method. This naming convention helps in differentiating the source and methodological background of each variable, ensuring clarity and precision in our analysis and discussions.

The Kolmogorov–Smirnov Test of Normality [20] confirmed the normal distribution of VO₂ and HR values at the AerT. The outcomes were as follows: VO₂_RQA (D = 0.09, p = 0.26) and HR_RQA (D = 0.07, p = 0.64).

Repeated measures analysis of variance (RM ANOVA) was utilized to assess the main effects of methods (RQA vs. Gex) and the interaction effects.

Agreement between the RQA-based and Gex methods was assessed using an Ordinary Least Products (OLP) regression analysis [21].

The regression’s coefficients of determination (R²), regression parameters (slope and intercept), and 95% confidence intervals (CI) were calculated to identify fixed and proportional biases. The following criteria were adopted to interpret the magnitude of the correlation r between the estimates: <0.1 trivial, 0.1–0.3 small, 0.3–0.5 moderate, 0.5–0.7 large, 0.7–0.9 very large, and 0.9–1.0 almost perfect [20]. The hypothesis of proportional and fixed bias was rejected when the 95% CI for the slope included 1 and the intercept included 0.

Percentage differences between the RQA-based and Gex methods at AerT for HR and VO₂ were calculated using Typical Percentage Error (TE(%)) and Mean Percentage Difference (Meandiff (%)).

Intraclass Correlation Coefficients (ICC) were employed as measures of criterion validity for the RQA-based method compared to the Gex method at AerT [22].

The magnitude of differences between the RQA-based and Gex methods was assessed using Cohen’s d effect size statistics. Effect sizes were classified as follows: d < 0.2 was considered trivial, 0.2 < d < 0.6 as small, 0.6 < d < 1.2 as moderate, 1.2 < d < 2.0 as large, and d > 2.0 as very large. Additionally, the interpretation of ICC ranges was defined as excellent (ICC > 0.9), moderate (0.75 < ICC < 0.90), fair (0.50 < ICC < 0.75), and poor (ICC < 0.5) [23].

Bland–Altman plot analyses were employed to visualize the agreement between RQA-based and Gex methods [24].

Statistical significance was set at p < 0.05 in both studies, and all statistical analyses were conducted using SPSS version 24.0 software (SPSS Inc., Chicago, IL, USA).

3. Results

3.1. Variability of Heart Rate and Oxygen Consumption at Aerobic Threshold across Participant Groups

The final number of time series considered in this study was 111, since 12 time series did not meet the inclusion criteria introduced in the Materials and Methods section.

Table 1. Overview of participant profiles and study parameters.

	Group A	Group O	Group C	Group Y	All
Name	Athletes	Obese	Cardiac Patients	Young athletes	Athlete, obese
					Young athletes
					cardiac patients
n	27	19	16	49	111
Age(years)	15.30 ± 1.87	42.90 ± 12.67	52.80 ± 10.09	15.20 ± 2.02	25.35 ±16.62
BMI (kg/m2)	22.06 ± 3.28	39.98 ± 4.95	25.92 ± 4.29	22.65 ± 4.97	25.78 ± 7.60
Range (years)	14–20	19–63	35–70	11–17	11–70
Gender	25/2	6/13	16/0	26/23	73/38
(M/F)
Wstatus (NW/OW)	27/0	0/19	13/3	32/17	72 NW/39 OW
Hstatus	Healthy (H)	Disease (D)	Disease (D)	Healthy (H)	76 H/35 D
Incremental	20 or	1 km/h every	10 (Watts/min)	10 (Watts/min)	10, 15, 20
protocol	15(Watts/min)	3 min			(Watts/min)
					1 km/h /3 min
Type	Cycle-ergometer (C)	Treadmill (T)	Cycle-ergometer (C)	Cycle-ergometer (C)	90 C/19 T

Wstatus is the weight status (NW = normal weight, OW = overweight); Hstatus is health status.

serves as a comprehensive overview of participant profiles and study parameters, allowing for a detailed comparison and analysis across different demographic and health-related categories within the study population.

Figure 3 presents HR and VO₂ values at the AerT for the study participants. It showcases a wide spectrum of HR values, ranging from 60 to 180 beats per minute (bpm), and VO₂ values, spanning from 600 to 2500 mL/min at AerT. Notably, young athletes typically display lower VO₂ values, which can be attributed to their smaller body sizes. Conversely, cardiac patients often exhibit both lower VO₂ and HR values, reflecting their distinct physiological condition.

Figure 4 illustrates two representative HR time series and outlines the key steps of the automatic detection process. These time series were analyzed using the Python script with our automatic detection algorithm. The series on the right pertains to a young athlete from group Y, while the series on the left is from a cardiac patient in group C. The blue graded lines vividly depict the incremental test’s workload. As the intensity of physical exercise increases, so does the HR, as observed in panels (a) and (d). Panel (d) demonstrates how beta blockers effectively regulate the increase in HR. Through the RQA of HR, the determinism for each epoch and its second derivative were computed, as shown in panels (b) and (e). The colored points mark relative minima identified by our automatic algorithm. In panels (c) and (f), the most convex minimum, which aligns with the AerT value, is highlighted by a red vertical line.

3.2. Comparative Analysis of Traditional and RQA-Based Methods in Assessing Aerobic Threshold

The concordance between the traditional Gex method and the RQA-based automated technique for determining the AerT has been examined for both HR and VO₂. An analysis of variance with repeated measures (ANOVA RM) showed no significant differences between the two methods for HR and VO₂ at AerT (p > 0.05), as outlined in

Table 2. Main Parameters at AerT by two methods (mean ± SD) in athletes (A), cardiac patients (C), obese (O), and young athletes (Y).

Group	HR_RQA (bpm)	HR_Gex (bpm)	VO2_RQA (mL/Min)	VO2_Gex (mL/Min)
A (27)	136.81 ± 15.46	137.34 ± 15.26	1776.56 ± 445.95	1767.57 ± 335.89
O (19)	110.22 ± 14.06	112.37 ± 14.09	1379.28 ± 335.19	1359.12 ± 316.4
C (16)	87.31 ± 12.53	84.69 ± 13.01	1051.09 ± 462.49	902.79 ± 242.67
Y (49)	126.36 ± 14.13	121.16 ± 14.18	1134.54 ± 343.67	1056.57 ± 260.37
All (111)	120.5 ± 21.36	118.62 ± 21.44	1320.57 ± 472.98	1263.70 ± 427.71
p method *	0.50		0.10
p method X group	0.04		0.38
p method X protocol	0.03		0.33
p method X Wstatus *	0.97		0.62
p method X Hstatus *	0.87		0.86
p method X gender *	0.98		0.61

Two factors (Gex vs. RQA) repeated measures for HR and VO₂ values across all participants, groups (A, C, O, Y), protocol (Cycle-ergometer or treadmill), Wstatus (obese and normal weight), Hstatus (healthy and disease) and gender (male, female). * Gex method vs. RQA-based method all p > 0.05; post-hoc power = 0.998, n = 111 subjects.

.

Additionally, the ANOVA RM revealed no significant differences in HR and VO₂ across various groups when considering factors such as weight status (72 Normal or 39 Overweight), health status (76 Healthy from groups Y and A; 35 Diseased from groups C and O), and gender (73 Male or 38 Female).

The substantial difference in participant numbers between the treadmill (19 time series) and cycle-ergometer (92 time series) protocols could influence the statistical outcomes and significance observed in the comparison methodXprotocol for HR values (p = 0.03). A larger sample size in the cycle-ergometer group likely enhances statistical power, increasing the ability to detect significant differences between the automated detection methods. Moreover, larger samples generally yield more precise estimates and may better represent the target population. However, it is important to acknowledge potential biases and variations in participant characteristics across these protocols, which could impact the generalizability and interpretation of results.

Moreover, the comparison methodXgroup for HR values (p = 0.04) is influenced by the difference between Group Y (coupled data t-test, p < 0.05) and the other groups (coupled data t test p > 0.1 for A, C, and O group). Young athletes typically exhibit distinct physiological characteristics compared to adults, including higher heart rates at rest and during exercise due to factors such as increased cardiac output, lower resting heart rates, and greater cardiovascular efficiency. These differences can influence HR measurements and may contribute to variations observed in the comparison of HR metrics between groups.

Table 3. Agreement with the gold standard (RQA-based vs. Gex method) for HR (bpm) and VO₂ (mL/min), Pearson correlation (R²), slope and intercept values, mean difference (Mean_diff(%)), t-test (p) and Typical Percentage Error (TE(%)) of the whole sample.

Parameters	R2	Slope	Intercept	Meandiff (%)	p	TE (%)
HR	0.68	1.0 (0.57 to 1.75)	2.31 (−86.77 to 53.11)	2.23	0.121	7.71
VO2	0.60	1.11 (0.72 to 1.69)	−76.86 (−815.94 to 406.55)	6.51	0.053	17.54

reports various parameters and statistics related to the comparison between two methods (RQA and Gex) for measuring HR and VO₂ in the total sample of 111 time series. The R-squared (R2) values (0.68 for HR and 0.6 for VO₂) suggest a moderate to good level of agreement between the two measurement methods. The slope indicates how measurements from one method relate to measurements from the other method. For both HR and VO₂ the slope value is around 1, which indicates perfect agreement. The intercept provides insight into any systematic bias or offset between the two methods. The values for HR and VO₂ are within the accepted range, suggesting that there is no bias. Meandiff(%) and TE(%) reflect the magnitude and variability of differences between measurements obtained by the two methods.

The mean percentage difference between the two assessment methods was below 2.5% for HR and under 8% for VO₂, with these discrepancies not proving statistically significant (p > 0.05) at the Aerobic Threshold. TE was considered acceptable, with values below 8% for HR and around 18% for VO₂, indicating a high level of precision in the measurements [24,25].

The Intraclass Correlation Coefficients (ICC) demonstrated excellent reliability, with scores of 0.90 for HR and 0.87 for VO₂, indicating strong agreement between methods [22]. The effect size, calculated using Cohen’s d, was found to be trivial for both HR and VO₂, indicating minimal differences between the methods, as detailed in

Table 4. Agreement with Gex method for HR (bpm) and VO₂ (mL/min) value.

Parameters	r Pearson	ICC	Effect Size (d)
HR (bpm)	0.82 **	0.90 (0.86 to 0.93)	0.10 *
VO2 (mL/min)	0.77 **	0.87 (0.81 to 0.91)	0.14 *

Pearson correlation coefficient (r Pearson), Intraclass Correlation Coefficient (ICC), and Cohen’s Effect Size (d); ** all large; * d < 0.2 trivial; (0.2 < d < 0.6 moderate); ICC > 0.9 excellent, 0.75 < ICC < 0.90 moderate.

.

The Ordinary Least Products (OLP) regression analysis plots for HR and VO₂ at the AerT between both methods are visually depicted in Figure 5. Here, the OLP line of best fit is drawn as a solid line, while the line of unity is depicted as a dotted line, facilitating a clear comparison between the traditional and automated techniques.

Furthermore, Bland-Altman plots, illustrated in Figure 6, showed no systematic bias for either HR or VO₂. The plots indicated that the dispersion of differences was more contained for HR than VO₂, highlighting VO_2′s greater variability. Only a few outliers were present, but they did not significantly affect the overall concordance between the methods, thereby reinforcing the robustness of the findings.

4. Discussion

The most recent published study introduced an automated method for identifying the aerobic threshold using Heart Rate (HR) data gathered during the Cardiopulmonary Exercise Test (CPET), specifically in adult athletes, including rowers and amateurs [6]. This current research aimed to expand the application of this innovative method by validating its reliability across a broader sample. Specifically, we tested the method’s effectiveness across different demographic groups, including variations in age, gender, and fitness levels, as well as among individuals with diverse health conditions, such as obesity and cardiovascular diseases, and those undergoing treatment with beta-blockers. By examining these variations, we aimed to establish the method’s robustness and adaptability in detecting the AerT under a wide range of physiological conditions, ensuring its relevance for both general health and specialized athletic training. Our study extended the automated detection method to 111 time series from individuals between 11 and 70 years old and including various health conditions (healthy children, obese individuals, and cardiac patients). The tests were administered on both treadmills and cycle-ergometers, necessitating different incremental protocols to achieve a consistent total exercise time and thereby yielding comparable HR time series [6,7,8].

The post-hoc power of the F-test increased from 0.88 for a 27-subject sample [6] to 0.999 for our expanded sample of 111 subjects. In the present research, the HR values derived from the automated RQA-based method were not significantly different from those obtained using the gold standard ventilatory equivalents-based method (Gex), as evidenced by a p-value greater than 0.05. Furthermore, a strong correlation (r > 82%) and an acceptable Technical Error (TE) of 8% were observed. Importantly, despite the potential effects of beta-blockers in cardiac patients altering the HR variability during exercise, the RQA-based method proved to be reliable for this group. However, it’s worth noting that in individuals with heart conditions, the use of heart rate as an exercise threshold marker is compromised by medication effects. For a more accurate exercise prescription in these cases, reliance on the Rating of Perceived Exertion (RPE) is advised [4]. General physiological considerations could be confirmed through intergroup comparisons. For instance, regular exercise enhances cardiovascular efficiency, reflected in higher HR and maximum oxygen consumption (VO_2max) values among athletes and young athletes when compared to obese or cardiac patients [1,3,8]. Age-related factors like lung volume also exhibited expected trends.

The RQA-based method shows promising future applications in the detection of AerT, a fundamental and complex challenge of exercise physiology [26]. By analyzing the recurrence patterns in physiological signals, RQA could provide deeper insights into the complex interplay of physiological processes and the critical thresholds demarcating different metabolic states. This enables personalized threshold identification, capturing unique patterns specific to each individual test. For the young athletes’ group, it’s crucial to consider the physiological uniqueness of adolescence, including growth spurts and changes in aerobic capacity. Tailored reference values and assessment methods are needed for this population [27]. Additional applications include identifying nonlinear markers related to fatigue and performance degradation, aiding in training optimization and performance prediction. As wearable technologies advance, the integration of RQA-based methods for real-time metabolic monitoring is increasingly feasible [28]. Machine learning algorithms could further refine threshold detection, particularly when combined with wearable device data or during tele-exercise sessions. Such methods would make fitness levels readily accessible and customizable, which are also useful for conditions like obesity [7] and cardiac issues [8]. Advanced machine learning techniques can efficiently explore functional network variations based on extensive data [28]. Recently, Kaufmann et al. 2023 compared established lactate (LT) and ventilatory thresholds (VT) with HR-derived thresholds for exercise intensity prescription in endurance athletes [29]. Their systematic review concluded that HR-derived thresholds offer a promising alternative, for example, can be found in the results of Schaffarczyk et al. 2023; they have shown good agreement between the heart rate and oxygen uptake at first and second ventilatory threshold by a nonlinear analysis named detrended fluctuation analysis (DFA) in a cohort of randomly selected women [30].

Advantages and Limitations

The main advantages of the RQA method are simplicity, non-invasiveness, and automation, which reduces subjectivity. The method is less expensive and more practical compared to traditional methods that require complex equipment and human intervention for interpretation. The issue of time consumption is likewise addressed by the automatic method. After 30 to 40 min (including the preparation, exercise, and recovery phases) of the CPET test, the time needed to determine metabolic thresholds drops from 10s of minutes for an operator to about 3 s. This is the amount of time needed for the computational time of the Python program, which includes all stages (preprocessing, data cleaning, and threshold calculation). This proves that it is more than possible to include this automatic system for real-time detection in wearable technology.

Limitations of the study include the number of subjects: a larger and more diverse sample could provide a more comprehensive understanding of the methodology’s applicability and accuracy. Group C is limited to a specific cardiovascular pathology; further investigations are needed to verify the applicability of the methodology to individuals with different health challenges. The within-subject variability could be studied if multiple tests on the same subject were conducted. The sample consisted of individuals from a predominantly European background. We acknowledge the importance of ethnic diversity and plan to conduct further studies on more heterogeneous populations to investigate any potential differences and further validate our method.

5. Conclusions

In summary, nonlinear heart rate variability analysis (RQA) during graded exercises holds tremendous promise for future research and practical applications in the identification of metabolic thresholds based on HR. This method can facilitate personalized training plans, performance monitoring, and real-time adjustments through wearable technology integration. Further research is likely to deepen our understanding of exercise physiology, thereby advancing the fields of sports science, exercise physiology, and healthcare.