Next Article in Journal
Macroinvertebrate Assemblages along the Longitudinal Gradient of an Urban Palmiet River in Durban, South Africa
Next Article in Special Issue
Associations between Objectively Determined Physical Activity and Cardiometabolic Health in Adult Women: A Systematic Review and Meta-Analysis
Previous Article in Journal
Environmental Factors Affecting Feather Taphonomy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 11
Issue 5
10.3390/biology11050704
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

What Mathematical Models Are Accurate for Prescribing Aerobic Exercise in Women with Fibromyalgia?

by
Santos Villafaina
1,2,*,
Clarissa Biehl-Printes
3,
José A. Parraca
2,4,
Fabiane de Oliveira Brauner
3 and
Pablo Tomas-Carus
2,4
1
Facultad de Ciencias del Deporte, Universidad de Extremadura, Avenida de la Universidad s/n, 10003 Caceres, Spain
2
Departamento de Desporto e Saúde, Escola de Saúde e Desenvolvimento Humano, Universidade de Évora, 7004-516 Evora, Portugal
3
Instituto de Geriatria e Gerontologia, Pontifícia Universidade Católica do Rio Grande do Sul, Porto Alegre 6681, Brazil
4
Comprehensive Health Research Centre (CHRC), University of Évora, 7005 Evora, Portugal
*
Author to whom correspondence should be addressed.
Biology 2022, 11(5), 704; https://doi.org/10.3390/biology11050704
Submission received: 8 April 2022 / Revised: 21 April 2022 / Accepted: 28 April 2022 / Published: 4 May 2022
(This article belongs to the Special Issue The Biological Underlying Processes of Physical Exercise Benefits)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Intensity prescription for cardiorespiratory exercises is crucial for achieving health/fitness benefits. However, not all of the population can access a cardiopulmonary exercise test, either for economic reasons or location resources, to determine their ventilatory thresholds. Therefore, different mathematical models can predict the intensity based on the maximum or reserve heart rate. Exercise prescription guidelines indicate that people with fibromyalgia should exercise at 60% of their VO2max. However, people with fibromyalgia suffer from dysautonomia, which could lead to chronotropic incompetence, the inability to increase heart rate with increasing exercise intensities. Therefore, this study aimed to investigate the relationship and level of agreement between different mathematical models and the heart rate obtained from a cardiopulmonary exercise test at their ventilatory threshold 1. The results showed that the well-known “220 − age” at 76% and the mathematical model designed for people with fibromyalgia “209 − 0.85 × age” at 76% showed a significant level of agreement. However, Tanaka and Karvonen’s formula did not show a significant level of agreement. Thus, the “220 − age” at 76% and “209 − 0.85 × age” at 76% can be used in people with FM to prescribe aerobic exercise.

Abstract

Objectives: This article aims to verify the agreement between the standard method to determine the heart rate achieved in the ventilatory threshold 1 in the cardiopulmonary exercise testing (VT1) and the mathematical models with exercise intensities suggested by the literature in order to check the most precise for fibromyalgia (FM) patients. Methods: Seventeen women with FM were included in this study. The VT1 was used as the standard method to compare four mathematical models applied in the literature to calculate the exercise intensity in FM patients: the well-known “220 − age” at 76%, Tanaka predictive equation “208 − 0.7 × age” at 76%, the FM model HRMax “209 – 0.85 × age” at 76%, and Karvonen Formula at 60%. Bland–Altman analysis and correlation analyses were used to explore agreement and correlation between the standard method and the mathematical models. Results: Significant correlations between the heart rate at the VT1 and the four mathematical estimation models were observed. However, the Bland-Altman analysis only showed agreement between VT1 and “220 − age” (bias = −114.83 + 0.868 × x; 95% LOA = −114.83 + 0.868 × x + 1.96 × 7.46 to −114.83 + 0.868 × x − 1.96 × 7.46, where x is the average between the heart rate obtained in the CPET at VT1 and “220 − age”, in this case 129.15; p = 0.519) and “209 − 0.85 × age”(bias = −129.58 + 1.024 × x; 95% LOA = −129.58 + 1.024 × x + 1.96 × 6.619 to −129.58 + 1.024 × x − 1.96 × 6.619, where x is the average between the heart rate obtained in the CPET at VT1 and “209 − 0.85 × age”, in this case 127.30; p = 0.403). Conclusions: The well-known predictive equation “220 − age” and the FM model HRMax (“209 − 0.85 × age”) showed agreement with the standard method (VT1), revealing that it is a precise model to calculate the exercise intensity in sedentary FM patients. However, proportional bias has been detected in all the mathematical models, with a higher heart rate obtained in CPET than obtained in the mathematical model. The chronotropic incompetence observed in people with FM (inability to increase heart rate with increasing exercise intensities) could explain why methods that tend to underestimate the HRmax in the general population fit better in this population.

