Machine learning and statistical methods for the prediction of maximal oxygen uptake: recent advances

Abstract

Maximal oxygen uptake (VO₂max) indicates how many milliliters of oxygen the body can consume in a state of intense exercise per minute. VO₂max plays an important role in both sport and medical sciences for different purposes, such as indicating the endurance capacity of athletes or serving as a metric in estimating the disease risk of a person. In general, the direct measurement of VO₂max provides the most accurate assessment of aerobic power. However, despite a high level of accuracy, practical limitations associated with the direct measurement of VO₂max, such as the requirement of expensive and sophisticated laboratory equipment or trained staff, have led to the development of various regression models for predicting VO₂max. Consequently, a lot of studies have been conducted in the last years to predict VO₂max of various target audiences, ranging from soccer athletes, nonexpert swimmers, cross-country skiers to healthy-fit adults, teenagers, and children. Numerous prediction models have been developed using different sets of predictor variables and a variety of machine learning and statistical methods, including support vector machine, multilayer perceptron, general regression neural network, and multiple linear regression. The purpose of this study is to give a detailed overview about the data-driven modeling studies for the prediction of VO₂max conducted in recent years and to compare the performance of various VO₂max prediction models reported in related literature in terms of two well-known metrics, namely, multiple correlation coefficient (R) and standard error of estimate. The survey results reveal that with respect to regression methods used to develop prediction models, support vector machine, in general, shows better performance than other methods, whereas multiple linear regression exhibits the worst performance.

Keywords:

• machine learning methods
• maximal oxygen consumption
• prediction models
• feature selection

Introduction

Maximal oxygen uptake (VO₂max) refers to the maximal amount of oxygen that an individual can utilize during intense or maximal exercise. VO₂max is expressed either as an absolute rate in liters of oxygen per minute (L min⁻¹) or as a relative rate in milliliters of oxygen per kilogram of body mass per minute (mL kg⁻¹ min⁻¹). VO₂max plays a significant role in both sport and medical sciences for different purposes. In sport sciences, it is often used as an indicator for endurance capacity of athletes, representing the upper limit of their aerobic fitness. Knowledge of VO₂max can assist coaches in designing efficient training programs to maximize performance, minimize injury and fatigue of athletes, too.^1,2 In medical sciences, it can serve as a metric to estimate the disease risk of a person, suggesting an unusually large amount of cholesterol, body fat, and blood pressure.³ Consequently, an accurate knowledge of VO₂max has gained an increasing attention over the last decades.

The direct measurement of VO₂max during a maximal graded exercise test (GXT) is accepted as the most accurate method for the assessment of aerobic power. However, despite a high level of accuracy, determining VO₂max during maximal GXTs also entails a number of limitations. For example, conducting such tests require expensive and elaborated laboratory equipment, such as oxygen and carbon dioxide gas analyzers, an expiratory air flow probe, an air mixing chamber, a dehumidifier, a vacuum pump, and a data acquisition system. Trained staff is needed in order to perform the VO₂max measurements and interpret the test results. In addition, maximal GXTs are unappealing to some individuals because the test requires exhausting exercise to the point of volitional exhaustion. Because of this, some older or higher risk individuals should not perform the test without medical supervision. Also, when the number of subjects is large, it is not feasible to apply VO₂max tests for all subjects.⁴

The practical limitations of direct measurement have led to the development of various regression models for predicting VO₂max that are principally based on one of the three different tests as follows: exercise tests, nonexercise tests, and hybrid tests.

Exercise tests in turn comprise maximal and submaximal tests. Maximal tests are expensive to administer and can be risky because the subject reaches the point of exhaustion in terms of heart rate (HR). The requirement of expensive gas analysis, ventilation equipment, and the need for medical attendance are the other major disadvantages of maximal tests.⁵ Because of these drawbacks of maximal tests, other methods that do not require maximal effort have been developed to determine VO₂max. These methods use either submaximal exercise tests or nonexercise tests to predict VO₂max.

Submaximal tests indirectly estimate VO₂max using submaximal exercise variables and are usually performed on a treadmill, ergometer, or track. Submaximal tests have certain advantages over maximal tests, such as the specialized laboratory equipment is unnecessary, the test administrators require less training, and the exercise intensity is realistic for most participants. In addition, periodic submaximal testing provides a convenient way to monitor progress throughout an exercise program and educates participants about the selection of an appropriate intensity of exercise to develop individualized exercise programs. Submaximal tests are faster, safer, and cheaper to administer than maximal tests.

Nonexercise prediction equations enable an appropriate prediction of VO₂max without the need to conduct a maximal or submaximal exercise test. Researchers may prefer nonexercise models instead of maximal and submaximal exercise tests due to their simplicity of administration to large populations, their independence from expensive laboratory equipment, and exercise testing. The main disadvantage of nonexercise tests is that they rely on the truthful self-report of activity and cannot be utilized when the individuals have an interest in the results, giving them a reason to falsify the self-reported activity levels.^5–7

