Skip to main content

Exploring spatial variations and factors associated with childhood stunting in Ethiopia: spatial and multilevel analysis



Stunting reflects a failure to receive adequate nutrition over a long period of time. Stunting is associated with adverse functional consequences including poor cognition, low educational performance, low adult wages, and poor reproductive outcomes. The objective of the study was to investigate spatial variations and factors associated with childhood stunting in Ethiopia.


This study is a secondary data analysis of the 2011 Ethiopian Demographic and Health Survey (EDHS). A total of 9893 children aged 0–59 months were included in the analysis. The Getis-Ord spatial statistical tool was used to identify high and low hotspots areas of stunting. A multilevel multivariable logistic regression was used to identify factors associated with stunting.


Statistically significant hotspots of stunting were found in northern parts of the country whereas low hotspots where there was less stunting than expected were found in the central, eastern, and western parts of the country. In the final model of multilevel logistic regression analysis, individual and community level factors accounted for 36.6 % of childhood stunting. Short birth interval [AOR = 1.68; 95%CI: (1.46–1.93)], being male [AOR = 1.20; 95%CI: (1.08–1.33)], and being from a male-headed household [AOR = 1.18; 95 % CI: (1.01–1.38)] were the factors that increased the odds of stunting at the individual level. Children in the age group between 24–35 months were more likely to be stunted than children whose age was less than one year [AOR = 6.61; 95 % CI: (5.17–8.44)]. The odds of stunting among children with severe anemia were higher than children with no anemia [AOR = 3.23; 95%CI: (2.35–4.43)]. Children with mothers who had completed higher education had lower odds of being stunted compared to children whose mothers had no formal education [AOR = 0.42; 95%CI: (0.18–0.94)]. The odds of being stunted were lower among children whose fathers completed higher education [AOR = 0.58; 95%CI: (0.38–0.89)] compared to children whose fathers had no formal education. Children whose mothers who had high a Body Mass Index (BMI) (≥25.0 kg/m2) were less likely to be stunted compared with children whose mothers had a normal BMI (18.5 kg/m2-24.9 kg/m2)[AOR = 0.69; 95%CI: (0.52–0.90)]. Children from the poorest wealth quintile had higher odds of being stunted compared to children from the richest wealth quintiles [AOR = 1.43; 95 % CI: (1.08–1.88)]. Unavailability of improved latrine facilities and living in the northern parts of the country (Tigray, Affar, Amhara and Benishangul-Gumuzregions) were factors associated with higher odds of stunting from the community-level factors.


Stunting in children under five years old is not random in Ethiopia, with hotspots of higher stunting in the northern part of Ethiopia. Both individual and community-level factors were significant determinants of childhood stunting. The regions with high hotspots of child stunting should be targeted with additional resources, and the identified factors should be considered for nutritional interventions.

Peer Review reports


Height-for-age Z-score is used as an indicator of linear growth retardation and cumulative growth deficits in children [1]. Stunting is defined as a child with a height-for-age Z-score (HAZ) less than minus two standard deviations (<−2 SD) below the median of a reference height-for-age standard. Stunting reflects a failure to receive adequate nutrition over a long period of time, and is affected by both recurrent and chronic illness [2]. Stunting in childhood is associated with adverse functional consequences later in life including poor cognition, poor educational performance, low adult wages and poor reproductive outcomes [3, 4]. Stunted children also have a higher risk of being overweight or obese later in life, putting them at risk of chronic disease in adulthood [3, 57].

Childhood stunting varies not only across various regions of the world but also within and between local authorities, regional space dimensions and/or countries [8]. Although Ethiopia has made steady progress in reducing stunting (from 2000–2011 stunting declined from 58 % to 44 %)[9, 10], the prevalence of stunting is still one of the highest in the world [10], and the issue is a national priority. In order to accelerate efforts in reducing stunting, along with other nutritional problems, Ethiopia revised its national nutrition program in 2013 [11]. Identifying hotspots, or areas where the prevalence of stunting is higher than the national average, would help the Ethiopian Government strategically intensify interventions in order to reduce the prevalenceof stunting in the country.

There aremany studies which determine the prevalence of stunting and analyze socioeconomic, demographic and cultural factors associated with childhood stunting in Ethiopia [12, 13]. However, no spatial analysis has identified the hotspotsof stunting in the country. Moreover, many of the previous studies used standard logistic regression to identify the independent predictors of stunting. Analyzing variables from different levels at one single common level using the standard binary logistic regression model leads to biased results (loss of power or Type I error). Households in the same geographic cluster have common characteristics such as seasonal variability, types of crops and housing which can have a similar impact on the nutritional status of children in the cluster. The assumptions of independence among individuals within the same cluster and of equal variance across clusters are violated in the case of grouped data. Hence, a multilevel analysis is the appropriate statistical analysis method for such a study. This study employed a multilevel logistic regression analysis which has a number of advantages over standard logistic regression, as described in detail by Guo and Zhao [14]. This study aimed to investigate the spatial variation of stunting and the factors associated with stunting in Ethiopia using spatial and multilevel analyses.


Study design and setting

An in-depth secondary data analysis was conducted using Ethiopian Demographic and Health Survey (EDHS) data from 2011. The EDHS is carried out every five years to provide health and health-related indicators at the national and regional levels in Ethiopia. Administratively, regions in Ethiopia are divided into zones, and zones, into administrative units called woredas. Each woreda is further subdivided into the lowest administrative unit, called kebeles. During the 2007 census each kebele was subdivided into census enumeration areas (EAs), which were convenient for the implementation of the census. The 2011 EDHS sample was selected using a stratified, two-stage cluster design where EAs were the sampling units for the first stage, and households for the second stage. The detailed sampling procedure is presented in the full EDHS report [10].