1. Introduction

Fibromyalgia (FM) is characterized by generalized chronic pain, fatigue, or sleep disorders, among other symptoms [1]. These symptoms are exacerbated by physical inactivity, negatively affecting the activities of daily living and the quality of life of people with FM [2,3,4]. However, according to the American Heart Association (AHA), most people with FM are below the average levels of physical conditioning. Consequently, the maximal oxygen uptake and the anaerobic threshold decrease while the pain or fatigue increases [5,6]. Thus, aerobic exercise and strength training have been strongly recommended to keep a healthy functional capacity for people with FM [7,8].
Intensity prescription for cardiorespiratory exercises is crucial for achieving health/fitness benefits. According to the overload principle of training, the exercise below a minimum intensity or threshold will not cause physiological adaptation [9]. However, this minimum threshold varies depending on different factors such as age, health status, genetics, or habitual physical activity, among others [9,10,11]. Previous articles suggested that people with FM should perform activities above the ventilatory threshold 1 (VT1) [6,12,13]. This threshold can be obtained by analyzing ventilatory parameters or lactate values in a cardiopulmonary exercise test (CPET) using a maximal or submaximal exercise incremental protocol. In this line, guidelines for exercise prescription in people with FM suggested aerobic training [7] between 40 to 60% heart rate reserve (HRR) or VO2 reserve [14]; or 64 to 76% maximal heart rate (HRmax) [9,14].
However, not all of the population can access a CEPT, either for economic reasons, knowledge, or location resources. Thus, from a pedagogical perspective, Garber, Blissmer, Deschenes, Franklin, Lamonte, Lee, Nieman, and Swain [9] recommended the application of mathematical models (formulas) to calculate the target heart rate zone and predict the maximum heart rate. Among those mathematical models are those that predict HRR and HRmax. Fox 3rd, et al. [15] proposed the well-known formula to predict the HRmax (HRmax = 220 − age) [15,16]. Although there is a discussion about this formula, it has been described as one of the most used due to its practical application and simplicity [17,18]. On the other hand, studies based on meta-analysis reported that “HRmax = 208 − 0.7 × age” is a more reliable option to predict HRmax [19,20]. However, people with FM show reduced cardiorespiratory fitness [6,12,21] as well as autonomic dysfunction (dysautonomia) [13,22,23]. Dysautonomia of people with FM could lead to chronotropic incompetence, the inability to increase heart rate with increasing exercise intensities [24]. Taking into account these problems, Lemos, Valim, Zandonade, and Natour [13] developed a specific formula to calculate the HRmax of people with FM (HRmax = 209 − 0.85 × age), showing a correlation to “220 − age” and “HRmax = 208 − 0.7 × age”. However, due to the dysautonomia in people with FM, another interesting formula is the Karvonen formula. This formula (HRB + % intensity desired [HRmax – HRB]) has been recommended for people with cardiac diseases and athletes due to the accurate characteristics of the formula (using the percentage of HRR to determine the exercise intensity) [25].
Several methods can be used to calculate exercise intensity. As previously mentioned, intensity for people with FM should be higher than VT1 in order to cause physiological adaptation. However, there are people with FM or physical trainers who could not conduct a CEPT prior to an exercise intervention. Therefore, a mathematical model to estimate the intensity of exercise is recommended. However, to the best of our knowledge, any previous studies have studied the agreement between these two recommendations [6,7,12,13,14]: (1) exercising above VT1; and (2) exercising at 60% HRR or 76% HRmax. Therefore, the aim of this study was to assess the level of agreement between the HR at VT1 and the predicted HR at which people with FM would be at 60% HRR or 76% HRmax. The results would be relevant due to the dysautonomia detected in this population that could cause predictive mathematical models to underestimate or underestimate the optimal intensity to produce physiological adaptations.

2. Materials and Methods

2.1. Participants