Representative examples for maximal variables include maximal heart rate (HRmax), rating of perceived exertion (RPE), and exercise time, whereas treadmill stage, treadmill speed, and exercise HR are popular examples for submaximal variables. On the other hand, questionnaire variables, such as perceived functional ability (PFA), physical activity rating (PAR), and activity code (AC), belong to the class of nonexercise variables. Figure 1 illustrates an overview of various types of VO₂max prediction models along with their representative examples of predictor variables.

Figure 1 Overview of various types of VO₂max prediction models.

Abbreviations: HRmax, maximal heart rate; RER, respiratory exchange ratio; PFA, perceived functional ability; PAR, physical activity rating; AC, activity code; VO₂max, maximal oxygen uptake.

The purpose of this paper is to give a detailed overview about the data-driven modeling studies, such as linear regression models and artificial neural network models, for the prediction of VO₂max conducted in recent years. Also, the performance and prediction accuracies of such various regression models reported in related literature are compared with each other in terms of multiple correlation coefficient (R) and standard error of estimate (SEE) in order to identify the best set of predictor variables along with the used machine learning and statistical methods. The R and SEE values are the most popular and frequently used metrics in the field of sport physiology, and also vast majority of the studies related to the prediction of VO₂max use R and SEE for performance and accuracy evaluations of prediction models.

The rest of the paper is structured as follows. The “Prediction models based on exercise tests” section presents the prediction models based on exercise tests that have used either maximal or submaximal tests for the prediction of VO₂max. The “Nonexercise prediction models” section introduces the nonexercise prediction models. The “Hybrid models” section gives details about the hybrid prediction models, which combine data obtained from both exercise and nonexercise tests. The “Overview of prediction methods and model evaluation techniques” section gives a brief overview of the machine learning and statistical methods as well as model evaluation techniques encountered within the scope of this study. Finally, the “Conclusion and future work” section concludes the paper and gives pointers to future work.

Prediction models based on exercise tests

Maximal prediction models

Chatterjee et al⁸ validated the relevance of the 20 m multistage shuttle run test (20 m MST) in Nepalese adult females. Forty female students, aged from 20 years to 24 years, from diverse colleges of Nepal were selected for the study. Direct estimation of VO₂max included treadmill exercise followed by expired gas analysis by scholander microgas analyzer, whereas VO₂max was indirectly predicted by the 20 m MST. Maximal shuttle run velocity was used to predict the VO₂max with the help of multiple linear regression (MLR). The R and SEE values were reported to be 0.94 and 0.53 mL kg⁻¹ min⁻¹, respectively.

Bandyopadhyay⁹ examined the effectivity of 20 m MST for the prediction of VO₂max in male university students using MLR. Eighty-four stationary male university students of the same socioeconomic background were selected from students of University of Calcutta, Kolkata, India, to test the efficiency of 20 m shuttle run test (SRT) for indirect estimation of VO₂max. Predictor variables that are used to develop the linear regression model incorporate age, height, weight, and velocity. The R and SEE values were reported to be 0.94 and 1.38 mL kg⁻¹ min⁻¹, respectively.

Mahar et al¹⁰ created and cross-validated regression models to predict VO₂max from progressive aerobic cardiovascular endurance run (PACER) 20 m shuttle run performance in 10–16-year-old boys and girls. PACER performance and VO₂max were evaluated using a dataset composed of 244 participants. The dataset was randomly separated in validation (n=174) and cross-validation (n=70) samples. The validation sample was utilized to build the regression models for predicting VO₂max from sex, age, body mass index (BMI), and PACER. Values of R and SEE for the best model were found as 0.75 and 6.17 mL kg⁻¹ min⁻¹, respectively.

Akay et al¹¹ developed VO₂max models by using different regression methods, such as support vector machine (SVM), multilayer perceptron (MLP), generalized regression neural network (GRNN), and MLR. The dataset included data of 439 subjects (211 males and 228 females) who are apparently healthy volunteers ranging in age from 20 years to 79 years. The predictor variables used to develop the models were sex, age, BMI, body fat percentage (BF%), respiratory exchange ratio (RER), RPE, HRmax, and time to exhaustion from treadmill test. For model testing, tenfold cross-validation has been conducted. The most accurate R and SEE values were obtained by MLP with 0.91 and 3.73 mL kg⁻¹ min⁻¹, respectively. In a follow-up study, Akay et al¹² combined SVM and MLR methods with the Relief-F feature selection algorithm to identify the discriminative predictors of VO₂max using the same dataset. The models included the predictor variables sex, age, BMI, HRmax, RPE, RER, and exercise time. According to the results obtained, the model made up of five predictor variables including sex, age, BMI, HRmax, and exercise time performed better than the model including the full set of predictor variables, independent of whether it was evaluated with SVM or MLR. SVM-based models, in turn, outperformed MLR-based models. The values of R and SEE for the most accurate SVM-based model were reported as 0.90 and 4.58 mL kg⁻¹ min⁻¹, respectively.

