Growth is an important indicator of a child’s health(1). To optimally monitor growth during childhood, the WHO established universal Child Growth Standards (WHO-CGS) for children aged 0–5 years for length/height-for-age, weight-for-age, weight-for-length/height and BMI-for-age(2,3) in 2006. These standards were derived from combined growth data of affluent populations from Brazil, Ghana, India, Norway, Oman and the USA, representing all continents. The populations were carefully selected for growth under ‘optimal conditions’. As the distribution of growth data were highly similar across the different ethnic groups between birth and 5 years, the growth standards are generally considered universal for this age range and represent how children ‘should grow’(2,3) ; as such, it was implemented in over 125 countries(Reference de Onis, Onyango and Borghi4). To detect abnormal growth and/or nutritional problems, the WHO defined different indicators based on cut-off centiles: stunting (height-for-age < –2 sd), wasting (weight-for-height < –2 sd) and underweight (weight-for-age < –2 sd)(5).
Several studies have raised doubts about the applicability of such normative child growth standards to all ethnic groups(Reference Wit and Hermanussen6–Reference Christesen, Pedersen and Pournara8). South Asia still has the highest rates of stunting (39 %), underweight (29 %) and wasting (15 %) of children under 5 in the world(9), despite having lower poverty levels and more favourable socioeconomic conditions than Sub-Saharan Africa(Reference Nube10,Reference Coffey11) Even in a high-income country such as the Netherlands, implausibly high thinness rates were found in affluent children of South Asian descent based on WHO cut-off criteria; however, this was not the case when using ethnic-specific BMI criteria(Reference de Wilde, Zandbergen-Harlaar and van Buuren12,Reference de Wilde, Dekker and Middelkoop13) . Consequently, the high prevalence found using WHO BMI criteria could be down to misclassification. As South Asians are characterised by a predisposition to a ‘thin-fat’ phenotype, essentially a muscle-thin but adipose body composition, from birth until adulthood(Reference D’Angelo, Yajnik and Kumaran14–Reference Van Steijn, Karamali and Kanhai16), ethnic-specific growth standards/references are likely to better estimate the nutritional status. However, normative weight-for-age and weight-for-height references are not yet available for this group.
With regard to height, Asian children are also among the shortest in the world(Reference Ulijaszek17), while Dutch children are the tallest(Reference Schonbeck, Talma and van Dommelen18). Due to these differences, Dutch growth curves for the children of Dutch descent are likely to be an insufficient representation of the growth of children of South Asian descent living in the Netherlands(Reference de Wilde, van Dommelen and van Buuren19); the universal WHO-CGS may not be applicable in this case either.
To our knowledge, the prevalence of stunting, wasting and underweight among South Asian children living in affluent countries has not yet been studied. Whether stunting, wasting and underweight rates in Dutch children aged 0–5 in the Netherlands based on ethnic-specific references differ from those derived from WHO-CGS is unknown. It is also unknown whether the largest population of South Asian descent in the Netherlands, which originates from the former Dutch colony Suriname, is similar in growth to children from parents born in India or Pakistan.
The main objective of this study was to identify potential misclassification of stunting, wasting and underweight based on WHO-CGS in 0–5-year-old children of Dutch descent and in two distinct populations of 0–5-year-olds of South Asian descent: one originating from the Dutch ex-colony Suriname and the other directly from South Asia (India and Pakistan).
Methods
Population
Most South Asians in the Netherlands are descendants of Asian Indians who migrated from India to Suriname, a former Dutch colony, between 1873 and 1916. When Suriname became independent from the Netherlands in 1975, many of these Surinamese Asian Indians migrated to the Netherlands, especially to larger cities. The total population of Asian Indians with Surinamese ancestry in the city of The Hague is estimated at around 35 000(Reference Oudhof, Harmsen and Loozen20). Almost 5000 originate directly from India and Pakistan(21).
Data collection
Youth HealthCare in the Netherlands routinely offers periodical health assessments to all children living in the Netherlands, of which at least ten are conducted from birth till the age of 5; the results are recorded. The current study was based on this routine data.
Firstly, growth data of a historic cohort of affluent Surinamese children of Asian Indian descent, which was previously used to construct normative ethnic-specific BMI-for-age curves(Reference de Wilde, van Dommelen and Middelkoop22), were used in the current study to develop ethnic-specific weight-for-age and weight-for-length/height curves for children 0–5 years of age. As children from this cohort were born between 1974 and 1976 (and as such represent the period before the obesity epidemic) and had no known growth constraints, their growth data were intended to set the norm for a ‘normal’ weight-for-age and weight-for-length/height in children of Asian Indian descent. Data extraction and handling were previously published(Reference de Wilde, van Dommelen and Middelkoop22).
The second dataset used for the current study was also based on a series of routine growth measurements in Youth HealthCare, but now from contemporary children. All length/height and weight data of children aged 0–5 years (0–60 months), measured between 1 January 2012 and 31 December 2015, were extracted from digital health records of the Youth HealthCare centre in the city of The Hague (the Netherlands). Additionally, demographic data, such as parental and child family names, date of birth, sex, parental country of birth and parental level of education, were obtained from the Youth HealthCare database. For this study, the level of education was divided into three subgroups: low (≤prevocational education), middle and high (≥undergraduate), according to parents with the highest level of education.
Inclusion criteria
Children were classified as Surinamese Asian Indian if they met two criteria: both parents had to be born in either Suriname (first generation) or the Netherlands (second generation) and both parents needed to have a Surinamese Asian Indian family name, which was determined by matching the names with a list of almost 2500 (Surinamese) Asian Indian family names. To determine an Indian or Pakistani ethnicity, both parents had to be born in either India or Pakistan. Dutch children were included if both parents were listed as born in the Netherlands and they had a typical Dutch family name.
Similar to the data selection of WHO-CGS, only measurements of children who were born singleton at term (≥37 and <42 weeks of gestation) were included. The criterion of exclusive breastfeeding for at least 4 months followed by partial breastfeeding up to 12 months could not be fulfilled as this information was unavailable or incompletely recorded. Children <2500 g were not excluded as this is generally regarded as normal, considering the variation within a normal distribution.
Anthropometric measurements, references and cut-offs
Length/height and/or weight were measured by trained Youth HealthCare physicians, nurses or medical assistants at each routine health assessment. Between birth and the age of 18 months, length was measured with a measuring board in recumbent position with head against the base of the board, fully stretched legs and soles against the end of the measuring board. The length was read to the nearest 0·1 cm. A calibrated baby scale was used to measure weight to the nearest 0·01 kg. In children >18 months, standing height was measured using a calibrated stadiometer that was read to the nearest 0·1 cm. A physician’s balance scale or calibrated mechanical weight scale was used to measure weight, rounding to the nearest 0·1 kg. All measurements were registered on paper health records or the digital record system of Youth HealthCare.
Several growth references were used to calculate length/height-for-age, weight-for-age and weight-for-length/height z-scores for all measurements of contemporary children (Table 1). Firstly, WHO-CGS from the WHO Multicenter Growth Reference Study(2) was used as a universal reference to measurements of the three ethnic groups. Secondly, length/height-for-age z-scores were calculated using the (descriptive) Dutch and Surinamese Asian Indian references of the fifth national growth study (measured 2007–2010)(Reference Schonbeck, Talma and van Dommelen18,Reference de Wilde, van Dommelen and van Buuren19) , of which the latter was also applied to measurements of children of South Asian descent (Pakistan/India). Finally, weight-for-age and weight-for-length/height z-scores for Dutch children were based on the third national growth study references (1978–1979)(Reference Roede and Van Wieringen23), which are considered normative as the included children were not or minimally influenced by the obesity epidemic. The weight-for-age and weight-for-length/height references created in the current study are also intended as normative references and were applied to all measurements of children of Surinamese Asian Indian and South Asian descent in the current study.
Sources of different ethnic-specific growth references
NGS, National Growth Study; SAI, Surinamese Asian Indian; SA, South Asian.
Stunting, wasting and underweight were defined as length/height-for-age, weight-for-length/height and weight-for-age <–2 sd(5). When considering a normal (Gaussian) distribution, a prevalence of 2·3 % is to be expected.
Statistical analyses
Based on the growth data from the Surinamese Asian Indian cohort between 1974 and 1976, ethnic-specific weight-for-age and weight-for-length/height references were calculated, separately for males and females, by employing the LMS method(Reference Cole24) in R (v 3.3.2) with the package for Generalised Additive Models for Location, Scale and Shape (GAMLSS)(Reference Stasinapoulos and Rigby25). This method results in a reference that can be described by three curves: skewness (L-curve), median (M-curve) and coefficient of variation (S-curve). The fit of the LMS model to the data and the necessary amount of smoothing of the curves were determined by assessing ‘worm plots’ (AGD package for R), which are essentially (diagnostic) tools for analysing the residuals(Reference van Buuren and Fredriks26). To achieve a stable fit of the LMS model at 5 years of age, measurements of children aged 5–9·99 were also added to the dataset before fitting. To establish a proper fit in periods of rapid growth, age was cube root-transformed for the weight-for-age references, and for the weight-for-length/height references, the length/height value was log-transformed. After an optimal fit was achieved, the transformed values were scaled back to the original values.
Descriptive statistics of the data of contemporary children (2012–2015) were calculated using Chi-square (categorical variables) and ANOVA (continuous variables). Ethnic differences in z-scores of length/height-for-age, weight-for-age and weight-for- length/height were determined using a linear model within Generalised Estimating Equations (GEE), which was set to account for the repeated measures (autoregressive model; AR1). The specific z-scores based on either WHO-CGS or ethnic-specific references were the dependent variable, and ethnicity was the independent variable, while adjusting for sex and parental level of education. Within this procedure, Estimated Marginal Means (EMM) by age were calculated, which for this study yielded z-score means by age that were adjusted for other factors in the model, also taking the correlation between repeated measurements into account.
To determine the prevalence of stunting, wasting and underweight, dichotomised variables (stunting v. no stunting, wasting v. no wasting, underweight v. no underweight) were calculated including the standard error. Statistical analyses were performed with IBM SPSS Statistics v24.
Results
First, ethnic-specific growth references for weight-for-age and weight-for-length/height were constructed based on the growth of a cohort of 954 Surinamese Asian Indian children (51·3 % boys) born between 1974 and 1976. For the construction of weight-for-age references, the total sample included 7342 measurements of weight, while for weight-for-length/height references, 4667 measurements of height and weight were included. The resulting growth references, including the LMS values, can be found in online Supplemental file S1.
The study of growth of contemporary children had a total of 59 495 measurements from 11 935 children (Table 2). Most children were of Dutch descent (85 %). Parents of the children of Surinamese Asian Indian descent had the lowest education level compared with parents of the children of Dutch and South Asian descent.
Study population characteristics by ethnic group
* Difference between ethnic groups, P < 0·001.
WHO Child Growth Standards
When applying WHO-CGS, a similar pattern for length/height was found in all three ethnic groups (Fig. 1). Compared with the (reference) WHO population, the length at 1 month of age was relatively short in all ethnic groups; this was especially true for Surinamese Asian Indian children. From that age onwards, length increased sharply to a stable trajectory from around 6 months of age in all ethnic groups to a mean z-score of –0·09 (95 % CI–0·14,–0·05) in children of Surinamese Asian Indian descent, and of 0·29 (95 % CI 0·27, 0·31) and 0·37 (95 % CI 0·22, 0·52) in children of Dutch and South Asian descent, respectively (data not shown).
Application of WHO Child Growth Standards for length/height-for-age, weight-for-age and weight-for-length/height: (estimated marginal) mean z-scores by age in months and by ethnic group, adjusted for sex and parental education level (
, Dutch; 
, Surinamese Asian Indian; 
, South Asian).
Weight-for-age also showed a typical trajectory that was quite similar in the three ethnic groups, although the starting level differed. The initial level was relatively low, especially in the Surinamese Asian Indian group, which was followed by an increase between 3 and 12 months, after which a plateau was reached. Surinamese Asian Indian children had the lowest weight-for-age over the whole age range, consistently below the mean of the WHO standard. From around 10 months of age, children of Dutch and South Asian descent had, on average, a higher mean weight-for-age than WHO-CGS (P < 0·05).
In contrast to height-for-age and weight-for-age, weight-for-length/height started at a relatively high level in all ethnic groups, declining in the first 3 months, and then increasing to plateau at an age of 14 months. Children of Dutch descent had the highest weight-for-length/height z-scores, with values up to 0·50 (95 % CI 0·49, 0·53) at an age of 14 months, while children of Surinamese Asian Indian descent had the lowest z-scores, up to –0·40 (95 % CI –0·48, –0·33) at an age of 36 months, which are consistently below WHO-CGS over the whole age range.
At most ages, length/height-for-age, weight-for-age and weight-for-length/height z-scores were significantly higher in South Asian children than in children of Surinamese Asian Indian descent (data not shown).
When stunting, wasting and underweight criteria (<–2 sd) of the WHO reference (Table 3) were applied, a relatively high prevalence of stunting at 1 month of age was evident in all ethnic groups, which declined to the range of expected values at an age of 6 months. The highest level of stunting was found in children of Surinamese Asian Indian descent (16·0 %; 95 % CI 13·6, 18·4 %); simultaneously they had a relatively high underweight prevalence at 1 month of age (7·2 %; 95 % CI 5·5, 8·8 %), which decreased to near-normal levels at 18 months after which it increased again to higher-than-expected rates, up to 4·9 % (95 % CI 3·1, 6·6 %), at 36 months. Also, children of Surinamese Asian Indian descent had higher-than-expected wasting rates over almost the whole age range, up to 6·7 % (95 % CI 4·7, 8·7 %), at 36 months.
Prevalence of stunting, wasting and underweight in Dutch, Surinamese Asian Indian and South Asian children by age in months, based on WHO and ethnic-specific criteria
For all indicators, Dutch children had consistently lower rates at most ages than their Surinamese Asian Indian and South Asian counterparts. Wasting in South Asian children was quite high, up to an age of 9 months, after which the prevalence fluctuated around near-normal levels.
Ethnic-specific growth references
With ethnic-specific normative references applied, length/height-for-age showed largely similar trajectories for the three ethnic groups, with more or less stable levels up to 18 months, after which length/height-for-age declined (Fig. 2). Surinamese Asian Indian and South Asian children also had similarly shaped curves for weight-for-age and weight-for-length/height, although the mean z-scores of South Asian children were consistently higher than that of Surinamese Asian Indian children over almost the whole age range.
Application of ethnic-specific references for length/height-for-age, weight-for-age and weight-for-length/height: (estimated marginal) mean z-scores by age in months and by ethnic group, adjusted for sex and parental education level (, Dutch; , Surinamese Asian Indian; , South Asian)
At 1 month after birth, the weight-for-length/height levels were high in all ethnic groups, even reaching up to +1·15 sd in South Asian children. Subsequently, the levels dropped quite steeply to the bottom at around 10 months of age in Surinamese Asian Indian and South Asian children, and at 18 months in children of Dutch descent. Then the levels rose in both South Asian populations, while the levels in Dutch children remained near the baseline of 0 sd.
For the three ethnic groups, stunting, wasting and underweight rates based on ethnic-specific criteria were not significantly lower and did not significantly differ from the 2·3 % level at most ages (Table 3). Rates of wasting and underweight were especially low in the South Asian group.
Discussion
This study showed that in the first 6–12 months of life, the growth pattern of children of all three ethnic groups deviated from WHO-CGS, after which it plateaued at a relatively stable level. The rates of stunting based on WHO-CGS were relatively high in the first month after birth in Dutch (4·4 %), Surinamese Asian Indian (16·0 %) and South Asian (5·3 %) children, but declined to the expected level of 2·3 % at 6 months of age in children of Surinamese Asian Indian and South Asian descent. However, in Dutch children of the same age, stunting, wasting and underweight rates were considerably lower than expected. The prevalence of wasting and underweight in children of Surinamese Asian Indian descent was high at most ages, up to 6·7 and 7·2 %, respectively, while in children of South Asian descent, these remained around the 2·3 % level.
At all ages in the three ethnic groups, the ethnic-specific criteria for stunting, underweight and wasting led to the expected 2·3 % (or even lower rates than expected), which suggests that the WHO-CGS criteria are likely to misclassify many children as stunted, underweight or wasted. Especially length-for-age and weight-for-age in the first 6–12 months showed large contrasts between WHO-CGS and ethnic-specific references. Where low but increasing values were found using WHO-CGS, the ethnic-specific references showed relatively high but declining values.
In addition, the generally high trajectories of all growth parameters in Dutch children after 12 months of age relative to WHO-CGS, together with very low rates of stunting, wasting and underweight, render WHO-CGS less suitable for this group. Equally, the opposite pattern found in (affluent) children of Surinamese Asian Indian descent, generally low trajectories of growth parameters together with high rates of underweight and wasting over the whole age range, showed the unsuitability of WHO-CGS for children of Surinamese Asian Indian descent.
Besides the Netherlands, other high-income countries also experience the same problem, with quite low values of length and weight shortly after birth when using WHO-CGS, which increase to normal or above-normal levels in the first year(Reference Rolland-Cachera and Peneau27), as well as generally high trajectories of mean z-scores based on WHO-CGS together with lower rates of stunting, underweight and wasting than expected. As children in those countries are generally taller and heavier than those in the WHO-CGS population, the distribution shifts to the right relative to WHO-CGS. Therefore, the detection of growth problems at the lower end of the scale – for instance, Turner syndrome – is delayed, while at the higher end, an overestimation of potential health problems may occur(Reference Christesen, Pedersen and Pournara8,Reference Natale and Rajagopalan28–Reference Heude, Scherdel and Chalumeau30) . At the same time, in children of shorter and thinner populations, such as Surinamese Asian Indian children in our study, wasting and underweight are likely to be overdetected(Reference de Wilde, Zandbergen-Harlaar and van Buuren12,Reference de Wilde, Dekker and Middelkoop13) . For all the above reasons, WHO-CGS was considered unsuitable for use in the studied populations(Reference Rolland-Cachera and Peneau27,Reference Natale and Rajagopalan28,Reference Hermanussen and Wit31) .
Explaining the results
The deviations of height and weight from WHO-CGS, consistently found in many European and other high-income countries, suggest a structural or systematic error in WHO-CGS. Selective dropout from the WHO Multicenter Growth Reference Study was suggested as a plausible explanation(Reference Ziegler and Nelson32), as the data of 50 % of all children enrolled in that study were not used due to non-compliance to the rules set out at the start of the study. For instance, all children had to be exclusively breastfed(3). As breastfed infants generally have a more rapid weight gain in the first 3–4 months of life than formula-fed infants, followed by a slower body growth between 3 and 12 months(Reference Heude, Scherdel and Chalumeau30), the weight patterns seen in other populations relative to WHO-CGS suggest that only very-well-growing babies might have been included.
The large discrepancies found in the current study among children of Surinamese Asian Indian descent are a magnification of this problem of WHO-CGS, as these children, despite their affluence, are generally lighter than other ethnic groups(Reference de Wilde, Dekker and Middelkoop13). The current study also shows that this finding is not restricted to infancy alone: weight-for-age and weight-for-length/height were considerably lower over the whole age range, 0–5 years, leading to an overestimation of underweight and wasting rates.
Furthermore, instead of low mean z-score levels for weight and weight-for-height in the first 6–12 months, the application of ethnic-specific references even showed relatively high mean z-score levels in children from all ethnic groups. The higher levels in Dutch children may be the result of a general increase of mean birth weight of term babies in the past 30 years(Reference Visser, Eilers and Elferink-Stinkens33), but in Surinamese Asian Indian children, the mean birth weight was found to be stable from 1974 to 2009(Reference De Wilde, van Buuren and Middelkoop34). The currently found higher level at 1 month of age may suggest a right shift of the normal curve, with less children below –2 sd (underweight and wasting) as a consequence.
Compared with people of European descent at similar weight or BMI levels, South Asian infants, children and adults generally have a smaller lean body mass and a larger (abdominal) fat mass. This typical ‘thin-fat’ body composition is described in South Asians living in South Asia, as well as South Asians living in Suriname and other countries(Reference D’Angelo, Yajnik and Kumaran14,Reference Van Steijn, Karamali and Kanhai16,Reference Deurenberg, Deurenberg-Yap and Guricci35–Reference Wickramasinghe38) . However, in the present study, South Asian children in the Netherlands who originated directly from India or Pakistan differed considerably from Surinamese Asian Indian children in all growth parameters, although the growth patterns seemed similar. The length/height-for-age was higher than that of Dutch children, while the weight-for-age was only slightly lower. The weight-for length/height was, on average, in-between Surinamese Asian Indian and Dutch children. In Pakistani babies in the UK, a similar growth pattern was found compared with white children(Reference Fairley, Petherick and Howe39).
Previously, the length/height-for-age references and BMI distribution of Surinamese Asian Indian children aged 5–20 in the Netherlands were shown to be quite similar to those based on the growth of affluent Asian Indian children(Reference de Wilde, van Dommelen and van Buuren19,Reference de Wilde, van Dommelen and Middelkoop22) . Nevertheless, the South Asian population aged 0–5 in the Netherlands seemingly grow differently. Therefore, the current South Asian-specific references in the Netherlands, which were based on the growth of a norm population of Surinamese Asian Indian children living in the Netherlands, may not be applicable to this specific group. Considering the high parental level of education of these children, it is more likely that the South Asian population in the city of The Hague is a highly selective sample of Pakistani and Asian Indian children, who may not be representative of Pakistani or Asian Indian children in general. In addition, it cannot be ruled out that the high z-score levels of weight-for-age and weight-for length/height may indicate that the population of Asian Indian and Pakistani children in the Netherlands are in fact overfat.
Practical implications
From a clinical perspective, growth references should be able to detect, with a reasonably high sensitivity and specificity, children at risk of disease or premature death. A universally applicable reference, such as WHO-CGS, should do this well, irrespective of ethnic origin. Because of the potential misclassifications when using WHO-CGS, many studies raised doubts about its suitability to detect such at-risk children, and recommended the use of local (or ethnic-specific or country-specific) references instead of WHO-CGS(Reference Wit and Hermanussen6,Reference Ziegler and Nelson32,Reference Kerac, Blencowe and Grijalva-Eternod40,Reference Milani, Buckler and Kelnar41) . As such references can likely reduce the risk of misclassification, these have the potential of elucidating true differences in growth and nutritional status between ethnic groups. However, as the creation of such references is not always feasible, the use of WHO-CGS is an alternative, the interpretation of which justifies a critical eye.
To be able to properly assess the growth of children of Dutch and Surinamese Asian Indian descent, we recommend using local, ethnic-specific references available in the Netherlands. The use of WHO-CGS may have potential adverse effects in Surinamese Asian Indian children, such as overfeeding, which may ultimately result in a larger (abdominal) fat mass. Ethnic-specific references for weight-for-age and weight-for-height are likely to better reflect the typical body composition of Surinamese Asian Indian children compared with WHO-CGS. For children who directly originated from Pakistan or India, further research on the applicability of ethnic-specific references is recommended, preferably in relation to body composition measures. The found differences between the subpopulations of the same ethnic group of ‘South Asians’ stress that heterogeneity within populations should be taken into account when interpreting results based on one reference.
Strengths and limitations
Strengths of this study are the large sample size, reliability of data and stringent criteria for determining ethnicity, which restricted the inclusion to only children whose both parents met the criteria for ethnicity. Moreover, analyses were adjusted for individual socioeconomic status (highest level of education of parents). A limitation of the current study is that normative references based on the 1974–1976 cohort could not be corrected for the correlation of longitudinal measurements, which may have influenced the precision of the centiles. Also, formula-fed children could not be excluded due to this information not being available. Although ethnic differences in breastfeeding in the Netherlands are known(Reference Quittner, Franssen and Steenkamer42), most infants in the Netherlands are still breastfed during the first 6 months(Reference Peeters, Lanting and Van Wouwe43). Therefore, the effects of formula feeding on the results are expected to be attenuated. Because the general prevalence of growth disorders and the use of growth-interfering medication is low and most of the involved children receive paediatric care instead of Youth HealthCare, the number of measurements from such children in our database are expected to be low. Therefore, we expect a negligible influence of it on the results of this study.
Conclusion
The use of WHO-CGS misclassifies a considerable proportion of children from all ethnic groups as stunted in the first 6 months of life, while underweight and wasting are considerably overestimated in children of Surinamese Asian Indian descent. To properly assess the growth of children aged 0–5 in the Netherlands, ethnic-specific growth references are recommended for use in Dutch and Surinamese Asian Indian children. Further research into the suitability of normative South Asian references for other groups originating from the South Asian subcontinent is needed.
Acknowledgements
Acknowledgements: The authors thank Ron Smit for his help with extracting data from the Youth HealthCare database. Financial support: This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors. Conflicts of interest: There are no conflicts of interest. Authorship: J.A.d.W. designed the study, conceptualised the paper, guided the analyses and wrote the manuscript. M.P.K drafted the first version of the paper. All co-authors contributed to the interpretation of results, reviewed/edited different drafts of the manuscript, and read and approved the final version. Ethics of human subject participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki. This study was based on routinely collected health data and, as such, does not fall under the Dutch Medical Research Involving Human Subjects Act(1). The use of such data is allowed without ethical approval as long as strict data management rules regarding privacy protection are followed. For that reason, personal information was kept separate from the research database, in which only derived data were included.
Supplementary material
For supplementary material accompanying this paper visit https://doi.org/10.1017/S1368980019004464
Open access
