search
for
 About Bioline  All Journals  Testimonials  Membership  News


Iranian Journal of Pediatrics
Tehran University of Medical Sciences Press
ISSN: 1018-4406 EISSN: 2008-2150
Vol. 21, Num. 1, 2011, pp. 88-94

Iranian Journal of Pediatrics, Vol. 20, No. 1, January-March, 2011, pp. 88-94

Standardized Percentile Curves of Body Mass Index of Northeast Iranian Children Aged 25 to 60 Months

Maryam  Emdadi1, MSc; Mohammad Safarian2, PhD, and Hassan Doosti1, PhD

  1. Department of Statistics, Ferdowsi University, Mashhad, Iran
  2. Department of Biochemistry and Nutrition, Mashhad University of Medical Sciences, Mashhad, Iran

* Corresponding Author; Address: Department of Biochemistry and Nutrition, Mashhad University of Medical Sciences, Azadi Square, Mashad, Iran E-mail: safarianm@mums.ac.ir

Received:  Feb 24, 2010; Final Revision: Jul 12, 2010; Accepted: Aug 30, 2010

Code Number: pe11015

Abstract

Objective: Growth charts are widely used to assess children’s growth status and can provide a trajectory of growth during early important months of life. Racial differences necessitate using local growth charts. This study aimed to provide standardized growth curves of body mass index (BMI) for children living in northeast Iran.
Methods: A total of 23730 apparently healthy boys and girls aged 25 to 60 months recruited for 20 days from those attending community clinics for routine health checks. Anthropometric measurements were done by trained health staff using WHO methodology. The LMSP method with maximum penalized likelihood, the Generalized Additive Models, the Box-Cox power exponential distribution distribution, Akaike Information Criteria and Generalized Akaike Criteria with penalty equal to 3 [GAIC(3)], and Worm plot and Q-tests as goodness of fit tests were used to construct the centile reference charts.
Findings: The BMI centile curves for boys and girls aged 25 to 60 months were drawn utilizing a population of children living in northeast Iran.
Conclusion: The results of the current study demonstrate the possibility of preparation of local growth charts and their importance in evaluating children's growth. Also their differences, relative to those prepared by global references, reflect the necessity of preparing local charts in future studies using longitudinal data.

Key Words: Body Mass Index; Growth Charts; Anthropometry; Children; Reference Growth Curves

Introduction

Age-related reference curves are widely used in pediatric clinics and infant health units for growth screening of children. These charts are used to determine the trajectory of growth as a marker of nutritional status at individual and population levels to assist clinicians and public health policy makers. Differences between populations with respect to growth patterns because of racial, environmental and genetic factors affecting growth beside nutrition, necessitates the use of curves obtained from local data. Using the National Center for Health Statistics (NCHS) and international norms for assessing growth in Iranian children may be misleading and therefore, we aimed to develop local growth charts from our standardized and representative local data[1-4].  

Several methods for developing age-related curves have been developed, some of them are based on distributional assumptions and others are not. One of the methods of growth chart drawing named LMS method is based on distributional assumptions. This method had been offered by Cole and Green in 1992[4,5]. The LMS method is widely used to fit a centile curve[6]. In this method, to fit a suitable model, one has the measurement in the reference sample following a standard normal distribution on all ages after a chosen Box-Cox transformation. In this research, we have used the LMSP model, a generalized model of LMS method. This model is proposed by Rigby and Stasinopoulos[7]. LMSP provides a model for a dependent variable which shows both skewness and kurtosis. On the other hand, LMSP summarizes the changing distribution with age according to 4 curves representing the median (M), the coefficient of variation (S), skewness (L), and kurtosis (P). In this method the maximum penalized likelihood is also used. The use of the maximum penalized likelihood approach allows us to provide smooth estimates of L, M, S, P curves directly[8]. In this paper, percentile curves of body mass index [BMI=weight/ (height)2] for Iranian boys and girls aged 25-60 months are presented based on data on children living in northeast Iran. The paper compares the results with the newest World Health Organization (WHO) and Center for Dieses Control (CDC) results.

