Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Observed and predicted ages at peak height velocity in soccer players

  • Robert M. Malina,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

    Affiliations Department of Kinesiology and Health Education, University of Texas, Austin, Texas, United States of America, University of Louisville, School of Public Health and Information Sciences, Louisville, Kentucky, United States of America

  • Manuel J. Coelho-e-Silva ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    mjcesilva@hotmail.com

    Affiliations University of Coimbra, FCDEF, Coimbra, Portugal, University of Coimbra, CIDAF (uid/dtp/04213/2020), Coimbra, Portugal

  • Diogo V. Martinho,

    Roles Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations University of Coimbra, FCDEF, Coimbra, Portugal, University of Coimbra, CIDAF (uid/dtp/04213/2020), Coimbra, Portugal

  • Paulo Sousa-e-Siva,

    Roles Investigation, Methodology, Validation, Visualization, Writing – original draft

    Affiliations University of Coimbra, FCDEF, Coimbra, Portugal, University of Coimbra, CIDAF (uid/dtp/04213/2020), Coimbra, Portugal

  • Antonio J. Figueiredo,

    Roles Data curation, Investigation, Validation, Visualization, Writing – original draft

    Affiliations University of Coimbra, FCDEF, Coimbra, Portugal, University of Coimbra, CIDAF (uid/dtp/04213/2020), Coimbra, Portugal

  • Sean P. Cumming,

    Roles Formal analysis, Supervision, Validation, Visualization

    Affiliation Department of Health, University of Bath, Bath, United Kingdom

  • Miroslav Králík,

    Roles Formal analysis, Methodology, Software, Validation

    Affiliation Faculty of Science, Department of Anthropology, Masaryk University, Brno, Czech Republic

  • Sławomir M. Kozieł

    Roles Formal analysis, Investigation, Methodology, Software, Validation, Writing – original draft

    Affiliation Department of Anthropology, Polish Academy of Sciences, Hirszfeld Institute of Immunology and Experimental Therapy, Wrocław, Poland

Abstract

The purpose of the study was to evaluate predicted maturity offset (time before age at PHV) and age at PHV (chronological age [CA] minus maturity offset) in a longitudinal sample of 58 under-13 club level soccer players in central Portugal for whom ages at PHV were estimated with the SITAR model. Two maturity offset prediction equations were applied: the original equation which requires CA sitting height, estimated leg length, height and weight, and a modified equation which requires CA and height. Predicted maturity offset increased, on average, with CA at prediction throughout the age range considered, while variation in predicted maturity offset and ages at PHV within CA groups was considerably reduced compared to variation in observed ages at offset and at PHV. Predicted maturity offset and ages at PHV were consistently later than observed maturity offset and age at PHV among early maturing players, and earlier than observed in late maturing players. Both predicted offset and ages at PHV with the two equations were, on average, later than observed among players maturing on time. Intra-individual variation in predicted ages at PHV with each equation was considerable. The results for soccer players were consistent with similar studies in the general population and two recent longitudinal studies of soccer players. The results question the utility of predicted maturity offset and age at PHV as valid indicators of maturity timing and status.

Introduction

Predicted maturity offset, defined as time before peak height velocity (PHV) [1,2], and estimated age at PHV, i.e., chronological age (CA) minus predicted offset, are widely used as estimates of maturity status (state of maturation at the time of observation) and/or timing (age at which a specific maturational event occurs) in studies of youth athletes and to a lesser extent in studies youth physical activity and fitness [36]. The original sex-specific equations require CA, sitting height, estimated leg length, height and weight [1], while the modified equations [2] require CA and height (both sexes) or sitting height (boys). Validation studies of the original equations in three independent longitudinal series, the Wroclaw Growth Study [7,8], the Fels Longitudinal Study [9] and the Cracow Growth Study [10], and of the modified equations in two of the samples [10,11] have indicated major limitations of the predictions in both males and females. The validity of the prediction equations has also been questioned in longitudinal samples of female artistic gymnasts [12] and soccer players [13,14], but sample sizes in longitudinal samples tend to be limited and to some extent select as they are limited to athletes who have persisted in the respective sports [15]. Similarly, cross-sectional studies of tennis [16] and soccer [17] players have questioned maturity status classifications based on the original prediction equations [1] relative to classifications based on skeletal age.

Current interest in the application of the maturity offset prediction equations in samples of youth athletes is considerable [36]. Predictions based on the original and modified equations, however, depend upon CA and body size at prediction, have reduced variation relative to observed ages at PHV, and have major limitations with early and late maturing youth as defined by observed ages at PHV [711]. The latter are problematic as advanced (early for CA) or delayed (late for CA) maturity status is often of major concern in developmental studies of youth athletes [15].

In the context of the preceding, the purpose of the present study is to evaluate predicted maturity offset and derived ages at PHV with the original [1] and modified [2] equations in a sample of 58 soccer players for whom age at PHV was determined from longitudinal height records. The study specifically considers variation in the predictions in three contexts: (i) relative to actual maturity offset and observed age at PHV at each observation, (ii) among players differing in the timing of observed age at PHV, and (iii) within individual players.

Materials and methods

Research design and procedures

The data set for the present study was extracted from the Coimbra Soccer Longitudinal Project [18]. This project was conducted according to the standards established by the declaration of Helsinki [19], and formal approval was obtained from the University of Coimbra Sports Sciences and Physical Education Board (FCDEFUC/AAC/2003; FCDEFUC/ADCA/2003; FCDEFUC/CFM/2003; FCDEFUC/CFUC/2003; FCDEFUC/GRVM/2003). Signed institutional agreements were also obtained from the Presidents of the respective clubs. All players were registered with the Portuguese Soccer Federation. Male players 11–14 years of age were recruited from five clubs in the Midlands of Portugal; the initial sample included 159 players [20]. Written consent was obtained from parents or legal guardians of the participants who were informed that contribution to the study was voluntary and that they could withdraw from the study at any time. All observations were completed at the Biokinetics Laboratory of the Coimbra University Stadium.