Eighteen women with FM participated in this cross-sectional study. Participants were recruited by a Local Fibromyalgia Association (Myos, Lisbon, Portugal). Participants were included if they fulfilled the following inclusion criteria: (1) were women between 35 and 55 years old and (2) had FM diagnosis according to the American College Rheumatology (ACR) criteria (Wolfe et al., 2010). However, participants were excluded if: (1) had a condition where exercise was contraindicated, (2) had been enrolled in physical therapy in the last six months and (3) if they had a body mass index equal to or greater than 30 (kg/m2). One was excluded from the eighteen women recruited after conducting the CPET since the values were not correctly registered. Thus, seventeen women with FM were analyzed.
Procedures were approved by the Biomedical Ethics Committee from Cordoba University (approval number: 10396/7269), following the updates of the Helsinki Declaration.

2.2. Procedures

The main socio-demographic characteristics and anthropometric data (weight and height) were collected. Furthermore, participants were assessed through a CPET to obtain the VT1. In addition, four mathematical models were used to calculate 60% HRR or 76% HRmax, which was the intensity recommended by exercise prescription guidelines for people with FM [7,14].
Cardiopulmonary exercise test(CPET). All participants performed a CPET using a maximal exercise incremental protocol in a cycle ergometer (JAEGER, 0 to 600 Watt) in order to determine the VT1. The protocol contained: (1) a 3 min warm-up; (2) the incremental test started at 20 watts; (3) each minute, the load was increased by 10 watts. HR and respiratory gases (Oxycon Pro-JAEGER, Hoechberg, Germany) were monitored during the protocol as well as during recovery (Medizintechnik System mit Erich Jaeger Gmbh®, Wuerzburg, Germany). Once the CPET was performed, VT1 was detected using Skinner’s three-phase model [26]. Participants were encouraged to report dizziness, confusion, or weakness during the CPET. If any of the participants manifested any of these symptoms, the CPET would be stopped, and blood pressure monitored. However, none of the participants reported any adverse events.
Mathematical Models. The following formulas were used to predict the HR at VT1: (1) “220 − age” at 76% HRmax [15,16]; (2) “208 − 0.7 × age” at 76% HRmax [20]; (3) “209 − 0.85 × age” at 76% HRmax) [13], and (4) “(HRB + 60% [HRMax – HRB])” at 60% HRR [25].

2.3. Statistical Analysis

Data normality was initially tested by Shapiro–Wilk test and Q-Q plots. Taking into account both the Shapiro–Wilk test and Q-Q plots, normality was assumed (see Supplementary File S1 for further details). The Bland–Altman analysis [27] was used to show agreement between the reference method VT1 obtained in a CPET and the mathematical models of training heart rate estimation. In this regard, the Bland–Altman analysis shows the dispersion graph between the difference of two variables (VT1–mathematical models) “in Y axis” and the mean of the two variables (VT1 and mathematical models) “in X axis”. To evaluate the existence of proportional bias (if the methods do not agree equally through the range of measurements), linear regressions were conducted. In this line, the differences between the methods were regressed on the average of the two methods. Following the procedure described by Bland and Altman [28], the bias (corresponding to the regression model) and the 95% limits of agreement (LOA), expressed as the regression model ± 1.96 * standard error of the residuals (the smaller range between these two limits = the better agreement) were depicted in Bland-Altman plots. The hypothesis of the bias to be or not equal to zero was tested with a t-test on paired samples. Pearson correlation tests were calculated between heart rate at VT1 and “220 − age” at 76%; heart rate at VT1 and “208 − 0.7 × age” at 76%; heart rate at VT1 and “209 − 0.85 × age” at 76%; and heart rate at VT1 and the Karvonen Formula at 60%. For all the tests, the significance level was defined at p < 0.05. All the analyses were conducted using SPSS (IBM, New York, NY, USA) (v21).

3. Results

