Evaluation of robust Bayesian mixed-effects models of longitudinal childhood BMI

This is one of the first studies to rigorously evaluate Fractional Polynomials (FP), Reed2 (R2), and Breakpoint (BP) models using data from three longitudinal cohorts, assessing performance based on bias, heteroskedasticity, and Bayesian goodness-of-fit.</p> <p>Results indicate that model choice matters, as FP and Reed2 models generally outperform BP models in capturing the complex non-linearity of BMI trajectories, with residual errors accounting for less than 10% of BMI variance.

Among five likelihoods considered, the Laplace distribution provided the most robust fit. The residual correlation structures were complex and non-stationary. A normal likelihood with a stationary correlation function improved fit, but less so than an independent Laplace likelihood.