Silva et al¹³ developed and tested two models to predict VO₂max in Portuguese youths, aged from 10 years to 18 years, using a 20 m SRT. Participants (60 males and 54 females) were divided into estimation and cross-validation groups, and their VO₂max was directly measured by wearing a portable gas analyzer during the SRT. MLP and MLR methods were applied by using sex, age, height, weight, BMI, and SR stage as predictors of VO₂max. The reported R and SEE values were 0.84 and 4.90 mL kg⁻¹ min⁻¹, respectively.

Daros et al¹⁴ developed a maximum aerobic power test and prediction equations to estimate VO₂max for soccer athletes using MLR. The proposed test composed of applying progressive, persistent, and maximal speed running that covers 80 m, organized in a square (20 m ×20 m), where the athletes performed until they became exhausted. The dataset was made up of 24 young soccer players aged 16.66±1.49 years, each of them was a practitioner of the sport for at least 4 years. Two VO₂max estimation equations were developed that incorporated total distance and maximum speed as predictor variables. The highest R and the lowest SEE were found as 0.76 and 4.29 mL kg⁻¹ min⁻¹, respectively.

Costa et al¹⁵ developed equations to predict the VO₂max of nonexpert adult swimmers. Participants were 22 male nonexpert swimmers (aged between 18 years and 30 years). A progressive swim test was conducted, in which there occurs an increase in the intensity of the swim every two laps. Three different MLR-based equations were developed using the predictor variables, such as weight, number of laps performed, and after HR. According to the performances obtained by the three models, best values of R and SEE were determined as 0.80 and 7.20 mL kg⁻¹ min⁻¹, respectively.

Machado and Denadai¹⁶ presented a specific predictive equation to predict the VO₂max from boys aged between 10 years and 16 years, with the help of MLR. Forty-two boys conducted a treadmill running ergospirometric test, with the introductory speed adjusted at 9 km h⁻¹, until voluntary exhaustion. The predictor variables weight and velocity have been utilized to establish the prediction equation. R and SEE values of the equation were reported as 0.93 and 4.10 mL kg⁻¹ min⁻¹, respectively.

Aktürk and Akay¹⁷ employed MLP combined with Relief-F feature selection algorithm to develop VO₂max prediction models for healthy subjects ranging in age from 18 years to 65 years. The dataset comprised 100 subjects (50 males and 50 females) and included the predictor variables, such as sex, age, BMI, HRmax, RPE, RER, and grade. For model testing, tenfold cross-validation was used. The results showed that the model, including six predictor variables (sex, age, BMI, HRmax, RER, and grade), yielded the highest R and the lowest SEE with 0.86 and 4.62 mL kg⁻¹ min⁻¹, respectively. The model including the full set of predictor variables, on the other hand, yielded the second highest R and lowest SEE with 0.86 and 4.65 mL kg⁻¹ min⁻¹, respectively.

Akay et al¹⁸ predicted VO₂max of cross-country skiers using different regression methods, including SVM using the radial basis function, linear SVM and MLP. The dataset contained information pertaining to 175 cross-country skiers (101 males and 74 females) aged between 15 years and 30 years. Predictor variables that were used to build the models were sex, age, BMI, HRmax, HR at lactate threshold, exercise time, and protocol. Tenfold cross-validation was applied on the dataset. The values of R and SEE for the most accurate model were attained by SVM using the radial basis function with 0.83 and 4.39 mL kg⁻¹ min⁻¹, respectively.

Table 1 gives an overview of recent studies in literature that developed maximal models for the prediction of VO₂max. The studies are sorted in chronological order. For each study, only the R and SEE values for the most accurate model along with its predictor variables are reported. For the case where more than one regression method is utilized in a study, the best values of R and SEE for a prediction model are reported for every regression method. The R values are the averages of cross-validation testing data, where cross-validation has been applied. In all other cases, R values belong to the testing data.

Table 1 Summary of recent studies in literature that developed maximal models for the prediction of VO₂max

Study	Year	Method	Predictor variables	R	SEE (mL kg^−1 min^−1)
Chatterjee et al^Citation8	2010	MLR	Velocity	0.94	0.53
Bandyopadhyay^Citation9	2011	MLR	Age, height, weight, velocity	0.93a	1.38
Mahar et al^Citation10	2011	MLR	Sex, age, BMI, PACER	0.74a	6.17
Akay et al^Citation11	2012	SVM	Sex, age, BMI, BF%, RER, RPE, HRmax, time	0.92	4.09
Akay et al^Citation11	2012	MLP	Sex, age, BMI, BF%, RER, RPE, time	0.91	3.73
Akay et al^Citation11	2012	GRNN	Sex, age, BMI, BF%, RER, time	0.82	5.03
Akay et al^Citation11	2012	MLR	Sex, age, BMI, BF%, RER, RPE, HRmax, time	0.81a	4.54
Silva et al^Citation13	2012	MLP	Sex, age, height, weight, BMI, SR stage	0.79	5.60
Silva et al^Citation13	2012	MLR	Sex, age, height, weight, BMI, SR stage	0.68a	4.90
Daros et al^Citation14	2012	MLR	Meter	0.76	4.29
Veronese da Costa et al^Citation15	2013	MLR	Weight, NLP, AHR	0.76a	7.21
Machado and Denadai^Citation16	2013	MLR	Weight, velocity	0.92a	4.10
Aktürk and Akay^Citation17	2014	MLP combined with Relief-F	Sex, age, BMI, HRmax, RER, grade	0.86	4.63
Akay et al^Citation12	2014	SVM combined with Relief-F	Sex, age, BMI, HRmax, time	0.90	4.58
Akay et al^Citation12	2014	MLR combined with Relief-F	Sex, age, BMI, HRmax, time	0.88a	4.77
Akay et al^Citation18	2014	SVM	Sex, age, BMI, HRmax, Pr	0.83	4.38
Akay et al^Citation18	2014	MLP	Sex, age, BMI, HRmax, Pr	0.81	4.64
Akay et al^Citation18	2014	MLR	Sex, age, BMI, HRmax, Pr	0.78a	4.73