Participants had a mean age of 47.5 ± 6.1 (years) and a body mass index of 25.7 ± 3.4 (kg/m2). Women with FM had 16.7 ± 2.0 numbers of tender points (scale 1–18) and a mean duration of symptoms of 11.5 ± 9.0 years (4.6 ± 1.5 diagnostics (years)). Participants took 2.5 ± 2.3 specific drugs (antidepressants, muscular relaxants, analgesics). Table 1 showed a poor to fair [14] aerobic fitness of people with FM (21.6 mL/min/kg). However, the physiological and cardiorespiratory parameters showed normal values (%VO2Max = 97.7 and %VO2VT1 = 69.5) [29].
Figure 1, Figure 2, Figure 3 and Figure 4 show the correlations (A) and the level of agreement (Bland–Altman plot) (B) between the gold-standard method (VT1 obtained in a CPET) and the mathematical models of HR measurement. Figure 1 shows a moderate correlation between heart rate at VT1 and “220 − age” at 76% HRmax (r = 0.563; p = 0.019) (Figure 1A). Regression showed a significant proportional bias (p-value = 0.002), and then the Bland–Altman plot was depicted according to regression results. The results from the t-test showed a significant agreement (Figure 1B) between heart rate at VT1 and the “220 − age” at 76% HRmax [15,16] (bias = −114.83 + 0.868 × x; 95% LOA = −114.83 + 0.868 × x + 1.96 × 7.46 to −114.83 + 0.868 × x − 1.96 × 7.46, where x is the average between the heart rate obtained in the CPET at VT1 and “220 − age”, in this case 129.15; p = 0.519). However, the proportional bias indicated that the heart rate obtained by CPET was higher than that obtained in the mathematical model. In addition, six measures exceeded the LOA.
Figure 2 shows a moderate correlation between the heart rate at VT1 and “208 − 0.7 × age” at 76% HRmax (r = 0.563; p = 0.019) (Figure 2A). Regression showed a significant proportional bias (p-value < 0.001), then, Bland-Altman plot was depicted according to the regression results. Results from t-test (Figure 2B) did not show a significant agreement between heart rate at VT1 and the “208 − 0.7 × age” formula [20] at 76% HRmax (bias = −157.888 + 1.024 × x; 95% LOA = −157.888 + 1.024 × x + 1.96 × 6.619 to −157.888 + 1.024 × x − 1.96 * 6.619, where x is the average between the heart rate obtained in the CPET at VT1 and “208 − 0.7 × age”, in this case 136.66; p < 0.001). The proportional bias indicated that heart rate obtained in CPET is higher than that obtained in mathematical model.
Figure 3 shows a moderate correlation between the heart rate at VT1 and “209 − 0.85 × age” at 76% HRmax (r = 0.563; p = 0.019) (Figure 3A). Regression showed a significant proportional bias (p-value < 0.001), and then the Bland–Altman plot was depicted according to the regression results. Results from the t-test (Figure 3B) showed significant agreement between the heart rate at VT1 and the “209 − 0.85 × age” formula at 76% HRmax [13] (bias = −129.58 + 1.024 × x; 95% LOA = −129.58 + 1.024 × x + 1.96 × 6.619 to −129.58 + 1.024 × x − 1.96 × 6.619, where x is the average between the heart rate obtained in the CPET at VT1 and “209 − 0.85 × age”, in this case 127.30; p = 0.403). However, the proportional bias indicated that the heart rate obtained by CPET was higher than that obtained in the mathematical model. In addition, one measure exceeded the LOA.
Figure 4 shows a strong correlation between heart rate at VT1 and the Karvonen Formula at 60% HRR (r = 0.801; p = <0.001) (Figure 4A). Regression showed a significant proportional bias (p-value = 0.001), and then the Bland–Altman plot was depicted according to the regression results. Results from the t-test (Figure 4B) did not show a significant agreement between heart rate at VT1 and the Karvonen Formula at 60% HRR (bias = −92.135 + 0.653 × x; 95% LOA = −92.135 + 0.653 × x + 1.96 × 5.518 to −92.135 + 0.653 × x − 1.96 × 5.518, where x is the average between the heart rate obtained in the CPET at VT1 and Karvonen Formula, in this case 131.04; p = 0.009) [25]. However, the proportional bias indicated that the heart rate obtained by CPET was higher than that obtained in the mathematical model.

4. Discussion