Subjects and Methods

A total of 23730 apparently healthy boys and girls aged 25 to 60 months recruited for 20 days from those attending community clinics for routine health check. This study was conducted in Grand Khorasan Province (northeast Iran) during 2005. Weight and height of the subjects were measured based on WHO methodology by trained health staff and their age was recorded in months. Data quality control showed that the data follow the world health organization standards with respect to the sample quantity in each age category, the mean weight standard deviation and data distribution.

In this study, BMI charts for boys and girls were constructed using LMSP method proposed by Rigby and Stasinopoulos, and implemented in the Generalized Additive Model for Location, Scale and Shape (GAMLSS) package in R software[9]. In LMSP method, GAMLSS parameters[8] and the Box-Cox power exponential distribution (BCPE) were used for model fitting to data[7]; Also Worm plot and Q-test were used for goodness of fit test[10,11]. The GAMLSS is a general class of statistical models for a univariate response variable. For maximizing the penalized likelihood of data under GAMLSS, Rigby and Stasinopoulos and Cole and Green algorithms were used[8]. The BCPE distribution has 4 parameters and is denoted as BCPE μ, σ, ν, t. The para-meters μ, σ, ν and t may be interpreted as relating to location (median), scale (approximate coefficient of variation), skewness (transformation to symmetry) and kurtosis (power exponential parameter), respectively.

Smooth centile curves are drawn by modeling each of the 4 parameters of the distribution as a smooth non-parametric function of an explanatory variable[7]. Now, we introduce the 2 tools for goodness of fit tests that was used in this research. 1) Worm plot: This is a graphical diagnostic tool for assessing goodness of fit and for analysis of residuals. The shapes of the Worm plot indicate how the data differ from the assumed underlying distribution, and suggest useful modifications of the model[10]; 2) Q-tests: Q-tests are described by Royston and Wright[12]. These tests include tests for location, scale, skewness, kurtosis and non-normality[11]. Comparison of the 2 tools for goodness of fit, shows that the Worm plot is the ideal tool for preliminary model building and the Q-test is useful to determine the final model[4]. In fitting  procedure,  before  fitting the model, an age transformation was needed to stretch the age scale for values close to zero. For this, a power transformation is applied to age, i.e. f(λ)=ageλ [4]. Therefore, at first we should search the best λ for the age power transformation and then determine the age power transformation and then determine the best degrees of freedom for the parameter curves. In order to draw centile curves by the LMSP method, the following sequence of steps should be implemented. At first, the M curve should be selected and then S curve, then L curve and finally P curve. The sequence of steps is sensible because the M curve describes the most important variation, while the influence of S, L and P is relatively small. Also in order to select the best model for μ parameter, we used the Akaike Information Criteria (AIC)[13], defined as: AIC=-21+ 2p, where l is the penalized maximum likelihood and p is the number of parameters (or the total number of degrees of freedom); to determine the best model for σ, υ and τ, we had used the Generalized Akaike Criteria with penalty equal 3 [GAIC (3)] as defined in Rigby and Stasino­poulos: GAIC=-21+3p. According to these 2 criteria, the best model was the one with smallest AIC or GAIC value.   

There were 11897(50.2%) and 11833 (49.8%) observations for boys and girls respectively. A total of 23730 records with both weight and height were available for constructing the BMI-for-age curves for girls and boys. The cross-sectional sample sizes by age and sex are given in Table 1.

Findings

Centile curves (3rd, 10th, 25th, 50th, 75th, 90th, 97th) were generated using the LMSP method for BMI-for-age. These are shown by sex in Fig. 1 and 2.

The fitted models of data are listed below:

BCPE [x=age, df (μ)=3, df (σ)=9, df(υ)=1, df(τ)=1]; for boys

BCPE [x=age, df (μ)=9, df (σ)=4, df(υ)=1, df(τ)=1]; for girls

Comparisons of BMI 50th percentile curves for northeast Iran and references from United States are shown in Fig. 3 and 4. In Tables 2 and 3, values of BMI 50th percentile of northeast Iran and WHO reference are given.

