Analyzing determinants from both compositional and contextual level impeding desired linear growth of children in Indian context

Abstract

Childhood stunting is recognized as significant public health concern in India. It is a form of malnutrition with impaired linear growth and creates a range of adversaries among children, including under-5 mortality, morbidity, and physical and cognitive growth. The purpose of the present study was to recognize the various leading determinants causing childhood stunting from both individual and contextual level in Indian context. Data were obtained from the India’s Demography and Health Survey (DHS) conducted in 2019–2021. A total of 1, 46,521 children aged 0–59 months were included in this present study. The study applied a multilevel mixed-effect logistic regression model in which individual factors nested within community based contextual-level factors estimating the likelihood of childhood stunting phenomena among Indian children. The variance explained in full model accounted for about 35.8% of the odds of stunting across the communities. The present study elucidates that the recognized factors from individual level characteristics have really increased the odds of childhood stunting: gender of child, multiple births, low birth weight, low BMI among mothers, less educational attainment by mothers, maternal anemic status, breast feeding duration longer than usual, < 4 antenatal care (ANC) visits during gestation period. Similarly, contextual-level factors like rural places of residence, western Indian children, and communities with high poverty rates, lower literacy rates, improper sanitation, and unsafe drinking water were also found to have a significant positive association with childhood stunting. The study finally concludes that cross level interaction between individual and contextual-level factors are identified as significant determinants of linear growth retardation among child in India. In order to reduce this type of malnutrition among the child one should more concentrate on both individual and contextual-level factors as a notable reasons.

Introduction

Childhood stunting is one of the three anthropometric faces of malnutrition and a serious public health concern. It is an indication of severe, irreversible physical, psychological, and cognitive impairment among children caused by early-life chronic malnutrition. Compared to the developed countries, the increasing trend of childhood stunting is worst in developing countries [1] and India is not an exception. The National Family and Health Survey (2019–21) estimated that approximately 36% of Indian under-five children are suffering from stunting (short height for age). However, the rate of childhood stunting has decreased by 2%, according to NFHS-5 estimates; the phenomenon of childhood stunting is still alarming in India. After the COVID-19 scenario, the wrath of pandemic’s impact on childhood stunting would become terrible and could persist for years, as estimated by UNICEF [2]. According to UNICEF, WHO, and the World Bank Group's joint programme 2021, there are 149.2 million under-five children who are stunted worldwide and constitute 18.95% out of all children. The World Health Organization set a goal in 2012 to reduce the prevalence of stunting by up to 40% by 2025 [3]. In the year 2015, the United Nations Sustainable Development Goals (SDGs) set the goal of “ending all forms of malnutrition globally by 2030” [4]. Despite setting the several worldwide goals for reducing stunting phenomena, India still bears about 30.90% of the total burden of childhood growth retardation as stated by UNICEF [5].

Numerous studies conducted in India at the national and sub-national levels have provided shreds of evidence that a number of factors were associated with childhood stunting. The underlying risk factors which are identified as dominant factors shaping childhood stunting phenomena are the age of the child, male gender [6], higher birth order [7, 8], severe anemia level [9], delayed introduction of breastfeeding [10], small child size at birth [11], unsafe stool disposal practice of child [12], lower maternal age [13], lesser maternal education [14, 15], maternal underweight [13, 16], not attending recommended antenatal care visits [17], maternal short height [18]. Concerning, socio-economic & environmental factors which are comprised with rural setting of residence [19], the Muslim contextual, backward castes [20], households with lower income [21], inadequate or no toilet facilities, unimproved drinking water sources [22] and are highly associated with childhood stunting. An optimum child's growth may be hindered and become irreversibly stunted if they do not receive appropriate nutrition during their first 1000 days of life [23]. Maternal height plays a significant role to avert the stunting phenomena because stunting is a cyclical process, and mothers who themselves were stunted in their childhood are more likely to have stunted children. This leads to a transmission of cycle of poverty from one generation to next generation and diminished human capital which is hard to overcome [24]. Notably, stunting has a strong association with socio-economic & environmental factors of children, as children from low-income families are more likely to be stunted due to lack of improved defecation places, taking unsafe drinking water, as compared to high-income families [22, 25, 26]. Children who have been failure to achieve the ideal height due to linear growth anomaly would ultimately have less working capacity in the near future and lose their capability to work in contrast to other people [15].

