Background Malnutrition is a major cause of child death, and many children suffer from acute and chronic malnutrition. Nigeria has the second-highest burden of stunting globally and a higher-than-average child wasting prevalence. Moreover, there is substantial spatial variation in the prevalence of stunting and wasting in Nigeria. This paper assessed the socioeconomic inequalities and determinants of the change in socioeconomic inequalities in child stunting and wasting in Nigeria between 2013 and 2018. Methods Data came from the 2013 and 2018 Nigeria Demographic and Health Survey. Socioeconomic inequalities in stunting and wasting were measured using the concentration curve and Erreygers’ corrected concentration index. A pro-poor concentration index is negative, meaning that the poor bear a disproportionately higher burden of stunting or wasting than the wealthy. A positive or pro-rich index is the opposite. Standard methodologies were applied to decompose the concentration index (C) while the Oaxaca-Blinder approach was used to decompose changes in the concentration indices (ΔC). Findings The socioeconomic inequalities in child stunting and wasting were pro-poor in 2013 and 2018. The concentration indices for stunting reduced from -0.298 (2013) to -0.330 (2018) (ΔC = -0.032). However, the concentration indices for wasting increased from -0.066 to -0.048 (ΔC = 0.018). The changes in the socioeconomic inequalities in stunting and wasting varied by geopolitical zones. Significant determinants of these changes for both stunting and wasting were changes in inequalities in wealth, maternal education and religion. Under-five dependency, access to improved toilet facilities and geopolitical zone significantly explained changes in only stunting inequality, while access to improved water facilities only significantly determined the change in inequality in wasting. Conclusion Addressing the socio-economic, spatial and demographic determinants of the changes in the socioeconomic inequalities in child stunting and wasting, especially wealth, maternal education and access to sanitation is critical for improving child stunting and wasting in Nigeria.
Data come from the 2013 and 2018 rounds of the nationally representative Nigeria Demographic and Health Survey (NDHS), the two most recent rounds of the NDHS. A stratified three-stage cluster design was used for the 2013 NDHS, comprising 904 primary sampling units (PSUs) which served as the clusters: 372 urban and 532 rural. The 2013 sample selected 40,680 households, where a minimum of 943 interviews were completed in each state. (Nigeria has 36 states and an autonomous Federal Capital Territory). A more detailed account of the dataset is also available elsewhere [10]. The 2018 NDHS was conducted between August and December 2018 via a two-stage stratified cluster sampling design, with each cluster or PSU defined based on enumeration areas from the 2006 population census frame. The survey, which was billed to be conducted in 1 400 clusters comprising 580 urban and 820 rural clusters, actually took place in 1 389 clusters due to feasibility issues. In total, 41 821 (13 311) eligible women (men) were interviewed as part of the 2018 NDHS. A more detailed description of the 2018 NDHS sampling design is available elsewhere [11]. Both datasets allow for the consistent measurement of indicators at the national, zonal (i.e. the six geopolitical zones) and state levels, while survey weights were appropriately calculated to ensure national representativity. This paper used the Birth Recode data file, which contains information on every child ever born to each interviewed woman (up to a maximum of 20 children). This data file also links children to their mothers, and it contains other relevant household information. For each NDHS round, the analysis was restricted to the sample of children aged 0–59 months. After data cleaning, the final sample sizes were 23,992 and 11 150 children in the 2013 and 2018 NDHS, respectively. The outcome variables are indicators of child stunting and wasting. Both follow their respective standard definitions. A child with a height-for-age z-score that is less than the negative of twice the standard deviation of the WHO Child Growth Standards median is considered stunted. A child is wasted if the weight-for-height z-score is less than the negative of twice the standard deviation of the WHO Child Growth Standards median [12]. Given that the NDHS does not contain household expenditure information, this paper used the wealth index created in the DHS dataset as an indicator of a child’s socioeconomic status [10]. The wealth index in the NDHS was obtained by applying principal components analysis on household assets [13], accounting for rural-urban differences. Variables used to decompose the change in stunting and wasting were the child’s characteristics (age, sex and religion), mother’s characteristics (education, marital and employment status), and household characteristics (wealth, location, sanitary condition, dependency ratio, and the characteristics of the household head). These variables are correlated with child nutritional status in prior empirical studies. For instance, a child’s age is associated with a child’s nutritional health outcome [14, 15]. Similarly, a child’s sex is associated with nutritional health outcomes, particularly stunting [16–18], while religion and culture are important correlates of a child’s nutritional outcomes [19, 20]. Furthermore, a child with an educated mother is less likely to be stunted or wasted [21–23]. Other mothers’ characteristics like marital status and employment are also associated with children’s nutritional status [24, 25]. Children’s nutritional health outcomes are also affected by household characteristics like wealth and sanitary conditions [15, 21, 23]. In this paper, sanitary conditions were proxied by a household’s water and toilet conditions. Improved water and toilet conditions were created based on the WHO/UNICEF Joint Monitoring Programme for Water Supply, Sanitation and Hygiene classification [26]. A child’s location is also associated with the child’s health status [27–29]. Moreover, the age and sex of the household head are correlated with children’s nutritional outcomes [30, 31], while high dependency ratios are detrimental to children’s health outcomes [32]. Basic descriptive statistics were computed and compared between 2013 and 2018. Percentages were computed for categorical variables while mean values were reported for continuous/count variables. Even though there may be skewness, the mean of continuous/count variables were reported to ease comparison with results in the literature. A concentration curve, suited to assess socioeconomic inequalities [33], depicts the cumulative shares of stunted or wasted children against the cumulative population shares, ranked by socioeconomic status. A line of equality is depicted by a 45-degree line, with a proportional concentration curve coinciding with this line. A pro-rich concentration curve—which depicts a disproportionate concentration of stunting or wasting on the rich—lies below the line of equality. Conversely, a pro-poor curve indicating the disproportionate concentration of stunting or wasting on the poor lies above the line of equality [34]. The concentration index, obtained from the concentration curve, is twice the covariance between the health outcome and the fractional rank in the wealth distribution, divided by the mean of the health outcome indicator [34–37]: where CH refers to the concentration index of the health outcome (H); μH refers to the mean of the health outcome (i.e. stunting or wasting), and r is the fractional rank of the individual/household in the wealth distribution. For continuous outcomes, the index lies in the [-1, +1] interval. A negative (positive) index typically indicates a pro-poor (pro-rich) distribution of stunted or wasted children, while a zero index indicates perfect equality. A zero concentration index may also result from a complex relationship where the concentration curve crosses the line of equality [34]. For categorical variables like indicators of stunting and wasting, the concentration index should be normalised as it may not lie between -1 and +1 [38]. The Erreygers’ normalisation approach was used in this paper [39, 40] to obtain the Erreygers’ normalised concentration index (EC) as follows [41]: where a and b refer to the lower and upper limits of the ordinal health indicator, respectively; and μH and CH remain as earlier defined. Changes in the concentration indices or curves between 2013 and 2018 were assessed using the general characterisation of a pro-poor or a pro-rich change or shift laid out by Ataguba [42]. Let Ct−1 and Ct represent the value of the concentration indices at time t − 1 and t, respectively. A pro-poor change in the concentration indices occurs when Ct−1 > Ct while a pro-rich change corresponds to Ct−1 < Ct. Only two indices are suitable for assessing relative socioeconomic inequalities in health—the slope index of inequality and the concentration index [33]. They are consistent with ranking individuals across socio-economic groupings, sensitive to changes in population distribution across socio-economic groups and consistent in the distribution of stunting and wasting across the distribution of socioeconomic status [33, 43]. The concentration index was used in this paper because it is decomposable to ascertain the factors that significantly explain socioeconomic inequalities in stunting and wasting, which is relevant to policymakers for evidence-based policy interventions. Socioeconomic inequalities in stunting and wasting were decomposed using the Wagstaff et al.’s [44] approach. The relationship between a child’s health outcomes (i.e. stunting or wasting) (H) and associated determinants (z) can be denoted as: where φ and β are parameters, and ε is the error term. Eq (3) was estimated using linear probability models because binary indicators of stunting and wasting were used, with each model appropriately weighted to the population while correcting for heteroscedasticity. For each year, the concentration index in Eq (1) can be decomposed as (the time subscripts are suppressed here to avoid notational clutter): where (βkz-kμH=ηk) denotes the elasticity of stunting or wasting to marginal changes in the k-th explanatory variable, while Ck refers to the concentration index of the k-th explanatory variable. GCε refers to the generalised concentration index of the error term, and (GCεμH) captures the unexplained/residual component. While Eq (4) enables us to decompose observed socioeconomic inequalities in stunting and wasting at each point in time, it does not explain temporal changes in such inequalities. Sequel to decomposing the health concentration index for the relevant health outcome in each year as shown in Eq (4), the change in the concentration index can be decomposed between both periods using the Oaxaca-Blinder decomposition as follows [34]: where t indicates the time period (with the convention that t − 1 and t denote the first (i.e. 2013) and second (i.e. 2018) periods respectively, while other terms in Eq (5) are as earlier defined). Eq (5) decomposes the temporal change in the socioeconomic inequality in stunting or wasting into a part that captures the temporal changes in socioeconomic inequality in the determinants of stunting or wasting: ∑k ηk,t(Ck,t − Ck,t−1), and another component that captures the changes in the elasticities of stunting or wasting with respect to these determinants: ∑k Ck,t−1(ηk,t − ηk,t−1). The final term in Eq (5), Δ(GCε,tμt), measures the temporal difference in the error/residual component [34]. Note that in Eq (5), the temporal changes in the concentration indices of the predictors/determinants are weighted by the elasticities of the second period, while the first-period concentration indices are used to weight the change in the elasticities. The Oaxaca-Blinder decomposition is not unique [34], as an alternative decomposition would be to weight the difference in the concentration indices by the first-period elasticities and the difference in the elasticities by the second-period concentration indices. Given that there are no analytical standard errors for the estimates derived from Eqs (4) and (5), bootstrap routines [45, 46] with 1,000 replications were used, accounting for the full sampling structure of the NDHS. All analyses were implemented using Stata® [47]. The conindex routine was employed in obtaining the Erreygers’ normalised or corrected concentration indices, while accounting for the survey design (clustering, stratification and the resulting survey weights) in each NDHS round [48].
N/A