In this research, the Worm plot and Q-test were used together for selecting the preliminary model but for refining the final model the Q-test was used. Group-specific Q-test statistics resulting in absolute values of z1, z2, z3 or z4 that were larger than 2 were interpreted to indicate a misfit of, mean, variance, skewness or kurtosis respectively. The overall Q-test statistics combining all groups were based on Chi-square distribution, which assumes that observations from different groups are independent. The results of Q-test show misfits in fitted models of 4 parameters (μ, σ, υ and τ) of distribution for a few age groups for boys and girls, although overall Q-test statistics for 4 parameters were non-significant or P-values were large than 0.05 (for boys: P-values for M, S, L and P curves were 0.1, 0.06, 0.9 and 0.2 respectively, and for girls: 0.07, 0.07, 0.6 and 0.1 respectively).

Discussion

The present study was conducted in Grand Khorasan  Province  (Iran) during 2005. The main purpose of this study was to obtain smoothed growth charts for children aged 25 to 60 months. The recruited sample was 11% of the total children living in Grand Khorasan province with above-mentioned age.  Children were apparently healthy   and   attended   local   community   health clinics for routine health checking. It was conducted as a cross sectional study with a representative sample of children living in northeast Iran. Racial similarities, population distribution characteristics and the extension of Grand Khorasan province probably make these curves applicable to children in Iran. The present study shows that the BMI of children aged 25 to 60 months living in northeast Iran were lower than WHO centiles for both genders.

Regarding data quality, WHO standards were followed by the obtained data[4]. Construction of local curves to measure growth pattern in pediatric clinics and to determine the population at risk of malnutrition may be much more helpful than using international norms or charts obtained from other populations. Comparing our charts with those of WHO and CDC showed significant differences between growth patterns of Iranian children and other populations. Similar differences have been highlighted in other studies of various countries[6].

Comparisons of Qatari children of less than 5 years old with NCHS showed a significant deviation from NCHS pattern[14]. In comparison with NCHS, 11% of children suffered from dwarfism and 4.4% were underweight and only 4.9% were above +2SD of NCHS indicating growth pattern differences. Similar study in Saudi Arabia for children under 3 years has also reported a significant difference between NCHS charts and Saudi children[15]. In Iran, though a study in Isfahan, in 1997, reported that the NCHS curve can be suitably used for growth evaluation in Iranian children[16], others have reported reverse. In a study in Shiraz in 1996, all Shirazi children curves were lower than the corresponding NCHS curves[17].

Recently, In Jahrom, the study carried out on infants growth charts showed that height and weight of children under two years were generally lower than those of NCHS[18]. Although the WHO or NCHS charts can give a good measure of children's growth status and may be used for risk assessment, they cannot reliably be applied to other populations such as Iranians and because of inconsistency with the growth pattern of Iranian children, they may even be misleading. Therefore, constructing local charts seem to be of high priority in health policies for child care. Because of the nature of cross-sectional studies as in this study, obtained charts may need to be reconstructed or be confirmed by longitudinal data with more representative samples across the country.

It is important to consider limitations of current study when using obtained charts. First of all, data regarding general health and nutritional practices were not available. However, the obtained charts are regarded as reliable considering the sample size and healthy condition of the recruited children. Secondly, there is not enough data regarding feeding practice of children younger than two years of age, namely, the age of commencing the complementary feeding or weaning age and whether they have been breast-fed or bottle fed.

Growth patterns of these two groups may be different and our study was unable to recognize these differences. Thirdly, the data used in this study are cross-sectional and obviously of less validity compared to longitudinal data.

Despite all these obstacles, the charts prepared by using well-known methods and based on a data set of a representative sample of Iranian children maybe preferred over those prepared by using data of other populations. Therefore, we recommend using our local charts instead of international ones until our longitudinal data has been prepared or our charts updated.

Conclusion

The current study demonstrates the possibility of local growth charts preparation and their importance for evaluating children's growth and risk assessment. The differences between our local charts and those prepared using international references, showed the necessity of preparing and using local charts for clinical and epidemiologic applications.