The length of children aged < 24 months was measured during the EDHS in a recumbent position to the nearest 0.1 cm using a locally made board with an upright wooden base and movable headpieces. Children ≥24 months were measured while standing upright. The height-for-age Z-score, an indicator of nutritional status, was compared with reference data from the WHO Multicentre Growth Reference Study Group, 2006 [15]. Children whose height-for-age Z-score is < −2 SD from the median of the WHO reference population are considered stunted (short for their age). A wealth index was constructed using principal components analysis on household asset data to categorize individuals into five wealth quintiles (poorest, poorer, medium, richer and richest). Variables included in the wealth index were ownership of selected household assets, size of agricultural land, quantity of livestock and materials used for house construction [16]. Three steps were used in the construction of the wealth index to permit greater adaptability of the wealth index to both urban and rural areas. In the first step, a subset of indicators common to urban and rural areas was used to create wealth scores for households in both areas. In the second step, separate factor scores were produced for households in urban and rural areas using area-specific indicators. The third step combined the separate area-specific factor scores to produce a nationally-applicable combined wealth index by adjusting area-specific scores through a regression on the common factor scores. A more detailed description of the wealth index is presented in the full EDHS report [10].

Explanatory variables

The individual- and community-level variables included in the study as explanatory variables are shown in Table 1, along with the coding and definitions. Individual-level variables include socio-demographic and economic characteristics (Level one). Community-level variables describe the cluster of the communities living in the same geographical living environment (Level two). These two hierarchal levels were used to create a multilevel analysis for this study. Communities were based on sharing a common primary sampling unit (cluster) within the EDHS data. A multilevel logistic regression model was applied for three reasons. First, in the EDHS sample, primary sampling unit (PSU) was used to define the clusters. Second, it has been shown that for most of the DHS data, the sample size per cluster meets the optimum size with a tolerable precision loss to do a multilevel analysis [17]. Third, multilevel modeling systematically analyzes how covariates at various levels of hierarchal structure affect the outcome variable and how the interactions among covariates measured at different levels affect the outcome variable. Moreover, multilevel modeling corrects for the biases in parameter estimates resulting from clustering and provides correct standard errors [14].

Table 1 Variables definition

Data analysis

Data analysis was carried out using STATA version 12(StataCorp, College Station, Texas, United States) statistical software and spatial analysis was done using ArcGIS software, version 10.0 (ESRI, Redlands, CA, USA). The authors used the “svy” command in STATA version 12 to weight the survey data as per recommendation of the EDHS. Sample weights were applied in order to compensate for the unequal probability of selection between the strata that were geographically defined, as well as for non-responses. A detailed explanation of the weighting procedure can be found in the EDHS methodology report [10]. Multilevel logistic regression was carried out using STATA version 12data analysis and statistical software.

Spatial analysis

Spatial analysis was done using GIS Getis-Ord statistics. The prevalence rates of stunting were exported into ArcGIS to visualize key estimations, and the excess risk of stunting of each region was calculated. Excess risk is defined as a ratio of the observed number over the expected number of cases. In this study, the Local Getis-Ord G index (LGi) was applied to do spatial statistical analysis. The LGi is described in detail in the literature [18]. The spatial heterogeneity of significant high prevalence/low prevalence areas of stunting were computed for each cluster using the Getis-Ord G-statistic tools in ArcGIS. The Local Getis-Ord G index helped to classify the autocorrelations into positive and negative correlations. If prevalence rates had similar high values or low values, they were defined as positive autocorrelation hotspots (represented as High-High or Low-Low autocorrelation). If the attributes held opposing high and low values, they were considered to have negative autocorrelation (represent as High-Low or Low-High autocorrelation). To determine the significance of these statistics, Z-scores and P-values were used. A Z-score near zero indicates no apparent clustering within the study area. A positive Z-score with P-value of <0.05 indicates statistical clustering of hotspots of childhood stunting whereas a negative Z-score with p-value of <0.05 indicates statistical clustering of children who are not stunted.

Multilevel logistic regression

Multivariable multilevel logistic regression was used to analyze factors associated with childhood stunting at two levels: individual and community (cluster) levels. Four models were constructed for this multilevel logistic regression analysis. The first model was an empty model without any explanatory variables, to evaluate the extent of the cluster variation on stunting. The second model adjusted for the individual-level variables, the third model adjusted for community-level variables while the fourth model adjusted for both the individual- and community-level variables simultaneously. A P-value of <0.05 was used to define statistical significance. Adjusted Odds Ratios (AOR) with their corresponding 95 % confidence intervals (CIs) were calculated to identify the independent predictors of stunting. Intra-cluster correlation (ICC), Median Odds Ratio (MOR) and proportional change in variance (PCV) statistics were calculated to measure the variation between clusters. ICC was used to explain cluster variation while MOR is a measure of unexplained cluster heterogeneity [19]. ICC is the measure of variation attributed to contextual neighborhood factors(residential level factors), and is often used to operationalise the concept of contextual phenomena [20]. The ICC was calculated using this formula [\( ICC=\frac{\uptau^2}{\left({\uptau}^2+\frac{\uppi^2}{3}\right)}, \) where τ2 is the estimated variance of clusters)] described elsewhere [21]. MOR is defined as the median value of the odds ratio between the area at highest risk and the area at lowest risk when randomly picking out two areas and it was calculated using the formula \( \left[\mathrm{MOR}= \exp \kern0.28em \left(\sqrt{2\mathrm{x}{\uptau}^2\times 0.6745}\right)\approx \exp \kern0.28em \left(0.95\uptau \right)\right]. \) In this study MOR shows the extent to which the individual probability of being stunted is determined by residential area [21]. The proportional change in variance (PCV) measures the total variation attributed by individual level factors and area level factors in the multilevel model. MOR and the formula for PCV have been described elsewhere [20, 22, 23].

