Department of Mathematics, The Hong Kong University of Science & Technology, Hong Kong SAR, China.
Department of Paediatrics, Faculty of Medicine, The Chinese University of Hong Kong, Hong Kong SAR, China; School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China.
Comput Methods Programs Biomed. 2018 Sep;163:79-85. doi: 10.1016/j.cmpb.2018.05.030. Epub 2018 Jun 1.
We compared three methods for estimating the individual peak height velocity (PHV) and age at peak height velocity (APHV) from the SuperImposition by Translation and Rotation (SITAR) model.
We fitted the SITAR model using simulated data and heights of 12 girls from the Chard Growth Study and obtained individual PHVs and APHVs from three methods: the model method, the quadratic function method and the numerical method, which are available in our newly developed R package"iapvbs". The mean, interquartile range, range of biases in estimated APHV and PHV as well as the rates of warning and unreasonable cases, i.e. estimated APHVs being outside the range of age measurements, from the three methods were presented and compared.
When the growth curves of all individuals were well fitted by the SITAR model, all three methods estimated individual APHVs with similarly small biases, with a few unreasonable cases (0.16%) observed when the model method was used while more computation time required for the numerical method. When the growth curves of some individuals were not very well fitted, the model method generated more unreasonable individual APHV (8.15%) and more bias in PHV and APHV, compared to those estimated by the numerical method and quadratic function method. In line with the observations from the simulated data, the real data analysis demonstrated that the numerical method generated more reliable PHV and APHV for individuals with growth curve not well fitted by the SITAR model.
The performance of different methods estimating individual APHV depends largely on how well the growth curves are fitted by the SITAR model. The quadratic function method is more superior when growth curves of all individuals are well fitted by the SITAR model; otherwise, the numerical method should be adopted for getting most robust estimates of PHV and APHV. The model method generates unreasonable APHV estimates, particularly when the growth curves are not well fitted.
我们比较了三种通过平移和旋转叠加(SITAR)模型估算个体峰值身高速度(PHV)和峰值身高速度年龄(APHV)的方法。
我们使用模拟数据和 Chard 生长研究中的 12 名女孩的身高拟合 SITAR 模型,并从三种方法中获得个体 PHV 和 APHV:模型法、二次函数法和数值法,这些方法均可在我们新开发的 R 包“iapvbs”中使用。展示并比较了三种方法的估计 APHV 和 PHV 的均值、四分位间距、偏差范围以及警告和不合理案例的比例,即估计的 APHV 超出年龄测量范围的情况。
当所有个体的生长曲线都被 SITAR 模型很好地拟合时,所有三种方法估算的个体 APHV 偏差都很小,使用模型法时观察到少数不合理案例(0.16%),而数值法需要更多的计算时间。当一些个体的生长曲线拟合不太好时,与数值法和二次函数法相比,模型法生成的个体 APHV 更不合理(8.15%),PHV 和 APHV 的偏差也更大。与模拟数据的观察结果一致,真实数据分析表明,对于 SITAR 模型拟合不佳的个体,数值法生成的 PHV 和 APHV 更可靠。
不同方法估算个体 APHV 的性能在很大程度上取决于 SITAR 模型对生长曲线的拟合程度。当所有个体的生长曲线都被 SITAR 模型很好地拟合时,二次函数法更优越;否则,应采用数值法以获得最稳健的 PHV 和 APHV 估计值。模型法生成的 APHV 估计值不合理,特别是当生长曲线拟合不佳时。