Obviously further studies and longitudinal data would be required to confirm these results.      Constructing similar charts using longitudinal studies and regular updating to account the secular trends is recommended.

Acknowledgment

The authors appreciate the Deputy of Health of Mashhad University of Medical Sciences for their financial support and give an especial thank to Dr. Abasalti and Dr. Shakeri for editorial assistance and preparing the data and finally the University Review Board for approval of this study.

Conflict of Interest: None

References

  1. Panpanich R, Garner P. Growth monitoring in children (Cochrane Review). In: The Cochrane Library, Issue 3, 2003. Oxford: Update Software.
  2. Hamill PV, Drizd TA, Johnson CL, et al. Physical growth: National Center for Health Statistics Percentiles. Am J Clin Nutr 1979;32(3):607-29.
  3. Ogden CL, Kuczmarski RJ, Flegal KM, et al. Centers for Disease Control and Prevention 2000 growth charts for the United States: improvements to the 1977 National Center for Health Statistics version. Pediatrics 2002;109(1):45-60.
  4. WHO Child Growth Standards. Geneva: World Health Organization. 2006.
  5. Cole TJ, Green PJ. Smoothing reference centile curves: the lms method a penalized likelihood. Stat Med 1992;11(10):1305-19.
  6. Hosseini M, Carpenter RG, Mohammad K, et al. Standardized percentile curves of body mass index of Iranian children compared to the US population reference. Int J Obes Relat Metab Disord 1999;23(8):783-7.
  7. Rigby RA, Stasinopoulos DM. Smooth centile curves for skew and kurtotic data modelled using the Box-Cox power exponential distribution. Stat Med 2004; 23(19):3053-76. 
  8. Rigby RA and Stasinopoulos. Generalized additive models for location, scale and shape. Appl Statist 2005;54(part 3):1-38.
  9. Stasinopoulos D, Rigby B, Akantziliotou C. Instructions on how to use the GAMLSS package in R. 2nd ed, Availabe at: http://studweb.north. londonmet.ac.uk/~stasinom/papers/gamlss-manual.pdf. Access date: Jan 11, 2008.   
  10. van Buuren S, Fredriks M. Worm plot: a simple diagnostic device for modeling growth reference curves. Stat Med 2001;20(8):1259-77.
  11. Pan H, Cole TJ. A comparison of goodness of fit tests for age-related reference ranges. Stat Med 2004; 23(11):17-22.
  12. Royston P, Wright EM. Goodness of fit statistics for age specific reference intervals. Stat Med 2000; 19(21):2943-62.
  13. Akaike, H. A new look at the statistical model identification. IEEE Trans Autom Control 1974;19:716-23.
  14. Kamal AA, Bener A, Kareem Al-Mulla AM. Growth pattern of Qatari preschool children. Croat Med J 2004;45(4):461-5.
  15. Bamgboye EA, Al-Nahedh N. Factors associated with growth faltering in children from rural Saudi Arabia. Afr J Med Med Sci 2003;32(4):343-7.
  16. Aminorroya A, Amini M, Naghdi H, Zadeh AH. Growth charts of heights and weights of male children and adolescents of Isfahan, Iran. J Health Popul Nutr 2003;21(4):341-6.
  17. Ayatollahi SMT, Ahmadi K. Infants'growth chart for southern Iran. Ann Hum Biol 2001;28(3):337-45.
  18. Heydari ST, Emamghoreishi F, Amini M. Infants' growth charts in Jahrom, Iran. Iran J Pediatr 2009;19(1):25-34.

Copyright 2011 - Iran Journal of Pediatrics


The following images related to this document are available:

Photo images

[pe11015f3.jpg] [pe11015t2.jpg] [pe11015f4.jpg] [pe11015t3.jpg] [pe11015t1.jpg] [pe11015f2.jpg] [pe11015f1.jpg]
Home Faq Resources Email Bioline
© Bioline International, 1989 - 2024, Site last up-dated on 01-Sep-2022.
Site created and maintained by the Reference Center on Environmental Information, CRIA, Brazil
System hosted by the Google Cloud Platform, GCP, Brazil