Note:

a Represents adjusted R (R_adj).

Abbreviations: VO₂max, maximal oxygen uptake; SEE, standard error of estimate; MLR, multiple linear regression; BMI, body mass index; PACER, progressive aerobic cardiovascular endurance run; SVM, support vector machine; BF%, body fat percentage; RER, respiratory exchange ratio from treadmill test; RPE, self-reported rating of perceived exertion from treadmill test; HRmax, maximal heart rate; MLP, multilayer perceptron; GRNN, general regression neural network; SR stage, shuttle run stage; NLP, number of laps performed; AHR, after heart rate; Pr, protocol.

Submaximal prediction models

Coquart et al¹⁹ presented an indirect method for predicting VO₂max from a submaximal GXT in obese women by using MLR. Thirty obese women performed maximal GXT to volitional exhaustion to determine the reference value for VO₂max. During GXT, the power at RPE 15 was measured, too. The predictor variables, age and power at RPE 15 elicited during a submaximal GXT, were found to produce a reasonably accurate prediction of VO₂max in obese women. The values of R and SEE were reported as 0.83 and 0.16 mL kg⁻¹ min⁻¹, respectively.

Akay et al²⁰ developed an accurate MLP model to predict VO₂max of fit adults from a single-stage submaximal treadmill jogging test. For the study, 126 participants (81 males and 45 females) aged from 17 years to 40 years successfully completed a maximal GXT to determine the VO₂max. The variables, such as sex, age, weight, HR, and jogging speed, were used to build the MLP-prediction model. The highest R and the lowest SEE values were reported as 0.95 and 1.80 mL kg⁻¹ min⁻¹, respectively.

Açıkkar et al²¹ built MLP models to predict VO₂max from submaximal endurance exercise at varying distances (0.5 mile, 1.0 mile, and 1.5 mile), involving walking, jogging, or running. In total, 185 (115 males and 70 females) “apparently healthy” college students, aged 18–26 years, successfully completed a submaximal 1.5-mile endurance test and a maximal GXT. Predictor variables utilized to create the prediction models involved sex, age, BMI, the elapsed exercise times (MIN1, MIN2, and MIN3), and HRs (HR1, HR2, and HR3) at the 0.5 mile, 1 mile, and 1.5 mile distances, respectively. Using tenfold cross-validation on the dataset, values of R and SEE for the most accurate prediction model were reported as 0.89 and 2.22 mL kg⁻¹ min⁻¹, respectively. In a follow-up study, Akay et al²² developed various VO₂max prediction models based on SVM using the same dataset. The predictor variables employed to build the models remained the same, with the exception of BMI, instead of which the height and weight variables have been utilized. For comparison purposes, VO₂max prediction models based on MLP and MLR were also developed. In conclusion, it was observed that SVM-based models showed better performance than the ones obtained by using MLP and MLR-based prediction models. The values of R and SEE for the most accurate SVM-based model were reported as 0.87 and 2.90 mL kg⁻¹ min⁻¹, respectively.

Billinger et al²³ inspected the capability of the Young Men’s Christian Association submaximal exercise test protocol using a total body recumbent stepper to predict VO₂max. One-hundred and ten 18–60-year-olds were screened for the study by fulfilling the participation criteria of low-to-moderate cardiovascular disease risk. The subjects were requested to keep a pace of 100 stages every minute, and the resistance rose in 3 minutes intervals as specified by the protocol until the subjects became tired, or 85% of HRmax was attained. The predictor variables, sex, age, weight, HR, and work rate expressed in Watts, were utilized to form the MLR-based prediction model. A cross-validation study was conducted to validate the accuracy of the prediction equation. The presented R and SEE values were 0.92 and 4.09 mL kg⁻¹ min⁻¹, respectively.