Sample

According to the Portuguese Soccer Federation, male soccer players were grouped as infantiles (aged 11–12 years, n = 87) and initiates (aged 13–14 years, n = 72). The analysis in the current study is limited to under-13 players (U13) who were measured annually in December for four or five years (n = 59). CAsat baseline ranged from 10.98 to 12.94 years. The heights of one player were not successfully modeled; consequently, the final sample was composed of 58 players. All players were of European ancestry, except one. Participants trained and competed September through May. They had a median of 3 years of soccer experience at baseline (range: 2–6 years). The clubs had 3–5 training sessions per week (each about 90–120 minutes) and usually one game, mainly on Saturdays.

Anthropometry

Participants wore shorts and a t-shirt; shoes were removed. Height and sitting height were measured to the nearest 0.1cm using, respectively, a stadiometer (Harpenden 98.603, Holtain Ltd, Croswell, UK) and a table (Harpenden sitting height table, model 98.607, Holtain Ltd, Crosswell, UK). Body weight was measured to the nearest 0.1 kg using a scale (SECA 770, Hanover, MD, USA). The heights, sitting heights and weights of players who continued at the respective clubs were subsequently measured on an annual basis. Measurements were made by a single observer (MJCS). Intra-observer technical errors of measurement for height, sitting height and weight were 0.27 cm, 0.31 cm and 0.47 kg, respectively [19].

Age at PHV

The longitudinal height records were fit with the Superimposition by Translation and Rotation (SITAR) model [2124] to derive an age at peak height velocity for each player. As noted, the heights of one player were not successfully modeled; the estimated age at PHV for the player was outside of the empirical data range and was inconsistent with his advanced skeletal maturity status at observations one and three. Mean age at PHV for the remaining 58 players was 13.60±0.85 years, with a range from 11.89 to 15.49 years [25].

Predicted maturity offset

Maturity offset, defined as time before PHV, was predicted at each observation for the 58 players with the original equation for boys [1]: Eq (1)

Leg length was estimated at each observation as standing height minus sitting height. The need to multiply the weight by height ratio by 100 was overlooked in the original report [1]; in some publications using the equation, it is not clear if the adjustment was applied. Maturity offset was also predicted at each observation with a modified equation for boys that incorporated age and height [2]: Eq (2)

The equation with age and height was selected for evaluation as it is increasingly used [2631]. Standard errors for the original [1] and modified [2] equations were, respectively, 0.592 and 0.542 year. Predicted maturity offset and predicted age at PHV with the respective equations are subsequently labelled in the text, tables and figures as Mirwald and Moore, respectively.

Predicted age at PHV

Predicted age at PHV was estimated as CA minus predicted maturity offset at each observation for individual players with the respective equations.

Observed maturity offset

Observed or actual maturity offset at each observation was estimated as CA at prediction minus observed age at PHV based on the SITAR model.

Analyses

Descriptive statistics were calculated at each observation for CA and actual offset, for predicted maturity offset and age at PHV, and for the difference of predicted age at PHV minus observed age at PHV with the two prediction equations. The players were also classified as advanced (early), average (on time) or delayed (late) maturing relative to the mean and standard deviation for age at PHV (SITAR) in the total sample. Average was defined as an age at PHV within ±1 SD of the mean age at PHV for the total sample of 58 players (13.60±0.85 years); delayed was a PHV >14.45 years and advanced was a PHV <12.75 years.

Linear mixed-effect models with the data grouped by subjects (random effects) were used with a maximum-likelihood estimator to evaluate the variance structure of the dependent variable, i.e., the differences between observed and predicted ages at PHV. Separate analyses were done for predictions with the Mirwald and Moore equations. The difference of the dependent variable from zero (so-called unconditional means model) was initially tested. In the second step, the effect of predictions at observations 1–5 (as a time variable) on the dependent variable at the population level (fixed effect) and at the intra-individual level (i.e., random slopes model, or unconditional growth model) were tested. Note that the term growth in the statistical context refers to the general change in a dependent variable with a time variable; it does not refer to growth in the biological sense. Finally, maturity status was added as a fixed factor to test the effect of variation in maturity timing on the predictions. The model was run separately for predictions with the Mirwald and Moore equations in the R-software [22] with the help of the nlme statistical package [32].

Weighted Deming regression using the jackknife procedure [33] was used to compare observed age at PHV estimated with the SITAR model (i.e., the reference) with predicted ages at PHV at each observation based on the two equations. The weighted Deming procedure considers both x (observed age at PHV) and y (predicted age at PHV) as subject to measurement error whereas simple regression permits only the y variable to have an associated error. With the weighted Deming regression, systematic differences between x and y are indicated by the intercept, while proportional differences are indicated by the slope.

Many applications of the equations use predicted maturity offset to classify youth as pre-PHV, at/circa/mid-PHV, or post-PHV using a band of -0.5 to +0.5 year to define the interval at PHV [29,3441]; a band of -1.0 to +1.0 year is used less often [4246]. On the other hand, some studies do not report the specific cut-offs that were used [4749]. The standard errors of the prediction equations, 0.592 and 0.542 year, also approximate the narrow cut-offs used in many studies. Thus, the number and percentage of predicted ages at PHV with each equation within ±0.5 year of observed age at PHV at each observation were estimated for players of contrasting maturing status and for the total sample.

Results

Descriptive statistics for CA at prediction, observed and predicted maturity offset, predicted ages at PHV and the difference of predicted age at PHV minus observed age at PHV (the reference for comparison) with the original (Mirwald) and modified (Moore) equations in the soccer players are summarized in Table 1. Observed and predicted maturity offset increase linearly, on average, across the five observations. Predicted offset with the Moore equation is similar to actual offset at observation one but is then less than actual offset at subsequent observations. Predicted maturity offset with the Mirwald equation is less than actual offset and predicted offset with the Moore equation across the five observations.