Ethical considerations

The data were downloaded and analyzed after the purpose of the analysis was communicated and approved by MEASURE DHS. The original EDHS data were collected in conformation with international and national ethical guidelines. Ethical clearance for the survey was provided by the Ethiopian Public Health Institute (EPHI), former Ethiopian Health and Nutrition Research Institute (EHNRI) Review Board, the National Research Ethics Review Committee (NRERC) at the Ministry of Science and Technology, the Institutional Review Board of ICF Macro International, and the United States Center for Disease Control and Prevention(CDC).


Characteristics of the study participants

Majority of the study subjects (87.1 %) were from rural area. Most of the mothers (69.3 %) had no formal education while 50.22 % of their partners had also no formal education. The poorest wealth quintile comprises about 22.8 % of the total population. A total of 11.43 % of the children were infants less than a year. Most of the children had a preceding birth interval of 24 months and above. Only 7.46 % of the households have improved latrine while 10.36 % used piped water as source of drinking water (Table 2).

Table 2 Socio-demographic and economic characteristics of respondents included in the analysis, 2011 EDHS

Spatial variation of stunting

The overall prevalence of childhood stunting in Ethiopia was 44.5 % (95 % CI: 43.6–45.5). There were regional variations, with Tigray, Amhara, Afar and Benishangul-Gumuzregions having a higher prevalence of stunting, while the lowest prevalence of stunting was found in Addis Ababa(Table 3). Rural areas had a higher prevalence of stunting than urban areas. The most excess cases of childhood stunting were found in Amhara, Tigray and Afar regions, while lower excess risk was found in Addis Ababa and Gambella Regions (Table 3).

Table 3 Regional variation of prevalence rate and risk of stunting in Ethiopia, DHS 2011

All zones in Addis Ababa, in Harari, four zones in Oromia region (Jimma, East Wollega, Arsi and West Showa) and two zones in Somali region (Jijiga and Shinile) were significantly clustered with low prevalence of stunting (negative Z-score and Gipvalue < 0.05). Five zones in Amhara region (Kamissie zone, Debub Gonder,SemienWello,Wag Himra and Misrak Gojjam), and Southern zone in Tigray region were significantly clustered with high prevalence of stunting (Positive Z-score and Gipvalue < 0.05). The rest of the zones were not significantly clustered with either low or high prevalence of stunting (Fig. 1).

Fig. 1
figure 1

Statistical significant hotspots of childhood stunting at zonal level, DHS 2011

Figure 2 shows the spatial variation of stunting at the cluster level (lower level). The spatial analysis at the cluster level shows that statistically significant high hotspots of stunting were found in northern parts of the country(Amhara, Benishangul-Gumuz, Tigray and Affar regions), whereas statistically significant low spots of stunting were found in the western (Gambella), central (Addis Ababa) and eastern (DireDawa) parts of the country (Table 4 and Fig. 2).

Fig. 2
figure 2

Statistical significant hotspots of childhood stunting at cluster level, DHS 2011

Table 4 Hotspot and coldspot analysis of stunting among enumeration areas (clusters) per regional state in Ethiopia, EDHS 2011

Multilevel analysis

The results of multilevel logistic regression models for individual and community level factors are displayed in Table 5. In the final full model where all individual and community level factors are included, child’s age; gender; birth interval; maternal body mass index (BMI); education status of mother; educational status of father; sex of the household head; child’s anemia status; low household wealth index; region; and lack of availability of an improved latrine were factors that were significantly associated with childhood stunting.

Table 5 Factors associated with childhood stunting in Ethiopia by multilevel logistic regression analysis, EDHS 2011

Individual level factors

Children 24–35 months old were 6.61 times (AOR = 6.61; 95 % CI: 5.17–8.44) more likely to be stunted than children less than one year old. The odds of stunting were increased by 20 % (AOR = 1.20; 95 % CI: 1.08–1.33) in male children compared to females. The odds of stunting among children with severe anemia were 3.32 times higher (AOR = 3.23; 95 % CI: 2.35–4.43) than in children with no anemia. Children with mothers who completed higher education were 58 % (AOR = 0.42; 95 % CI: 0.18–0.94) less likely to be stunted compared to those children whose mothers had no formal education. The odds of stunting were42% lower among children with fathers who completed higher education (AOR = 0.58; 95 % CI: 0.38–0.89) compared to children whose fathers had no formal education.

Children with mothers who had a high BMI (≥25.0 kg/m2) (AOR = 0.69; 95 % CI: 0.52–0.90) were less likely to be stunted compared with children whose mothers had normal BMI (18.5 kg/m2-24.9 kg/m2). The odds of being stunting were 68 % higher among children with a shorter (<24 months) birth interval compared to children with longer birth interval (≥24 months) (AOR = 1.68; 95 % CI: 1.46–1.93). Children from the poorest wealth quintile had 43 % higher odds of being stunted compared with children from the richest wealth quintile (AOR = 1.43; 95 % CI: 1.08–1.88). Children from male-headed households were 18 % more likely to be stunted compared to children from female-headed households (AOR = 1.18; 95 % CI: 1.01–1.38).