Tönis et al²⁴ determined parameters that can be utilized for predicting VO₂max, from a combination of HR and accelerometer data. Data were obtained from 41 (23 males and 18 females) healthy subjects, ranging in age from 20 years to 29 years. The measurement protocol consisted of a submaximal single-stage treadmill walking test for VO₂max estimation followed by a walking test at two different speeds (4 km h⁻¹ and 5.5 km h⁻¹) for parameter determination. MLR analysis was performed using general subject measures such as sex, age, height, weight, BMI, and intercept and slope of the relation between HR and accelerometer output during walking as independent variables to predict the VO₂max. The R and SEE values were reported as 0.90 and 2.05 mL kg⁻¹ min⁻¹, respectively.

Cao et al²⁵ built new VO₂max prediction models using a perceptually regulated 3-minute walk distance (3MWD) test. VO₂max was measured with a maximal incremental cycle test in 283 Japanese adults. A 3MWD test was performed at a self-regulated intensity corresponding to RPE 13. Three prediction models were developed by using MLR for predicting VO₂max using data on sex, age, 3MWD, BMI, waist circumference, and BF%. The R and SEE values of the best model were reported as 0.84 and 4.57 mL kg⁻¹ min⁻¹, respectively.

Table 2 gives an overview of recent studies in literature that used submaximal models for the prediction of VO₂max in chronological order.

Table 2 Summary of recent studies in literature that developed submaximal models for the prediction of VO₂max

Study	Year	Method	Predictor variables	R	SEE (mL kg^−1 min^−1)
Coquart et al^Citation19	2010	MLR	Age, power	0.81a	0.16
Akay et al^Citation20	2011	MLP	Sex, age, weight, HR, jogging speed	0.95	1.80
Açıkkar et al^Citation21	2012	MLP	Sex, age, BMI, MIN3, HR3	0.89	2.22
Billinger et al^Citation23	2012	MLR	Sex, age, weight, WR, HR	0.91a	4.09
Tönis et al^Citation24	2012	MLR	Sex, age, height, weight, HR, AD	0.87a	2.05
Cao et al^Citation25	2013	MLR	Sex, age, 3MWD, BF%	0.83a	4.57
Akay et al^Citation22	2014	SVM	Sex, age, height, weight, MIN3, HR3	0.87	2.90
Akay et al^Citation22	2014	MLP	Sex, age, height, weight, MIN3, HR3	0.85	3.13
Akay et al^Citation22	2014	MLR	Sex, age, height, weight, MIN3, HR3	0.83a	3.25

Note:

a Represents adjusted R (R_adj).

Abbreviations: VO₂max, maximal oxygen uptake; SEE, standard error of estimate; MLR, multiple linear regression; MLP, multilayer perceptron; HR, heart rate; BMI, body mass index; MIN3, elapsed exercise time at the end of 1.5 mile; HR3, exercise heart rate at the end of 1.5 mile; WR, work rate; AD, accelerometer data; 3MWD, 3-minute walk distance; BF%, body fat percentage; SVM, support vector machine.

Nonexercise prediction models

Akay et al²⁶ developed nonexercise VO₂max prediction models by using SVM and MLP. VO₂max values of 100 subjects (50 males and 50 females) were measured by using a maximal GXT. The variables, such as sex, age, BMI, PFA and PAR, were used to build two nonexercise prediction models. The PFA question is designed to determine how fast participants believe they could walk, jog, or run one mile (PFA-1) and three mile (PFA-2) distance without becoming overly breathless or fatigued. The sum of the responses to the PFA-1 and PFA-2 questions is the PFA score. The PAR question allows participants to self-assess their level of PA within the last past 6 months. For model testing, tenfold cross-validation was utilized. It has been concluded that SVM-based nonexercise prediction models performed better than the MLP-based models and can be valid predictors of VO₂max. The highest R and the lowest SEE values were obtained by the SVM-based model with 0.91 and 3.23 mL kg⁻¹ min⁻¹, respectively.

Cao et al²⁷ examined the use of the accelerometer-determined PA variables as the objective PA variables for predicting VO₂max in Japanese adult men. In total, 127 Japanese adult men aged from 20 years to 69 years were selected as subjects for the study. The nonexercise prediction models were developed and cross-validated with the help of MLR and using the predictor variables, such as age, BMI, step count, and vigorous PA (VPA). The reported values of R and SEE are 0.86 and 3.90 mL kg⁻¹ min⁻¹, respectively.

Schembre and Riebe²⁸ established a VO₂max estimation equation derived from the International Physical Activity Questionnaire (IPAQ)–Short Form. To conduct the regression analysis, 80 college-aged males and females completed the IPAQ–Short Form and completed a maximal exercise test. The estimation equation was created using MLR along with predictor variables sex and VPA. It has been shown that the IPAQ–Short Form can be used in successfully estimating VO₂max. The most accurate values of R and SEE were reported as 0.65 and 5.45 mL kg⁻¹ min⁻¹, respectively.

