Cole Tim J, Cortina-Borja Mario, Sandhu Jat, Kelly Frank P, Pan Huiqi
Centre for Paediatric Epidemiology and Biostatistics, UCL Institute of Child Health, London WC1N 1EH, UK.
Biostatistics. 2008 Jan;9(1):159-71. doi: 10.1093/biostatistics/kxm020. Epub 2007 Jun 16.
Higher moments of the frequency distribution of child height and weight change with age, particularly during puberty, though why is not known. Our aims were to confirm that height skewness and kurtosis change with age during puberty, to devise a model to explain why, and to test the model by analyzing the data longitudinally. Heights of 3245 Christ's Hospital School boys born during 1927-1956 were measured twice termly from 9 to 20 years (n=129508). Treating the data as independent, the mean, standard deviation (SD), skewness, and kurtosis were calculated in 40 age groups and plotted as functions of age t. The data were also analyzed longitudinally using the nonlinear random-effects growth model H(t)=h(t-epsilon )+alpha, with H(t) the cross-sectional data, h(t) the individual mean curve, and epsilon and alpha subject-specific random effects reflecting variability in age and height at peak height velocity (PHV). Mean height increased monotonically with age, while the SD, skewness, and kurtosis changed cyclically with, respectively, 1, 2, and 3 turning points. Surprisingly, their age curves corresponded closely in shape to the first, second, and third derivatives of the mean height curve. The growth model expanded as a Taylor series in epsilon predicted such a pattern, and the longitudinal analysis showed that adjusting for age at PHV on a multiplicative scale largely removed the trends in the higher moments. A nonlinear growth process where subjects grow at different rates, such as in puberty, generates cyclical changes in the higher moments of the frequency distribution.
儿童身高和体重频率分布的高阶矩随年龄变化,尤其是在青春期,但其原因尚不清楚。我们的目的是确认青春期期间身高偏度和峰度随年龄变化,设计一个模型来解释原因,并通过纵向分析数据来检验该模型。对1927年至1956年出生的3245名基督医院学校男生的身高进行测量,从9岁到20岁每学期测量两次(n = 129508)。将数据视为独立数据,计算40个年龄组的均值、标准差(SD)、偏度和峰度,并将其绘制为年龄t的函数。还使用非线性随机效应生长模型H(t)=h(t - ε)+α对数据进行纵向分析,其中H(t)为横断面数据,h(t)为个体平均曲线,ε和α为反映身高增长峰值速度(PHV)时年龄和身高变异性的个体特异性随机效应。平均身高随年龄单调增加,而标准差、偏度和峰度分别有1个、2个和3个转折点,呈周期性变化。令人惊讶的是,它们的年龄曲线在形状上与平均身高曲线的一阶、二阶和三阶导数密切对应。以ε的泰勒级数展开的生长模型预测了这种模式,纵向分析表明,在乘法尺度上对PHV时的年龄进行调整,在很大程度上消除了高阶矩的趋势。在青春期等受试者以不同速度生长的非线性生长过程中,会在频率分布的高阶矩中产生周期性变化。