Application of multilevel linear spline models for analysis of growth trajectories in a cohort with repeat antenatal and postnatal measures of growth: a prospective cohort study ================================================================================================================================================================================= * Linda M. O’Keeffe * Cara A. Yelverton * Helena C. Bartels * Kate N. O’Neill * Ciara McDonnell * Fionnuala M. McAuliffe ## Abstract **Introduction** Antenatal and postnatal growth are important indicators of fetal and child health and development. Studies frequently have repeat antenatal and postnatal measures of growth available and require approaches that can maximise the use of these measures to examine growth trajectories. We demonstrate the use of multilevel linear spline modelling to model growth trajectories with repeated antenatal and postnatal measures of growth from 20 weeks gestation to five years in a cohort at high risk of macrosomia. **Methods** Prospective follow-up data from 720-759 mother-child pairs from the ROLO study (initially a randomized controlled trial of a low glycemic index diet in pregnancy to prevent recurrence of macrosomia [birthweight > 4K]) were analysed. Fetal measurements were obtained from ultrasound scans performed on mothers at 20-and 34-weeks gestation, including abdominal circumference (AC) and head circumference (HC). An estimated fetal weight was obtained at 20-and 34-weeks gestation, calculated using the Hadlock 4-parameter formula. At delivery, AC, HC, weight and length were recorded. Follow-up anthropometry assessments (AC, HC, weight and length/height) were also obtained in childhood at six months, two years and five years. Linear spline multilevel models were used to examine trajectories of AC, HC and weight from 20 weeks gestation to five years and length/height from birth to five years. **Results** 754, 756 and 759 participants were included in analyses of AC, HC and weight respectively, while 720 participants were included in analysis of length/height. Over 50% of women had 3rd level education and over 90% were of White ethnicity. Women were a mean (SD) age of 32 (4.2) at recruitment. Following exploration of a series of different models for each growth measure, the best fitting model for AC, HC and weight included a model with knots at each measurement occasion giving rise to five linear spline periods from: 20 weeks to 34 weeks gestation, 34 weeks gestation to birth, birth to six months, six months to two years and two years to five years. The best fitting models for length/height included a model with three linear spline periods from birth to six months, six months to two years and two years to five years. Comparison of observed and predicted values for each model demonstrated good model fit. For all growth measures, fetal growth rates were generally fastest in pregnancy or immediately postpartum (for length/height), with rates of growth slowing after birth and becoming slower still as infancy and childhood progressed. We found little difference in growth trajectories between the intervention and control group. There was some evidence of slightly lower HC, weight and length among females compared with males at birth which appeared to widen by age five years due to slower postnatal growth rates among females. **Conclusion** We demonstrate the application of multilevel linear spline models for examining growth trajectories when both antenatal and postnatal measures of growth are available. The approach may be useful for cohort studies or randomised controlled trials with repeat prospective assessments of growth spanning pregnancy and childhood. ## Introduction Antenatal and childhood growth are important indicators of fetal and child health and development and are associated with health in adult life (1, 2). Consequently, modelling of growth trajectories, identifying causes and predictors of different growth trajectories and relating growth trajectories in the early life course to later life health is important for informing a life course approach to disease prevention (3-5). A key aspect of understanding growth patterns, their causes, predictors and outcomes includes appropriate modelling of longitudinal growth data (3). Since repeated measures of growth within individuals are not independent of each other and the scale and variance of growth measures often changes over time, traditional approaches to analysis of growth data, such as Z-score based methods analysed using multiple regression, do not take account of the clustering of repeated measures within individuals (3). Moreover, the true shape of growth trajectories cannot be modelled using such approaches. While appropriate methods for the study of longitudinal growth data have been applied to antenatal and childhood growth measures in many cohort studies, most studies to date have examined antenatal growth (6, 7) or postnatal growth as separate processes/trajectories (8-14). Appropriate modelling of growth data as a continuum from antenatal to postnatal life is important to accurately characterise the shape of growth from early gestation into childhood to better understand it’s aetiology. In addition, it also allows such trajectories to be examined as outcomes for pre-conception or early pregnancy exposures or to be examined themselves as exposures for later health outcomes (3). Using data from the prospective follow-up of a randomised controlled trial of a low glycaemic index diet in pregnancy (ROLO study), we demonstrate the application of linear spline multilevel models for modelling antenatal and postnatal growth trajectories using four measures of anthropometry (abdominal circumference [AC], head circumference [HC], weight and length/height) from 20 weeks’ gestation to age five years. ## Methods ### Study population The ROLO study is a randomised control trial of a low glycaemic index diet in pregnancy that recruited 800 secundigravid women who had previously given birth to a baby weighing over 4kg between 2007-2011 at the National Maternity Hospital, Dublin, Ireland (15). Women were recruited at first antenatal consultation. Women with any underlying medical disorders, including a previous history of gestational diabetes, those on any drugs, those unable to give full informed consent, aged less than 18 years, of gestation greater than 18 weeks, and having multiple pregnancies were excluded. Women were randomised to either the intervention group which received dietary advice on a low glycaemic diet, or the control group who received routine antenatal care. ### Measurement of anthropometry #### Antenatal measures Fetal measurements were obtained from ultrasound scans performed on mothers at medians of 20 + 6 (Interquartile Range [IQR]: 20 + 1 to 21 + 5) and 34 + 1 (IQR: 33 + 5 to 34 + 5) weeks’ gestation, including AC and HC. An estimated fetal weight (EFW) at 20- and 34-weeks’ gestation was calculated using the Hadlock 4-parameter formula. Ultrasound measurements were taken by two ultrasonographers using a Voluson 730 Expert (GE Medical Systems, Germany) using standard procedures. #### Postnatal measures At delivery, infants’ AC, HC, weight and length were recorded. Follow-up anthropometry assessments were also obtained in childhood at six months, two years and five years (15-17). All measurements were obtained and calculated by a trained member of the research team. At six months, two years and five years, weight (kg) of the child was measured using a calibrated stand on digital weighing scale (SECA 813) to the nearest 0.1 kg by a trained research team member. Children were measured in light clothing without shoes. Standing height was measured, without shoes, with head aligned in the Frankfort plain, using a free-standing stadiometer (SECA 217) and measurements recorded to the nearest 0.1cm. The child’s head and abdominal circumferences were measured using a SECA ergonomic circumference measuring tape, to the nearest 0.1cm. All measurements were recorded three times and the average calculated to improve reliability. ### Statistical analysis We used multilevel models to examine trajectories of change in AC, HC, weight and length/height from 20 weeks gestation to age five years (18, 19). Multilevel models estimate mean trajectories of the outcome while accounting for the non-independence (i.e. clustering) of repeated measurements within individuals, change in scale and variance of measures over time and differences in the number and timing of measurements between individuals (using all available data from all eligible participants under a Missing at Random [MAR] assumption) (3, 20). Change in all four growth measures was estimated here using linear spline multilevel models (two levels: measurement occasion and individual) (3). Linear splines allow knot points to be fit at different ages to derive periods in which change is approximately linear. The optimal linear spline model for each growth measure was selected by examining observed data for each growth measure and comparing model fit statistics of different models including models that assumed linear change over time to models with knot points at different ages (strategies for selection of knot points are described elsewhere in detail (3)). Model fit statistics examined included Akaike’s Information Criterion and observed and predicted values of each growth measure across the age range of the model. Following exploration of a series of models, the best fitting model for AC, HC and weight included a model with knots at each measurement occasion giving rise to five linear spline periods from 20 weeks’ to 34 weeks’ gestation, 34 weeks’ gestation to birth, birth to six months, six months to two years and two years to five years while the best fitting models for length/height included a model with three linear spline periods from birth to six months, six months to two years and two years to five years. All outcomes were normally distributed at each measurement occasion. Except for length/height which did not include antenatal measures, trajectories were centred on the first available measure (20 weeks gestation) for AC, HC and weight. Length/height trajectories were centred at birth. For all models we placed no restrictions on the variance-covariance matrices of level two (individual level) random effects. Given the substantial change in scale and variance of growth from antenatal to postnatal life, we also aimed to allow occasion level measurement error to vary with age (level one random effects for the slope). Therefore, all models included a level one random effect for the slope while the HC model also included a level one random effect for the intercept. The final models for growth trajectories from 20 weeks gestation took the following form: ACij / HCij /weightij= β + u0j + (β1+ u1j)sij1 + (β2+ u2j)sij2 + (β3 + u3j)sij3 + (β4 + u4j)sij4 + (β5 + u5j)sij5 + eij where for person j at measurement occasion i; β represents the fixed effect coefficient for the average intercept, β1 to β5 represent fixed effect coefficients for the average linear slopes of each linear spline, u0j to u5j indicate person-specific random effects for the intercept and slopes respectively, and eij represents the occasion-specific residuals or measurement error which were allowed to vary with age. The final model for length took a similar form but with only three linear spline periods due to the absence of measures prior to birth. Code for the application of these models using the “runmlwin” command from MlWin (21) within Stata 16 (22) is included in Supplementary Material. ## Results 754, 756 and 759 offspring were included in analyses of AC, HC, and weight respectively while 720 offspring were included in analyses of length/height. Table 1 includes the number of measures of each growth measure at each measurement occasion with number of measures available broadly similar across growth measures; for example, weight measures available on each occasion included 655 measures at 20 weeks gestation, 730 at 34 weeks gestation, 756 at birth, 280 at six months, 339 at two years and 387 at five years. View this table: [Table 1](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T1) Table 1 N repeated measures included in analyses for each growth measure Of participants included in analyses (Table 2), over 50% had completed third level education and a majority (>90%) were of White ethnicity. Among mothers of male babies, mean age (standard deviation (SD)) at delivery was approximately 32.3 (4.2) years, mean (SD) BMI at delivery was 27.1 (5.2) kg/m2, mean (SD) birthweight at delivery was 4.1 (0.5) kg and median (interquartile range (IQR)) gestational age was 40.4 (39.6, 41.1) weeks. Mothers of male babies had relatively low levels of deprivation as indicated by the mean (SD) Pobal HP index of 5.3 (10.8). Characteristics were broadly similar for mothers of female babies though mothers of female babies had somewhat higher levels of third level education (∼60%). Model fit as judged by differences between observed growth measures and those predicted by the models for AC, HC, weight and length are shown in Tables 3-6. Overall, our models have good model fit as all reference ranges for the difference between observed and predicted are less than the SD of the observed or less than 10% of the observed value which can be used as a rule of thumb for the assessment of model fit. View this table: [Table 2](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T2) Table 2 Characteristics of ROLO participants included in the analysis of length/height, by sex View this table: [Table 3](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T3) Table 3 Model details for abdominal circumference View this table: [Table 4](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T4) Table 4 Model details for head circumference View this table: [Table 5](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T5) Table 5 Model details for weight View this table: [Table 6](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T6) Table 6 Model details for length Trajectories of AC, HC and weight from 20 weeks’ gestation to five years and trajectories of length/height from birth to five years by intervention status and sex are shown in Table 7 and Figures 1-4. AC and HC had the fastest rates of growth from 20 to 34 weeks’ gestation with growth rates continuing to slow thereafter up to age five years. Weight had the fastest growth rate from 34 weeks’ gestation to birth with growth rates slowing somewhat from birth to six months and continuing to slow thereafter until five years. Length/height had the fastest growth rates from birth to two years, with the growth rate decreasing thereafter and slowing further from two years to five years. View this table: [Table 7](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/T7) Table 7 Mean trajectories of anthropometry and mean difference in trajectories by intervention status and sex in the ROLO cohort ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/05/2022.05.03.22274407/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/F1) Figure 1 Trajectories of abdominal circumference from 20 weeks gestation to age five years by intervention status and sex. **Legend:** Note that X axis displays time in months because trajectory spans the antenatal and postnatal period. ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/05/2022.05.03.22274407/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/F2) Figure 2 Trajectories of head circumference from 20 weeks gestation to age five years by intervention status and sex. **Legend:** Note that X axis displays time in months because trajectory spans the antenatal and postnatal period. ![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/05/2022.05.03.22274407/F3.medium.gif) [Figure 3](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/F3) Figure 3 Trajectories of weight from 20 weeks gestation to age five years by intervention status and sex. **Legend:** Note that X axis displays time in months because trajectory spans the antenatal and postnatal period. ![Figure 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/05/2022.05.03.22274407/F4.medium.gif) [Figure 4](http://medrxiv.org/content/early/2022/05/05/2022.05.03.22274407/F4) Figure 4 Trajectories of length/height from birth to age five years by intervention status and sex. **Legend:** Note that X axis displays time in months. We found no strong evidence of differences in trajectories of AC, weight and length/height between the intervention and control group, but we found some evidence of slightly greater HC (difference 0.27 cm (95% Confidence Interval [CI] = 0.03, 0.51) emerging among the control group at five years. AC trajectories did not differ between males and females, though we found some evidence of modest differences in HC, weight and length/height trajectories between males and females. Females had lower HC at 20 weeks gestation with this difference widening at birth and persisting at age five years (difference at five years: - 0.91cm, 95% CI = -1.14, -0.68). Females had -0.15kg (95% CI = -0.21, -0.08) lower birth weight and slower postnatal growth rates in weight leading to -0.50 kg (95% CI = -0.96, - 0.05) lower weight among females at five years. Similarly, females were -0.83 cm (95% CI = - 1.17, -0.48) shorter in length at birth and had slower postnatal growth rates in length/height leading to -1.22 cm (95% CI = -2.01, -0.43) shorter height among females at five years. ## Discussion In this prospective follow-up of a randomised control trial of approximately 750 infants at high risk of macrosomia, we demonstrated the use of linear spline multilevel modelling to examine trajectories of AC, HC, weight and length/height from 20 weeks’ gestation to age five years. We showed their applicability to data with repeated measures of growth which span the antenatal and postnatal period, even when as few as four repeat assessments are available (in the case of length/height) and measures are sparse. This work may be of value to other studies including randomised control trials with follow-up data such as ours in demonstrating the application of a multilevel modelling approach to examine growth trajectories which can subsequently be used as exposures or outcomes to better understand determinants and outcomes of growth in early life. There are several strengths to the approach used here including ability to maximise sample sizes for analyses and reduce selection bias compared with traditional Z-score approaches since multilevel models can include all participants with at least one growth period under a MAR assumption (3). This is particularly advantageous where attrition rates from cohorts are high. Further advantages include more precise standard errors which consider the non-independence of repeated measures and here we have shown that the approach is implementable with as little as four repeated measures and with repeated measures that span antenatal and postnatal life. Limitations of this work include an inability to explore other non-linear growth patterns such as fractional polynomials due to the sparsity of measures which did not allow a range of possible shapes of growth trajectories to be explored (3). In cohorts with greater numbers of repeated measures and density of repeats, linear spline multilevel modelling can be implemented and compared to other possible shapes include fractional polynomials which have been shown to provide a more biologically intuitive shape of change(3). However, the linear spline approach demonstrated here provides many practical advantages including being more easily interpretable, allowing analysts to split trajectories into distinct periods of change that can then be easily related to exposures and outcomes. It should be noted that this cohort are unlikely to represent the growth rates or trajectories of a general population since their development is above average compared to what would be expected from an age and gender matched general population (the cohort is roughly approximated to the 75th centile based on a crude comparison of means and SDs on the UK-WHO (Ireland) chart) (23). ## Conclusion We demonstrate the application of multilevel linear spline models for examining growth trajectories when both antenatal and postnatal measures of growth are available. The approach may be useful for cohort studies or randomised control trials with repeat prospective assessments of fetal growth spanning pregnancy and childhood. ## Data Availability All data produced in the present study are available upon reasonable request to the authors. ## Supplementary Material ### Sample code for implementing linear spline multilevel models using “runmlwin” command This syntax utilises the user-written command ‘runmlwin’ which must be installed prior to use. The most recent version of MLwiN must be installed to be able to use this command and this package is available for use within Stata and R. Below we demonstrate the basic steps involved in implementing linear spline multilevel modelling suing “runmlwin” in Stata. Code below assumes data are in long format and that a variable called “occasion” exists identifying the ordering of observations within individuals. Sample code below applies to length/height from birth to five years. ### Generate the spline variable First, three new variables are created: s1 (spline 1 from birth to 6 months), s2 (6 months to 2 year), s3 (2 years to 5 years). mkspline s1\_birth\_6m 27 s2\_6m\_2 107 s3\_2\_max = age_lw ### Generate a constant term MLwiN does not automatically include a constant term, so this must be generated and included in models. gen cons=1 ### Identify the location of MLwiN global MLwiN_path “C:\Program Files\MLwiN v3.05\mlwin.exe” ### Run the multilevel model, sorting the data by person and occasion/age first sort study_id age runmlwin length cons s1\_birth\_6m s2\_6m\_2 s3\_2\_max /// level2 (study\_id: cons s1\_birth\_6m s2_6m_2 s3_2_max, reset(var) residuals (res, var)) /// level1 (occ: age_lw, reset(var) diag) nopause maxiterations(150) ### Adding covariates The following assumes covariates are binary and coded 0 and 1 or for covariates with multiple categories, dummy variables have been created. The addition of continuous covariates should be undertaken in the same manner as for categorical covariates but continuous covariates should be centred on the mean so that the baseline trajectory in the model is for the individuals with the mean level of the continuous covariate. Here we demonstrate the steps required for addition of sex as a covariate. ### Multiply covariate by splines Once the covariate is coded in the format of 0/1 representing 0 for the baseline category, we multiply the covariate by the splines, creating interaction terms for inclusion in our model. gen s1\_birth\_6m\_fem = s1\_birth_6m*female gen s2\_6m\_2\_fem = s2\_6m_2*female gen s3\_2\_max\_fem = s3\_2_max*female ### Run model now including covariate terms The model is then ran as before but this time including a term for the covariate in question, here “female” and each of the above female*spline interaction terms generated. This allows the mean trajectory to differ for females and males. Because in this example the variable female is coded 0 for male and 1 for female the baseline trajectory is now for males with coefficients for “female”, s1\_birth\_6m\_fem, s2\_6m\_2\_fem, s3\_2\_max_fem representing the difference in the intercept, spline 1 and spline 2 and spline 3 in females compared with males. sort study_id age runmlwin length cons s1\_birth\_6m s2\_6m\_2 s3\_2\_max female2*, /// level2 (study\_id: cons s1\_birth\_6m s2_6m_2 s3_2_max, reset(var) residuals (res, var)) /// level1 (occ: age_lw, reset(var) diag) nopause maxiterations(150) ## Acknowledgements The authors would like to thank all the ROLO participants for their involvement and all the staff of the National Maternity Hospital and the Perinatal Research Centre. ## Footnotes * **Funding:** LMOK and KNON are supported by a Health Research Board of Ireland Emerging Investigator Award (EIA-FA-2019-007 SCaRLeT). The ROLO study was funded by the Health Research Board of Ireland, Health Research Centre for Health and Diet Research, and the European Union’s Seventh Framework Programme (FP7/2007-2013), project Early Nutrition under grant agreement no. 289346. * **Disclosures:** None of the authors have any conflicts of interest to declare. * Received May 3, 2022. * Revision received May 3, 2022. * Accepted May 5, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Kelishadi R, Haghdoost AA, Jamshidi F, Aliramezany M, Moosazadeh M. Low birthweight or rapid catch-up growth: which is more associated with cardiovascular disease and its risk factors in later life? A systematic review and cryptanalysis. Paediatrics and international child health. 2015;35(2):110–23. 2. 2.Baird J, Fisher D, Lucas P, Kleijnen J, Roberts H, Law C. Being big or growing fast: systematic review of size and growth in infancy and later obesity. Bmj. 2005;331(7522):929. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjEyOiIzMzEvNzUyMi85MjkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNS8wNS8yMDIyLjA1LjAzLjIyMjc0NDA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 3. 3.Howe LD, Tilling K, Matijasevich A, Petherick ES, Santos AC, Fairley L, et al. Linear spline multilevel models for summarising childhood growth trajectories: a guide to their application using examples from five birth cohorts. Statistical methods in medical research. 2013:0962280213503925. 4. 4.Ben-Shlomo Y, Kuh D. A life course approach to chronic disease epidemiology: conceptual models, empirical challenges and interdisciplinary perspectives. Oxford University Press; 2002. 5. 5.Kuh D, Ben-Shlomo Y, Lynch J, Hallqvist J, Power C. Life course epidemiology. Journal of epidemiology and community health. 2003;57(10):778. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiamVjaCI7czo1OiJyZXNpZCI7czo5OiI1Ny8xMC83NzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNS8wNS8yMDIyLjA1LjAzLjIyMjc0NDA3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 6. 6.Villar J, Gunier RB, Tshivuila-Matala CO, Rauch SA, Nosten F, Ochieng R, et al. Fetal cranial growth trajectories are associated with growth and neurodevelopment at 2 years of age: INTERBIO-21st Fetal Study. Nature Medicine. 2021;27(4):647–52. 7. 7.El Marroun H, Tiemeier H, Steegers EA, Jaddoe VW, Hofman A, Verhulst FC, et al. Intrauterine cannabis exposure affects fetal growth trajectories: the Generation R Study. Journal of the American Academy of Child & Adolescent Psychiatry. 2009;48(12):1173–81. 8. 8.O’Keeffe L, Simpkin A, Tilling K, Anderson E, Hughes A, Lawlor DA, et al. Sex-specific trajectories of cardiometabolic risk factors during childhood and adolescence: a prospective cohort study Atherosclerosis. Atherosclerosis. 2018;278:190–6. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F05%2F2022.05.03.22274407.atom) 9. 9.O’Keeffe L, Simpkin A, Tilling K, Anderson E, Hughes A, Lawlor DA, et al. Data on trajectories of measures of cardiovascular health in the Avon Longitudinal Study of Parents and Children (ALSPAC). Data in Brief. 2019;23(103687). 10. 10.Howe LD, Tilling K, Galobardes B, Smith GD, Gunnell D, Lawlor DA. Socioeconomic differences in childhood growth trajectories: at what age do height inequalities emerge? Journal of epidemiology and community health. 2012;66(2):143–8. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiamVjaCI7czo1OiJyZXNpZCI7czo4OiI2Ni8yLzE0MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA1LzA1LzIwMjIuMDUuMDMuMjIyNzQ0MDcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 11. 11.O’keeffe LM, Frysz M, Bell JA, Howe LD, Fraser AJB. Puberty timing and adiposity change across childhood and adolescence: disentangling cause and consequence. 2019:578005. 12. 12.O’Keeffe LM, Tilling K, Bell J, Lee MA, Lawlor DA, Smith GD, et al. Sex-specific trajectories of molecular cardiometabolic trait concentrations through childhood, adolescence and young adulthood: a cohort study. medRxiv. 2021. 13. 13.Bright HD, Howe LD, Khouja JN, Simpkin AJ, Suderman M, O’Keeffe LM. Epigenetic gestational age and trajectories of weight and height during childhood: a prospective cohort study. Clinical epigenetics. 2019;11(1):1–8. 14. 14.O’Keeffe LM, Fraser A, Howe LD. Accounting for height in indices of body composition during childhood and adolescence. Wellcome Open Research. 2019;4(105):105. 15. 15.Walsh JM, McGowan CA, Mahony R, Foley ME, McAuliffe FM. Low glycaemic index diet in pregnancy to prevent macrosomia (ROLO study): randomised control trial. Bmj. 2012;345. 16. 16.Donnelly JM, Walsh JM, Byrne J, Molloy E, McAuliffe F. Impact of maternal diet on neonatal anthropometry: a randomized controlled trial. Pediatric Obesity. 2015;10(1):52–6. 17. 17.Callanan S, Yelverton CA, Geraghty AA, O’Brien EC, Donnelly JM, Larkin E, et al. The association of a low glycaemic index diet in pregnancy with child body composition at 5 years of age: A secondary analysis of the ROLO study. Pediatric Obesity. 2021:e12820. 18. 18.Goldstein H. Multilevel statistical models; 2nd edition ed. London: : Edward Arnold; 1995. 19. 19.O’Keeffe LM, Howe LD, Fraser A, Hughes AD, Wade KH, Anderson EL, et al. Associations of Y chromosomal haplogroups with cardiometabolic risk factors and subclinical vascular measures in males during childhood and adolescence. Atherosclerosis. 2018;274:94–103. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F05%2F2022.05.03.22274407.atom) 20. 20.Tilling K, Macdonald-Wallis C, Lawlor DA, Hughes RA, Howe LD. Modelling childhood growth using fractional polynomials and linear splines. Annals of Nutrition and Metabolism. 2014;65(2-3):129–38. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000362695&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F05%2F2022.05.03.22274407.atom) 21. 21.University of Bristol Centre for Multilevel Modelling. MLwiN Version 3.04. 2019. 22. 22.Statistical Software: Release 16.0 [Internet]. 2019. 23. 23.Health Services Executive. Growth Monitoring Resources 2022 [Available from: [https://www.hse.ie/eng/health/child/growthmonitoring/](https://www.hse.ie/eng/health/child/growthmonitoring/).