Statistical analysis
All analyses were done using STATA version 16.0 (Stata-Corp LLC). We calculated the means (The primary analysis included the comparison of neurodevelopment outcomes at 12 and 24 months of age between the three study groups and was based on the intention-to-treat principle. Effect sizes (difference in means and 95 % CI) for the continuous outcomes were calculated using generalised linear models of the Gaussian family with an identity-link function. The primary analysis was unadjusted as there were no significant differences in the baseline characteristics among children in the three groups. Additionally, we also conducted an adjusted analysis after including variables in the models that have been shown to influence neurodevelopment outcomes, based on previously published studies