This study aimed to assess the level of agreement between the heart rate at VT1 obtained by CPET and the heart rate obtained by well-known mathematical models. Results showed an agreement between the “220 − age” equation and the FM model HRMax (“209 − 0.85 × age”) and the standard method (VT1) to predict the heart rate at VT1. The other two mathematical models, the Karvonen Formula at 60% of HRR and “208 − 0.7 × age”, did not show agreement, although significant correlations were observed. Thus, “220 − age” and “209 − 0.85 × age” could be considered the most precise mathematical model to determine the VT1 in people with FM. However, proportional bias has been detected in all the mathematical models, with a higher heart rate obtained in CPET than obtained in the mathematical model. Therefore, results should be taken with caution, and methods cannot be used interchangeably.
Our results are in line with those of a previous study, which reported correlations between different mathematical models and heart rate at VT1 in people with FM [13]. However, our study focused not only on the correlation analysis but also on analyzing the level of agreement between mathematical models and heart rate at VT1. In this regard, the Bland and Altman plot shows if an agreement exists when the measures of two methods report a bias around or equal to zero and an error that does not cause significant clinical impact to be replaced [27]. Thus, this statistical procedure would help to identify the agreement between the gold standard and alternative procedures, such as in our case between CPET and the mathematical models to predict VT1 [30]. According to our results, two mathematical models (“209 − 0.85 × age” and “220 − age”) showed agreement between CPET and formulas to predict VT1.
Regarding the “209 − 0.85 × age” predictive equation, it was specifically designed for people with FM [13]. According to our results, it was the mathematical model which showed better agreement and a lower bias. However, this method also showed a proportional bias that overestimated the CPET results. This could be due to the dysautonomia suffered by people with FM [13,22,23]. In this regard, dysautonomia could make people with FM unable to increase their heart rate while increasing exercise intensity [24]. This could explain why methods that underestimate the heart rate, such as HRmax predictive equations [20,31], fit better in people with FM than others based on the HRR. In the same line, Lemos, Valim, Zandonade, and Natour [13] recommended using the 76.2% of the HRmax (using the “220 − age” equation) and 52 to 60% of HRR (using Karvonen Formula) in order to achieve strong correlations between the mathematical models and the heart rate at VT1.
As in other studies, people with FM showed proximity between VT1 and HRmax in the CPET (Table 1), which was not observed in a healthy population [6,13]. Therefore, considering that people with FM did not achieve the maximum effort, the VT1 might be considered a better fitness index than VO2max [6,32]. This can be due to the dysautonomia observed in people with FM that alter the sympathetic response to different stress factors such as orthostasis [22,23]. This can also explain the constant fatigue and other symptoms associated with low blood pressure, such as dizziness, confusion, and weakness [33].
This study has some limitations that should be acknowledged. One limitation is the relatively small size that can contribute to the decrease in statistical power to detect agreement between the methods studied. Nevertheless, this study adds relevant information in order to prescribe exercise to women with FM. However, results must be taken with caution in men with FM as well as in other people with FM who present other physical fitness levels.

5. Conclusions

The well-known predictive equation “220-age” and the FM model HRMax (“209 − 0.85 × age”) showed agreement between the standard method (VT1 obtained in a CPET), revealing that it is a precise model to calculate the exercise intensity in sedentary FM patients. The other mathematical models using HRmax (“208 − 0.7 × age”) or HRR, which despite showing a significant correlation, did not show agreement between the heart rate at VT1 and mathematical models as well as overestimated the heart rate at VT1. However, proportional bias has been detected in all the mathematical models, with a higher heart rate obtained in CPET than obtained in the mathematical model.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11050704/s1, Figure S1. Normal Q-Q plot of heart rate at ventilatory threshold 1, Figure S2. Normal Q-Q plot of “220 – age” formula, Figure S3. Normal Q-Q plot of “220 – age” formula Figure S4. Normal Q-Q plot of “208 – 0.7 × age” formula, Figure S5. Normal Q-Q plot of Karvonen formula, Table S1. Shapiro-Wilk results of heart rate achieved in the cardiopulmonary exercise test and mathematical models.

Author Contributions

Conceptualization, C.B.-P. and P.T.-C.; Data curation, C.B.-P.; Formal analysis, S.V. and P.T.-C.; Funding acquisition, S.V. and J.A.P.; Investigation, C.B.-P.; Methodology, P.T.-C.; Resources, S.V.; Software, F.d.O.B.; Validation, C.B.-P., F.d.O.B. and P.T.-C.; Visualization, J.A.P.; Writing—original draft, S.V. and C.B.-P.; Writing—review & editing, J.A.P., F.d.O.B. and P.T.-C. All authors have read and agreed to the published version of the manuscript.

Funding

This study benefited from the support of the UÉvora—UniverCIDADE VII program. Portuguese Institute for Sportand Youth—I.P., Support for Sport Activity 2022, Sport Development Program Agreement, CP/DD/2022. The author SV was supported by a grant from the Universities Ministry of Spain and the European Union (NextGenerationUE) “Ayuda del Programa de Recualificación del Sistema Universitario Español, Modalidad de ayudas Margarita Salas para la formación de jóvenes doctores” (MS-03).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Biomedical Ethics Committee from the Cordoba University (protocol code: 10396/7269).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data will be available upon reasonable request to corresponding author.