thumbnail
Table 1. Sample sizes and descriptive statistics for chronological age (CA) at prediction, observed maturity offset and predicted maturity offset and ages at PHV, and the difference of predicted age at PHV minus observed ages at PHV (criterion) with the original (Mirwald) and modified (Moore) equations in soccer players at each observation.

https://doi.org/10.1371/journal.pone.0254659.t001

Corresponding trends in predicted ages at PHV and the difference of predicted minus observed ages at PHV parallel those for maturity offset. Standard deviations for predicted maturity offset and ages at PHV with both equations are consistently lower than corresponding standard deviations for observed offset and age at PHV across the five observations. Variability is reduced more so with the Moore compared to the Mirwald equation.

Results of the mixed effects model indicate significant differences between observed and predicted ages at PHV with the Mirwald equation (F = 40.95; p<0.001) and also with the Moore equation (F = 9.39; p<0.01). Details of the analytical protocol and results are summarized in S1 Table. The differences between observed and predicted ages at PHV with the respective equations at each observation increase significantly with subsequent observations with the Mirwald (F = 22.81; p<0.001) and the Moore (F = 172.97; p<0.001) equations, although the 95% confidence intervals indicate that the difference at observation 1 for the Moore equation is not different from zero.

Intercepts based on weighted Deming regressions for the two prediction equations are well above zero at each observation, indicating that the methods differ significantly by a constant error (Table 2). Estimated slopes for each regression are <1.0, indicating significant proportional differences between predicted ages at PHV with each equation and observed age at PHV. Overall, the results indicate systematic error for predicted ages at PHV.

thumbnail
Table 2. Intercepts and slopes, and respective standard errors (SE) and 95% confidence limits based on the weighted Deming regression of predicted ages at PHV (y-axis) and observed (actual) ages at PHV (x-axis) for the Mirwald and Moore prediction equations at each observation (Obs) in youth soccer players*.

https://doi.org/10.1371/journal.pone.0254659.t002

Descriptive statistics for the three maturity groups are summarized in Table 3. Sample sizes and ages at PHV of players in each maturity timing group were as follows: advanced, n = 8, 12.30±0.27 years; average, n = 38, 13.50±0.51 years; delayed, n = 12, 14.76±0.27 years.

thumbnail
Table 3. Sample sizes and descriptive statistics (mean ± standard deviation) for chronological age (CA) at prediction, observed maturity offset and predicted maturity offset, predicted ages at PHV and the difference of predicted age at PHV minus observed ages at PHV (criterion) with the original (Mirwald) and modified (Moore) equations at each observation in players classified as advanced, average and delayed based on observed ages at PHV.

https://doi.org/10.1371/journal.pone.0254659.t003

Results of the mixed-effects model comparing the three maturity groups indicate that maturity status as a fixed factor has a significant effect on the dependent variable (predicted ages at PHV) with both the Mirwald (F = 36.85, p<0.001) and the Moore (F = 51.28, p<0.001) equations. By inference, predicted ages at PHV differ relative to observed age at PHV in each group. However, the interaction between observation and maturity group is significant for the Mirwald equation (F = 5.39, p = 0.005) and indicates different slopes of change in predicted ages at PHV with consecutive observations in the maturity groups. The latter reflects the trend in differences between predicted and observed ages at PHV for the Mirwald equation which are not significant except at observation 1. In contrast, the interaction term is not significant for the Moore equation.

In the context of the results of the mixed-effects model, differences between predicted and observed ages at PHV are significant and positive across the five observations among players advanced in maturity timing (i.e., early ages at PHV). The predicted ages at PHV are consistently later observed ages at PHV.

Results are similar for players maturing on time (average), i.e., predicted ages at PHV with the two equations are later than observed age at PHV. Across the five observations, the differences between predicted and observed ages at PHV with each equation are significant, although predicted age at PHV with the Moore equation at observation one approaches that for observed age at PHV.

Among late maturing players, in contrast, differences between predicted and observed ages at PHV for the Mirwald equation are not significant except at observation one, while differences between predicted and observed ages at PHV for the Moore equation are significant at each observation. The differences between predicted and observed ages at PHV with the Mirwald equation are also smaller than corresponding differences with the Moore equation.

Predicted ages at PHV (y-axis) for individual players with the Mirwald and Moore equations are illustrated relative to their respective observed ages at PHV (x-axis) in Figs 1 and 2, respectively. Intra-individual variation in predicted ages at PHV is considerable and ranges of predicted ages are reduced with the Moore equation. Relatively few predicted ages approximate observed ages at PHV in early and late maturing players.

thumbnail
Fig 1. Predicted ages at PHV with the Mirwald equation plotted relative of observed age at PHV at each observation for individual soccer players classified as early, average and late maturing.

The diagonal corresponds to the line of identity (x = y).

https://doi.org/10.1371/journal.pone.0254659.g001

thumbnail
Fig 2. Predicted ages at PHV with the Moore equation plotted relative of observed age at PHV at each observation for individual soccer players classified as early, average and late maturing.

The diagonal corresponds to the line of identity (x = y).

https://doi.org/10.1371/journal.pone.0254659.g002

Across the five observations (Table 4), no predicted ages at PHV with the Mirwald equation are within ±0.5 year of observed age at PHV (SITAR) among the eight early maturing players (0 of 39). Corresponding estimates for predicted ages at PHV with the Mirwald equation within ±0.5 year of observed age at PHV across the five observations are 35 of 53 (66%) among late and 76 of 177 (43%) among average maturing players. For the total sample, 111 of 269 (41%) of predicted ages at PHV with the Mirwald equation are within ±0.5 year of observed age at PHV.

thumbnail
Table 4. Number of participants by maturity status according to predicted ages at PHV with the Mirwald equation and, separately, with the Moore equation who were within ±0.50 year of observed age at PHV (SITAR model) at each observation in youth soccer players.

https://doi.org/10.1371/journal.pone.0254659.t004