Jang et al²⁹ determined whether work-related PA (WPA) is a potential predictor of VO₂max and developed a VO₂max equation using a MLR-based nonexercise regression model for the cardiorespiratory fitness test in Korean adult workers. In total, 217 (113 males and 104 females) adult subjects aged from 21 years to 63 years were recruited for the study. Self-report questionnaire survey was conducted on study subjects, and VO₂max of each subject was measured with the exercise test. The variables used for the prediction of VO₂max include sex, age, BMI, smoking, leisure-time PA (LPA), and the factors representing WPA. In conclusion, the WPA was revealed to be a predictor of VO₂max. The R and SEE values were presented as 0.89 and 3.36 mL kg⁻¹ min⁻¹, respectively.

Shenoy et al³⁰ employed MLR to produce an equation-based model to predict treadmill VO₂max scores based on nonexercise data. In this study, 120 college-aged participants (60 males and 60 females) voluntarily participated and successfully completed a maximal GXT on a motorized treadmill to assess VO₂max. Relevant nonexercise data included PFA and PAR scores as well as the body surface area. The accuracy of the model was evaluated by conducting a cross-validation analysis. The presented R and SEE values were 0.90 and 0.42 mL kg⁻¹ min⁻¹, respectively.

Akay et al³¹ employed SVM to develop new nonexercise VO₂max prediction models for healthy subjects ranging in age from 18 years to 65 years. The dataset included VO₂max values of 126 subjects as well as the predictor variables such as age, sex, BMI, BF%, lean body mass, and finally, the questionnaire variable AC that describes the level of PA. In order to compare the performance of SVM-based VO₂max prediction models, MLP and MLR prediction models were developed, too. Tenfold cross-validation was used on the dataset. The highest R and the lowest SEE values for the most accurate SVM-based model were presented as 0.89 and 4.68 mL kg⁻¹ min⁻¹, respectively.

A summary of recent studies that developed nonexercise models for the prediction of VO₂max is presented in Table 3 in chronological order.

Table 3 Summary of recent studies in literature that developed nonexercise models for the prediction of VO₂max

Study	Year	Method	Predictor variables	R	SEE (mL kg^−1 min^−1)
Akay et al^Citation26	2009	SVM	Sex, age, BMI, PFA, PAR	0.93	3.41
Akay et al^Citation26	2009	MLP	Sex, age, BMI, PFA, PAR	0.91	3.23
Cao et al^Citation27	2010	MLR	Age, BMI, SC, VPA	0.85a	3.90
Schembre and Riebe^Citation28	2011	MLR	Sex, VPA	0.63a	5.45
Jang et al^Citation29	2012	MLR	Sex, age, BMI, smoking, LPA, WPA	0.88a	3.36
Shenoy et al^Citation30	2012	MLR	Sex, PFA, BSA	0.89a	0.42
Akay et al^Citation31	2014	SVM	Sex, age, BMI, BF%, AC	0.89	4.68
Akay et al^Citation31	2014	MLP	Sex, age, BMI, BF%, LBM, AC	0.87	5.00
Akay et al^Citation31	2014	MLR	Sex, age, BMI, BF%, AC	0.84a	5.12

Note:

a Represents adjusted R (R_adj).

Abbreviations: VO₂max, maximal oxygen uptake; SEE, standard error of estimate; SVM, support vector machine; BMI, body mass index; PFA, perceived functional ability; PAR, physical activity rating; MLP, multilayer perceptron; MLR, multiple linear regression; SC, step count; VPA, vigorous physical activity; LPA, leisure-time physical activity; WPA, work-related physical activity; BSA, body source area; BF%, body fat percentage; AC, activity code; LBM, lean body mass.

Hybrid models

Nielson et al³² built a MLR model for predicting treadmill VO₂max scores using both exercise and nonexercise variables. For the study, 105 college-aged participants (53 males and 52 females) successfully accomplished a submaximal cycle ergometer test and a maximal GXT on a motorized treadmill. Relevant data incorporated weight, ending steady-state HR, ending work rate, PFA, and PAR. The values of R and SEE were reported as 0.91 and 3.36 mL kg⁻¹ min⁻¹, respectively.

Akay et al³³ built a MLP model using submaximal and nonexercise variables to predict VO₂max in healthy volunteers. The dataset included data of 100 subjects (50 males and 50 females) ranging in age from 18 years to 65 years. The predictor variables of the dataset were sex, age, BMI, PAR, HR, and treadmill stage, and tenfold cross-validation has been carried out on the dataset. The presented SEE and R values were 0.96 and 2.01 mL kg⁻¹ min⁻¹, respectively. In a follow-up work, Akay et al³⁴ developed VO₂max prediction models based on maximal and nonexercise variables, with the help of MLP and MLR using the same dataset. Predictor variables used to form the prediction models incorporated the physiological variables sex, age, and BMI; the maximal variables grade, RPE and HRmax; and finally the questionnaire variables PFA and PAR. The highest R and the lowest SEE values were reported for the MLP-based model with 0.94 and 2.28 mL kg⁻¹ min⁻¹, respectively.

