- 1III Klinika Chorób Wewnętrznych i Kardiologii, Warszawski Uniwersytet Medyczny (WUM), Warsaw, Poland
- 2Department of Physical Education and Health in Biala Podlaska, Jozef Pilsudski University of Physical Education in Warsaw Faculty in Biala Podlaska, Biala Podlaska, Poland
- 3Public Health School Centrum Medyczne Kształcenia Podyplomowego (CMKP), Warsaw, Poland
- 4Wydział Matematyki i Nauk Informacyjnych, Politechnika Warszawska, Warsaw, Poland
- 5Department of Experimental and Clinical Pharmacology, Center for Preclinical Research and Technology (CEPT), Medical University of Warsaw, Warsaw, Poland
- 6Institute of Primary Care, University of Zurich, Zurich, Switzerland
- 7Medbase St. Gallen Am Vadianplatz, St. Gallen, Switzerland
Maximal heart rate (HRmax) is associated mostly with age, but age alone explains the variance in HRmax to a limited degree and may not be adequate to predict HRmax in certain groups. The present study was carried out on 3374 healthy Caucasian, Polish men and women, clients of a sports clinic, mostly sportspeople, with a mean age of 36.57 years, body mass 74.54 kg, maximum oxygen uptake (VO2max, ml∗kg–1 ∗min–1) 50.07. Cardiopulmonary exercise tests (CPET) were carried out on treadmills or cycle ergometers to evaluate HRmax and VO2max. Linear, multiple linear, stepwise, Ridge and LASSO regression modeling were applied to establish the relationship between HRmax, age, fitness level, VO2max, body mass, age, testing modality and body mass index (BMI). Mean HRmax predictions calculated with 5 previously published formulae were evaluated in subgroups created according to all variables. HRmax was univariately explained by a 202.5–0.53∗age formula (R2 = 19.18). The weak relationship may be explained by the similar age with small standard deviation (SD). Multiple linear regression, stepwise and LASSO yielded an R2 of 0.224, while Ridge yielded R2 0.20. Previously published formulae were less precise in the more outlying groups of the studied population, overestimating HRmax in older age groups and underestimating in younger. The 202.5–0.53∗age formula developed in the present study was the best in the studied population, yielding lowest mean errors in most groups, suggesting it could be used in more active individuals. Tanaka’s formula offers the second best overall prediction, while the 220-age formula yields remarkably high mean errors of up to 9 bpm. In conclusion, adding the studied variables in multiple regression models improves the accuracy of prediction only slightly over age alone and is unlikely to be useful in clinical practice.
Introduction
Heart rate (HR) is a commonly measured parameter, often used in clinical practice, sports and scientific research; it is easy to reliably measure with very little equipment (Robergs and Landwehr, 2002). HR increases in a linear way with increasing physical exertion, until a maximum (HRmax) is reached by the individual at maximal workload (Kostis et al., 1982). HRmax is useful to prescribe exertion levels in sports training, or to carry out electrocardiogram (ECG) or ECHO exercise stress tests (Robergs and Landwehr, 2002). HRmax is measured in a graded exercise test, often along with other parameters like maximum oxygen uptake (VO2max) or respiratory exchange ratio (RER), usually with a treadmill or cycle ergometry, until maximum exertion is achieved (Pescatello et al., 2006; Beltz et al., 2016). Several test protocols exist which are commonly used, depending on the clinic’s experience, or the type of patient (e.g., sportspeople or patients with cardiovascular disease). These tests may use a stepwise increase in speed/Watts, or a ramped increase and they may also vary in length. All of these differences make accurate comparison of the results of different studies on maximum exertion parameters difficult (Beltz et al., 2016). Since determining an individual’s actual HRmax (with the use of a maximal exercise test) is difficult and not always possible or advisable, it is usually estimated with the use of several formulae (Robergs and Landwehr, 2002).
The most simple and widely used is the 220-age formula, the origin of which is unclear, but first appeared in scientific writing in a review by Fox and Haskell (1970). The downside to its simplicity is the high standard error of estimate (SEE) of ∼7–12 beats per minute (bpm) (Robergs and Landwehr, 2002). Tanaka et al. (2001) described a new formula (208–0.7∗age) in 2002, calculated from a meta-analysis of 351 studies involving 18,712 subjects, which was then validated on a group of over 500 subjects. This new formula had a SEE of ∼10 bpm. The Tanaka equation is currently quite often applied, alone or with the Fox equation (Nikolaidis et al., 2018). There are also less frequently used equations, including Inbar’s formula (Inbar et al., 1994), or the Londeree and Moeschberger formula (Londeree and Moeschberger, 1982). One of the more recent formulae is that published by Nes et al. (2013).
The HRmax is difficult to predict exactly (Robergs and Landwehr, 2002). This is partly due to the fact that the rate of decline of HRmax is non-linear, as demonstrated by Zhu et al. (2010) the rate of decline is significantly different in various age groups, lower in the younger population, and higher in the older population, actually being curvilinear. Londeree and Moeschberger (1982) conducted a meta-analysis and established with the use of multivariate analysis that age accounts for 70–75% of the HRmax variance, the other factors being sex, level of fitness, type of ergometer used, continent of residence and race. The study did not, however, take body composition or BMI into account and the fitness level was not a continuous variable based on experimental protocols, but simply a categorization as sedentary, active, or athlete. The accuracy of commonly used HRmax prediction equations varies strongly in different groups, often yielding errors higher than 10 bpm, which then impact on the adequacy of training and clinical stress tests (Robergs and Landwehr, 2002; Nes et al., 2013).
Few studies evaluate the accuracy of currently available formulae for HR max calculation in physically active individuals. Also, many factors other than age have very rarely been considered, for example training modality or body composition. The aim of this study was to compare the accuracy of currently used equations and assess how the addition of other variables, such as BMI, body composition or maximum oxygen uptake (VO2max) impacts the precision of the calculations. We also set out to design an optimal formula using multiple linear regression. We also attempted to verify the accuracy of commonly used formulae in a physically active population. We had the following hypotheses:
Body composition could perhaps impact HRmax due to differences in fat/muscle proportions, as recent studies demonstrated for VO2max and respiratory compensation point (Maciejczyk et al., 2014a,b).
• Higher aerobic capacity may impact HRmax, perhaps due to adaptations in the heart muscle to training (Whyte et al., 2008).
• The accuracy of certain tests could be worse in more active, younger individuals according to a recent study (Shookster et al., 2020).
Materials and Methods
Participants
Participants included clients of the Sportslab clinic1 in years 2013–2019 who had commercial cardiopulmonary exercise tests (CPET) performed. They were recruited via internet and social media advertisements, or via recommendation from trainers or other clients. The tests were carried out on personal request of the participants as part of training optimization and diagnostics. Most participants were runners, cyclists or triathletes of various levels. A detailed breakdown of sports disciplines trained and the level of performance are not available. Inclusion criteria for the database were age over 18 years, training for at least 3 months and meeting the maximum exertion criteria described below. Exclusion criteria were any chronic or acute medical conditions (including musculoskeletal system disorders like new fractures and sprains, as well as addiction to nicotine, alcohol or other substances) or ongoing intake of any medication. Over 4,000 tests were performed, out of which N = 3374 met the inclusion criteria for the study. During the entire period, the same methods and procedures were used. Participants received information via e-mail on how to prepare for the test. Participants received information via e-mail on how to prepare for the test. Participants were advised to prepare for the test by avoiding any exercises 48 h prior to the test, as well as eating a light carbohydrate meal 2–3 h before the test and keeping good hydration by drinking isotonic beverages. They were also instructed to avoid medicines, caffeine and cigarettes before the test.
Populational data was calculated as means for men, women and the entire population. There was a total of n = 2,816 treadmill tests (n = 510 female) and n = 958 cycle ergometer tests (n = 110 female). Basic anthropometric data as well as VO2max are presented in Table 1.
Incremental Exercise Test
Body mass and fat mass were determined with a body composition analyzer (Tanita, MC 718, Japan) before every test using the method of multifrequency 5 kHz/50 kHz/250 kHz electrical bioimpedance. The body composition tests were conducted directly prior to exercise testing. All measurements (body composition and CPET) took place under similar conditions in the medical clinic Sportslab (Warsaw, Poland). The conditions were: 40 m2 of indoor, air-conditioned space, altitude 100 m, MSL temperature 20–22 degrees Centigrade, 40–60% humidity. The cardiopulmonary exercise tests (CPET) were a part of commercial endurance testing offered in the clinic. Apart from medical contraindications, there were no other exclusion criteria to perform tests.
Exercise tests were performed on a cycle ergometer Cyclus-2 (RBM elektronik-automation GmbH, Leipzig, Germany) or on a mechanical treadmill (h/p/Cosmos quasar, Germany), depending on the dominant physical activity of each individual (running or cycling). During all tests cardio-pulmonary indices like HR and oxygen uptake (VO2) were recorded using a Cosmed Quark CPET device (Rome, Italy), calibrated before each test according to the manufacturer’s instructions. HR was measured using the ANT + chest strap which is part of the Cosmed Quark CPET device (declared accuracy to ECG, ± 1 bpm.).
For each person, the initial power (Watt) or speed (km/h) were determined based on an interview carried out before the test. The initial power for cycle ergometer tests was the lowest power at which the participant subjectively felt resistance. For treadmill tests the start speed was a slow running pace, also based on subjective experience, between 5 and 12 km/h based on interview. The test began after 5 min warm up (walking or pedaling with no resistance). The speed was increased by 1 km/h or the power was increased by 20–30 W every 2 min. For treadmill tests, 1% incline was applied.
In order to assess the maximum level of aerobic fitness, participants were instructed to keep the effort as long as possible. They could terminate the test at any moment if they felt they had reached their limit of exertion. Participants were under cardiopulmonary monitoring during the test. The test was terminated by the operator if either VO2 or HR showed no further increase with increasing speed/power.
The results of the body composition analysis and CPET were saved as an Excel (Microsoft corporation, Washington, United States) spreadsheet. The raw data were anonymized and processed with the use of a custom program created in Python in order to identify data like VO2max, HRmax, or anaerobic threshold. CPET data were recorded breath by breath and then averaged across 10-s intervals to limit file sizes. The highest HR in the interval was recoded and HR values were not averaged.
For statistical evaluation we included only cases where 3 of 4 following criteria were met: respiratory exchange ratio (RER) during test reaching > 1.10, VO2 plateau (an increase in VO2 with increasing speed/power lower than 100 ml/min), respiratory frequency over 45/min, perceived exertion over 18 in Borg scale (Kaminsky et al., 2015).
Data Analysis
Statistical analysis was carried out with the use of R language in Rstudio (Rstudio, PBC, United States). All obtained data were cleaned using Rstudio. Entries with erroneously high or low results were removed based on the following adopted range constraints: age (years)—18–100; body mass (kg)—40–200; VO2max (ml/min)—1,000–8,000; BMI—10–50; body fat percentage—4–60%; HRmax—30–250; VO2max (ml⋅kg–1⋅min–1)—10–100; finally, results were also manually verified.
Apart from the measured variables, an additional variable was created, which was labeled as fitness level and was calculated as follows. The population was divided into four age groups according to age percentiles (0–10, 10–50, 50–90, 90–100), separately for males and females, yielding eight groups. The distribution of VO2max in each group was analyzed and the difference of the VO2max result from the average for the given age group was defined as fitness level (Figures 1, 2). The results were smoothed using the generalized additive model method. The higher the fitness level, the higher the VO2max in relation to the age group’s average.
Figure 1. Fitness level distribution in males (A) and females (B), the gray dotted lines represent the 10th, 50th, and 90th quantile in the population.
Figure 2. Univariate linear regression of maximum heart rate (HRmax) and age (202.5–0.53*age; R2 = 0.198).
Several regression models (multiple linear regression, stepwise regression, LASSO regression and Ridge regression) were applied to establish the relation between HRmax and six variables which were considered: age, VO2max, body mass, body mass index (BMI), body fat and fitness level. Two data subsets were randomly created: a training subset of 70% for building the regression model and a test subset of 30% in order to avoid overfit and provide test data and unseen data for the model. Normality of residuals was tested, and the residuals were normally distributed.
Previously published formulae by Fox and Haskell (1970), Tanaka et al., Nes et al., Londeree et al., and Inbar et al. were evaluated by calculating the mean error of the equations on the training data subset. The training subset was then further divided according to age (every 10 years), sex, fitness level (3 equally sized groups selected with the use of a histogram), test used (cycle ergometer or treadmill), BMI (<20, 20–25, >25) and body fat percentage separately for males and females according to normal ranges (American Council on Exercise [ACE], n.d.).
Results
The fitness levels in males and females are presented in the form of graphs in Figure 1. The gray dotted lines represent the 10th, 50th, and 90th quantile in the population.
Starting speeds and power are shown in Table 1.
HRmax was linearly and inversely related to age and explained by the regression formula 202.5–0.53∗age (Figure 2). This formula’s R2 was 19.18, thus explaining 19.18% of the observed variability. Results of the multiple regression models are presented in Table 2. In the multivariate formulae, body mass was significantly (p < 0.01) inversely linked with HRmax but contributed only 0.01 to the R2. Body fat percentage and fitness level were positively and significantly (p < 0.05) related to HRmax, contributing 0.00003 and 0.02 to the overall R2, accordingly. Sex and testing device were also significantly related to HRmax contributed to the R2 (0.01, p < 0.01 and 0.002, p = 0.01, accordingly). BMI and VO2max were not related to HRmax. The multivariate formula was:
229–0.64∗age –0.23∗body mass + 0.02∗BMI –0.38∗ VO2max + 0.33∗body fat + 0.02∗fitness level + 8.74∗ sex + 0.97∗testing modality.
Where VO2max is expressed in ml∗kg–1*min–1, age in years, body mass in kg, sex is a binary variable: 1 for male, 0 for female and testing modality is binary: 1 for treadmill, 0 for cycle ergometer. The R2 for the formula was 0.224, MAE = 7.04, ME = 0.22.
Stepwise regression yielded a formula with 7 variables, excluding only BMI. The formula was:
229–0.64∗age–0.23∗body mass −0.38∗ VO2max + 0.33∗body fat + 0.44∗fitness level + 8.74∗sex + 0.97∗testing modality.
The R2 for the formula was also 0.2224, ME = 0.22 and MAE = 7.04, while using 1 variable fewer, therefore this model was selected as the most adequate and further discussed. Results of the regression are presented in Table 2. A graph illustrating the accuracy of prediction using this formula is presented in Figure 3.
Figure 3. Graph of observed (y-axis) vs. predicted (x-axis) HRmax (bpm) in the stepwise regression model.
The LASSO regression formula was:
225–0.62∗age–0.23∗body mass + 0.02∗BMI −0.32∗ VO2max + 0.33∗body fat + 0.37∗fitness level + 8.32∗ sex + 0.97∗testing modality.
R2 = 0.22 and lambda = 0.002, MAE = 7.04, ME = 0.22.
The Ridge regression formula was:
202–0.46∗age −0.05∗body mass + 0.03∗ VO2max.
R2 = 0.20 and lambda = 0.45, MAE = 7.11, ME = 0.21.
The results of the comparison of different formulae in several subgroups are presented in Table 3. Fitness levels 1, 2, and 3 correspond to three equally sized groups divided according to a histogram, where level 1 is the lowest tertile.
Discussion
The results of this paper demonstrate that adding additional variables other than age to HRmax prediction formulae brings only a very small (although statistically significant) improvement in R2.
The studied population represents healthy males and females, including sportspeople of various levels, due to the activity profile of the Sportslab clinic. Despite no activity level questionnaire having been recorded, the relatively high fitness level can be observed by comparing VO2max in the present study with Kaminsky et al.’s (2015) reference ranges. Few other studies analyzed VO2max as a factor potentially influencing HR max (Nes et al., 2013). Nes et al. analyzed VO2max as two separate variables, one continuous and the other expressed as VO2max tertiles in the studied population and neither was found to significantly impact the HRmax estimation formula’s accuracy. In the present study a variable labeled fitness level was created based on VO2max adjusted for average VO2max in the participant’s age group. VO2max is known to decrease with age and therefore is not sufficient to express an individual’s fitness level relative to people of different age (Betik and Hepple, 2008). For example an older athlete who is very fit relative to his peers would have a VO2max which may be no higher than the average in a younger population. The new variable expresses the fitness level of the individual independent of age and in a continuous way. To the best of our knowledge, such a variable has not been used in previous studies and therefore direct comparison is impossible at present, but we believe it is more useful for the assessment of fitness than VO2max alone or self-assessed fitness.
The R2 of the 202.5–0.53∗age formula was 19.18, thus explaining 19.18% of the observed variability. This is considerably lower than values recorded in the studies of Tanaka et al. (80%) or Londeree et al. (70–75%) (Londeree and Moeschberger, 1982; Tanaka et al., 2001). The main factor influencing this difference is most likely the narrower age distribution in our group, the SD being only 8.43. Neither Londeree nor Tanaka published the average age and its standard deviation in their papers, but the age range was 5–81 in Londeree’s paper and 18–81 in that of Tanaka et al. and 18–74 in this paper. The average age was 36 (± 8) years, therefore only approximately 5% of the population were over 52, representing a rather young group with small age variance. In this case, age naturally explains much less of the variation, as it is more similar across the studied population. Nes et al. studied a group (N = 3320) with a mean age of 46 (± 13) years, a standard deviation still much higher than in the present study and explained 36% of the variability with their regression formula, supporting this line of thought. The more similar the age is in the group, the lower the proportion of the variance that can be explained by age. Therefore, studies with a normal age distribution across a wide age range demonstrate a high correlation of age and maximum heart rate, but do not address the issue of differences in maximum heart rate in individuals of the same age. These differences must be explained by other factors, which remain largely unknown, as demonstrated by the low R2 in studies on populations with a narrow age range. Another factor influencing the lower R2 might be that the more fit one is, the slower the HRmax decline with age, which further limits the amount of variance explained by age alone. It is difficult to otherwise explain the above-mentioned high values demonstrated by Tanaka et al. and Londeree et al. due to the little data they offer on the statistical methods and results, and little population data. This lack of sufficient detail was also noted by Robergs and Landwehr in their hallmark review (Robergs and Landwehr, 2002).
The multivariate models revealed only a small rise in R2 after adding body weight, BMI, body fat percentage, fitness level, sex and testing modality into the equations. Most of the added variables apart from BMI affected the model significantly, but to a small extent, together explaining approximately only additional 3% of the observed variance. This does not seem to justify practical usage of a much more complex multivariate formula. These observations are similar to those of previous researchers who attempted to provide a multivariate model (Londeree and Moeschberger, 1982; Nes et al., 2013), but the impact of body composition and fitness level on HRmax had not previously been properly studied. It has been observed in several studies that HRmax decreases with aerobic training and increases with detraining in endurance athletes (Zavorsky, 2000). However, these studies had a mean of only 10 participants (ranging from 5 to 24) and results were not consistent. The present study shows that fitness level is responsible for only a small part of the variance in HRmax. Few other papers exist which would explain the variance in HRmax with any other variables than age. Graettinger et al. found that relative wall thickness in echocardiography explained 9% of the variance in hypertensive individuals, but not in healthy ones (Graettinger et al., 1995). Londeree and Moeschberger (1982) concluded that race is not an important factor, but their paper lacks demographic data and the number of people of different races is unknown.
We found that the apart from the presently developed formula, the second most accurate formula in terms of lower ME across the subgroups was Tanaka’s. These observations support the results of previous smaller studies on highly active individuals such as amateur marathon runners or physically active young males, which also revealed that HRmax predicted with Tanaka’s formula was similar to observed HRmax (Barboza et al., 2016; Nikolaidis et al., 2018). Previously authors observed that the 220-age formula tends to overestimate the HRmax in the younger population and underestimate it in the older population (Londeree and Moeschberger, 1982; Nes et al., 2013). This is also visible in the present study, where the 220-age formula underestimates HRmax by a mean of 7 bpm in the oldest group, and underestimates by a mean of 9 bpm in the youngest group, while being relatively accurate only in the age group 30–39 years. This is a serious consideration especially when prescribing exercise stress tests in older individuals based on their HRmax, as such large errors may lead to invalid results or excessive cardiac stress (Lauer et al., 2005). Also when prescribing exercise intensity to young sportspeople, training with a significantly underestimated HRmax may be inefficient (Jain et al., 2011). We observed that many well established organizations still recommend the use of the 220-age formula, despite a significant body of work advising its discontinuation and replacement with Tanaka’s formula (Tanaka et al., 2001; Robergs and Landwehr, 2002; Nes et al., 2013; Roy and McCrory, 2015; Nikolaidis et al., 2018). When searching the phrase “maximum heart rate training” in the Google search engine (as of 16.11.2020), the first two results are from mayoclinic.com and heart.org (the American Heart Association’s website) and both suggest the use of 220-age (Mayo Clinic Staff [MCS], 2019; Target Heart Rates Chart [THRC], n.d.). We observed that Tanaka’s formula (and all the other formulae except that of Nes et al.) also overestimated the maximum heart rate in the 40–49 age group by several bpm, perhaps due to the population in this paper being relatively fit, and thus demonstrating a lower maximum heart rate (Zavorsky, 2000). The 202.5–0.53∗age formula displayed a lower mean error than Tanaka’s and other formulae in most subgroups, especially in the 30–49 age range, but its accuracy was lower in older subgroups, probably due to the low number of subjects aged 50 or above. The formula developed by Nes et al. underestimated the results by a ME of 4.32 in the entire group (Nes et al., 2013). This discrepancy is difficult to explain. The population studied by Nes et al. does not seem to be physically more active than the one in the present study. Criteria for a maximal heart rate test were also well defined in their study, therefore the possibility of their subjects not reaching HRmax in unlikely.
Robergs et al. described Inbar’s formula to be the most precise. However, this conclusion was based entirely on the error described by the authors from their study on 1,424 healthy subjects (Inbar et al., 1994). In our study group, Inbar’s formula tended to overestimate HRmax in most subgroups and was overall less precise than Tanaka’s formula (Inbar et al., 1994; Tanaka et al., 2001).
Tanaka’s formula also appears to slightly overestimate HRmax in males, while slightly underestimating in females, but the differences are approximately 1 bpm, therefore using targeted male or female equations is likely to be unnecessary (Tanaka et al., 2001). Tanaka, Haskell and Nes underestimate the HRmax in individuals with low BMI (<20). It has previously been shown that the Tanaka formula is accurate in overweight individuals (BMI > 25), but little data exists on HRmax in individuals with lower BMI (Franckowiak et al., 2011).
Study Limitations
The study was carried out in a single center. Activity level questionnaires were not available, but the activity profile of the group was analyzed based on VO2max comparison with reference ranges relative to age. Two different modes of testing were used. The initial loads were selected subjectively. Bioimpedance was used for body composition measurements, which may be influenced by the hydration status. Due to the demographics in Poland, the study does not present racial differences. The group included mostly endurance athletes with little representation of other sports. However, a more homogenous group allows for drawing more valid conclusions for this group.
Conclusion
The present study demonstrates that factors other than age, including body composition, gender, fitness level, VO2max, BMI, or testing modality, add little to the accuracy of HRmax estimation. We demonstrated that multiple regression is only slightly superior to linear regression, which does probably not justify the use of multiple regression in practice.
Clinical Implications
The 202.5–0.53∗age formula developed in the present study was the best in the studied population, yielding lowest mean errors in most groups, suggesting it could be used in more active individuals.
A simple formula predicting HRmax based on age only may be used when the exact HRmax is not needed or is difficult to obtain. Tanaka’s formula proved to be fairly accurate in the studied population, although slightly underperforming the calculated formulae. All the formulae have a MAE of approximately 7–8 bpm and when possible, HRmax should be tested for precise results.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Komisja Bioetyczna przy Warszawskim Uniwersytecie Medycznym. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JL and SW: conceptualization and resources. JL, SW, DŚ, and AM: investigation. MZ, SP, SW, and IC: statistical analysis. SP: writing—original draft preparation. JL, SW, DŚ, BK, and MZ: writing—review and editing. MP, BK, and AM: supervision. All authors have read and agreed to the published version of the manuscript.
Funding
The study was funded by the affiliated university (Medical University of Warsaw).
Conflict of Interest
SW was the owner of the Sportslab clinic.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Katarzyna Wiecha and Paweł Witecki for their help in creating the database for this study. We also thank I. Cieśliński for his assistance in reviewing the paper.
Footnotes
References
[ACE] American Council on Exercise [ACE]. (n.d.). Percent Body Fat Norms for Men and Women. Retrieved December 1, 2020, from https://www.acefitness.org/education-and-resources/lifestyle/tools-calculators/percent-body-fat-calculator/
Barboza, J., Nogueira, F. S., Paulucio, D., and Pompeu, F. (2016). The validity of maximum heart rate predictions based on age for young men. Gazz. Med. Ital. 175, 154–157.
Beltz, N. M., Gibson, A. L., Janot, J. M., Kravitz, L., Mermier, C. M., and Dalleck, L. C. (2016). Graded Exercise Testing Protocols for the Determination of VO(2)max: historical Perspectives, Progress, and Future Considerations. J. Sports Med. 2016:3968393. doi: 10.1155/2016/3968393
Betik, A. C., and Hepple, R. T. (2008). Determinants of VO2 max decline with aging: an integrated perspective. Appl. Physiol. Nutr. Metab. 33, 130–140. doi: 10.1139/H07-174
Fox, S. M., and Haskell, W. L. (1970). “The exercise stress test: needs for standardization,” In Cardiology: Current topics and progress, eds M. Eliakim and H. N. Neufeld (New York: Academic Press), 6, 149–154.
Franckowiak, S. C., Dobrosielski, D. A., Reilley, S. M., Walston, J. D., and Andersen, R. E. (2011). Maximal heart rate prediction in adults that are overweight or obese. J. Strength Cond. Res. 25, 1407–1412. doi: 10.1519/JSC.0b013e3181d682d2
Graettinger, W. F., Smith, D. H., Neutel, J. M., Myers, J., Froelicher, V. F., and Weber, M. A. (1995). Relationship of left ventricular structure to maximal heart rate during exercise. Chest 107, 341–345. doi: 10.1378/chest.107.2.341
Inbar, O., Oren, A., Scheinowitz, M., Rotstein, A., Dlin, R., and Casaburi, R. (1994). Normal cardiopulmonary responses during incremental exercise in 20- to 70-yr-old men. Med. Sci. Sports Exerc. 26, 538–546.
Jain, M., Nkonde, C., Lin, B. A., Walker, A., and Wackers, F. J. T. (2011). 85% of maximal age-predicted heart rate is not a valid endpoint for exercise treadmill testing. J. Nucl. Cardiol. 18, 1026–1035. doi: 10.1007/s12350-011-9454-0
Kaminsky, L. A., Arena, R., and Myers, J. (2015). Reference Standards for Cardiorespiratory Fitness Measured With Cardiopulmonary Exercise Testing: data From the Fitness Registry and the Importance of Exercise National Database. Mayo Clin. Proc. 90, 1515–1523. doi: 10.1016/j.mayocp.2015.07.026
Kostis, J. B., Moreyra, A. E., Amendo, M. T., Di Pietro, J., Cosgrove, N., and Kuo, P. T. (1982). The effect of age on heart rate in subjects free of heart disease. Studies by ambulatory electrocardiography and maximal exercise stress test. Circulation 65, 141–145. doi: 10.1161/01.cir.65.1.141
Lauer, M., Froelicher, E. S., Williams, M., and Kligfield, P. (2005). Exercise testing in asymptomatic adults: a statement for professionals from the American Heart Association Council on Clinical Cardiology, Subcommittee on Exercise, Cardiac Rehabilitation, and Prevention. Circulation 112, 771–776. doi: 10.1161/CIRCULATIONAHA.105.166543
Londeree, B. R., and Moeschberger, M. L. (1982). Effect of Age and Other Factors on Maximal Heart Rate. Res. Q. Exerc. Sport 53, 297–304. doi: 10.1080/02701367.1982.10605252
Maciejczyk, M., Wiecek, M., Szymura, J., Cempla, J., Wiecha, S., Szygula, Z., et al. (2014a). Effect of body composition on respiratory compensation point during an incremental test. J. Strength Cond. Res. 28, 2071–2077. doi: 10.1519/JSC.0000000000000347
Maciejczyk, M., Więcek, M., Szymura, J., Szyguła, Z., Wiecha, S., and Cempla, J. (2014b). The influence of increased body fat or lean body mass on aerobic performance. PloS One 9:e95797. doi: 10.1371/journal.pone.0095797
Mayo Clinic Staff [MCS]. (2019). Exercise Intensity: How to Measure it. Retrieved November 16, 2020, from Mayo Clinic website: Available online at: https://www.mayoclinic.org/healthy-lifestyle/fitness/in-depth/exercise-intensity/art-20046887#:~:text=You can calculate your maximum, beat per minute during exercise (accessed April 11, 2021).
Nes, B. M., Janszky, I., Wisløff, U., Støylen, A., and Karlsen, T. (2013). Age-predicted maximal heart rate in healthy subjects: the HUNT fitness study. Scand. J. Med. Sci. Sports 23, 697–704. doi: 10.1111/j.1600-0838.2012.01445.x
Nikolaidis, P. T., Rosemann, T., and Knechtle, B. (2018). Age-Predicted Maximal Heart Rate in Recreational Marathon Runners: a Cross-Sectional Study on Fox’s and Tanaka’s Equations. Front. Physiol. 9:226. doi: 10.3389/fphys.2018.00226
Pescatello, L., Arena, R., Riebe, D., and Thompson, P. D. (2006). ACSM’s Guidelines for Exercise Testing and Prescription. Philadelphia: Lippincott Williams & Wilkins. 60–113.
Robergs, R., and Landwehr, R. (2002). The surprising history of the “HRmax=220-age” equation. J. Exerc. Physiol. 5, 1–10.
Roy, S., and McCrory, J. (2015). Validation of Maximal Heart Rate Prediction Equations Based on Sex and Physical Activity Status. Int. J. Exerci. Sci. 8, 318–330.
Shookster, D., Lindsey, B., Cortes, N., and Martin, J. R. (2020). Accuracy of Commonly Used Age-Predicted Maximal Heart Rate Equations. Int. J. Exerc. Sci. 13, 1242–1250.
Tanaka, H., Monahan, K. D., and Seals, D. R. (2001). Age-predicted maximal heart rate revisited. J. Am. Coll. Cardiol. 37, 153–156. doi: 10.1016/s0735-1097(00)01054-8
Target Heart Rates Chart [THRC]. (n.d.). Target Heart Rates Chart. Retrieved November 16, 2020, from Heart.org website: https://www.heart.org/en/healthy-living/fitness/fitness-basics/target-heart-rates
Whyte, G. P., George, K., Shave, R., Middleton, N., and Nevill, A. M. (2008). Training induced changes in maximum heart rate. Int. J. Sports Med. 29, 129–133. doi: 10.1055/s-2007-965783
Zavorsky, G. S. (2000). Evidence and possible mechanisms of altered maximum heart rate with endurance training and tapering. Sports Med. 29, 13–26. doi: 10.2165/00007256-200029010-00002
Keywords: hrmax, 220-age, cardiopulmonary testing, formulae, body composition, aerobic performance, treadmill ambulation, cycle ergometer
Citation: Lach J, Wiecha S, Śliż D, Price S, Zaborski M, Cieśliński I, Postuła M, Knechtle B and Mamcarz A (2021) HR Max Prediction Based on Age, Body Composition, Fitness Level, Testing Modality and Sex in Physically Active Population. Front. Physiol. 12:695950. doi: 10.3389/fphys.2021.695950
Received: 15 April 2021; Accepted: 14 June 2021;
Published: 30 July 2021.
Edited by:
Anthony S. Leicht, James Cook University, AustraliaReviewed by:
Daniel Boullosa, Federal University of Mato Grosso do Sul, BrazilTrine Karlsen, Nord University, Norway
Copyright © 2021 Lach, Wiecha, Śliż, Price, Zaborski, Cieśliński, Postuła, Knechtle and Mamcarz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel Śliż, ZGFuaWVsLnNsaXpAd3VtLmVkdS5wbA==; c2xpei5kYW5pZWxAZ21haWwuY29t