For the Moore equation, only 1 of 39 predicted ages at PHV (3%) is within ±0.5 year of observed age at PHV among early maturing players. On the other hand, 88 of 177 predicted ages at PHV (50%) among average and 23 of 53 predicted ages at PHV (43%) among late maturing players are within ±0.5 year of observed age at PHV. For the total sample, 112 of 269 (42%) of predicted ages at PHV with the Moore equation are within ±0.5 year of observed age at PHV.

Discussion

The 58 players comprising the present study was larger than samples in five other longitudinal studies of European youth soccer players, 8 to 33 players [25]. Results of the application of the maturity offset prediction equations in the longitudinal series of Portuguese youth soccer players were consistent with recent studies of English [13] and Dutch [14] soccer players. The study of English players was limited to the Mirwald et al. [1] equation, while that of the Dutch players considered the original and modified [2] equations in addition to an equation which predicted a maturity ratio [50]. Although the three studies varied in design, scope and focus, the results were consistent in highlighting major limitations of predicted maturity offset and predicted age at PHV in longitudinal samples of soccer players.

Applications of the original and modified equations in the longitudinal series of Portuguese youth soccer players were also consistent with validation studies of the maturity offset prediction equations in three longitudinal series of youth spanning late childhood through adolescence, one in the U.S. [9] and two in Poland [10,11]. The three studies and the present study of soccer players used similar analytical methods and noted several major limitations of the prediction equations.

First, predicted maturity offset increased, on average, with CA at prediction throughout the age range considered in each study. In the study of soccer players, mean predicted ages at PHV based on the Moore equation increased, on average, with CA, while those based on the Mirwald equation increased from observation one to two, changed negligibly through observation four and then increased to observation five. The age-related trend probably reflects the predictors comprising each equation which increase, on average, with CA. The preceding is apparent in the correlations between predicted maturity offset and predicted age at PHV with CA, height, sitting height, estimated leg length and weight at each observation for the Mirwald equation and with CA and height with the Moore equation (Table 5). For predicted maturity offset, correlations within each CA group range from moderately high to high; correlations are highest for sitting height and tend to be lowest for estimated leg length. For predicted age at PHV, correlations are relatively low and positive for CA, but are negative and moderate to high for the anthropometric variables. Within an age group, taller and heavier players tended to have an earlier predicted age at PHV.

thumbnail
Table 5. Correlations at each observation between predicted maturity offset and predicted APHV (CA–predicted maturity offset) with chronological age (CA), height (Ht), sitting height (SitHt), estimated leg length (LegLt) and body weight (Wt) for the Mirwald equation and with CA and Ht for the Moore equation.

https://doi.org/10.1371/journal.pone.0254659.t005

Second, variation in predicted maturity offset and ages at PHV within CA groups was consistently reduced compared to variation in observed ages at offset and at PHV. Variation was reduced more so with the Moore compared to the Mirwald equation.

Third, predictions varied with maturity status defined by observed age at PHV. Predicted maturity offset and ages at PHV were consistently later than observed maturity offset and age at PHV among early maturing boys, and earlier than observed in late maturing boys. By inference, maturity status defined by observed age at PHV influenced predicted ages at PHV in both early and late maturing boys. It should be noted, however, that Moore et al. [2, p. 1761] cautioned that "Our sample was not large enough to rigorously assess variation in prediction error due to early- and late-maturing children". This caution, however, is overlooked in applications of the equations. In contrast to early and late maturing youth, predicted ages at PHV appeared to be reasonably accurate for average maturing boys within approximately ±1 year of observed PHV. Unfortunately, the maturity status and/or timing of individuals is not ordinarily known in studies applying the prediction protocols.

Fourth, intra-individual variation in predicted ages at PHV with each equation was considerable in the present sample of soccer players and in each of the longitudinal studies. Ranges of predicted ages at PHV were reduced with the Moore compared to the Mirwald equation.

The dependency of predicted maturity offset upon CA and body size at prediction merits attention. Means for predicted maturity offset are plotted relative to means for CA and height in the present study and in samples of male soccer players extracted from the literature [13,38,43,46,5169] are illustrated in Figs 3 and 4. The means plotted in the figures were limited to studies using the Mirwald equation, as it was more widely used in studies of soccer players. Predicted maturity offset increased linearly with CA and with height at prediction. The plotted means were largely based on one or two year CA groups, although several were based on players spanning age ranges of three or more years.

thumbnail
Fig 3. Means for predicted maturity offset (MO) plotted relative to means for chronological age for each year of observation in the current sample of players (filled diamonds) and in samples of male soccer players extracted from the literature (open squares; references: 13, 38, 43, 46, 51–69).

https://doi.org/10.1371/journal.pone.0254659.g003

thumbnail
Fig 4. Means for predicted maturity offset (MO) plotted relative to means for height for each year of observation in the current sample of players (filled diamonds) and in samples of male soccer players extracted from the literature (open squares; references: 13, 38, 43, 46, 51–69).

https://doi.org/10.1371/journal.pone.0254659.g004

Studies reporting maturity offset by relative age categories within an age group and studies classifying players across variable age ranges as pre-, at/circa- or post-PHV based on predicted offset were not included in the graphs. In the former, players born early in the year were, on average, older and taller than those born later in the year, while in the latter, CA and height systematically increased from pre-, to at/circa- to post-PHV groups (see above). Given the trends, studies applying predicted maturity offset as a maturity indicator beg the following question: Is predicted maturity offset an indicator of time before age at PHV or is it an indicator of CA and size at prediction? By inference, the validity of predicted maturity offset as an indicator of the time before or after PHV can be questioned.

The preceding has implications for studies using predicted maturity offset as an indicator of maturity status among youth athletes in soccer and other sports. Predicted maturity offset is used most often to classify youth as pre-, at-/circa- and post- PHV, although mean CAs, heights and weights show, on average, a clear gradient across the respective maturity groups. Many studies simply compare the three groups with analysis of variance without controlling for the variation in CA and body size among groups [31,35,36,40,70]. It is also unclear as to how CA-related variation in predicted offset or ages at PHV was addressed in studies applying the prediction equations in short-term longitudinal studies [7173].

