In 2016, the American Heart Association launched a series of publications promoting the clinical evaluation of cardiorespiratory fitness (CRF) with the overall aim of improving the prevention and treatment of cardiovascular disease (CVD; Ross et al., 2016). Furthermore, the association urged the US federal government to compile a registered CRF database (Kaminsky et al., 2013); this highlights the importance of CRF. CRF is generally defined as the integrated ability to transport oxygen from the atmosphere to the mitochondria for physical activity. Notably, CRF involves the respiratory, circulatory, and neuromuscular systems and has a clear and direct relationship with the functions of various systems. Individuals with weak CRF have an up to 70% all-cause mortality rate and 56% cardiovascular mortality rate (Kodama et al., 2009). Similarly, every 1-MET increase in athletic ability reduces all-cause mortality and cardiovascular mortality rates by 15% and 13%, respectively (Kodama et al., 2009). Numerous studies have suggested that CRF and CVD are related to all-cause mortality and cancer mortality (Blair et al., 1989; Laukkanen et al., 2004; Sui, LaMonte & Blair, 2007; Sawada et al., 2014; Sui, LaMonte & Blair, 2007). A recent meta-analysis reported CRF to be a predictor of the risk of sudden cardiac death (Jiménez-Pavón, Lavie & Blair, 2019). Therefore, CRF assessment is crucial for diagnosing health risks and assessing interventions.
CRF can be measured using the respiratory data of exercising participants. Specifically, these data are used to calculate maximal oxygen uptake (V̇O2 max), the gold standard for CRF measurement; in the measurement, participants either run on a treadmill or use an ergometer at an exercise intensity that increases progressively until a given maximum is reached. Although submaximal exercise models and nonexercise models (without an exercise test) are alternatives for estimating V̇O2 max in measuring CRF (Abut, Akay & George, 2016), the direct measurement of V̇O2 max yields more objective and accurate results. However, such measurement is inconvenient because it requires expensive equipment and well-trained experimenters. In addition, participants perceive such measurement tests to be exhausting, time-consuming, and relatively risky and are thus less willing to participate. Accordingly, researchers have developed various submaximal exercise tests to indirectly estimate V̇O2 max; moreover, retrospective studies conducted by the American Heart Association have demonstrated that CRF indicators, whether directly measured or indirectly estimated, are robust indicators of health (Ross et al., 2016).
Submaximal exercise is a common method for estimating V̇O2 max, particularly in epidemiological research and large-scale physical fitness testing that involve numerous participants. The field tests in these measurement procedures include running, shuttle running, and the step test, with the step test being the most common method for evaluating cardiovascular function (Grant, Joseph & Campagna, 1999). In particular, the YMCA step test is widely used to predict V̇O2 max (Beutner et al., 2015). Currently, the Sports Administration of Taiwan’s Ministry of Education uses the 3-min Harvard step test for its National Physical Fitness and Cardiovascular Test. Specifically, three heart rate measurements are used to calculate the step-up index. However, previous studies have reported considerable differences in the validity of using the step test index to evaluate V̇O2 max, with the corresponding correlation coefficient (R) being 0.35–0.94 (Buckley et al., 2004; Chang & Lin, 1995; Mazic et al., 2001; Su, Lin & Hsieh, 2006; Chang & Lin, 1995; Yoopat, Vanwonterghem & Louhevaara, 2002). Furthermore, step tests require the use of step-up boxes, and the overall test time must be at least 6 min to allow for heart rate recovery. Participants who are less physically fit or who have knee conditions may find it difficult to complete the test and may also fall in the process of going up and down the stairs. A team of Japanese researchers developed a new 3-min walking test (Cao et al., 2013). Specifically, their main evaluation criteria comprised participant characteristics such as age, sex, and BMI as well as participants’ RPE during exercise. These criteria were determined to be effective predictors of V̇O2 max, and participants thought that this method was quicker and easier.
Tests of general CRF are crucial to the clinical evaluation of CVD. Additionally, the advantages and disadvantages, such as venue size, participant willingness, and the instruments, of various past field tests should be considered during the formulation of new methods, as done in the present study. Accordingly, we conducted the present study with the aim of developing a rapid, convenient, and low-risk model that can predict V̇O2 max in Taiwanese adults. Additionally, our model accords with the principle that physical exercise ought to be progressive. We investigated the feasibility of using a 3-min progressive knee-ups and step (3MPKS) test to predict V̇O2 max.
Materials and Methods
Prospective participants were excluded if they (1) had cardiovascular, pulmonary, or metabolic diseases; (2) had neurological, muscular, or skeletal disorders that affected their athletic ability; (3) had other health conditions that made them unsuited for moderate or intense exercise; or (4) were taking medications that could affect the outcome of this study. In total, among 200 participants recruited for this experiment, 166 completed the test. The data of the 166 participants were included in the analysis (age: 20–64 years; 65 men, 101 women). Among the 34 participants excluded, one participant withdrew from the experiment after experiencing suspected symptoms of arrhythmia during exercise; 11 were excluded because they failed to complete the step test within the requisite time (3 min); 12 were excluded because they could not attain the requisite step frequency and knee height for 20 consecutive seconds; nine were excluded because they had missing or improperly measured heart rate data; and one was excluded for having a “0” in their heart rate data. All participants signed an informed consent form after understanding their rights, the risks when participating in this study, and the purpose and method of our research. Our research plan was approved by the Institutional Review Boards (IRBs) of the Industrial Technology Research Institute and of Taipei Medical University (IRB No: N201808055). Participant characteristics are detailed in Table 1.
|Total||Training dataset||Testing dataset|
|Age(years)||41.9 ± 9.6||42.2 ± 9.4||40.8 ± 10.2|
|Height||164.83 ± 8.35||164.33 ± 8.07||166.30 ± 9.06|
|Weight||65.63 ± 13.60||65.22 ± 14.08||66.85 ± 12.15|
|Body fat (%)||27.81 ± 7.92||27.79 ± 7.65||27.86 ± 8.75|
|V̇O2 max (ml kg−1 min−1)||34.45 ± 8.69||34.06 ± 8.14||35.61 ± 10.15|
|HR0||86.04 ± 12.78||86.04 ± 12.99||86.02 ± 12.29|
|ΔHR3- HR0||71.00 ± 13.24||71.10 ± 13.41||70.69 ± 12.87|
|ΔHR3-HR4||14.64 ± 13.72||14.14 ± 13.94||16.65 ± 14.09|
Data are presented as mean ± standard deviation.
heart rate at the beginning
difference between third minute heart rate and beginning heart rate
difference between third minute and fourth minute heart rates
The anthropometric and body composition measures were height, weight, and body fat percentage (BF%). BF% was measured using bioelectrical impedance analysis (InBody 720, Biospace, USA; McLester et al., 2020), and body mass index (BMI, in kg/m2) was calculated as the quotient that is weight (in kilograms) divided by the squared height (in meters).
We conducted two exercise tests in a counterbalanced design. The second test was conducted exactly 1 week after the first and at the same time of the day to ensure that the participants recovered adequately from the first exercise. The participants underwent 5–10 min of dynamic warm-up prior to both exercise tests; to mitigate extraneous influence on the results, the participants were also asked not to engage in moderate or intense exercise 48 h before both exercise tests.
To measure the V̇O2 max of the participants, we used a bicycle ergometer (839E, Monark, Varberg, Sweden) for a maximal graded exercise test. After participants sat still for 2 min, they sat on the stationary bicycle and started cycling at the speed of 70 ± 10 rpm. The participants began the exercise with a 2-min warm-up at 25 W loading, where the loading was increased by 15 W every 2 min. The testing was terminated when the participants could no longer continue the exercise due to bradypnea or fatigue, although the bicycle speed was maintained at 70 rpm. Subsequently, the participants rested for 3 min at a loading of 0 W (no resistance). Throughout the exercise testing, the participants wore a watch to monitor their heart rate and a mask to monitor their breathing. Breath-by-breath analysis was conducted on the participant data through a cardiopulmonary testing system (MetaMax 3B, Cortex, Germany). V̇O2 max was defined as the maximum average oxygen uptake for 20 consecutive seconds. To ensure that every participant reached V̇O2 max, we defined V̇O2 max as being reached if two of the three following conditions were met: (1) V̇O2 plateaus with increases in work rate; (2) the maximum respiratory exchange ratio is ≥1.10; and (3) 90% of the expected maximal heart rate, obtained by subtracting the participant’s age from 220, is reached (American College of Sports Medicine, 2009). Nearly all participants satisfied the criteria for an acceptable V̇O2 max, with only one participant excluded from the V̇O2 max test due to suspected symptoms of arrhythmia observed in the step test.
Prior to the 3MPKS test, the participants wore a sports watch with heart rate (Polar V800, USA) and stride sensors (Polar S3 BlueTooth Stride Sensor, USA). The heart rate sensor was placed at the center of each participant’s chest using a heart rate belt (Polar H10), and the step sensor was fixed on a pair of shoes, with shoelaces, to monitor their heartbeat and number of steps taken. After the devices were worn, we measured the midpoint of the line connecting the anterior epicondyle to the midpoint of the sacrum. We marked the midpoint on the wall using colored tape as a reference for the height at which the knee should be lifted to when stepping. After the test started, the participants followed the appropriate rhythm and were required to lift their knee to the marked height at each step. The participants began the test at a pace of 80 spm (steps per minute), which increased by 16 spm every 30 s in six stages. The participants walked in stages 1 to 4 and had to perform stationary running in stages 5 and 6 (Fig. 1). We stopped the exercise if the participants could not achieve the requisite knee height or rhythm for 30 s. For their safety, the participants were asked to relax at a step rate of 80 spm in the first 30 s before resting in a standing position. We recorded the participants’ heart rate during the exercise, at the end of the exercise, and 1 min after the end of the exercise. Thirty-four participants were excluded because (1) their heart rate data were missing, (2) their heart rate was 0, (3) they did not maintain the requisite step frequency or knee height for 20 consecutive seconds, (4) they failed to complete the step test within the requisite duration, and (5) they were suspected of having heart arrhythmia. Potential predictor variables for the results of the 3MPKS test were based on per-second heart rate data collected during the test. The data included heart rate at the beginning as well as at the first, second, third, and fourth minutes, denoted by HR0, HR1, HR2, HR3, and HR4, respectively, and were used for subsequent analysis.
To construct and subsequently evaluate a model for estimating relative oxygen uptake, we divided the full sample set (n = 166) into a 75% training sample set (n = 124) and 25% test sample set through simple random sampling. We analyzed the descriptive statistics for the main parameters, for the whole sample, and for the two subsamples.
Development of prediction model
Using Pearson correlation coefficients, we examined the relationship between the predicted and actual relative oxygen uptakes. Multiple regression analysis was used to construct a method for selecting which variables to include in the model for predicting relative oxygen uptake. Through a backward-selection regression approach, the initial model included all possible predictors, including sex (men = 1, women = 0), age, BMI, BF%, HR0, HR1, HR2, HR3, HR4, △HR0 − HR1, △HR1 − HR2, △HR2 − HR3, △HR3 − HR0, and △HR3 − HR4. Additionally, we constructed a BMI model and BF% model to predict body composition. The goodness of fit and precision of the regression equations were evaluated using the multiple coefficient of determination (R2), absolute standard error of estimate (SEE), and relative SEE (%SEE).
To construct an accurate regression model, the regression assumptions were verified. We conducted a Kolmogorov–Smirnov test to examine the normality of the residuals, and we calculated the variation inflation factor (VIF) to check for multicollinearity.
All statistical analyses were performed using SPSS version 20 (IBM, USA). Statistical significance was indicated by an alpha level of 0.05.
The 166 participants had an average age of 41.9 ±9.6 years (range: 22–64 years), and 40% of them were men. Their mean relative oxygen uptake was 34.45 ±8.69 mL/kg/min. The training sample and test sample did not differ significantly with respect to their parameter values (p > 0.05) Table 1.
The test–retest reliability of the 3MPKS test, as evaluated in our laboratory, was excellent: the intraclass correlation coefficient (ICC) was 0.88 (95% confidence interval [CI]: 0.77–0.94), and 60 Taiwanese adults tested 1 week apart participated in this evaluation. In general, good, moderate, and poor reliability levels are indicated by ICC values of >0.75, 0.5–0.75, and <0.5, respectively.
According to the correlation matrix, V̇O2 max had the strongest correlation with BF% among all variables (R = −0.662; training data set, n = 124). In addition, V̇O2 max was significantly correlated with the heart rate parameters (HR0, HR2, HR3, and HR4), whose data were collected in the step test. V̇O2 max was most and least correlated with HR4 (R = −0.442) and HR3 (R = −0.289), respectively. Despite the high correlation between V̇O2 max and the heart rate parameters at different stages, the heart rates of the participants were expected to increase continuously from the first to third minutes of stepping, if performed properly. An individual’s heart rate typically reaches its peak immediately after exercise, and it either decreases at 1 min after exercise or does not decrease at all depending on whether the individual recovers quickly or poorly. Because heart rate is dynamic, to establish a regression model, we used combinations of heart rate parameters and adopted the difference between predicted and measured heart rate data at each stage as inputs (Table 2).
The results of our other cross-validation analyses are presented in terms of CE (Constant error) values. The absolute CE values for subgroups stratified by sex and age were <1.00 for the two models (both in training and testing data sets, n = 124 and 42). Regarding the subgroups stratified by V̇O2 max, the CE values were negative in low-fitness, middle-fitness subgroups in training data set and low-fitness in testing data set. On the other hand, the CE values were positive in high-fitness in all two data sets (Table 3).
|Subgroup||n(%)||BF% model(%)||BMI model(kg m−2)|
|Training set(n = 124)|
|Testing set(n = 42)|
Figures 2 and 3 present the Bland–Altman plots produced by the BF% and BMI models based on the testing data set (n = 42). As evident in the plots, the differences between the predicted and measured data were within an acceptable range. The mean error of the BF% model was −0.36 mL/kg/min (95% CI [−12.38–11.98]). For the BMI model, the mean error was 0.4 mL/kg/min (95% CI [−12.35–13.58]). In the BF% and BMI models, the errors for three and two participants, respectively, fell outside the 95% CI.
We constructed a model to predict relative oxygen uptake by using multiple regression analysis. The parameters selected for the BF% model were sex, age, BF%, HR0, ΔHR3 − HR0, and ΔHR3 − HR4; R2 = 0.624 and SEE = 4.982 (training data set, n = 124) (Fig. 4). The parameters selected for the BMI model were sex, age, BMI, initial heart rate, ΔHR3 − HR0, and ΔHR3 − HR4; R2 = 0.567 and SEE = 5.153 (training data set, n = 124) (Fig. 5). We used BF% as a predictor of body composition; it is more accurate relative to BMI, which is calculated using only height and weight (Table 4). Table 4 presents the cross-validation results for the predicted residual sum of squares (PRESS) statistics (R2p = 0.64 and SEE p = 4.84), which demonstrated minimal shrinkage in the accuracy of the regression model.
All regression assumptions were satisfied in our V̇O2 max prediction models. Specifically, the Kolmogorov–Smirnov test indicated normality in the residuals (p > 0.05). No pattern was determined in the scatter plot between the residuals and predicted V̇O2 max. Multicollinearity was absent among the predictor variables: the VIF ranges for the BF% and BMI models were 1.09–1.49 and 1.10–1.40, respectively; multicollinearity is absent if VIF ≤ 10 (O’brien, 2007).
This study developed a practical and easy-to-use model for predicting V̇O2 max in Taiwanese people. We recruited 166 Taiwanese adults and constructed and then evaluated a prediction model. Our results suggest that age, sex, and BF% as well as heart rate during the step test are excellent predictors of V̇O2 max. We also developed a novel 3MPKS test.
|BF% model (%)||BMI model (kg m−2)|
|V̇O2max (ml kg−1 min −1)||Coefficients||β||p value||Coefficients||β||p value|
|Sex (0=women, 1=men)||4.366||0.258||.000||9.338||0.551||.000|
body mass index
body fat percentage
standardized regression weights
standard error of estimate
SEE / mean of measured V̇O2 max ×100.
predicted residual error sum of squares
PRESS standard error of estimate
PRESS squared multiple correlation coefficient
Nes et al. (2011) conducted large-scale V̇O2 max tests on 4,260 participants. They developed a nonexercise model and determined four variables (age, waist circumference, physical activity, and resting heart rate) to be excellent predictors of V̇O2 max; for their model, R2 was 0.61 and SEE was 5.70 mL/kg/min for men, and R2 was 0.56 and SEE was 5.14 mL/kg/min for women. Jackson et al. (2012) conducted a 27-year study that examined the V̇O2 max of 11,365 people and used variables such as age, sex, BMI, waist circumference, resting heart rate, physical activity, and smoking habits to estimate CRF; for their model, R was 0.78–0.81 and SEE was 5.3–5.6 mL/kg/min. Although the nonexercise model is an excellent predictor of V̇O2 max, its SEE is generally higher than those of submaximal exercise models; compared with nonexercise models, our developed BF% model had better predictive performance and a lower standard error of estimate (R2 = 0.624 and SEE = 4.982). Abut, Akay & George (2016) reported that (1) when perceived functional ability (PFA) was used as the sole predictor of V̇O2 max, an R value of 0.73 and a higher RMSE of 6.08 mL/kg/min could be obtained; (2) when submaximal ending speed (SM-ES) of a treadmill was used as the sole predictor, the R value increased to 0.82 and the RMSE was relatively low at 4.99 mL/kg/min; and (3) when both PFA and SM-ES were used as predictors, the R value was 0.89 and RMSE was 4.14 mL/kg/min. These findings indicate that predicted values of V̇O2 max that are based only on participant self-reports are likely to deviate from their measured values. Although predictive performance is ostensibly improved when motion is added to the prediction model, the cost of exercise tests due to the use of this method restricts its application in large-scale tests.
Several studies have developed simple models involving submaximal motion. Lee et al. (2019) investigated 568 adults and used sex, age, height, and weight and inverse recovery heart rate during a YMCA step test to predict V̇O2 max; for their model, R was 0.78 and SEE was 4.74 mL/kg/min. The duration of their exercise test plus recovery time was only 4 min, and they used exercise-induced heart rate as a predictor; their results are similar to ours. Their study provided a simple and practical method for simultaneously estimating CRF in many Korean adults. Cao et al. (2013) used age, sex, and physical composition as well as stepping distance over a 3-min period to develop a set of prediction methods. They determined that BF% (a measure of body composition) was a better predictor than BMI (R2 = 0.83 vs. 0.80, SEE = 4.565 vs. 5.037 mL/kg/min). In contrast to our method, their method has the considerable advantages of a shorter testing time of 3 min and the fact that participants need not wear a heart rate monitor. However, their test is limited by its need for a 20-m open space. Similarly, we found that sex, age, and BF% as well as heart rate during the 3MPKS test yielded the best prediction performance (R = 0.79, SEE = 4.982 mL/kg/min). Because BMI is based on only height and weight and may not accurately represent the body characteristics of participants, BMI is a less accurate predictor than BF%.
Most submaximal exercise models proposed by previous studies involve a fixed-height step test. However, the height and leg length of participants when standing may affect their physiological response in the step test (Culpepper & Francis, 1987). Relative to their European counterparts, Asian adults have shorter heights and leg lengths when standing (Stanfield et al., 2012). Therefore, differences in heart rate and oxygen consumption potentially affect the model’s prediction. The 3MPKS test employs the knee-ups and step test to measure the physical fitness and cardiopulmonary endurance of older adults (Rikli & Jones, 2001). In the test, participants must execute tasks at various knee heights based on their thigh length, and individualized exercise testing goals are provided. Moreover, most field tests involve average speed tests, such as step tests and running. In running tests specifically, if the distance is used as the capacity index but the speed or frequency of exercise is not progressively increased, participants may exercise intensely at the beginning of the test (i.e., run at a higher speed). However, due to the lack of appropriate speed allocation, decremental loading occurs in participants as their physical strength decreases. The difficulty of diagnosing potential heart diseases in advance increases the risk of sudden death during running tests. To the best of our knowledge, research has not been conducted on the ethics of running tests. Most previous studies have investigated the rate of sudden death among athletes in long-distance competitions. However, cases of sudden cardiac death occur frequently worldwide during running tests, and the principle that physical activities ought to be progressive must be adhered to in physical fitness tests. Our research method used body composition and heart rate as variables. The advantages of the 3MPKS test are that it does not require a step-up box and is not subject to venue restrictions. These make the 3MPKS test accord with the principle that physical activities ought to be progressive, thus making it safer.
Considering the immediacy of heart rate measurement and that of confounding factors, we used a chest-worn heart rate monitor in the experiment. Although the requirement of heart rate monitoring constitutes a disadvantage for the 3MPKS test, it is ameliorated by the prevalence of low-cost wearable devices. More comfortable than the chest-worn heart rate belt, products that combine running clothes with heart rate belts have also appeared on the market. Research has also suggested a high correlation between the heart rate measurements of various types of optical devices and chest-worn heart rate belts (Stahl et al., 2016). Therefore, when conducting a large-scale cardiorespiratory general test, the use of easily wearable optical heart rate monitors can be considered. The whole-range monitoring of heart rate can also considerably improve test safety in a field study. Notably, through whole-range monitoring, we found that one research participant was likely to have an unknown heart disease. We then terminated the experiment for the participant and recommended that the participant seek medical treatment. This example illustrates a side benefit of CRF tests.
In our research model, heart rate during stepping at each stage was used as the main variable. Therefore, the test may be unsuitable for individuals who have psychological sensitivity or dysautonomia or who are taking medication. Furthermore, because our participants were adults between 20 and 64 years old, it was unclear whether our 3MKPS test is appropriate as a physical fitness and cardiorespiratory test for students (7–23 years old) and older adults (≥65 years old). Future research must include samples with greater diversity in age and ethnicity to assess whether our 3MKPS test can be applied to the wider global population.
This study, involving Taiwanese adults, constructed and verified a model for predicting V̇O2 max, which is used to measure CRF. This model comprises the predictors sex, age, and body composition as well as heart rate changes during a step test. Our 3MKPS test has three advantages: it has a short testing time of 4 min, it has no venue limitations, and it does not require a step box. Furthermore, measurements can be taken for many participants simultaneously by asking them to wear a heart rate monitor and move according to a beat. Our model can also be applied to large-scale epidemiological research. In future applications, the model can be combined with smartwatches or used to develop health and well-being apps, helping users to track their V̇O2 max. Future research can further explore the correlation between various diseases and V̇O2 max, as predicted using our simple and reliable method for measuring CRF.