Acknowledgments

The authors thank local Fibromyalgia association—Myos for helping us to make contact with patients. We are also very grateful to all participants of this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wolfe, F.; Smythe, H.A.; Yunus, M.B.; Bennett, R.M.; Bombardier, C.; Goldenberg, D.L.; Tugwell, P.; Campbell, S.M.; Abeles, M.; Clark, P. The American College of Rheumatology 1990 criteria for the classification of fibromyalgia. Arthritis Rheum. Off. J. Am. Coll. Rheumatol. 1990, 33, 160–172. [Google Scholar] [CrossRef] [PubMed]
  2. Segura-Jiménez, V.; Álvarez-Gallardo, I.C.; Carbonell-Baeza, A.; Aparicio, V.A.; Ortega, F.B.; Casimiro, A.J.; Delgado-Fernández, M. Fibromyalgia has a larger impact on physical health than on psychological health, yet both are markedly affected: The al-Ándalus project. Semin. Arthritis Rheum. 2015, 44, 563–570. [Google Scholar] [CrossRef] [PubMed]
  3. Carbonell-Baeza, A.; Aparicio, V.A.; Sjöström, M.; Ruiz, J.R.; Delgado-Fernández, M. Pain and functional capacity in female fibromyalgia patients. Pain Med. 2011, 12, 1667–1675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Huijnen, I.P.J.; Verbunt, J.A.; Meeus, M.; Smeets, R.J.E.M. Energy expenditure during functional daily life performances in patients with fibromyalgia. Pain Pract. 2015, 15, 748–756. [Google Scholar] [CrossRef] [PubMed]
  5. American College of Cardiology; American Heart Association. Clinical competence statement on stress testing. J. Am. Coll. Cardiol. 2000, 36, 1441–1453. [Google Scholar] [CrossRef] [Green Version]
  6. Valim, V.; Oliveira, L.M.; Suda, A.L.; Silva, L.E.; Faro, M.; Neto, T.L.B.; Feldman, D.; Natour, J. Peak oxygen uptake and ventilatory anaerobic threshold in fibromyalgia. J. Rheumatol. 2002, 29, 353–357. [Google Scholar]
  7. Bidonde, J.; Busch, A.J.; Schachter, C.L.; Overend, T.J.; Kim, S.Y.; Góes, S.M.; Boden, C.; Foulds, H.J.A. Aerobic exercise training for adults with fibromyalgia. Cochrane Database Syst. Rev. 2017, 6, CD012700. [Google Scholar] [CrossRef]
  8. Busch, A.J.; Webber, S.C.; Brachaniec, M.; Bidonde, J.; Bello-Haas, V.D.; Danyliw, A.D.; Overend, T.J.; Richards, R.S.; Sawant, A.; Schachter, C.L. Exercise therapy for fibromyalgia. Curr. Pain Headache Rep. 2011, 15, 358–367. [Google Scholar] [CrossRef] [Green Version]
  9. Garber, C.E.; Blissmer, B.; Deschenes, M.R.; Franklin, B.A.; Lamonte, M.J.; Lee, I.M.; Nieman, D.C.; Swain, D.P. Quantity and quality of exercise for developing and maintaining cardiorespiratory, musculoskeletal, and neuromotor fitness in apparently healthy adults: Guidance for prescribing exercise. Med. Sci. Sports Exerc. 2011, 43, 1334–1359. [Google Scholar] [CrossRef]
  10. Gill, T.M.; DiPietro, L.; Krumholz, H.M. Role of exercise stress testing and safety monitoring for older persons starting an exercise program. JAMA 2000, 284, 342–349. [Google Scholar] [CrossRef]
  11. Gonzalez, R.R.; Cheuvront, S.N.; Montain, S.J.; Goodman, D.A.; Blanchard, L.A.; Berglund, L.G.; Sawka, M.N. Expanded prediction equations of human sweat loss and water needs. J. Appl. Physiol. 2009, 107, 379–388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Valim, V.; Oliveira, L.; Suda, A.; Silva, L.; de Assis, M.; Neto, T.B.; Feldman, D.; Natour, J. Aerobic fitness effects in fibromyalgia. J. Rheumatol. 2003, 30, 1060–1069. [Google Scholar]
  13. Lemos, M.C.D.; Valim, V.; Zandonade, E.; Natour, J. Intensity level for exercise training in fibromyalgia by using mathematical models. BMC Musculoskelet. Disord. 2010, 11, 54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Pescatello, L.S.; Riebe, D.; Thompson, P.D. ACSM’s Guidelines for Exercise Testing and Prescription; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2014. [Google Scholar]
  15. Fox, S.M., III; Naughton, J.P.; Gorman, P.A. Physical activity and cardiovascular health. II. The exercise prescription: Intensity and duration. Mod. Concepts Cardiovasc. Dis. 1972, 41, 21–24. [Google Scholar]
  16. Fox, S.M., III; Haskell, W.L. Physical activity and the prevention of coronary heart disease. Bull. N. Y. Acad. Med. 1968, 44, 950. [Google Scholar]
  17. Goldberg, L.; Elliot, D.L.; Kuehl, K.S. Assessment of exercise intensity formulas by use of ventilatory threshold. Chest 1988, 94, 95–98. [Google Scholar] [CrossRef] [PubMed]
  18. Mesquita, A.; Trabulo, M.; Mendes, M.; Viana, J.F.; Seabra-Gomes, R. The maximum heart rate in the exercise test: The 220-age formula or Sheffield’s table? Port. J. Cardiol. 1996, 15, 139–144. [Google Scholar]
  19. Robergs, R.A.; Landwehr, R. The surprising history of the” HRmax = 220-age” equation. J. Exerc. Physiol. Online 2002, 5, 1–10. [Google Scholar]
  20. Tanaka, H.; Monahan, K.D.; Seals, D.R. Age-predicted maximal heart rate revisited. J. Am. Coll. Cardiol. 2001, 37, 153–156. [Google Scholar] [CrossRef] [Green Version]
  21. García-Hermoso, A.; Saavedra, J.M.; Escalante, Y. Effects of exercise on functional aerobic capacity in adults with fibromyalgia syndrome: A systematic review of randomized controlled trials. J. Back Musculoskelet. Rehabil. 2015, 28, 609–619. [Google Scholar] [CrossRef]
  22. Clauw, D.J.; Radulovic, D.; Antonetti, D.; Bagati, R.; Baraniuk, J.; Barbey, J.T. Tilt table testing in fibromyalgia. Arthritis Rheum 1996, 39, S276. [Google Scholar]
  23. Martínez-Lavín, M.; Hermosillo, A.G. Autonomic nervous system dysfunction may explain the multisystem features of fibromyalgia. Semin. Arthritis Rheum. 2000, 29, 197–199. [Google Scholar] [CrossRef]
  24. Da Cunha Ribeiro, R.P.; Roschel, H.; Artioli, G.G.; Dassouki, T.; Perandini, L.A.; Calich, A.L.; de Sá Pinto, A.L.; Lima, F.R.; Bonfá, E.; Gualano, B. Cardiac autonomic impairment and chronotropic incompetence in fibromyalgia. Arthritis Res. Ther. 2011, 13, R190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Karvonen, M.J. The effects of training on heart rate: A longitudinal study. Ann. Med. Exp. Biol. Fenn. 1957, 35, 307–315. [Google Scholar]
  26. Skinner, J.S.; McLellan, T.H. The transition from aerobic to anaerobic metabolism. Res. Q. Exerc. Sport 1980, 51, 234–248. [Google Scholar] [CrossRef]
  27. Bland, J.M.; Altman, D. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1986, 327, 307–310. [Google Scholar] [CrossRef]
  28. Bland, J.M.; Altman, D.G. Measuring agreement in method comparison studies. Stat. Methods Med. Res. 1999, 8, 135–160. [Google Scholar] [CrossRef]
  29. Saltin, B.; Rowell, L.B. Functional adaptations to physical activity and inactivity. Fed. Proc. 1980, 39, 1506–1513. [Google Scholar]
  30. Hirakata, V.N.; Camey, S.A. Bland-Altman analysis of agreement between methods. Rev. HCPA 2009, 29, 261–268. [Google Scholar]
  31. Sporis, G.; Vucetic, V.; Jukic, I.; Omrcen, D.; Bok, D.; Custonja, Z. How reliable are the equations for predicting maximal heart rate values in military personnel? Mil. Med. 2011, 176, 347–351. [Google Scholar] [CrossRef] [Green Version]
  32. Denis, C.; Fouquet, R.; Poty, P.; Geyssant, A.; Lacour, J.R. Effect of 40 weeks of endurance training on the anaerobic threshold. Int. J. Sports Med. 1982, 3, 208–214. [Google Scholar] [CrossRef] [PubMed]
  33. Jacomini, L.C.L.; da Silva, N.A. Dysautonomia: An emerging concept in fibromyalgia syndrome. Rev. Bras. Reumatol. 2007, 47, 354–361. [Google Scholar] [CrossRef] [Green Version]