There are many earlier researches which have determined the prevalence of childhood stunting and analyzed socio-economic, demographic and cultural factors associated with childhood stunting in India [27,28,29,30]. However, most of the existing studies employed a standard binary logistic regression model to identify the individual predictors of childhood stunting. Sometime analysis of individual factors from different levels at one single common level by applying a standard binary logistic regression model leads to biased results (Loss of power or type I error). Furthermore, with the presence of differences in regional settings and other contextual differences in the study of malnutrition becomes complicated. Household and individual level characteristics within same geographic region share common contextual characteristics such as seasonal or climatic variability, types of crops and housing which can have a similar impact on the nutritional status of children within the same cluster. The assumptions of independence among individuals within the similar homogeneous cluster are maintained in hierarchical data structure and sometime equal variance across clusters are violated in the case of cluster arrangement of data or nested data. Hence, the suitable statistical technique for such kind of study must be selected as hierarchical modeling or multilevel analysis. According to Guo and Zhao, multilevel logistic regression analysis has several benefits over traditional logistic regression analysis [30].

Unfortunately, none of the previous studies in the existing literature has highlighted the possible contextual effect of impact of location, homogeneous cluster of socio-economic characteristics and environmental hygiene on childhood stunting in Indian context. In developing countries, the likelihood of the impact of contextual effect on malnutrition is quiet high [31]. From this backdrop, this study aims to identify the significant factors of childhood stunting from both individual and contextual level hierarchical dataset.

Study design

Data

This study utilized data from the 5^th round of the National Family Health Survey (NFHS-5), conducted in 2019–2021. This population-based, large-scale, nationally representative sample survey, collected data from all 707 districts in 29 states and 7 Union Territories and gathered information from 636,699 households with a response rate of 98%, while724, 115 women aged 15–49 years with a response rate of 97%, and 101, 839 men aged 15–54 years with a response rate 92%. The information related to child (0–59 months) health indicators were taken from the mothers. The survey was conducted by the International Institute for Population Science (IIPS) under the stewardship of the Ministry of Health and Family Welfare (MoHFW), Government of India. The multi-stage sampling design was adopted to collect information from the respondents. Firstly, the country was divided into rural and urban parts. The rural sample was selected through a two-stage sample design where the villages were treated as the primary sampling units (PSUs) at the first stage (Selected with probability proportional to size), followed by a random selection of 22 households from each PSU at the second stage. In the case of urban areas, there was also a two-stage sampling design with Census Enumeration Blocks (CEB) selected at the first stage and a random selection of 22 households from each CEB at the second stage. In the second stage in both rural and urban areas, households were selected after conducting a complete mapping and household listing operation that was selected at first-stage units.

Study participants

NFHS-5 collected information for 2, 32,920 children aged 0–59 months. Child-related information was recorded from the response of their mother. Children’s information paired with the mother’s information as individual level factors and location, literacy and poverty rate at a glance and neighbor hood’s hygiene from both PSU and CEB as community level contextual factors was taken in the backdrop of childhood stunting for every child. However, there was missing information for 8,702 children. After removing all missing information and pairing children’s information with mothers' and including contextual information the full information was available for 146,521 children. That’s why the final study participants for the present study were 1, 46,521.

Outcome variable

Childhood stunting (0–59 months) was considered as the only outcome of interest for the present study. Linear growth retardation of children usually called childhood stunting was measured in height for age during the survey. According to the World Health Organization (WHO) stunting is defined as height for age < -2 standard deviation and severe stunting defined as height for age < -3 standard deviation from the median age (both fall left tail of standard normal distribution) as recommended by WHO in 2006 as child growth standards [32]. In this study children are coded as ‘1’ if they were stunted while not stunted children are coded as’0’. In this study both the stunted and severe stunted children are considered as stunted. Though the indicator was simplistic and the definition reflected on the assessment of the dynamics of growth and could be used consistently over the whole growth trajectory.