Community level factors

Children from households without access to an improved latrine (defined as private and hygienic with a cleanable slab) had 26 % higher odds of stunting compared to children from households that reported access to an improved latrine(AOR = 1.26; 95 % CI:1.01–1.59). The odds of childhood stunting in Tigray (AOR =1.58; 95 % CI: 1.09–2.29), Amhara (AOR 1.50; 95 % CI 1.04–2.17) and Benishangul-Gumuz(AOR =1.71; 95 % CI: 1.16–2.52) were higher compared to Dire Dawa respectively. However the odds of childhood stunting at Gambella region were lower by 43 % as compared to Dire Dawa (AOR = 0.57; 95 % CI: 0.37–0.86) (Table 5).

As shown in Table 6, the empty model (the null model) revealed that childhood stunting was not random across the communities (τ2 = 0.363, P < 0.001). About 9.9 % of the variance in the odds of childhood stunting could be attributed to the community level factors as calculated by the ICC based on estimated intercept component variance. The full model, after adjusting for individual and community level factors, shows that the variation in childhood stunting across communities remained statistically significant. About 36.6 % of the odds of childhood stunting variation across communities was observed in the full model.

Table 6 Results from random intercept model (measure of variation) for childhood stunting at cluster level by multilevel logistic regression analysis

Moreover, the MOR confirmed that childhood stunting was attributed to community level factors. The MOR for stunting was 2.01 in the empty model; this indicated that there is variation between communities (clustering) since MOR is two times higher than the reference (MOR = 1). The unexplained community variation in stunting decreased to MOR of 1.75 when all factors were added to the null model (empty model). This shows that when all factors are considered, the effect of clustering is still statistically significant in the full model.


This study found that the prevalence of stunting was above the national average of 44 % in six out of the eleven regions of the country. Excess cases of stunting were found in the northern parts of the country characterized by highlands and midlands. This finding is consistent with the study done by Hagos et al. which found that a higher prevalence of stunting was found in the highlands and midlands compared to the lowlands [24]. That study concluded that rainfall and temperature are the main predictors of stunting variation across the country. Other studies also showed that climate change could indeed increase stunting rates in areas of the country dependent on rain-fed agriculture [25, 26]. The spatial variation of childhood stunting is also determined by environmental or geographical factors (e.g. population density, climate and disease environment) in addition to the individual and household level factors [27]. In this study Getis-Ord spatial statistics showed spatial variation of childhood stunting at the cluster level. Statistically significant hotspot areas of childhood stunting were found particularly in the northern parts of the country, in Benishangul-Gumuz, Amhara, Tigray and Affar regions, when we compare to clusters found in other regions. Exploring spatial variation is important to identify aggregations of cases in order to target nutritional interventions [8]. The identified clusters might be the areas where childhood stunting prevention and control interventions should be given priority [28].

This study confirmed that the variation in childhood stunting can be attributed to both individual and community level factors. In the final model, individual and community-level factors accounted for about 36.6 % of the variation observed for childhood stunting. The national average for stunting in children under five years of age in Ethiopia is 44.4 %. Children aged >12 months were at higher odds of being stunted compared to infants of less than one year. The trend of stunting increased for each age group of children age up to 35 months, and then declined to just under 50 % for children aged 48–59 months. This finding is consistent with many previous studies [2934]. These studies reported that there is a rapid fall in children’s height-for-age Z-score from birth to 24 months, particularly during the time when children are being weaned off of exclusive breastfeeding and also becoming more mobile and crawling. These activities expose the child to contaminants in water and food, as well as soil and contaminants picked up on their hands that then go into their mouths. There continues to be an increase in stunting after 24 months, but the rate of increase is much slower. The cumulative effect in older age children might be one possible justification for this pattern. Studies report that children living in most developing countries are introduced directly to the regular household diet made of cereal or starchy root crops, which is the major cause for the high incidence of child malnutrition and morbidity [3537]. Another article found that linear growth failure occurs within a complex interplay of more distal community and societal factors, such as access to healthcare and education, political stability, urbanisation, population density and social support networks [38].

This study found that male children had a slightly higher chance of being stunted than females. This result is consistent with previous studies [31, 3941]. One possible explanation for this could be that childhood morbidity is higher among males than females in early life, even after adjusting for gestational age and body size [42]. Moreover, the proportion of male preterm births is higher than female preterm births which could also contribute to childhood stunting [39, 4244]. However there are also studies which showed that stunting was not significantly associated with sex [32, 45], and programmatically this difference in stunting rates may not be useful for targeting interventions.

Maternal education was found to be significantly negatively associated with childhood stunting. Many previous studies have found that maternal education has a positive effect on reducing childhood stunting [29, 31, 4651]. The knowledge that mothers acquire from formal education could help them to adopt essential nutrition and hygiene behaviors that prevent childhood stunting. Another possible reason might be that educated mothers have better health-seeking behavior for childhood illnesses as compared to uneducated mothers which can help prevent stunting [52, 53]. Besides maternal education, father’s education was also significantly associated with childhood stunting. Those children whose fathers had attended formal education had less chance of being stunted. Similar findings have been reported by other studies [46, 54]. Those fathers with formal education might be more knowledgeable on proper child feeding and hygiene practices, which have a positive effect on preventing childhood stunting. Father’s higher education has also been associated with good health seeking behavior for childhood illness [55]. In Ethiopia lack of education is a widespread problem, about 59 % and 33 % of the women and men were illiterate, respectively which shows how this factor is important in prevention of stunting in Ethiopia [10].