Although studies of youth athletes do not ordinarily indicate the ethnic composition of samples, the issue of ethnic variation is relevant as the maturity offset prediction equations were developed and validated on samples of European ancestry. The original Mirwald equation [1] requires sitting height and estimated leg length, while one of the Moore equations [2] for boys requires age and sitting height. Of potential relevance, population variability in the proportions of sitting height and estimated leg length to standing height is reasonably well established [74,75]. American youth of European (White), African (Black) and Hispanic ancestry, for example, vary in the proportions of sitting height and estimated leg length [76,77]. The proportions of the Portuguese youth soccer players in the present study, as reflected in the sitting height-height ratio, were, on average, generally similar to those for American White and Hispanic youth, but different from American Black youth who have proportionally longer lower extremities. This trend was also noted in a recent study of soccer players in which players of non-European ancestry were taller with a lower sitting height/height ratio, i.e., proportionally longer legs, than players of European ancestry [13].

Results of the current study also have practical implications for those working with youth athletes. The interval of PHV is central to the Long Term Athlete Development (LTAD) model for youth athletes [78,79], which calls for identifying youth of contrasting maturity status, i.e., early, average or late maturing. The LTAD, however, does not specify the method for doing so other than suggesting the monitoring of estimated growth velocities; the latter, however, have limitations over the short term. Estimated increments over short intervals (3–4 months), however, must be interpreted with care; they must be adjusted for the interval between measurements and evaluated relative to factors which influence short term estimates of growth rate–specifically measurement errors (both inter- and intra-observer), diurnal variation and perhaps seasonal variation [80]. Nevertheless, application of predicted maturity offset in this context has the potential for misclassification and thus implications for player development. Those using predicted maturity offset per se or variations of the method to identify when players enter and exit the interval of the adolescent growth spurt should employ these methods with caution. If predicted offset is used to inform training design and prescription, it is essential that variation in chronological age at prediction and error associated with the prediction equations be considered. Perhaps additional or alternative methods might be used as a complement, for example, percentage of predicted adult stature attained at the time of observation. As noted above, the utility of estimated velocities of growth in height based on short term height increments has limitations.

The inability of the maturity offset prediction methods to effectively differentiate between early and late maturing youth implies that they should not be used to group players by maturity status as in bio-banding [81], or to adjust fitness and performance scores to accommodate individual differences in maturation [70,82]. As age at PHV is over-estimated in early and under-estimated in late maturing youth, the majority of these players will likely be categorised as being on time and some will be grouped in equivalent bands. Similarly, maturity associated adjustments to performance or fitness scores in early and late maturing boys will, by virtue of these biases, be attenuated and regress towards a common mean.

Conclusions

Results of the study of applying maturity offset prediction equations to the sample of Portuguese soccer players were consistent with similar studies of soccer players and of the general population. Predicted maturity offset increased, on average, with CA at prediction throughout the age range considered. Variation in predicted maturity offset and ages at PHV within CA groups was consistently reduced compared to variation in observed ages at offset and at PHV. Predictions also varied with maturity status defined by observed age at PHV; predicted maturity offset and ages at PHV were consistently later than observed maturity offset and age at PHV among early maturing boys, and earlier than observed in late maturing boys. And, intra-individual variation in predicted ages at PHV with each equation was considerable and ranges of predicted ages at PHV were reduced more with the Moore compared to the Mirwald equations.

Supporting information

S1 Table.

Mixed effects linear models analysis of differences between predicted ages at PHV with the Mirwald (S1A Table) and Moore equations (S1B Table) and observed age at PHV based on the SITAR model, and likelihood ratio tests for differences between two consecutive models (S1C Table).

https://doi.org/10.1371/journal.pone.0254659.s001

(DOCX)