Explanatory variables

Several explanatory variables were included in the present study and considered as potential predictors of childhood stunting. The variables were taken from three different levels i.e. (i) Individual level (compositional variable) (ii) Household level and (iii) Contextual level (regional level). The detailed definition and coding for each explanatory variable have given in Table 1. Here the study used the term contextual level factor to represent the homogeneous clusters within the same geographical living environment and sharing the similar socio-economic conditions. Contextual factors (location, socio-economic condition and environmental hygiene) were chosen based on sharing a common PSU within the DHS data. The clusters were taken as a unit of analysis because of two reasons. First, in the DHS sample, PSUs were used to define the clusters based on homogeneous characteristics. Secondly, it has been shown that most of the DHS is conducted by taking the sample size per PSU that can meet the optimum size with a tolerable precision loss.

Table 1 Selection and description of the variables

Statistical analysis

At first univariate analysis was carried out to describe about the participant’s characteristics and after that the bivariate analysis along with Chi-square test statistics was performed to test the association of the sample distribution between the dependent and independent variables. The survey weights were applied throughout the analysis due to the hierarchal nature of the survey and to represent the results in accordance with the country’s population. Finally, a bi-variable and multivariable multilevel logistic regression model was performed to identify the factors that were significantly associated with childhood stunting. For the multilevel analysis, a total of four models were specified in the current study. The first model, commonly known as the empty or null model, was considered without any explanatory variable i.e. simple component or variance analysis. The first model also evaluates the extent of the cluster variation on stunting. The second model namely, the fixed effect model controlled for the individual level variables. The third model also known as the fixed effect model and controlled for the contextual-level variable without considering any individual-level variable and the final model namely, the random slope or random coefficient model controlled for both individual-level variable and contextual-level variables simultaneously. The results of the random slope model were reported as an adjusted odd ratio (aOR) at a 95% confidence interval (95% CI). The intra-class correlation coefficient (ICC) was used to report the variation of stunting at higher contextual level.

Ethical statement

The ethical approval of the NFHS-5 (2019–2021) was obtained from the ethics review board of the International Institute for Population Science (IIPS), Mumbai. The survey was also reviewed and approved by the ICF and International Review Board (IRB). The study also confirmed that all methods were performed in accordance with the relevant guidelines and regulations. For participation in this survey, informed written and verbal consent was obtained from the participants at the time of the survey. Each participant’s approval was sought, and after that, the interview was conducted. The study involved secondary analysis of publicly available data. Thus independent ethical approval was not needed. However, the corresponding author received permission from www.dhsprogram.com to use the data for analysis.

Results

Participant’s characteristics

A total of 146,521 children aged 0–59 months were included in the present study. Of whom, 53.75% were male and 46.25% are females. Among the total, more than 70% of children were single-born babies in a family. Most of the children (84.09%) had a normal birth weight that is more than or equal to 2500 g at birth. More than half of (58.46%) the children’s mother belongs to the age bracket between 25 to 34 years. Nearly 70% of mothers had a formal school education up to the secondary level and above (51.94 + 16.91) secondary level. The breastfeeding duration of the children varied by month, for example, one in every fourth of children had breastfeeding duration of at least 1 year while 23.87% of children had breastfeeding duration of more than two years. More than 60% of mothers had reported their BMI as between 18.5–24.9 kg/m². About 60% of mothers had visited the antenatal clinic at least 4 times or more than that and the remaining were taken antenatal care less than 4 times. Among all mothers, 40% reported that they are not anemic and the rests of the mothers were mild to severely anemic. The majority of the participants (88.25%) were from Scheduled Caste (SC) communities and male-headed families (85.01%). Participants are near about equally distributed in different wealth quintile categories. The majority of the study participants were from rural areas (72.76%). The participant hails from different Indian geographical regions like North (13.77%), South (16.52%), East (26.44%), West (12.83%), and North-East (3.96%) etc. About 23% respondents were from below poverty level and 32.32% mothers were from low literacy rate. About 16.43% respondents use to consume unsafe drinking water. The details of the participant’s characteristics have given in Table 2.