Yücel et al³⁵ investigated the effects of PFA-1, PFA-2, and PAR variables on the prediction of VO₂max by using MLP and MLR. Rest of the predictor variables used to form the prediction models were the physiological variables, such as sex, age, and BMI, and the submaximal variables ending speed, HR, and stage of the exercise test. According to the results obtained, two conclusions were made. First, inclusion of PFA and PAR scores in VO₂max models yielded 18.26%–28% and 17.92%–26.99% reductions in SEEs for MLR- and MLP-based models, respectively. Second, MLP-based models showed much better performance than MLR-based models. Among the results obtained by MLP-based models, the highest values for the R and SEE were found as 0.93 and 3.15 mL kg⁻¹ min⁻¹, respectively. A similar study on the same dataset was conducted by Akay et al,³⁶ with the intention to examine the effects of the PFA-1, PFA2, and PAR variables on the prediction of VO₂max using SVM. Predictor variables included the physiological variables sex, age, weight, and height; the submaximal variables ending speed and HR; and the three questionnaire variables PFA-1, PFA-2, and PAR. The inclusion of PFA and PAR scores in SVM-based VO₂max models yielded an average of 8.9% and 10.51% reductions in SEEs for SVM- and MLR-based models, respectively. The values for the highest R and the lowest SEE were found as 0.93 and 3.25 mL kg⁻¹ min⁻¹, respectively.

Table 4 lists recent studies in literature that developed hybrid models for the prediction of VO₂max in chronological order.

Table 4 Summary of recent studies in literature that developed hybrid models for the prediction of VO₂max

Study	Year	Method	Predictor variables	R	SEE (mL kg^−1 min^−1)
Nielson et al^Citation32	2010	MLR	Sex, weight, HR, WR, PFA	0.90a	3.36
Akay et al^Citation33	2010	MLP	Sex, age, BMI, PAR, HR, stage	0.96	2.01
Akay et al^Citation34	2011	MLP	Sex, age, BMI, HRmax, PFA, PAR	0.94	2.23
Akay et al^Citation34	2011	MLR	Sex, age, BMI, grade, PFA, PAR	0.85a	3.73
Yücel et al^Citation35	2013	MLP	Sex, age, BMI, MPH, PAR, PFA	0.93	3.15
Yücel et al^Citation35	2013	MLR	Sex, age, BMI, MPH, PAR, PFA	0.92a	3.27
Akay et al^Citation36	2014	SVM	Sex, age, weight, height, ES, PFA-1, PAR	0.93	3.25
Akay et al^Citation36	2014	MLR	Sex, age, weight, height, ES, PFA-1, PAR	0.87a	4.20

Note:

a Represents adjusted R (R_adj).

Abbreviations: VO₂max, maximal oxygen uptake; SEE, standard error of estimate; MLR, multiple linear regression; HR, heart rate; WR, ending work rate; PFA, perceived functional ability; MLP, multilayer perceptron; BMI, body mass index; PAR, self-reported level of physical activity; HR, heart rate; HRmax, maximal heart rate; MPH, miles per hour; SVM, support vector machine; RPE, self-reported rating of perceived exertion from treadmill test; ES, ending speed; PFA-1, PFA for one mile distance.

Overview of prediction methods and model evaluation techniques

Machine learning and statistical methods

SVM is considered as one of the most promising regression methods, which is widely utilized in many application areas due to its high prediction accuracy. Although SVMs were initially designed to solve classification problems, they have been lately reconstructed to also cope up with regression problems. In order to perform a regression analysis, SVM constructs a hyperplane or set of hyperplanes in a high- or infinite-dimensional space. The effectiveness of SVM-based models depends on the selection of the kernel function, the kernel function’s parameters, and the regularization parameter C. A common choice is the radial basis function kernel, which has a single-optimization parameter gamma (γ). The optimal values of C and γ are often chosen by a grid search method.³⁷

MLP belongs to the class of feed-forward artificial neural networks that can map sets of input data onto a set of convenient outputs. The basic constituents of a MLP are the neurons, which imitate the processing of human brain neurons. The neurons are organized into layers, and in any two contiguous layers, the neurons are linked in pairs with weighted edges. A practical MLP network possesses one input layer, one or more hidden layers, and one output layer. Usually, the back propagation algorithm is employed for training this network. The activation functions of the hidden and output layers, the number of neurons in the hidden layer, the number of epochs, the learning rate, and the momentum determine the performance of MLP-based models.³⁸

GRNN falls into the category of probabilistic neural networks and does not necessitate an iterative training procedure as back propagation networks. It approximates any arbitrary function between input and output vectors, drawing the function estimate directly from the training data. A GRNN has four layers, which are input layer, pattern layer, summation layer, and output layer. Minimum and maximum values for sigma, search step, and the type of kernel function influences the performance of GRNN-based models.³⁹

Finally, MLR is almost a standard technique frequently used in sport physiology to investigate the linear relationships between one or more predictor variables and a single-target variable. The case of regression analysis with one predictor variable is also called as simple linear regression. MLR can be considered as an extension of simple linear regression in such a way that it can use two or more predictor variables in an equation to predict a desired target variable.

