Bayesian Estimation of Correlation between Measures of Blood Pressure Indices, Aerobic Capacity and Resting Heart Rate Variability Using Markov Chain Monte Carlo Simulation and 95% High Density Interval in Female School Teachers

Background: Several explanations regarding the disparity observed in the literature with regard to heart rate variability (HRV) and its association with performance parameters have been proposed: the time of day when the recording was conducted, the condition (i.e., rest, active, post activity) and the mathematical and physiological relationships that could have influenced the results. A notable observation about early studies is that they all followed the frequentist approach to data analyses. Therefore, in an attempt to explain the disparity observed in the literature, the primary purpose of this study was to estimate the association between measures of HRV indices, aerobic performance parameters and blood pressure indices using the Bayesian estimation of correlation on simulated data using Markov Chain Monte Carlo (MCMC) and the equal probability of the 95% high density interval (95% HDI). Methods: The within-subjects with a one-group pretest experimental design was chosen to investigate the relationship between baseline measures of HRV (rest; independent variable), myocardial work (rate–pressure product (RPP)), mean arterial pressure (MAP) and aerobic performance parameters. The study participants were eight local female schoolteachers aged 54.1 ± 6.5 years (mean ± SD), with a body mass of 70.6 ± 11.5 kg and a height of 164.5 ± 6.5 cm. Their HRV data were analyzed in R package, and the Bayesian estimation of correlation was calculated employing the Bayesian hierarchical model that uses MCMC simulation integrated in the JAGS package. Results: The Bayesian estimation of correlation using MCMC simulation reproduced and supported the findings reported regarding norms and the within-HRV-indices associations. The results of the Bayesian estimation showed a possible association (regardless of the strength) between pNN50% and MAP (rho = 0.671; 95% HDI = 0.928–0.004), MeanRR (ms) and RPP (rho = −0.68; 95% HDI = −0.064–−0.935), SDNN (ms) and RPP (rho = 0.672; 95% HDI = 0.918–0.001), LF (ms(2)) and RPP (rho = 0.733; 95% HDI = 0.935–0.118) and SD2 and RPP (rho = 0.692; 95% HDI = 0.939–0.055). Conclusions: The Bayesian estimation of correlation with 95% HDI on MCMC simulated data is a new technique for data analysis in sport science and seems to provide a more robust approach to allocating credibility through a meaningful mathematical model. However, the 95% HDI found in this study, accompanied by the theoretical explanations regarding the dynamics between the parasympathetic nervous system and the sympathetic nervous system in relation to different recording conditions (supine, reactivation, rest), recording systems, time of day (morning, evening, sleep etc.) and age of participants, suggests that the association between measures of HRV indices and aerobic performance parameters has yet to be explicated.