To model cross-sectional growth data the LMS method is widely applied. In this method the distribution is summarized by three parameters: the Box-Cox power that converts outcome to normality (L); the median (M); and the coeficient of variation (S). Here, we propose an alternative approach based on fitting finite mixture models with several components which may perform better than the LMS method in case the data show an unusual distribution. Further, we explore fixing the weights of the mixture components in contrast to the standard approach where weights are estimated. Having fixed weights improves the speed of computation and the stability of the solution. In addition, fixing the weights provides almost as good a fit as when the weights are estimated. Our methodology combines Gaussian mixture modelling and spline smoothing. The estimation of the parameters is based on the joint modelling of mean and dispersion. We illustrate the methodology on the Fourth Dutch Growth Study, which is a cross-sectional study that contains information on the growth of 7303 boys as a function of age. This information is used to construct centile curves, so-called growth curves, which describe the distribution of height as a smooth function of age. Further, we analyse simulated data showing a bimodal structure at some time point. In its full generality, this approach permits the replacement of the Gaussian components by any parametric density. Further, different components of the mixture can have a diferent probabilistic (multivariate) structure, allowing for censoring and truncation.
@article{bwmeta1.element.doi-10_2478_bile-2013-0008,
author = {Marek Molas and Emmanuel Lesaffre},
title = {Finite mixture models with fixed weights applied to growth data},
journal = {Biometrical Letters},
volume = {49},
year = {2012},
pages = {103-119},
language = {en},
url = {http://dml.mathdoc.fr/item/bwmeta1.element.doi-10_2478_bile-2013-0008}
}
Marek Molas; Emmanuel Lesaffre. Finite mixture models with fixed weights applied to growth data. Biometrical Letters, Tome 49 (2012) pp. 103-119. http://gdmltest.u-ga.fr/item/bwmeta1.element.doi-10_2478_bile-2013-0008/
Cole T.J., Green P.J. (1992): Smoothing reference centile curves: The LMS method and penalized likelihood. Statistics in Medicine 11: 1305-1319.[Crossref]
Dempster A.P., Laird N.M., Rubin D.B. (1977): Maximum likelihood from incomplete data via the EM algorithm (with discussion). Journal of the Royal Statistical Society Series B 39: 1-38. | Zbl 0364.62022
Eilers P., Marx B. (1996): Flexible smoothing with b-splines and penalties. Statistical Science 11: 89-121.[Crossref] | Zbl 0955.62562
Ghidey W., Lesa re E., Eilers P. (2004): Smooth random e ects distribution in a linear mixed model. Biometrics 60: 945-953.[Crossref] | Zbl 1274.62238
Harrell F.E. (2001): Regression Modelling Strategies. Springer-Verlag, New York. | Zbl 0982.62063
Komarek A., Lesa re E., Hilton J.F. (2005): Accelerated failure time model for arbitrarily censored data with smoothed error distribution. Journal of Computational and Graphical Statistics 14: 726-745.[Crossref]
Lee Y., Nelder J.A., Pawitan Y. (2006): Generalized Linear Models with Random E ects. Chapman & Hall / CRC: Boca Raton. | Zbl 1110.62092
McLachlan G.J., Peel D. (2000): Finite Mixture Models. John Wiley and Sons, New York. | Zbl 0963.62061
Muthen B., Brown H.C. (2009): Estimating drug e ects in the presence of placebo response: Casual inference using growth mixture modelling. Statistics in Medicine 28: 3363-3385.[WoS][Crossref]
Nelder J.A., Pregibon D. (1987): An extended quasi-likelihood function. Biometrika 74: 221-232.[Crossref] | Zbl 0621.62078
Nelder J.A., Wedderburn R.W.M. (1972): Generalized linear models. Journal of Royal Statistical Society A 135: 370-384.
Ramsay J.O. (1988): Monotone regression splines in action. Statistical Science 3: 425-461.[Crossref]
van Buuren S., Fredriks M. (2001): Worm plot: a simple diagnostic device for modelling growth reference curves. Statistics in Medicine 20: 1259-1277. [Crossref]