The socioeconomic status of households is associated with access to nutritious foods at the household level, which in turn determines the growth and development of children [56]. Findings of this study revealed that children from households in the poorest wealth quintile were more likely to be stunted than children from the richest quintile, which is consistent with the findings of previous studies carried out in different developing countries [29, 4649, 5759]. This could be due to the fact that increased income improves dietary diversity [60, 61], which in turn improves the adequacy of nutrient intake and nutritional status, and underscores the importance of linking income generating activities with other nutritional interventions. In Ethiopia, above 45 % of the population was found in poor and poorest wealth quintile, which shows the proportion of children who are at risk for stunting due to low socioeconomic status [10].

The 2014 Mini EDHS data show that only 9.7 % of the population had access to an improved and shared latrine facilities [62]. This study found that the availability of an improved latrine was one of the significant factors associated with childhood stunting, which is consistent with other studies conducted in developing countries [41, 63]. One explanation could be a reduction in the pathogen load in the environment from correct and consistent use of improved sanitation, since exposure to feces has been linked to environmental enteropathy and stunting [6466].

A short birth interval is common in Ethiopia and a known risk factor for stunting [67, 68]. Twenty percent of the births have an interval of less than two years, and 9 percent of births are less than 18 months apart. This implies that 20 % of the Ethiopian children were at risk of stunting attributable to short birth interval [10]. This study found that having a birth interval ≥ 24 months decreased the odds of being stunted. This is consistent with other studies [68, 69].A short interval between births can have an adverse effect on child nutrition by causing intrauterine growth retardation, and/or undermining the quality of child care [70].

Birth spacing might also influence childhood under-nutrition through its association with preterm births and low birth weight. If a pregnancy occurs too soon after the previous birth, the mother may not recover her nutritional status, which can contribute to preterm birth and low birth weight [68]. Maternal nutritional status was associated with childhood stunting in this study. The prevalence of under nutrition (BMI < 18.5 kg/m) among women in Ethiopia was 27 % [10]. Mothers who had a BMI <18.5 kg/m2 were marginally associated with higher odds of childhood stunting compared to mothers with normal BMI (18.5 kg/m2 -24.9 kg/m2). Mothers who had BMI ≥25 kg/m2wereless likely to have child who is stunted. A study from Brazil also suggested that maternal nutritional status was associated with child nutritional status [71]. A study from Cambodia also found that maternal BMI was associated with childhood stunting [69].

The representativeness of the data with a large sample size, as well as a nationally representative population-based study with a high response rate are strengths of this study that give high statistical power to infer the characteristics of the study population. Another important strength of this study is the use of multilevel logistic regression analysis, which was able to identify other factors beyond individual-level factors that would not be identified by using standard logistic regression analysis. The use of a combination of methods (spatial and regression statistics) was a strength which allowed validation of the identified hotspots areas due to the assumptions of the statistical methods. One limitation is that the cross-sectional nature of the study means the results cannot be used to establish cause and effect relationship, and is a limitation of the study. Another limitation is the absence of other important variables like behavioral factors and quantitative dietary consumption to substantiate the findings, since the EDHS survey was designed to produce health and health related indices, but was not specifically administered for the purposes of this study.


This study found that childhood stunting was not random in Ethiopia. Statistically significant high hotspots of stunting were found in northern parts of the country whereas low hotspots of stunting were found in central, eastern, western parts of the country. Individual and community level factors accounted for 36.6 % of the variation in childhood stunting across the communities. Both individual- and community-level factors were significant determinants of childhood stunting. Being male, age above 11 months, short birth interval, having anemia, no formal education of mother and/or father, being from male-headed household, and being from a household in the lowest wealth quintile were the factors that increased the odds of stunting at the individual level, whereas a lack of access to an improved latrine and residence in the northern regions of the country were the factors associated with stunting at the community level.

Thus, the improvement of nutritional status of children requires multi-factorial interventions such as reducing poverty and ensuring food security, ensuring adequate birth interval, educating mothers and their partners, a healthy environment and good hygienic practices enabled by access to hygienic latrines. The regions identified as having high hotspots of childhood stunting should be prioritized for nutritional interventions.



adjusted odds ratio


confidence intervals


Ethiopian demographic health survey


geographic information system




intra-cluster correlation


median odds ratio


proportional change in variance


standard deviation