Table 2 Socio-demographic and economic characteristics of the respondents

Prevalence of childhood stunting

A total of 36% of children aged between 0–59 months were identified stunted from different parts of India. A significantly higher prevalence of childhood stunting was observed among male children (male: 35.50%; female: 32.13%; P: < 0.05) in comparison to female children. However, the prevalence of childhood stunting was found equally distributed in different age groups continuing with a certain interval but a significant difference was there. Childhood stunting was slightly higher in the case of multiple births (multiple births: 40.23%; Single birth: 33.88%; P < 0.05) and with low birth weight (BW) (BW < 2500:41.65%; BW ≥ 2500: 32.57; P < 0.05) babies. The proportion of childhood stunting is significantly highest in children born to mothers with no formal education (42.79%); low Body Mass Index (BMI) (41.52%); surviving with severe anemia (39%); taken less than 4 antenatal care treatment (36.58%) during pregnancy. The prevalence of childhood stunting was significantly highest among children whose parents were from the poorest (42.68%) background with respect to other wealth quintile categories. The proportion of childhood stunting was found higher in rural areas in comparison to urban (rural: 35.32%; urban: 28.47%; P < 0.05) counterparts. Prevalence of childhood stunting was also significantly higher in children born to communities with low poverty rate; (42.68%) low literacy rate (41.59%); no proper sanitation (36.56%) and consume unsafe drinking water (34.39%). The detail description about the prevalence of childhood stunting has been given in Table 3.

Table 3 Childhood stunting by compositional and contextual level

Compositional and contextual level factors of childhood stunting

Table 4 represents the result of Multi-level model for both individual level and contextual level factors. As mentioned above Model 1 in Table 4 was an empty model while Model 2 & 3 was step wise fixed effect model for individual level and contextual level factors and finally Model 4 represents a complete model adjusted for both individual level and contextual level factors at a time. Children’s age was statically significantly associated with the odds of childhood stunting. Male children were more likely to be stunted compared to the female children (aOR: 1.14; 95% CI: 1.11 – 1.17; P < 0.001). Children who were the products of multiple births were 26% more likely to be stunted (aOR: 1.26; 95% CI: 1.14 – 1.41; P < 0.001). Children born with low birth weight (< 2500 g) were 44% more likely (aOR: 1.44; 95% CI: 1.39 – 1.54; P < 0.001) to be stunted in respect to normal birth weight children. Increases with maternal age and level of education the probability of childhood stunting has decreased gradually. Result revealed that the odds of a child being stunted increases with increasing breastfeeding duration after 6 months. Children of mothers with BMI ≥ 25.0 kg/m²were 32% less likely (aOR: 0.68; 95% CI: 0.65 – 0.70; P < 0.001) to be stunted compare to the children of mothers with low BMI (< 18.5 kg/m²). Children born to mother with more than 4 ANC visits during pregnancy were 10% less likely to be stunted (aOR: 0.90; 95% CI: 0.87 – 0.95; P < 0.001). Children born to anemic mothers were more likely to be stunted than the children born to non-anemic mothers. Likewise children born to SC and ST families were more likely to be stunted than the children born to non-SC& ST families. The odds of being stunted decreases with increasing wealth index level. Although these are the factors discussed above were from individual or single level and household context but there has some contextual effects among them because the odd ratios were not identical in fixed effect models and random slope model or adjusted model.

Table 4 Model without and with cross level interaction; compositional and contextual level factors with childhood stunting (aged 0–59 Months) in India

When all other factors were controlled for in random slope model, contextual-level factors like residential place, geographical location, literacy rate, poverty condition, accompanied with sanitation and drinking water facilities were significantly associated with childhood stunting. Children from rural communities (aOR: 0.90; 95% CI: 0.87 – 0.95; P < 0.001) were less likely to be associated with childhood stunting compared to the children from urban communities but the result was different before adjusting with the Individual variables. The odds of being stunted were higher among the children from Western India (aOR: 1.11; 95% CI: 1.04 – 1.17; P < 0.001) compared to the children from the South Indian region. Children from communities with high poverty rates were 68% more likely to be associated with childhood stunting. The odds of being stunted were lower (aOR: 0.53; 95% CI: 0.39 – 0.61; P < 0.001) among the children who were from the communities with high literacy rates compared to the children belonging to the communities with low literacy rates. The existence of improper sanitation and unsafe drinking water facilities within the community were also significantly and positively associated with childhood stunting in India.

