Use of the Model
The fitness–fatigue model has been used in order to relate TL to changes in performance in different sports such as swimming,3–5 cycling,6–8 running,9,10 and triathlon.11 In these studies, a good to excellent model fit has often been found, suggesting that the model is able to relate TL with performance in an adequate manner. Based on these results, the model and derivatives of the model were integrated in several coaching platforms, making it accessible to a large group of athletes and coaches. This integration was made to support coaches in assessing and predicting the performance of their athletes without having to do extensive testing. It has also enabled coaches to make the right selection of TL in the planning and periodization phase, and to understand how the training program, with the corresponding TL, influences performance over time. However, in a lot of platforms, a simplification of the original model is used where the k parameter is removed, which then leads to a more simple (exponentially) weighted average. This means that the conversion to a real performance measure is no longer possible, and arbitrary units are used. This also implies that conclusions based on one model, or the other are not per se directly transferable to each other.
However, mostly average values for the τ and k parameters are used to establish these models since an individual model fit is hard to obtain in a real-world setting. For example, the τ parameter is often estimated to be 42 days for fitness and 7 days for fatigue. It is known, however, that an individual model fit is necessary when trying to relate TL and performance.12 In literature, values are found ranging from 4 to 51, and from 4 to 74 for τ1, and τ2, respectively.2,4,12–15 Although some coaching platforms enable coaches to adjust these values to the user’s preference, it is futile to estimate what parameter values will best suit the individual athlete without fitting the model.
To fit the model, TL is needed as an input. However, many different methods to quantify TL exist nowadays with the external TL methods that quantify the work completed (power, distance, speed, etc) and internal TL methods, which represent the psychophysiological response to training (perceived exertion, heart rate, lactate, etc). This leads to different TL scores for the same training sessions, depending on the method used. Moreover, studies have shown that different input variables (ie, TL calculations) can lead to substantially different results in the parameter values. The study of Mitchell et al13 and the study of Vermeire et al12 have both shown that the input of different TL quantification methods will lead to different τ values on the same data set. The parameter values of these studies differ considerably from the original papers, which gives rise to question the validity of the commonly used parameter values. Therefore, we want to point out 3 possible pitfalls when using performance models: (1) the interpretation of the model parameters, (2) the selection of the modeling technique, and (3) the input of the model.
Interpretation of the Model Parameters
When the model was originally constructed, it was stated that τ1 and τ2 are time constants which represent the number of days that the fitness and the fatigue, generated by a training session, will be apparent in the body.1 However, interpretation of this parameter is more complex and should be done with caution. To demonstrate this, a result is shown from a global optimization (see below) of the model to data of one of the subjects in an earlier study in Figure 1.12 As this optimization scanned a series of values for each combination of parameters in the fitness–fatigue model (a range of 0–3 for k and 0–60 for τ), a minimum absolute model error (residual sum of squares) could be plotted for each parameter value combination in Figure 1. From this figure, we can deduct that different subsets of parameter values can lead to similar model errors, which makes it hard to select the right parameter values. As previously shown by Hellard et al,4 k1 and k2 but also τ1 and τ2, are highly correlated. This ill-conditioning of the model could therefore explain the plethora of possible parameter values that give similar model errors.
In other words, this means that we cannot interpret the τ parameter as being a specific number of days, and that the k parameter does not solely serve as a magnitude and conversion factor, but rather that the fitness–fatigue model should be treated as one entity. Several authors have also previously indicated that the practical interpretation of the parameters might be difficult. For example, a positive correlation between the fatigue function and testosterone concentration was found, where a negative relation was expected.16 The model parameters were also criticized based on the fact that most of the models were poorly supported by physiological mechanisms.17
Selection of the Modeling Technique
The model deficiencies as explained above are noticeable since global optimization techniques were used to fit the model. However, global optimization techniques are technically advanced; and therefore, coaches in the field do not always have the possibility, or the means, to infer such a model fit. Nevertheless, there are more accessible ways to fit a model by means of local optimization. As proposed by Clarke and Skiba,2 one can use an Excel spreadsheet and fit the model to the individual by means of the Solver function in Microsoft Excel. The pitfall here is that local optimization values will depend greatly on the starting values. In Table 1, the locally optimized parameter values are presented, using the Solver function (generalized reduced gradient [GRG] nonlinear) in Excel (MS Excel, version 2111) without constraints, for 4 different sets of starting parameter values, taken from literature,2,14,15 using the same training and performance data. The weakness of using local optimization techniques is reflected in this table. The 4 sets of starting values lead to 3 different local optima (sets 2 and 4 are similar), which are determined to a great extent by the starting values.
Local Optimization Values for Different Sets of Starting Values
Set no. | Starting values | Local optima | Model error (RSS) | |
---|---|---|---|---|
Set 1 | k1 | 1 | 0.6407 | 378.47 |
k2 | 2 | 0.7146 | ||
τ1 | 42 | 10.5313 | ||
τ2 | 7 | 8.0583 | ||
Set 2 | k1 | 0.18 | 0.0649 | 323.06 |
k2 | 0.23 | 0.1057 | ||
τ1 | 36 | 24.1714 | ||
τ2 | 21 | 3.5542 | ||
Set 3 | k1 | 0.0048 | 0.037352 | 367.26 |
k2 | 0.3860 | 0.055929 | ||
τ1 | 49 | 39.94221 | ||
τ2 | 4.3 | 0.05 | ||
Set 4 | k1 | 0.0193 | 0.065958 | 323.05 |
k2 | 0.0148 | 0.106808 | ||
τ1 | 40.8 | 23.89291 | ||
τ2 | 9.0 | 3.619263 |
Abbreviation: RSS, residual sum of squared errors. Note: The different starting values are selected from literature.
In short, local optimization techniques search for the smallest model error in a smaller parameter space. The optimization will run as long as the model error keeps getting smaller. As soon as the model error starts rising again, the optimization will stop running at this local optimum (Figure 2). Initiating the optimization from different parameter values could lead to a smaller error and thus the starting set will determine the result. Global optimization techniques work around this problem in different ways, resulting in a true minimum error and thus more correct parameter values. Therefore, we advise, as Clarke and Skiba2 already suggested, when choosing for a local optimization method, different sets of starting parameter values should be used so to minimize the chance of missing the true optimum.
Input of the Model
A last point to consider when interpreting the model parameters is the input (TL) of the model. Different methods to quantify the TL exist nowadays, either based on internal measures (eg, heart rate, rate of perceived exertion, lactate), or external measures (eg, power, speed). At first glance, it seems that all these methods result in reasonably small model errors, suggesting it does not matter which TL method is used. However, research has shown that different TL quantification methods lead to different k and τ values.12,13 In an ideal situation, the τ values would not be different, depending on the input, and only the k values would be influenced since this is the factor that makes it possible to convert the input to the required output. Again, this implies that the absolute parameter values are to be interpreted with caution. Often in practice, different methods are used to quantify TL due to data registration restrictions or data capture failure. However, these data clearly show that fitting the fitness–fatigue model requires that only one quantification method is used throughout the entire period.
Also, since most of the existing TL methods are a combination of volume and intensity, the specificity of training is not represented in this score. For instance, a training session of 2 hours at low intensity would result in a TL score of 120 (using Lucia TL),18 but a high intensity bout of 40 minutes can result in an identical score. Since the training adaptations performing such different training sessions with a similar TL are totally different,19,20 the relationship with performance improvement will always be distorted.
Value of the Model
Despite the weaknesses of the model presented here and the clear message that the model should not be followed blindly to drive the training program, it is still the belief of the authors that the fitness–fatigue model can be of value to sport scientists and coaches. The true value of the model is in understanding not only how a certain training period influences performance, but also in interpreting possible other influencing factors during that period that may account for the unexpected values derived from the model (psychological issues, environmental conditions. . .). The model should be seen as an objective method to quantify this relationship, and as a means to highlight performances or training periods where this relationship seems off. It is the job of the coach and sport scientist in that case to interpret these data in light of their scientific knowledge on the physiology of training adaptations. It goes without saying that input from the athlete, and formal and informal communication between athlete, coach, and sport scientist are of utmost importance here. Also, fitting the model on historical data and altering the TL of that training period may lead to theoretical improvements in performance, which may in turn lead to insights for individual periodization strategies.
However, since the model parameters cannot be interpreted separately and since they will differ strongly between athletes and over different training periods,14 the conclusions taken from these models can only be used to support coaches. Nonetheless, if the pitfalls are accounted for, individualization of the parameters could aid in individualizing training programs.21
Practical Applications
For both coaches and sport scientists, we give some short recommendations when using the fitness–fatigue model:
- •Use only one TL metric to fit the fitness–fatigue model. However, apart from the model, different metrics should be collected and compared to each other. When the external internal TL ratio dissociates, this could indicate training (mal-)adaptations.22
- •When using local optimization techniques, use different sets of starting parameter values to infer the smallest error.
- •Only compare parameter values intrasubject, and track changes over different training periods.
- •Always interpret the parameter values using physiological knowledge and training experience.
Conclusion
Caution is needed when interpreting the fitness–fatigue model since the parameter values are influenced by the starting parameter values, the modeling technique, and the input of the model. The use of general parameter values should be avoided since they do not account for interindividual differences and differences between TL methods. Therefore, the authors advise coaches and sport sicentists to use the model as a way to work more data-informed rather than working data-driven.
References
- 1.↑
Banister EW, Calvert TW, Savage MV, Bach T. A systems model of training for athletic performance. Aust J Sports Med. 1975;7(3):57–61.
- 2.↑
Clarke DC, Skiba PF. Rationale and resources for teaching the mathematical modeling of athletic training and performance. Adv Physiol Educ. 2013;37(2):134–152. PubMed ID: 23728131 doi:10.1152/advan.00078.2011
- 3.↑
Mujika I, Busso T, Lacoste L, Barale F, Geyssant A, Chatard JC. Modeled responses to training and taper in competitive swimmers. Med Sci Sports Exerc. 1996;28(2):251–258. PubMed ID: 8775162
- 4.↑
Hellard P, Avalos M, Lacoste L, Barale F, Chatard JC, Millet GP. Assessing the limitations of the Banister model in monitoring training. J Sports Sci. 2006;24(5):509–520. PubMed ID: 16608765 doi:10.1080/02640410500244697
- 5.↑
Hellard P, Avalos M, Millet G, Lacoste L, Barale F, Chatard JC. Modeling the residual effects and threshold saturation of training: a case study of Olympic swimmers. J Strength Cond Res. 2005;19(1):67–75. PubMed ID: 15705048
- 6.↑
Busso T, Carasso C, Lacour JR. Adequacy of a systems structure in the modeling of training effects on performance. J Appl Physiol. 1991;71(5):2044–2049. PubMed ID: 1761506
- 7.
Busso T, Benoit H, Bonnefoy R, Feasson L, Lacour JR. Effects of training frequency on the dynamics of performance response to a single training bout. J Appl Physiol. 2002;92(2):572–580. PubMed ID: 11796666 doi:10.1152/japplphysiol.00429.2001
- 8.↑
Busso T, Denis C, Bonnefoy R, Geyssant A, Lacour JR. Modeling of adaptations to physical training by using a recursive least squares algorithm. J Appl Physiol. 1997;82(5):1685–1693. PubMed ID: 9134920
- 9.↑
Morton RH, Fitz-Clarke JR, Banister EW. Modeling human performance in running. J Appl Physiol. 1990;69(3):1171–1177. PubMed ID: 2246166 doi:10.1152/jappl.1990.69.3.1171
- 10.↑
Wood RE, Hayter S, Rowbottom D, Stewart I. Applying a mathematical model to training adaptation in a distance runner. Eur J Appl Physiol. 2005;94(3):310–316. PubMed ID: 15765236 doi:10.1007/s00421-005-1319-2
- 11.↑
Millet GP, Candau RB, Barbler B, Busso T, Rouillon JD, Chatard JC. Modelling the transfers of training effects on performance in elite triathletes. Int J Sports Med. 2002;23(1):55–63. PubMed ID: 11774068 doi:10.1055/s-2002-19276
- 12.↑
Vermeire KM, van de Casteele F, Gosseries M, Bourgois JG, Ghijs M, Boone J. The influence of different training load quantification methods on the fitness-fatigue model. Int J Sports Physiol Perform. 2021;16(9):1261–1269. doi:10.1123/ijspp.2020-0662
- 13.↑
Mitchell LJG, Rattray B, Fowlie J, Saunders PU, Pyne DB. The impact of different training load quantification and modelling methodologies on performance predictions in elite swimmers. Eur J Sport Sci. 2020;20(10):1329–1338. PubMed ID: 31955652 doi:10.1080/17461391.2020.1719211
- 14.↑
Busso T. Variable dose-response relationship between exercise training and performance. Med Sci Sports Exerc. 2003;35(7):1188–1195. PubMed ID: 12840641 doi:10.1249/01.MSS.0000074465.13621.37
- 15.↑
Chalencon S, Pichot V, Roche F, et al. Modeling of performance and ANS activity for predicting future responses to training. Eur J Appl Physiol. 2015;115(3):589–596. PubMed ID: 25359446 doi:10.1007/s00421-014-3035-2
- 16.↑
Busso T, Häkkinen K, Pakarinen A, et al. A systems model of training responses and its relationship to hormonal responses in elite weight-lifters. Eur J Appl Physiol. 1990;61:48–54. doi:10.5694/j.1326-5377.1950.tb106955.x
- 17.↑
Taha T, Thomas SG. Systems modelling of the relationship between training and performance. Sports Med. 2003;33(14):1061–1073. PubMed ID: 14599233 doi:10.2165/00007256-200333140-00003
- 18.↑
Lucía A, Hoyos J, Santalla A, Earnest C, Chicharro JL. Tour de France versus Vuelta a España: which is harder? Med Sci Sports Exerc. 2003;35(5):872–878. PubMed ID: 12750600 doi:10.1249/01.MSS.0000064999.82036.B4
- 19.↑
Vermeire KM, Vandewiele G, Caen K, Lievens M, Bourgois JG, Boone J. Training progression in recreational cyclists: no linear dose-response relationship with training load. J Strength Cond Res. 2021;35(12):3500–3505. PubMed ID: 31498226 doi:10.1519/jsc.0000000000003340
- 20.↑
Laursen PB. Training for intense exercise performance: high-intensity or high-volume training? Scand J Med Sci Sports. 2010;20(suppl 2):1–10. doi:10.1111/j.1600-0838.2010.01184.x
- 21.↑
Kellmann M, Bertollo M, Bosquet L, et al. Recovery and performance in sport: consensus statement. Int J Sports Physiol Perform. 2018;13(2):240–245. PubMed ID: 29345524 doi:10.1123/ijspp.2017-0759
- 22.↑
Bourdon PC, Cardinale M, Murray A, et al. Monitoring athlete training loads: consensus statement. Int J Sports Physiol Perform. 2017;12:161–170. doi:10.1123/IJSPP.2017-0208