world health organization


  1. WHO(World Health Organization). Physical status: the use and interpretation of anthropometry. Geneva: World Health Organization, 1995; 2006. p. 36. Technical report no.854.

    Google Scholar 

  2. Berkman DS, Lescano AG, Gilman RH, Lopez SL, Black MM. Effects of stunting, diarrhoeal disease, and parasitic infection during infancy on cogni-tion in late childhood:a follow-up study. Lancet. 2002;359:564–71.

    Article  PubMed  Google Scholar 

  3. Victora CG, Adair L, Fall C, Hallal PC, Martorell R, Richter L, Sachdev HS, Maternal and Child Undernutrition Study Group. Maternal and child undernutrition: consequences for adult health and human capital. Lancet. 2008;371(9609):340–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Black RE, Allen LH, Bhutta ZA, Caulfield LE, deOnis M, Ezzati M, et al. Maternal and child undernutrition: global and regional exposures and health consequences. Lancet. 2008;371:243–60.

    Article  PubMed  Google Scholar 

  5. Gluckman PD, Hanson MA, Beedle AS. Early life events and their consequences for later disease: a life history and evolutionary perspective. Am J Hum Biol. 2007;19:1–19.

    Article  PubMed  Google Scholar 

  6. Uauy R, Kain J, Mericq V, Rojas J, Corvalan C. Nutrition, child growth, and chronic disease prevention. Ann Med. 2008;40:11–20.

    Article  CAS  PubMed  Google Scholar 

  7. Black RE, Victora CG, Walker SP, Bhutta ZA, Christian P, de Onis M, et al. Maternal and child undernutrition and overweight in low-income and middle-income countries. Lancet. 2013;382:427–51.

    Article  PubMed  Google Scholar 

  8. Adekanmbi VT, Uthman OA, Mudasiru OM. Exploring variations in childhood stunting in Nigeria using league table, control chart and spatial analysis. BMC Public Health. 2013;13:361.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Central statistical agency -Ethiopia. Ethiopia Demographic and Health Survey. 2000.

    Google Scholar 

  10. Central Statistical Agency (CSA)Ethiopia. Ethiopia Demographic and Health Survey 2011. Addis Ababa, Ethiopia and Calverton, Maryland: CSA and ORC Macro; 2012.

    Google Scholar 

  11. Government of the Federal Democratic Republic of Ethiopia. National Nutrition Programme, June 2013 – June 2015. 2013.

    Google Scholar 

  12. Mengistu K, Alemu K, Destaw B. Prevalence of Malnutrition and Associated Factors Among Children Aged 6–59 Months at Hidabu Abote District, North Shewa, Oromia Regional State. Nutritional Disorders Therapy. 2013;T1:001.

    Google Scholar 

  13. Tariku B, Mulugeta A, Tsadik M, Azene G. Prevalence and Risk Factors of Child Malnutrition in Community Based Nutrition Program Implementing and Nonimplementing Districts from South East Amhara, Ethiopia. Open Access Library J. 2014;1:1–17.

    Google Scholar 

  14. Guo G, Zhao H. Multilevel Modeling for Binary Data. Annual Review Sociology. 2000;26:441–62.

    Article  Google Scholar 

  15. WHO/Multi Growth Resaerch Center Study Group (MGRSG. Complementary feeding in the WHO Multicentre Growth Reference Study. Acta Pediatrica. 2006;450:27–37.

    Google Scholar 

  16. Rutstein SO: The DHS Wealth Index: Approaches for Rural and Urban Areas. DHS Working Papers. Demographic And Health Research 2008(60). Available at Accessed 20 Oct 2013.

  17. Kravdal O. A simulation-based assessment of the bias produced when using averages from small DHS clusters as contextual variables in multilevel models. Demographic Res. 2006;15:1–20.

    Article  Google Scholar 

  18. Ord JK, Getis A. Local Spatial Autocorrelation Statistics: Distributional Issues and an Application. Geographical Analysis. 1995;27(4):286–306.

    Article  Google Scholar 

  19. Halonen JI, Kivimaki M, Pentti J, Kawachi I, Virtanen M, Martikainen P, Subramanian SV, Vahtera J. Quantifying neighbourhood socioeconomic effects in clustering of behaviour-related risk factors: a multilevel analysis. PLoS One. 2012;7(3):e32937.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Merlo J, Chaix B, Yang M, Lynch J, Rastam L. A brief conceptual tutorial of multilevel analysis in social epidemiology: linking the statistical concept of clustering to the idea of contextual phenomenon. J Epidemiol Community Health. 2005;59(6):443–9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, Rastam L, Larsen K. A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health. 2006;60(4):290-297.

  22. Larsen K, Merlo J. Appropriate assessment of neighborhood effects on individual health: integrating random and fixed effects in multilevel logistic regression. Am J Epidemiol. 2005;161(1):81–8.

    Article  PubMed  Google Scholar 

  23. Merlo J, Yang M, Chaix B, Lynch J, Rastam L. A brief conceptual tutorial on multilevel analysis in social epidemiology: investigating contextual phenomena in different groups of people. J Epidemiol Community Health. 2005;59(9):729–36.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Hagos S, Lunde T, Mariam DH, Woldehanna T, Lindtjorn B. Climate change, crop production and child under nutrition in Ethiopia; a longitudinal panel study. BMC Public Health. 2014;14:884.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Grace K, Davenport F, Funk C, Lerner AM. Child malnutrition and climate in Sub-Saharan Africa: An analysis of recent trends in Kenya. Applied Geography. 2012;35:405–13.

    Article  Google Scholar 

  26. Jankowskaa MM, Lopez-Carr D, Funk C, Husak GJ, Chafe ZA. Climate change and human health: Spatial modeling of water availability, malnutrition, and livelihoods in Mali, Africa. Applied Geography. 2011;1–12.

  27. Balk D, Pullum T, Storeygard A, Greenwell F, Neuman M. A spatial analysis of childhood mortality in West Africa. Population Space Place. 2004;10(3):175–216.

    Article  Google Scholar 

  28. Khamis FG, El-Refae GA. Association between Spatial Patterns of Acute Malnutrition and Household Income in Iraq-2004. Int J Statistics Probability. 2012;1(1):43–52.

    Article  Google Scholar 

  29. Kabubo-Mariara J, Ndenge GK, Mwabu DK. Determinants of children’s nutritional status in Kenya: evidence from demographic and health surveys. J African Economies. 2009;18:363–87.

    Article  Google Scholar 

  30. Shrimpton R, Victora CG, de Onis M, Lima RC, Blossner M, Clugston G. Worldwide timing of growth faltering: implications for nutritional interventions. Pediatrics. 2001;107(5):E75.

    Article  CAS  PubMed  Google Scholar 

  31. Adekanmbi VT, Kayode GA, Uthman OA. Individual and contextual factors associated with childhood stunting in Nigeria: a multilevel analysis. Matern Child Nutr. 2013;9(2):244–59.

    Article  PubMed  Google Scholar 

  32. Mushtaq MU, Gull S, Khurshid U, Shahid U, Shad MA, Siddiqui AM. Prevalence and socio-demographic correlates of stunting and thinness among Pakistani primary school children. BMC Public Health. 2011;11:790.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Wand H, Lote N, Semos I, Siba P. Investigating the spatial variations of high prevalences of severe malnutrition among children in Papua New Guinea: results from geoadditive models. BMC Res Notes. 2012;5:228.

    Article  Google Scholar 

  34. Fenske N, Burns J, Hothorn T, Rehfuess EA. Understanding Child Stunting in India: A Comprehensive Analysis of Socio-Economic, Nutritional and Environmental Determinants Using Additive Quantile Regression. PLoS One. 2013;8(11):78692.

    Article  Google Scholar 

  35. Dewey K, Brown K. Update on Technical Issues Concerning Complementary Feeding of Young Children in Developing Countries and Implications for Intervention Programs. Food Nutr Bull. 2003;24:5–28.

    Article  PubMed  Google Scholar 

  36. Krebs N, Westcott J. Zinc and Breastfed Infants: If and When Is There a Risk of Deficiency? Adv Exp Med Biol. 2002;503:69–75.

    Article  CAS  PubMed  Google Scholar 

  37. WHO/UNICEF. Complementary Feeding of Young Children in Developing Countries: A Review of Current Scientific Knowledge. Geneva: WHO/NUT/981 World Health Organization; 1998. p. 1–228.

    Google Scholar 

  38. Stewart CP, Iannotti L, Dewey KG, Michaelsen KF, Onyango AW. Contextualising complementary feeding in a broader framework for stunting prevention. Matern Child Nutr. 2013;9 Suppl 2:27–45.

    Article  PubMed  Google Scholar 

  39. Ingemarsson I. Gender aspects of preterm birth. BJOG. 2003;110 Suppl 20:34–8.

    Article  PubMed  Google Scholar 

  40. Wamani H, Astrom AN, Peterson S, Tumwine JK, Tylleskar T. Boys are more stunted than girls in sub-Saharan Africa: a meta-analysis of 16 demographic and health surveys. BMC Pediatr. 2007;7:17.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Keino S, Plasqui G, Ettyang G, van den Borne B. Determinants of stunting and overweight among young children and adolescents in sub-Saharan Africa. Food Nutr Bull. 2014;35(2):167–78.

    Article  PubMed  Google Scholar 

  42. Elsmen E, Hansen Pupp I, Hellstrom-Westas L. Preterm male infants need more initial respiratory and circulatory support than female infants. Acta Paediatrica. 2004;93(4):529–33.

    Article  CAS  PubMed  Google Scholar 

  43. Escobar GJ, Clark RH, Greene JD. Short-term outcomes of infants born at 35 and 36 weeks gestation: we need to ask more questions. Semin Perinatol. 2006;30(1):28–33.

    Article  PubMed  Google Scholar 

  44. Kilbride HW, Daily DK. Survival and subsequent outcome to five years of age for infants with birth weights less than 801 grams born from 1983 to 1989. J Perinatology. 1998;18(2):102–6.

    CAS  Google Scholar 

  45. Bloss E, Wainaina F, Bailey RC. Prevalence and Predictors of Underweight, Stunting, and Wasting among Children Aged 5 and Under in Western Kenya. J Trop Pediatr. 2004;50(5):264.

    Article  Google Scholar 

  46. Semba RD, de Pee S, Sun K, Sari M, Akhter N, Bloem MW. Effect of parental formal education on risk of child stunting in Indonesia and Bangladesh: a cross-sectional study. Lancet. 2008;371(9609):322–8.

    Article  PubMed  Google Scholar 

  47. Monteiro CA, Benicio MH, Konno SC, Silva AC, Lima AL, Conde WL. Causes for the decline in child under-nutrition in Brazil, 1996–2007. Rev Saude Publica. 2009;43(1):35–43.

    Article  PubMed  Google Scholar 

  48. Van de Poel E, Hosseinpoor AR, Jehu-Appiah C, Vega J, Speybroeck N. Malnutrition and the disproportional burden on the poor: the case of Ghana. Int J Equity Health. 2007;6:21.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Hong R, Hong R. Economic inequality and undernutrition in women: multilevel analysis of individual, household, and community levels in Cambodia. Food Nutr Bull. 2007;28(1):59–66.

    Article  PubMed  Google Scholar 

  50. Makoka D. The Impact of Maternal Education on Child Nutrition: Evidence from Malawi, Tanzania, and Zimbabwe. DHS working Papers. 2013;84.

  51. Urke HB, Bull T, Mittelmark MB. Socioeconomic status and chronic child malnutrition: Wealth and maternal education matter more in the Peruvian Andes than nationally. Nutrition Res. 2011;31(10):741–7.

    Article  CAS  Google Scholar 

  52. Tessema F, Asefa M, Ayele F. Mothers’ Health Services Utilization and Health Care Seeking Behavior During Infant Rearing: A Longitudinal Community Based Study, South West Ethiopia. Ethiopian J Health Development. 2002;16:51–8.

    Article  Google Scholar 

  53. Assefa T, Belachew T, Ayalew T, Deribew A. Mothers’ Health Care Seeking Behavior For Childhood Illnesses In Derra District, Northshoa Zone, Oromia Regional State, Ethiopia. Ethiopian J Health Development. 2008;18(3):90–3.

    Google Scholar 

  54. Moestue H, Huttly S. Adult education and child nutrition: the role of family and community. J Epidemiol Community Health. 2008;62(2):153–9.

    Article  CAS  PubMed  Google Scholar 

  55. Najnin N, Bennett CM, Luby SP. Inequalities in Care-seeking for Febrile Illness of Under-five Children in Urban Dhaka, Bangladesh. J Health Population Nutrition. 2011;29(5):523–31.

    Article  Google Scholar 

  56. Nelson M. Childhood nutrition and poverty. Proc Nutr Soc. 2000;59(2):307–15.

    Article  CAS  PubMed  Google Scholar 

  57. Hong R, Banta JE, Betancourt JA. Relationship between household wealth inequality and chronic childhood under-nutrition in Bangladesh. Int J Equity Health. 2006;5:15.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Hong R, Mishra V. Effect of wealth inequality on chronic under-nutrition in Cambodian children. J Health Popul Nutr. 2006;24(1):89–99.

    PubMed  Google Scholar 

  59. Hong R. Effect of economic inequality on chronic childhood undernutrition in Ghana. Public Health Nutr. 2007;10(4):371–8.

    Article  PubMed  Google Scholar 

  60. Taruvinga A, Muchenje V, Mushunje A: Determinants of rural household dietary diversity: The case of Amatole and Nyandeni districts, South Africa. International Journal of Development and Sustainability. 2013;2(4):1-15.

  61. Doan D. Does income growth improve diet diversity in China?Selected Paper prepared for presentation at the 58 the AARES Annual Conference, Port Macquarie, New South Wales. 2014.

    Google Scholar 

  62. Central Statistical Agency [Ethiopia]: Ethiopia Mini Demographic and Health Survey 2014. 2014. p. 56-58. Available at Accessed 20 Sept 2015.

  63. Masibo PK, Makoka D. Trends and determinants of undernutrition among young Kenyan children: Kenya Demographic and Health Survey; 1993, 1998, 2003 and 2008–2009. Public Health Nutr. 2012;15(9):1715–27.

    Article  PubMed  Google Scholar 

  64. Prendergast A, Kelly P. Enteropathies in the developing world: neglected effects on global health. Am J Trop Med Hyg. 2012;86(5):756–63.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Lin A, Arnold BF, Afreen S, Goto R, Huda TM, Haque R, Raqib R, Unicomb L, Ahmed T, Colford JM Jr., et al. Household environmental conditions are associated with enteropathy and impaired growth in rural Bangladesh. Am J Trop Med Hyg. 2013;89(1):130–7.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Spears D. How Much International Variation in Child Height Can Sanitation Explain? Policy Research Working Paper 6351.World Bank. 2013.

    Google Scholar 

  67. Rutstein SO. Effects of preceding birth intervals on neonatal, infant and under-five years mortality and nutritional status in developing countries: evidence from the demographic and health surveys. Int J Gynaecology Obstetrics. 2005;89 Suppl 1:S7–24.

    Article  Google Scholar 

  68. Gribble JN, Murray NJ, Menott EP. Reconsidering childhood undernutrition: can birth spacing make a difference? An analysis of the 2002–2003 El Salvador National Family Health Survey. Matern Child Nutr. 2008;5(1):49–63.

    Article  Google Scholar 

  69. Ikeda N, Irie Y, Shibuya K. Determinants of reduced child stunting in Cambodia: analysis of pooled data from three Demographic and Health Surveys. Bull World Health Organ. 2013;91:341–9.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Dewey KG, Cohen RJ. Does birth spacing affect maternal or child nutritional status? A systematic literature review. Maternal Child Nutrition. 2007;3:151–73.

    Article  PubMed  Google Scholar 

  71. Felisbino-Mendes M, Villamor E, Velasquez-Melendez G. Association of maternal and child nutritional status in Brazil: a population based cross-sectional study. PLoS One. 2014;24(9):87486.

    Article  Google Scholar 

Download references


The data used in this study were obtained from MEASURE DHS Archive and the authors acknowledge MEASURE DHS for allowing us to use the data. The analysis and the interpretation of the findings in this study are the responsibility of authors and do not necessarily reflect the views of MEASURE DHS. We thank Mr. Tesfaye Sorresa for editing the language and typological errors of this manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Muluken Azage.

Additional information

Competing interest

The authors declare that they have no competing interest.

Authors’ contributions

DH conceived the study, assisted the data analysis and interpretation, drafted the manuscript and critically reviewed the manuscript. MA conceived the study, conducted the data analysis and interpretation, and critically reviewed the manuscript. TM assisted the data analysis, critically reviewed the manuscript. RR critically reviewed the manuscript and did final editing. All of the authors read and approved the final version of the manuscript

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Haile, D., Azage, M., Mola, T. et al. Exploring spatial variations and factors associated with childhood stunting in Ethiopia: spatial and multilevel analysis. BMC Pediatr 16, 49 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Stunting
  • Children
  • Ethiopia