Proportion of total variance of stunting phenomena experienced on contextual level

As shown in Table 5, with respect to the empty model there was a significant variation in the odds of childhood stunting (variance: 0.416; P: < 0.001). The ICC pointed out that the estimated intercept component variance showed that an 11.5% variance in the odds of childhood stunting could be attributed to contextual-level determinants. However, after inclusion of contextual level determinants (11.5–5.2) 6.3% variation still exists due to any other contextual level determinants. The variation of childhood stunting across clusters remained statistically significant, even after controlling for individual-level and contextual-level factors (random slope model of the final model). As determined by the percentage change in variance of childhood stunting, the full model or random slope model accounts for about 35.8% of the odds of childhood stunting across the communities. The stepwise decreasing nature of DIC (-2 Loglikelihood) represents that the entire model has improved subsequently and represents a statistically significant model (Model 1 DIC: 92,826.87; Model 4 DIC: 86,153.18).

Table 5 Measures of variation and model fitting statistics

Discussion

Using data from nationally representative surveys, the study analyzed individual or single level factor from children, mothers and household level and contextual-level factors from the homogeneous areal PSU unit considering location, literacy and poverty as well as neighborhood’s hygiene as significant predictors estimating childhood stunting. It elaborates the significance of contextual-level variability in childhood stunting with respect to individual-level variability. About 36% of the variations (35.8%) in stunting were explained by individual and contextual-level characteristics in the final model. Children between the ages of 12 months and 23 months bear the highest risk of childhood stunting, which usually decline after 23 months. This finding has been validated by previous findings from other developing countries [32,33,34] including India. It indicates that children's height rapidly decreased from the time of birth to 23 months; however, stunting processes continued after 23 months, but at a considerably slower rate. This might be due to weaning and decreased breast milk intake and not enough supplementary dietary intakes in terms of frequency and variety which predisposes them to childhood stunting [15]. According to the study's findings, stunting in children under the age of five is strongly influenced according to child's sex. Although some of the earlier research demonstrated that, there is no valid link between a child's sex and childhood stunting. While some research revealed that stunting was more prevalent in male children; on the other hand few research also discovered that stunting was more common in female children. However, a recent systematic review and meta-analysis from South-East Asian countries revealed that in the majority of South-East Asian countries, male children were more likely to be stunted compared to female children [6]. Results from this recent study are consistent with those from earlier research [8, 25], which found that compared to single births, children belonging to the multiple birth cohort group are more likely to be stunted. Such kind of linear growth retardation in accordance with age are usually linked to poor nursing, low birth weight, and due to lack of sufficient nutritional intake, all of which are usually happened to multiple-birth children than in single-birth children. The results of this study corroborate those of related studies that found that mother’s educational attainment plays a crucial role as protective factor against childhood stunting [14, 15]. Mothers with formal education may be aware of what to door to avoid childhood stunting or to minimize its severity if it already exists. Desired antenatal care visits during the gestation period are really essential for both women and child and it was really consistently associated with childhood stunting. The risk of stunting was higher in children whose mothers had taken fewer than at least four professional prenatal care visits during pregnancy. This result of our study is consistent with previous finding [17, 35]. Children from the women were less likely to be malnourished if they had taken at least four ANC visits, and was able to discuss pregnancy and childbirth-related matters as well as optimum growth of foetus with professionals and their husband [35]. In contrary, those women who hardly followed any antenatal care visit because of lesser educational attainment and low household income, the foetus is more likely to be malnourished and the newborn offspring is more likely to be stunted as they are deprived of all the prenatal maternal precaution. This study found a substantial correlation between childhood stunting and mother's BMI, indicating that children under five whose mothers had low BMI are more likely to be prone to childhood stunting. A similar finding has been documented in some studies [14, 15]. To the best of our knowledge, we were unable to discover any literature to support the claim that the mother's height may selected as crucial measurement for evaluating this association. In this study, prolonged breastfeeding, especially after infancy (> 12 months), was found as pertinent one of childhood stunting. This conclusion may be explained by the fact that most households in developing nations are underprivileged and unable to provide their children with a sufficient diet of complementary foods of high quality. As a result, they continue to breastfeed exclusively for longer than 12 months without adding any additional foods [36,37,38]. Another theory is that some infants who are breastfed for longer than usual periods of time reject diets other than breast milk, as demonstrated by a Ghanaian study [39]. Additionally, this study's findings also revealed that under-five children from less affluent households have a higher likelihood of being stunted than under-five children belonging from wealthy households. Similar findings from earlier studies conducted in other developing nations provide more evidence that poverty is a strong predictor of childhood stunting [15, 18, 31, 40]. The socioeconomic condition of a household is important since it frequently affects the availability of healthy foods for children's growth and development. Additionally, under-fives from lower-income homes in the majority of developing nations, like India, where a publically funded health care system is not practiced, would not have access to quality and essential health care when they became unwell. The study revealed that maternal anemia significantly associated with childhood stunting. This could be due to children of anemic mothers experience slower growth because of their moms' severe reduction in working ability, which makes it difficult for them to care for their children and do household duties [22]. Therefore, when anemic moms are unable to provide adequate attention and care owing to own poor health, the harmful effects of maternal anemia that also impact on the infant's linear growth retardation expand and intensify.