Model validation and testing techniques

Both model validation and testing evaluate a model’s performance, with the difference that they are applied at different stages of a model’s lifetime. Validation is performed to evaluate the goodness of fit of a learned model, which in turn reflects the quality and performance of hyperparameters optimized for a learning algorithm during model training. Validation of a prediction model is carried out prior to real testing or deployment. Once the prediction model is trained and deployed, it can be continued with the stage of model testing, where the performance of a validated model is evaluated on a test dataset. This kind of testing does not necessitate retraining the model. For validation and testing purposes, one and the same dataset can be used as long as the test set does not contain the validation set.

For model validation and testing, there exists in principle two techniques, which are referred to as hold-out validation and cross-validation. Hold-out validation is straightforward in concept, and the dataset is partitioned into two disjoint training and testing sets. A prediction model is trained using the training set only. Then, the trained model is evaluated against the previously unseen data in the testing set to predict the desired target values. Normally, less than one-third of the initial dataset is utilized for the validation of data. The advantage of the hold-out validation is that it is executed on independent training and test sets, which may increase the reliability of model validation and testing, also it takes no longer to compute. However, its evaluation can have a high variance. The performance of this technique may rely heavily on the fact at which point the separation between training and test set is carried out. Another drawback refers to the fact that the number of training examples is reduced due to the strict separation of the original dataset in training and test sets. Cross-validation, as being the alternative to hold-out validation, is utilized in practice with the expectation of being more accurate than the hold-out estimate without reducing the number of training examples within a dataset. Although there exist many variants of cross-validation, the most commonly known variant is k-fold cross-validation. The k-fold procedure partitions the data into k equally sized folds. For a given a hyperparameter setting, every k-fold is handled similar to a hold-out set. The model is trained on the remaining k-1 folds, and its quality and performance are validated on the held-out fold. This principle is reiterated for each fold, and the overall performance is calculated as the average of the performance on all k-folds. In data mining and machine learning, tenfold cross-validation (k=10) is the most common variant. If k equals the sample size, this is called leave-one-out cross-validation. Cross-validation is especially useful when the dataset comprises a limited number of entries and cannot be reduced at the expense of validation purposes.

Conclusion and future work

This study presented an overview about the data-driven modeling studies for the prediction of VO₂max conducted in recent years. Numerous models have been presented to predict VO₂max of various target audiences, ranging from soccer athletes, nonexpert swimmers, cross-country skiers to healthy fit adults, teenagers, and children using various machine learning and statistical methods and different sets of predictor variables. Representative examples of variables for the prediction of VO₂max range from physiological variables, such as sex, age, weight, height, and BMI, to maximal or submaximal exercise variables, such as HRmax, RPE, treadmill stage, and treadmill speed, to nonexercise variables, including PFA, PAR, and AC. The performance of the presented prediction models has been compared with respect to their achieved R and SEE values.

Considering the survey results obtained, several conclusions can be reached. First of all, at the beginning of research efforts for the prediction of VO₂max, almost all of the studies used MLR for developing various VO₂max prediction equations. However, with the employment of intelligent data processing techniques, such as SVM and MLP, it has been shown that VO₂max prediction models based on these techniques led to higher Rs and lower SEEs than MLR-based prediction models. In more detail, SVM-based models, in general, have shown better performance than other regression methods for the prediction of VO₂max. In contrast, MLR-based models, in general, have yielded the lowest Rs and the highest SEEs, suggesting that VO₂max has nonlinear characteristics.

Second, when all prediction models presented in this study are investigated, it is observed that, on average, submaximal exercise-based models perform better than nonexercise-based models, whereas nonexercise-based models, in turn, outperform maximal exercise-based models for the prediction of VO₂max. Hybrid models, on the other hand, that combine nonexercise-based models, with maximal- or submaximal-based models, have the potential to improve the accuracy of VO₂max prediction over single models.

Third, feature selection algorithms have been applied on different datasets to reveal the relevant and irrelevant features for the prediction of VO₂max. The results show that using the full set of predictor variables does not always lead to most accurate results. In contrast, it has been reported that VO₂max can be predicted with higher Rs and lower SEEs using lesser number of predictor variables.

Fourth, the physiological variables, such as sex, age, height, weight, and BMI, in general, have been found to be the most important variables in building accurate prediction models, and they accordingly appear in almost each model giving predictions within acceptable limits of accuracy.

Future work can be performed in a number of different areas. Other promising regression methods as proposed in studies of Evrendilek and Karakaya^40,41 can be utilized that have never been applied for the development of VO₂max prediction models. Different further feature selection algorithms such as minimum redundancy maximum relevance⁴² and correlation-based feature selection⁴³ can be combined with these regression methods to identify the discriminative features for the prediction of VO₂max.

Also, ensemble methods can be utilized in developing new VO₂max models that combine the predictions of several base estimators built with a given learning algorithm in order to improve generalizability and robustness over a single estimator.

Disclosure

The authors report no conflicts of interest in this study.