Biology 11 00704 g001 550
Figure 1. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the mathematical model to predict the maximum heart rate “220 − age” at 76%.
Figure 1. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the mathematical model to predict the maximum heart rate “220 − age” at 76%.
Biology 11 00704 g001
Biology 11 00704 g002 550
Figure 2. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the mathematical model to predict the maximum heart rate “208 − 0.7 × age” at 76%.
Figure 2. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the mathematical model to predict the maximum heart rate “208 − 0.7 × age” at 76%.
Biology 11 00704 g002
Biology 11 00704 g003 550
Figure 3. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the mathematical model to predict the maximum heart rate “209 − 0.85 × age” at 76%.
Figure 3. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the mathematical model to predict the maximum heart rate “209 − 0.85 × age” at 76%.
Biology 11 00704 g003
Biology 11 00704 g004 550
Figure 4. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the Karvonen Formula at 60%.
Figure 4. Correlations (A) and Bland–Altman plots (B) between the heart rate achieved in the VT1 in the cardiopulmonary exercise testing (CPET) and the Karvonen Formula at 60%.
Biology 11 00704 g004
Table
Table 1. Cardiorespiratory parameters obtained in the CPET and predictive measures of heart rate of training by mathematical models of women with fibromyalgia.
Table 1. Cardiorespiratory parameters obtained in the CPET and predictive measures of heart rate of training by mathematical models of women with fibromyalgia.
Cardiorespiratory Parameters Obtained in the CPET
Basal Heart Rate (bpm)78.4 ± 11.6
Heart Rate at VT1 (bpm)127.2 ± 15.4
Maximum attained Heart Rate (bpm)145.3 ± 17.6
Maximum predicted Heart Rate 220-idade (bpm)173 ± 6.6
VO2VT1 (mL/min)1033.5 ± 175.0
VO2Max (mL/min)1363.4 ± 170.0
VO2Max predicted (mL/min)1486.7 ± 148.6
VO2Max (mL/min/kg)21.6 ± 4.3
VO2VT1% predicted69.5%
VO2Max% predicted97.7%
Mathematical Models and Predictive Measures of Heart Rate of Training
HRB + 60% (HRMax–HRB) [bpm] 134.8 ± 6.7
Heart Rate (220 − age) × 76% [bpm] 131 ± 5.7
Heart Rate (208 − 0.7 × age) × 76% [bpm]146 ± 4.9
Heart Rate (209 − 0.85 × age) × 76% [bpm]127.3 ± 4.9
Data expressed in Mean ± Standard deviation. VO2Max: maximum oxygen uptake; VT1: Ventilatory threshold 1; HRB: basal heart rate; HRMax: maximum heart rate; n.a.: not applicable.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Villafaina, S.; Biehl-Printes, C.; Parraca, J.A.; de Oliveira Brauner, F.; Tomas-Carus, P. What Mathematical Models Are Accurate for Prescribing Aerobic Exercise in Women with Fibromyalgia? Biology 2022, 11, 704. https://doi.org/10.3390/biology11050704

AMA Style

Villafaina S, Biehl-Printes C, Parraca JA, de Oliveira Brauner F, Tomas-Carus P. What Mathematical Models Are Accurate for Prescribing Aerobic Exercise in Women with Fibromyalgia? Biology. 2022; 11(5):704. https://doi.org/10.3390/biology11050704

Chicago/Turabian Style

Villafaina, Santos, Clarissa Biehl-Printes, José A. Parraca, Fabiane de Oliveira Brauner, and Pablo Tomas-Carus. 2022. "What Mathematical Models Are Accurate for Prescribing Aerobic Exercise in Women with Fibromyalgia?" Biology 11, no. 5: 704. https://doi.org/10.3390/biology11050704

APA Style

Villafaina, S., Biehl-Printes, C., Parraca, J. A., de Oliveira Brauner, F., & Tomas-Carus, P. (2022). What Mathematical Models Are Accurate for Prescribing Aerobic Exercise in Women with Fibromyalgia? Biology, 11(5), 704. https://doi.org/10.3390/biology11050704

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