The multilevel technique used in this study creates an opportunity to identify also the contextual level effects (area or neighborhood factors) rather than only compositional level effects (individual characteristics effects) entirely to account for variations in childhood stunting among different communities or clusters. Despite this, it has been discovered that differences in individual or single level characteristics account for the majority of the differences in childhood stunting rather than difference between various contextual characteristics. Moreover, the study also declared that contextual or regional characteristics also play a pivotal role in growing childhood stunting phenomena. These results are in line with the majority of research efforts to separate compositional effects from contextual effects [40,41,42]. The approach of multilevel analysis can help to avoid the traditional problem in single level analysis of likelihood of childhood stunting because it has aggregated individual data to the group level to conduct a group level analysis or disaggregated the group level data to conduct the individual analysis. Therefore, the analysis of likelihood of outcome of childhood stunting might be rational through compromise of both single or individual level and contextual level factors. Hence, the changing odds in the final model also signify the contextual-level factors, on their own, exert a significant influence on individual-level factors.

Conclusion

The study's conclusions have some important policy ramifications. Intervention is urgently required in lessening economic disparities, which will ultimately lead to poverty among the population. In this regard both short and long term approaches are required. Adult literacy programmes should be expanded across the communities, prioritizing the objectives to enhance the awareness regarding adequate nutritional requirement for both maternal and children’s health. As a long-term strategy, governments of developing nations should also make sure that women from illiterate communities receive a suitable formal education. All these footsteps may consider as an effective effort to enhance the awareness of women from different groups that they can better nourish their children during childhood period. Places, marked with a high rate of childhood stunting should get the priority to initiate any governmental and non-governmental organization’s project that aims to improve the nutritional status of children under-five years ultimately would brought the entire regional level within its ambit. The establishment of nutritional intervention programmes, such as the therapeutic feeding and nursery school feeding programmes, should be prioritized with a focus on higher risk groups, such as male children under the age of five, children with uneducated or less educated mothers, children of multiple births, and children who live in areas with a high illiteracy rate. Those children who are currently experiencing childhood stunting, the necessary interventions are required to minimize the severity of linear growth retardation. Considering and minimization of public service fees, particularly for health care, at the behest of less affluent members of society who are unable to pay for essential health care services will facilitate in reducing the severity of childhood stunting in the patients. Moreover, provision of free health care for pregnant mothers, new mothers and young children need to be carried out uninterruptedly. Finally, it is important to perform future studies that will take into account morbidity variables as well as causal linkages.