References

  1. 1. Mirwald RL, Baxter-Jones ADG, Bailey DA, Beunen GP. An assessment of maturity from anthropometric measurements. Med Sci Exerc Sports. 2002; 34:689–994. pmid:11932580
  2. 2. Moore SA, McKay HA, Macdonald H, Nettlefold L, Baxter-Jones ADG, Cameron N, et al. Enhancing a somatic maturity prediction model. Med Sci Sports Exerc. 2015; 47:1755–1764. pmid:25423445
  3. 3. Malina RM. Top 10 research questions related to growth and maturation of relevance to physical activity, performance, and fitness. Res Q Exerc Sport. 2014; 85:157–173. pmid:25098012
  4. 4. Moran J, Sandercock GR, Ramirez-Campillo R, Meylan C, Collison J, Parry DA. A meta-analysis of maturation-related variation in adolescent boy athletes’ adaptations to short-term resistance training. J Sports Sci. 2017; 35:1041–1051. pmid:27454545
  5. 5. Moran J, Sandercock G, Rumpf MC, Parry DA. Variation in responses to sprint training in male youth athletes: A meta-analysis. Int J Sports Med. 2017; 38:1–11. pmid:27793062
  6. 6. Towlson C, Salter J, Ade JD, Enright K, Harper LD, Page RM, et al. Maturity-associated considerations for training load, injury risk, and physical performance in youth soccer: One size does not fit all. J Sport Hlth Sci. 2020; S2095-2546(20)30119-8. pmid:32961300
  7. 7. Malina RM, Kozieł SM. Validation of maturity offset in a longitudinal sample of Polish boys. J Sports Sci. 2014; 32:424–437. pmid:24016098
  8. 8. Malina RM, Kozieł SM. Validation of maturity offset in a longitudinal sample of Polish girls. J Sports Sci. 2014; 32:1374–1382. pmid:24892233
  9. 9. Malina RM, Choh AC, Czerwinski SA, Chumlea WC. Validation of maturity offset in the Fels Longitudinal Study. Pediat Exerc Sci. 2016; 28:439–455. pmid:26757350
  10. 10. Malina RM, Kozieł SM, Králík M, Chrzanowska M, Suder A. Prediction of maturity offset and age at peak height velocity in a longitudinal series of boys and girls. Am J Hum Biol. 2020; Dec 11, e23551. pmid:33314450
  11. 11. Kozieł SM, Malina RM. Modified maturity offset prediction equations: Validation in independent longitudinal samples of boys and girls. Sports Med. 2018; 48:221–236. pmid:28608181
  12. 12. Malina RM, Claessens AL, Van Aken K, Thomis M, Lefevre J, Philippaerts P, et al. Maturity offset in gymnasts: Application of a prediction equation. Med Sci Sports Exerc. 2006; 38:1342–1347. pmid:16826033
  13. 13. Parr J, Winwood K, Hodson-Tole E, Deconick FJA, Parry L, Hill JP, et al. Predicting the timing of the peak of the pubertal growth spurt in elite youth soccer players: Evaluation of methods. Ann Hum Biol. 2020; 47:400–408. pmid:32543933
  14. 14. Teunissen JW, Rommers N, Pion J, Cumming SP, Rössler R, D’Hondt E, et al. Accuracy of maturity prediction equations in individual elite football players. Ann Hum Biol. 2020; 47:409–416. pmid:32996814
  15. 15. Malina RM, Rogol AD, Cumming SP, Coelho-e-Silva MJ, Figueiredo AJ. Biological maturation of youth athletes: Assessment and implications. Br J Sports Med 2015; 49:852–859. pmid:26084525
  16. 16. Myburgh GK, Cumming SP, Malina RM. Cross-sectional analysis investigating the concordance of maturity status classifications in elite Caucasian youth tennis players. Sports Med Open. 2019; 5:27. pmid:31264052
  17. 17. Malina RM, Coelho-e-Silva MJ, Figueiredo AJ, Carling C, Beunen GP. Interrelationships among invasive and non-invasive indicators of biological maturation in adolescent male soccer players. J Sports Sci. 2012; 30:1705–1717. pmid:22304621
  18. 18. Valente-dos-Santos J, Coelho-e-Silva MJ, Simões F, Figueiredo AJ, Leite N, Elferink-Gemser MT, et al. Modeling developmental changes in functional capacities and soccer-specific skills in male players aged 11–17 years. Pediatr Exerc Sci. 2012; 24: 603–621. pmid:23196767
  19. 19. Harris DJ, MacSween A, Atkinson G. Ethical standards in sport and exercise science research: 2020 update. Int J Sports Med. 2019; 40: 813–817. pmid:31614381
  20. 20. Figueiredo AJ, Gonçalves CE, Coelho e Silva MJ, Malina RM. Youth soccer players, 11–14 years: Maturity, size, function, skill and goal orientation. Ann Hum Biol. 2009; 36:60–73. pmid:19085511
  21. 21. Cao Z, Hui LL, Wong MY. IAPVBS: Individual Age at Peak Velocity Based on SITAR. R package version 0.0.2; 2018. https://rdrr.io/github/Zhiqiangcao/iapvbs/.
  22. 22. Cole T. SITAR: Super Imposition by Translation and Rotation Growth Curve Analysis. R package version 1.1.2; 2020. https://CRAN.R-project.org/package=sitar.
  23. 23. Cole TJ, Donaldson MDC, Ben-Shlomo Y. SITAR–a useful instrument for growth curve analysis. Int J Epidem. 2010; 39:1558–1566.
  24. 24. R Core Team. R: A language and environment for statistical computing. Vienna, R Foundation for Statistical Computing; 2020. https://www.R-project.org.
  25. 25. Malina RM, Coelho-e-Silva MJ, Sousa-e-Silva P, Figueiredo AJ, Cumming SP, Kralik M, et al (under review) Age at peak height velocity in soccer players.
  26. 26. Franco-Márquez F, Rodríguez-Rosell D, González-Suárez JM, Pareja-Blanco F, Mora-Custodio R, Yáñez-García , et al. Effects of combined resistance training and plyometrics on physical performance in young soccer players. Int J Sports Med. 2015; 36:906–914. pmid:26180903
  27. 27. Hammami R, Granacher U, Marthlouf I, Behm DG, Chaouachi A. Sequencing effects of balance and plyometric training on physical performance in youth soccer athletes. J Strength Cond Res. 2016; 30:3278–3289. pmid:27144955
  28. 28. John C, Rahlf AL, Hamacher D, Zech A. Influence of biological maturity on static and dynamic postural control among male youth soccer players. Gait Posture. 2019; 68:18–22. pmid:30439683
  29. 29. Rodríguez-Rosell D, Franco-Márquez F, Pareja-Blanco F, Mora-Custodio R, Yáñez-Garciá JM, González-Suárez JM, et al. Effects of 6 weeks resistance training combined with plyometric and speed exercises on physical performance of pre-peak height velocity soccer players. Int J Sport Phys Perf. 2016; 11:240–246. pmid:26218231
  30. 30. Krolo A, Gilic B, Foretic N, Pojskic H, Hammami R, Spasic M, et al. Agility testing in youth football (soccer) players: evaluating reliability, validity and correlates of newly developed testing protocols. Int J Env Res Pub Hlth. 2020; 17:294, 1–15, pmid:31906269
  31. 31. Zago M, Moorhead AP, Bertozzi F, Sforza C, Tarabini M, Galli M. Maturity offset affects standing postural control in youth male soccer players. J Biomech. 2020; 99: 109523, 1–5, pmid:31767282
  32. 32. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team (2020). nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–149. https://CRAN.R-project.org/package=nlme.
  33. 33. Cornbleet PJ, Gochman N. Incorrect least-squares regression coefficients in method—comparison analysis. Clin Chem. 1979; 25:432–438. pmid:262186
  34. 34. Lopes Machada DR, Bonfim MR, Costa LT. Pico de velocidade de crescimento como alternative para classificação maturacional associada ao desempenho motor. Rev Brasi Cineantrop Des Hum. 2009; 11:14–21.
  35. 35. Jakovljevic S, Macura M, Radivoj M, Jankovic N, Pajic Z, Erculj F. Biological maturity status and motor performance in fourteen-year-old basketball players. Int J Morphol. 2016; 34:637–643.
  36. 36. Lopez-Plaza D, Alacid F, Muyor JM, Lopez-Minarro PA. Sprint kayaking and canoeing performance based on the relationship between maturity status, anthropometry and physical fitness in young elite paddlers. J Sports Sci. 2017; 35:1083–1090. pmid:27433884
  37. 37. Read PJ, Oliver JL, de Ste Croix MBA, Myer GD, Lloyd RS. Hopping and landing performance in male youth soccer players: effects of age and maturation. Int J Sports Med. 2017; 38:902–908. pmid:28931173
  38. 38. Read PJ, Oliver JL, de Ste Croix MBA, Myer GD, Lloyd RS. Landing kinematics in elite male youth soccer players of different chronological ages and stages of maturation. J Ath Train. 2018; 53:372–378.
  39. 39. Read PJ, Oliver JL, Myer GD, de Ste Croix MBA, Lloyd RS. The effects of maturation on measures of asymmetry during neuromuscular control tests in elite male youth soccer players. Pediat Exerc Sci. 2018; 30:170–177. pmid:28787266
  40. 40. Peña-González I, Fernández-Fernández J, Cervelló E, Moya-Ramón M. Effect of biological maturation on strength-related adaptations in young soccer players. PloS One, 2019; 14: e0219355. pmid:31276566
  41. 41. Živković M, Stojiljković N, Antić V, Pavlović L, Stanković N, Jorgić B. The motor abilities of handball players of different biological maturation. Facta Universitatis, series: Physical Education and Sport. 2019; 17:125–133.
  42. 42. Mendez-Villanueva A, Buchheit M, Kuitunen S, Poon TK, Simpson B, Peltola E. Is the relationship between sprinting and maximal aerobic speeds in young soccer players affected by maturation? Pediat Exerc Sci. 2010; 22:497–510. pmid:21242600
  43. 43. Buchheit M, Mendez-Villanueva A. Reliability and stability of anthropometric and performance measures in highly-trained young soccer players: effect of age and maturation. J Sports Sci. 2013; 31:1332–1343. pmid:23656211
  44. 44. Cripps AJ, Hopper L, Joyce C. Maturity, physical ability, technical skill and coaches’ perception of semi-elite adolescent Australian footballers. Pediat Exerc Sci. 2016; 28:535–541. pmid:27046936
  45. 45. Hammami R, Chaouachi A, Makhlouf I, Granacher U, Behm DG. Associations between balance and muscle strength, power performance in male youth athletes of different maturity status. Pediat Exerc Sci. 2016; 28:521–534. pmid:27046937
  46. 46. Morris RO, Jones B, Myers T, Lake J, Emmonds S, Clarke ND, et al. Isometric mid-thigh pull characteristics in elite youth male soccer players: Comparisons by age and maturity offset. J Strength Cond Res. 2020; 34:2947–2955. pmid:29985220
  47. 47. Asadi A, Ramirez-Campillo R, Arazi H, Saez de Villareal E. The effects of maturation on jumping ability and sprint adaptations to plyometric training in youth soccer players. J Sports Sci. 2018; 36:2405–2411. pmid:29611771
  48. 48. Brownstein CG, Ball D, Micklewright D, Gibson NV. The effect of maturation on performance during repeated sprints with self-selected versus standardized recovery intervals in youth footballers. Pediat Exerc Sci. 2018; 30:500–505. pmid:30033816
  49. 49. Lloyd RS, Oliver JL, Myer GD, de Ste Croix M, Wass J, Read PJ. Comparison of drop jump and tuck jump knee joint kinematics in elite male youth soccer players: Implications for injury risk screening. J Sports Rehab. 2019; 29:760–765. pmid:31629336
  50. 50. Fransen J, Bush S, Woodcock S, Novak A, Deprez D, Baxter-Jones ADG, et al. Improving the prediction of maturity from anthropometric variables using a maturity ratio. Pediat Exerc Sci. 2018; 30:296–307.
  51. 51. Agostinete RR, Fernandes RA, Narciso PH, Maillane-Vanegas S, Werneck AO, Vlachopolous D. (2020) Categorizing 10 sports according to bone and soft tissue profiles in adolescents. Medicine and Science in Sports 52:2673–2681. pmid:32735110
  52. 52. Aquino R, Alves IS, Padilha MB, Casanova F, Puggina ER, Maia J. (2017) Multivariate profiles of selected versus non-selected elite youth Brazilian soccer players. Journal of Human Kinetics 60:113–121. pmid:29339991
  53. 53. Borges PH, Cumming S, Ronque ERV, Cardoso F, Avelar A, Rechenchosky L, et al. (2018) Relationship between tactical performance, somatic maturity and functional capabilities in young soccer players. Journal of Human Kinetics 64:160–169. pmid:30429908
  54. 54. Campa F, Silva AM, Iannuzzi V, Mascherini G, Benedetti L, Toselli S. (2019) The role of somatic maturation on bioimpedance patterns and body composition in male elite youth soccer players. International Journal of Environmental Research and Public Health 16: 4711, 1–11, pmid:31779215
  55. 55. Deprez DN, Fransen J, Boone J, Lenoir M, Philippaerts RM, Vaeyens R. (2015) Characteristics of high-level youth soccer players: variation by playing position. Journal of Sports Sciences 33:243–254. pmid:24998472
  56. 56. Deprez DN, Fransen J, Lenoir M, Philippaerts RM, Vaeyens R. (2015) A retrospective study of anthropometrical, physical fitness, and motor coordination characteristics than influence dropout, contract status, and first-team playing time in high-level soccer players aged eight to eighteen years. Journal of Strength and Conditioning Research 29:1692–1704. pmid:26010800
  57. 57. Doncaster G, Iga J, Unnithan V. (2018) Assessing differences in cardiorespiratory fitness with respect to maturity status in highly trained youth soccer players. Pediatric Exercise Science 30:216–228. pmid:29276855
  58. 58. Gibson NV, Henning G, Twist C. (2019) Movement characteristics, physiological and perceptual responses of elite standard youth football players to different high intensity running drills. Science and Medicine in Football, pp 1–7
  59. 59. Gibson NV, McCunn R, MacNay SA, Mullen T, Twist C. (2018) Playing exposure does not affect movement characteristics or physiological responses of elite youth footballers during an intensified period of competition. Science and Medicine in Football, pp 1–6,
  60. 60. Gill SM, Zabala-Lili , Bidaurrazaga-Leton I, Aduna B, Lekue JA, Santos-Concejero J, et al. (2014) Talent identification and selection process of outfield players and goalkeepers in a professional soccer club. Journal of Sports Sciences 32:1931–1939. pmid:25429718
  61. 61. Gonçalves RR, Severino V, Coelho-e-Silva MJ, Figueiredo AJ. (2011) Age-related variation in anthropometric and maturity characteristics of soccer goalkeepers 11–14 years. Annals of Research in Sport and Physical Activity (University of Coimbra) 1:71–81,
  62. 62. Hernandez Camacho JD, Huelva Leal AB, Martinez-Sanz JM, Lahoz Ruano MD, Vázquez Carrión J. (2018) Peak height velocity and muscle mass in young soccer players. Revista Española de Nutrición Humana y Dietética 22:219–226.
  63. 63. Laas MM, Wright MD, McLaren SJ, Eaves DL, Parkin G, Portas MD. (2020) Motion tracking in young male football players: A preliminary study of within-session movement reliability. Science and Medicine in Football, pp 1–8,
  64. 64. Lloyd RS, Oliver JL, Radnor JM, Rhodes BC, Faigenbaum AD, Myer GD. (2015) Relationships between functional movement screen scores, maturation and physical performance in young soccer players. Journal of Sports Sciences 33:11–19. pmid:24857046
  65. 65. Lovell TWJ, Bocking CJ, Fransen J, Coutts AJ. (2017) A multidimensional approach to factors influencing playing level and position in a school-based soccer programme. Science and Medicine in Football, pp 1–9,
  66. 66. Lovell TWJ, Bocking CJ, Fransen J, Kempton T, Coutts AJ. (2017) Factors affecting physical match activity and skill involvement in youth soccer. Science and Medicine in Football, pp 1–8,
  67. 67. Moreira A, Mortatti A, Aoki M, Arruda A, Freitas C, Carling C. (2013) Role of free testosterone in interpreting physical performance in elite young Brazilian soccer players. Pediatric Exercise Science 25:186–197. pmid:23504910
  68. 68. Seabra A, Marques E, Brito J, Krustrup P, Abreu S, Oliveira J, et al. (2012) Muscle strength and soccer practice as major determinants of bone mineral density in adolescents. Joint Bone Spine 79:403–408. pmid:22071408
  69. 69. Trecoci A, Longo S, Perri E, Iaia FM, Alberti G. (2019) Field-based physical performance of elite and sub-elite middle-adolescent soccer players. Research in Sports Medicine 33:1232–1236.
  70. 70. Till K, Jones B. Monitoring anthropometry and fitness using maturity groups within youth rugby league. J Strength Cond Res. 2015; 29:730–736. pmid:25226333
  71. 71. Mathys SPH, Vaeyens R, Fransen J, Deprez D, Pion J, Vandendriessche J, et al. A longitudinal study of multidimensional performance characteristics related to physical capacities in youth handball. J Sports Sci. 2013; 31:325–334. pmid:23078540
  72. 72. Till K, Cobley S, O’Hara J, Chapman C, Cooke C. A longitudinal evaluation of anthropometric and fitness characteristics in junior rugby league players considering playing position and selection level. J Sci Med Sport. 2013; 16:438–443. pmid:23072898
  73. 73. Zuber C, Zibung M, Conzelmann A. Holistic patterns as an instrument for predicting the performance of promising young soccer players–a 3 year longitudinal study. Front Psychol. 2016; 7:1088. pmid:27512378
  74. 74. Eveleth PB, Tanner JM. Worldwide Variation in Human Growth. Cambridge: Cambridge University Press, 1976.
  75. 75. Eveleth PB, Tanner JM. Worldwide Variation in Human Growth, 2nd edition. Cambridge: Cambridge University Press, 1990.
  76. 76. Malina RM, Brown KH, Zavaleta AN. Relative lower extremity length in Mexican American and in American Black and White youth. Am J Phys Anthropol. 1987; 72:89–94. pmid:3826332
  77. 77. Martorell R, Malina RM, Castillo RO, Mendoza FS, Pawson IG. Body proportions in three ethnic groups: children and youth 2–17 years in NHANES II and HHANES. Hum Biol. 1988; 60:205–222. pmid:3371962
  78. 78. Balyi I, Hamilton A. Long-Term Athlete Development: Trainability in Childhood and Adolescence–Windows of Opportunity, Optimal Trainability. Victoria, BC: National Coaching Institute British Columbia and Advanced Training and Performance Ltd, 2004.
  79. 79. Balyi I, Cardinal C, Higgs C, Norris S, Way R. Canadian Sport for Life: Long-term athlete development resource paper V2. Vancouver, BC: Canadian Sport Centres, 2005; http://www.canadiansportforlife.ca/default.aspx?PageID=1076&LangID=en.
  80. 80. Malina RM, Bouchard C, Bar-Or O. Growth, Maturation, and Physical Activity, 2nd edition. Champaign, IL: Human Kinetics; 2004.
  81. 81. Malina RM, Cumming SP, Rogol AD, Coelho-e-Silva MJ, Figueiredo AJ, Konarski JM, et al. Bio-banding in youth sport: Background, concept, and application. Sports Med. 2019; 49:1671–1685. pmid:31429034
  82. 82. Abbot S, Hogan C, Castiglioni MT, Yamauchi G, Mitchell LJG, Salter J, et al. Maturity-related developmental inequalities in age-group swimming: The testing of the ‘Mat-CAPS’ for their removal. J Sci Med Sport. 2020; S1440-2440(20)30782-9, pmid:33172611