Discussion
The formula proposed by Tanner in 1970 has been widely used by paediatric endocrinologists worldwide to estimate TH in patients for many decades. However, certain factors can potentially affect the accuracy of applying this formula to different demographic groups. Tanner’s formula was developed in the 1960s, and over the past century, FAH has increased in most populations of both sexes.7 These changes are widely observed across generations and countries and are thought to be mainly influenced by nutrition and health improvements.8–11 Therefore, a high possibility of underestimation is present when Tanner’s formula is used. Furthermore, differences between ethnicities may exist when Tanner’s formula is used to predict FAH. In addition, if the difference between both parents’ heights is too large or small, the accuracy in predicting the FAH of the offspring using Tanner’s formula decreases.12
Our data also showed that the TH calculated using Tanner’s formula underestimated the FAH in both sexes; therefore, the usefulness of this formula should be re-evaluated. However, factors that may influence FAH are quite different from those identified in previous studies. Economic status has influenced FAH directly and indirectly in previous studies.13–15 However, in this study, economic status was not associated with FAH. Nevertheless, we believe that educational differences reflected some aspects of economic influence on FAH because individuals with high economic status in Korea tend to be highly educated. When Tanner’s formula was applied, economic factors, nutritional factors, BMI and obesity were not associated with FAH or a positive height gap.
Many studies have shown that chronic disorders,16 such as inflammatory bowel disease,17–19 hepatitis,20 21 pulmonary tuberculosis,22 23 systemic lupus erythematosus,24 25 allergic disease (including atopic dermatitis),26 27 asthma,28 29 nephrotic syndrome,30 chronic kidney disease or hypothyroidism,31 which require high-dose steroid treatment, affect FAH. We compared these variables in our study but found no association between calculated height differences and FAH. Using data representing the entire Korean population, this study shows a realistic pattern of growth and relationship with various factors, suggesting that these factors may not significantly influence FAH.
To account for sex differences, Tanner’s formula was added to a constant value of 13 cm.
There have been varying opinions among studies regarding whether the sex difference is exactly 13 cm or 14 cm. However, on examining previous studies2 32 33 and our own research within the same era, the sex height difference was approximately 13 cm. According to the 2017 Korean paediatric growth chart, the 50th percentile height for 18-year-old men is 174.5 cm, while that for women is 161.1 cm, resulting in a sex-based height difference of 13.4 cm.34 Our study also revealed that parents and offspring had a sex-based height difference of 13 cm. In the newly proposed formula, the sex correction remains the same as that in Tanner’s formula at 13.013, with only the slope reduced.
In this study, the linear regression slope decreased from 0.5 to 0.369 in Tanner’s formula compared with our new formula (figure 2). Though it may not be relevant to compare directly because Tanner’s formula is not regression-based, there were many efforts to reduce the gap of 0.5 in clinical settings to reduce the error of predicted potential height of the offspring. This contrasts with findings of the previous Indian study, which had a slope of 0.615 for sons and daughters, as well as studies from Australia32 and Sweden,12 showing slopes of 0.78 in sons and 0.75 in daughters, respectively. The reduction in parental height’s influence on Korean offspring’s height suggests other factors (including environmental factors) exert an important influence on FAH. Unfortunately, our attempt to identify these other factors, including nutritional status, chronic illness and age at menarche, failed to reveal significant associations. Although genetic factors, such as parental height, still play a substantial role in children’s height, further research is required to identify other influential factors, as indicated by the new formula.
Our study encompassed data representative of the entire South Korean population and benefited from the participation of skilled experts in the research survey. In addition, this study has the advantage that KNHANES was conducted in an ethnically and culturally homogeneous region, lending credibility to the measured values. However, our study has some limitations. Owing to the cross-sectional nature of the study, it evaluates only participants’ immediate status, lacking a comprehensive view of the growth process from birth. Furthermore, medical history data such as chronic illness and age at menarche rely on participant’s recall, leading to the possibility for errors. Furthermore, participants enrolled in this study were those who had complete data of both parents’ height, and only 27.8% of initial candidates finally underwent analysis. The reason of missing data is various; usually due to participants’ preference in the process of consent and participants’ time affordability would be the major reason. These processes were purely random, however, the percentage of missing data is quite large, and this also should be considered in predicting this result. Therefore, further research should be proceeded not only explore factors known to impact FAH but also other environmental factors in the sequential growth in the well-controlled cohort. Additionally, there should be an assessment of the impact of endocrine-disrupting chemicals on growth, with a focus on regional differences among participants.
Finally, KNHANES is a nationally representative large-scale survey sample data in Korea with high variation on the dependent variable. Therefore, the low goodness of fit of our model with one or two independent variables is a limitation. In addition, the normality assumption was also violated, and transforming the outcome to normality was unsuccessful. We also analysed quantile regression (data not shown). It was found that the estimation of the regression parameters for sex and sum of the parental heights was significant at level α=0.05 for all quantiles.
In the 0.25th quantile regression model, sex and sum of the parental heights showed significant influences (p<0.001). The regression equation based on 0.25th quantile regression model for boys was Y=46.029+0.383X, and for girls was Y=33.091+0.383X. In the 0.50th quantile regression model, sex and sum of the parental heights showed significant influences (p<0.001). The regression equation based on 0.50th quantile regression model for boys was Y=50.217+0.381X, and for girls was Y=37.304+0.381X. In the 0.75th quantile regression model, sex and the sum of the parental heights showed significant influences (p<0.001). The regression equation based on 0.75th quantile regression model for boys was Y=62.254+0.354X, and for girls was Y=49.274+0.354X.
However, the aim of this study is to model outcomes with a simple formula and compare with Tanner’s formula directly. Therefore, results should be interpreted with caution and focus on estimation of linear tendency and compare with Tanner’s formula rather than prediction of each data point.
We also could not validate the new formula out of the KNHANES cohort. In the KNHANES cohort the new formula fits better than Tanner’s model, but further validation process also need to be used in clinical settings.