Abstract
Background Unsafe disposal of children’s stool makes children susceptible to fecal-oral diseases and children remain vulnerable till the stools of all children are disposed of safely. There is a paucity of data on spatial distribution and factors associated with unsafe child stool disposal in Ethiopia. Previous estimates, however, do not include information regarding individual and community-level factors associated with unsafe child stool disposal. Hence, the current study aimed (i) to explore the spatial distribution and (ii) to identify factors associated with unsafe child stool disposal in Ethiopia.
Methods A secondary data analysis was conducted using the recent 2016 Ethiopian demographic and health survey data. A total of 4145 children aged 0–23 months with their mother were included in this analysis. The Getis-Ord spatial statistical tool was used to identify high and low hotspots areas of unsafe child stool disposal. The Bernoulli model was applied using Kilduff SaTScan version 9.6 software to identify significant spatial clusters. A multilevel binary logistic regression model was fitted to identify factors associated with unsafe child stool disposal.
Results Unsafe child stool disposal was spatially clustered in Ethiopia (Moran’s Index = 0.211, p-value < 0.0001), and significant spatial SaTScan clusters of areas with a high rate of unsafe child stool disposal were detected. The most likely primary SaTScan cluster was detected in Tigray, Amhara, Afar (north), and Benishangul-Gumuz (north) regions (LLR: 41.62, p< 0.0001). Unsafe child stool disposal is more prevalent among households that had unimproved toilet facility (AOR = 1.54, 95%CI: 1.17-2.02), and those with high community poorer level (AOR: 1.74, 95%CI: 1.23-2.46). Higher prevalence of unsafe child stool disposal was also found in households with poor wealth quintiles. Children belong to agrarian regions (AOR: 0.62, 95%CI 0.42-0.91), children 6-11 months of age (AOR: 0.66, 95%CI: 0.52-0.83), 12-17 months of age (AOR: 0.68, 95%CI: 0.54-0.86), and 18-23 months of age (AOR: 0.58, 95%CI: 0.45-0.74) had lower odds of unsafe child stool disposal.
Conclusions Unsafe child stool disposal was spatially clustered. Higher odds of unsafe child stool disposal were found in households with high community poverty level, poor, unimproved toilet facility, and with the youngest children. Hence, the health authorities could tailor effective child stool management programs to mitigate the inequalities identified in this study. It is also better to consider child stool management intervention in existing sanitation activities considering the identified factors.
Background
Disposal of child stool in the open field, thrown into the garbage, put into open drains, burying, or left on the ground are considered unsafe [1,2]. Unsafe disposal of children’s stool makes children susceptible to many fecal-oral diseases [1,3–5]. A systematic review by Gil et al. reported that unsafe child stool disposal associated with a 23% (RR: 1.23, 95% CI: 1.15-1.32) increase in the risk of diarrheal diseases in children [6]. Another review of case-control studies showed that the disposal of a child’s feces into a latrine decreases the odds of diarrhea by about 25% in children under five years of age (OR: 0.73, 95% CI: 0.62-0.85) [7]. Furthermore, studies reported the effect of young children’s stool disposal and increased risk of stunting [8,9].
In Ethiopia, unsafe child stool disposal is a huge challenge and Ethiopia ranked number 26 for the percentage of children whose feces are safely disposed of, putting the country among the worst third of 38 African countries with available Multiple Indicator Cluster Survey (MICS) [2]. And diarrhea is one of the major contributors to deaths for under age 5 children in the country [10]. Based on the WHO/CHERG estimates, diarrhea contributes to more than one in every ten (13%) child deaths in Ethiopia [11]. According to a recent pooled data analysis using DHS surveys in Ethiopia, 77 percent of children feces disposed of unsafely; either throw outside the yard or not disposed of [12].
Previous studies identified multiple factors that contribute to the occurrence of unsafe child stool disposal [12–20]. Socio-demographic (age of the child, sex of the child, age of the mother, maternal educational status, place of residence) [12–14,21], socio-economic (household wealth index) [14,17], and access to sanitation facility [12–14,17,19] were associated with unsafe child stool disposal practices. Moreover, given the fact that unsafe child stool disposal is associated with open defecation, the prevalence of unsafe child stool disposal practices in rural areas is higher compared to urban areas [12,22]. In the context of this research, higher rates of unsafe child stool disposal were found in poor, rural households with the youngest children and where other household members defecate in the open [14–22].
So far few studies were conducted on child stool disposal in Ethiopia [12–14], and previous estimates, however, identified the determinants of child stool disposal using a standard logistic regression model that does not include information regarding individual and community-level factors [12–14]. Thus, a multilevel regression model is required, which considers the hierarchal and cluster nature of the Ethiopian Demographic and Health Survey (EDHS) data and enhances the accuracy of estimates. And to date, studies on child stool disposal in Ethiopia have not assessed the spatial distribution of unsafe child stool disposal. Therefore, the current study aimed (i) to explore the spatial distribution and (ii) to identify factors associated with unsafe child stool disposal in Ethiopia using a multilevel regression model.
Methods
Data source and sampling
A secondary data analysis of the Ethiopia Demographic and Health Survey (EDHS) conducted in the year 2016 was used. The EDHS-2016 is the recent survey implemented by the Central Statistical Agency (CSA). The survey was conducted from January 18, 2016, to June 27, 2016, based on a nationally representative sample that provides estimates at the national and regional levels and for urban and rural areas.
The surveys used a stratified two-stage cluster sampling technique. The sampling frame used for the EDHS-2016 is the Ethiopia Population and Housing Census (PHC), which was conducted in 2007 by the Ethiopia Central Statistical Agency. The census frame is a complete list of 84,915 enumeration areas (EAs) created for the 2007 PHC. An EA is a geographic area covering on average 181 households. In the first stage, a total of 645 EAs (202 in urban areas and 443 in rural areas) were selected with probability proportional to EA size and with independent selection in each sampling stratum. In the second stage of selection, a fixed number of 28 households per cluster were selected with an equal probability systematic selection from the newly created household listing [10]. The present study included the youngest child under age 2 living with the mother.
Data source and extraction
The EDHS-2016 data sets were downloaded in STATA format with permission from the Measure DHS website (http://www.dhsprogram.com). Unsafe child stool disposal and its potential predictor variables at the individual and community levels were extracted accordingly.
Study variables
Outcome variable
The outcome variable for this study was child stool disposal, which is dichotomized as unsafe and safe. A child’s stool was considered to be disposed of “safely” when the child used latrine/ toilet or child’s stool was put/rinsed into a toilet/latrine, whereas other methods were considered “unsafe” (i.e. put, rinsed in a drain, ditch, thrown in the garbage, left or buried in the open). The survey collected these data from mothers’ verbal reports on whether the child’s last stools were put in or rinsed into a toilet or latrine, buried, or the child used a toilet or latrine.
Independent variables
The independent variables for this study were classified as individual and community level factors. The individual-level variables of this study were (age of the child, sex of the child, presence of diarrhea in the last two weeks, the source of drinking water, sanitation facilities, mother educational level, mother occupation, and household wealth quintile). The contextual region, place of residence, and community poverty were identified as community-level variables (Table 1). The choice of independent variables was guided by the previous works of literature [12–15]. Community poverty level (proportion of women in the poorest and poorer quintile derived from data on wealth index which is categorized as low and high poverty community). The interest of the current study was not in the regions delineated for administrative purposes, which might not necessarily be related to child stool disposal of the population. Accordingly, in the current study, the regions were categorized into agrarian, pastoralist, and city. The regions of Tigray, Amhara, Oromiya, SNNP, Gambella, and Benshangul Gumuz were recorded as agrarian. The Somali and Afar regions were combined to form the pastoralist region and the city administrations-Addis Ababa, Dire Dawa, and Harar were combined as the city.
Data management and analysis
Descriptive measures were used to illustrate the overall characteristics of the study participants. A sampling weight was done to adjust for the non-proportional allocation of the sample to different regions and the possible differences in response rates. A detailed explanation of the weighting procedure can be found in the EDHS methodology parts [10].
In Ethiopian DHS data, children within a cluster more similar to each other than between clusters. When analyzing such datasets, a multilevel model is generally more appropriate than the standard regression model because it enables one to deal with the hierarchical structure of variables. For this reason, a multilevel model was used to identify factors associated with unsafe child stool disposal. As the response variable was dichotomous (safe, unsafe), multilevel binary logistic regression was fitted. The model goodness of fit was checked using deviance and Akakie Information Criteria (AIC). The model with the lowest deviance and AIC was chosen. The Proportional Change in Variance (PCV) was computed for each model with respect to the empty model to show the power of the factors in the model to explain unsafe child stool disposal. Accordingly, the PCV was calculated by the following formula [PCV = (Ve-Vmi)/Ve], where Ve is variance in unsafe child stool disposal in the empty model and Vmi is variance in successive models. [Median Odds Ratio ], where V is the estimated variance of clusters. The MOR measure is always greater than or equal to 1. If the MOR is 1, there is no variation between clusters. Variables with p-value < 0.25 in the bi-variable analysis were fitted in the multivariable model. Adjusted Odds Ratio (AOR) with a 95% Confidence Interval (CI) and p-value < 0.05 in the multivariable model were used to declare significant association with unsafe child stool disposal. Variance inflator factor (VIF) was employed for checking multicollinearity among the independent variables. STATA version 14 software was used for all statistical analyses.
Spatial autocorrelation analysis
In this study, the spatial analysis was performed using the spatial statistics tool (ArcGIS Version 10.3; Redlands, California, United States). The spatial autocorrelation (Global Moran’s I) statistic measures were used to evaluate whether unsafe child stool disposal was dispersed, clustered, or randomly distributed [23]. Spatial heterogeneity of high /low areas of unsafe child stool disposal was examined using the Getis-Ord Gi* statistics and associated Z-scores. Moreover, the spatial interpolation technique was applied (using the ordinary kriging interpolation technique) to predict the unsampled /unmeasured value from sampled measurements.
Spatial scan statistical analysis
Spatial scan statistical analysis was employed to identify the geographical locations of statistically significant spatial clusters of unsafe child stool disposal in Ethiopia using SaTScan™ version 9.6 software. Unsafe child stool was taken as cases and those with safe child stool disposal as controls to fit the Bernoulli model [24]. The default maximum spatial cluster size of < 50% of the population was used. A Likelihood ratio test statistic was used to determine whether the number of observed unsafe stool disposal cases within the potential cluster was significantly higher than the expected or not. Primary and secondary clusters were identified using p-values and log-likelihood ratio tests.
Results
Socio-demographic characteristics of participants
A total of 4145 children aged 0-23 months with their mother were included in this analysis. Of these, 2164 (52.2%) were female with a mean age of 10.66 months (SD ± 0.11). The majority of 3647 (88.0%) of the respondents were rural residents. About 2500 (60.3%) of the children’s mother had no formal education and about one-fifth were in the poorest wealth quintile (Table 2).
Unsafe child stool disposal status
The prevalence of unsafe child stool disposal was 63.10% (95%CI: 59.5-66.6%). Unsafe child stool disposal is varied across urban-rural areas, age of the child, and household wealth index (Table 3).
Spatial distribution of unsafe child stool disposal in Ethiopia
The analysis of spatial autocorrelation indicated that the spatial distribution of unsafe child stool disposal was clustered in Ethiopia. The Global Moran’s I value 0.211 (p-value < 0.0001) indicated that there was significant clustering of unsafe child stool disposal in Ethiopia (Additional File 1).
Hot-spot areas were found in Tigray (Central, and northeast), Amhara (Central, North, and Southeast), Afar (West, and South), Gambela (West), Oromia (South and East), North and Southeast parts of Somali regions, while cold-spot areas were found in SNNP (North, West, and East), Benishangul-Gumuz (Southwest), Addis Ababa, Harari and Dire Dawa (Figure 1).
Ordinary kriging interpolation analysis was conducted to predict child stool disposal in Ethiopia. High unsafe child stool disposal areas were found in Tigray, Amhara, Afar, Gambela, Southern Somali, and Southeastern parts of Oromia regions. In contrast, low unsafe child stool disposal areas were predicted in SNNP, Southern parts of Benishangul-Gumuz, Northern Somali, Western Oromia, and some parts of Amhara regions (Figure 2).
Spatial scan statistical analysis
Table 4 show significant spatial clusters of unsafe child stool disposal in Ethiopia. Most likely (primary clusters) and secondary clusters of unsafe child stool disposal were identified. A total of 270 significant clusters were identified at which 201 were most likely (primary) and 69 secondary clusters. The primary clusters were located in Tigray, Amhara, and Afar regions. The primary clusters were centered at 13.351814 N, 38.353591 E with 471.07 km radius, a relative risk (RR) of 1.26, and the Log-Likelihood Ratio (LLR) of 41.62, at p< 0.0001. The bright pink colors indicate that the most statistically significant spatial windows contain primary clusters of unsafe child stool disposal in Ethiopia. There was high unsafe child stool disposal within the cluster than outside the cluster (Figure 3).
Measures of variation (random-effects) and model fit statistics
Table 5 shows the measures of variation (random intercept models) and model fit statistics. Model comparison was done using deviance. The comparison was done among model with no independent variables (the null model), model 1 (a model with only individual-level factors), model 2 (a model with only community-level factors), and model 3 (a model with both individual and community level independent variables simultaneously). A model with the lowest deviance (model 3) was selected. According to the multilevel logistic regression model, about 39.61 % of the variance in the odds of unsafe child stool disposal could be attributed to community-level factors. And the between-cluster variability declined over successive models, from 39.61% in the empty model to 29.62% in the combined model. In the final model (model 3), individual and community-level factors accounted for about 35.87% of the variation observed for unsafe child stool disposal.
Factors associated with unsafe child stool disposal
The calculated value intra-cluster correlation (ICC) was 39.61%, which indicated that the assumption of independent observation was violated (Table 5). Thus, we used a multilevel logistic regression model to account for the cluster effect. In the multivariable multilevel binary logistic regression model, the age of the child, wealth index, type of toilet facility, region, and community poverty level were identified as associated factors with unsafe child stool disposal. The age of the child was an important predictor of unsafe child stool disposal. Children in the age group of 6-11 months, 12-17 months, and 18-23 months were about 34% (AOR: 0.66, 95%CI: 0.52-0.83), 32% (AOR: 0.68, 95%CI: 0.54-0.86), and 42% (AOR:0.58, 95%CI: 0.45-0.74) less likely to have unsafe child stool disposal than children younger than 0-5 months, respectively. Wealth index also showed a strong statistical association with unsafe child stool. Children belonging to the poorest wealth quintiles had a four times higher chance of unsafe child stool disposal (AOR: 4.62, 95%CI: 2.98-7.16) than the richest children. Similarly, children belonging to the poorer (AOR = 2.77, 95%CI: 1.82-4.23), middle AOR: 2.13, 95%CI: 1.41-3.22), and richer wealth quintiles (AOR: 1.56, 95%CI: 1.05-2.32) had higher odds of unsafe child stool disposal than the richest children. Children belong to households who had unimproved toilet facilities were about 54% (AOR: 1.54, 95%CI: 1.17-2.02) more likely to had unsafe child stool disposal than children in households with improved toilet facilities. The likelihood of unsafe child stool disposal among children belong to agrarian regions was about 38% (AOR: 0.62, 95%CI 0.42-0.91) lower compared with city dwellers. Unsafe child stool disposal is more prevalent among households that are high community poorer level (AOR: 1.74, 95%CI: 1.23-2.46) than those with younger children live in low community poverty level (Table 6).
Discussion
This study aimed to explore the spatial variations in unsafe child stool disposal in Ethiopia. It also aimed to identify the individual- and community-level factors associated with unsafe child stool disposal practice among women using the demographic health survey (DHS) conducting in Ethiopia. Unsafe child stool disposal was found to be a spatial problem in Ethiopia. It was also found that unsafe child stool disposal was significantly associated with different individual and community level factors. At the individual level, the age of the child, wealth index, and types of toilet facilities were associated with unsafe stool disposal. At the community level, region and community poverty were the identified factors associated with unsafe child stool disposal.
In Ethiopia, the proportion of unsafe child stool disposal was 63.10%; the highest proportion was reported in rural areas (p< 0.001). This is almost consistent with a previous similar study, which found unsafe child feces disposal reported by 67% of households in Ethiopia [14]. However, this figure was lower than unsafe child faces disposal reported from India 79% [3], 72.4% [21], and Bangladesh 80% [20,25], 84% [26]. This discrepancy may due to the operational definition used in classifying unsafe child stool disposal; as the study from Bangladesh includes buried as safe child feces disposal [20]. The other possible reasons for this disparity may be due to study participants; previous studies include under-5 children and we included the youngest child under age two. On the other hand, a much lower prevalence of unsafe child stool disposal was reported from Malawi, 14.5% [15]. This could be attributed to continued efforts in raising awareness, and ODF campaigns programs in Malawi [15].
In the global spatial autocorrelation analysis of this study, a clustering pattern of unsafe child stool disposal across the study areas was observed (Global Moran’s I = 0.211, p-value < 0.0001). This indicates that unsafe child stool disposal in Ethiopia was aggregated in specific areas. Accordingly, the hot-spot areas were found in Tigray (Central, and northeast), Afar (West, and South), Amhara (Central, North, and Southeast), Gambela (West), Oromia (East), North and some parts of Somali regions. The possible explanation for geographic variation in the prevalence of unsafe child stool disposal might be due to high open defecation practice in these identified hot spot areas. According to EDHS-2016 report, one in three households in Ethiopia have no toilet facility (39% in rural areas and 7% in urban areas) and open defecation was practices in 32.9% of the households (37.7% in rural areas and 6.8% in urban areas). Closer looks in these hot spot areas showed that unsafe child stool disposal is relatively aggregated in rural areas. Consistent with this affirmation, unsafe child feces disposal is more prevalent among households that defecate in the open and those in rural areas; over three fourth of the rural households in Ethiopia (81.2%) had unsafe child feces disposal while that is true only for (45.8%) of the urban households [12]. In the community-level factors (model 2), our finding also suggested that the odds of unsafe child stool disposal were two times higher among children residing in rural areas. Additionally, the high proportion of unsafe child stool disposal in this area might be due to disparity in access to improved sanitation facilities.
Consistent with previous studies in Ethiopia [12], Malawi [15], and in rural Bangladesh [20,27], women with older children were less likely to have unsafe child stool dispose of compared with those with younger children. This could be explained by the widespread belief and perception that the stool of young children is considered harmless in many communities, including Ethiopia [2,17]. Also, this strong association can be satisfactorily explained by the fact that a shift in safe disposal practices is usually seen as children grow: children being more likely to use potties and a toilet themselves as they get older and their behaviors tend to change as they get older [2,15,28].
In this study, the wealth index was a strong predictor of statistical association with unsafe child stool. Children belonging to the poorest and poorer wealth quintiles had a higher chance of unsafe child stool disposal than the children in households with the richest wealth quintiles. This finding was consistent with other related studies [12–15,19]. This association can be easily explained by the fact that households with a high wealth index are more likely aware of the negative effects of unsafe child stool disposal, and have access to improved sanitation, and therefore practice safe child stool disposal [15]. In connection, there is also evidence in the current study; unsafe child feces disposal is more prevalent among households with high community poverty levels.
In this study, children belong to households who had unimproved toilet facilities had a statistically significant association with high odds of unsafe child stool disposal, which is in line with the findings of some previous studies [12,14,29]. Yet reports also showed that among households with improved toilets or latrines, some unsafe child feces disposal behavior was reported [22].
At the community level, children belong to agrarian regions (like SNNP and Beneshangul Gumuz regions) were less likely to have unsafe child stool disposal than city dwellers. This finding highlighted the need for strong sanitation programs to strengthen their efforts in the city administration in Ethiopia. So far, the largest Community-Led Total Sanitation and Hygiene (CLTSH) campaign and other related efforts to end open defecation have mainly targeted rural communities in Ethiopia, with only a limited focus on the management of child stool among city dwellers [30,31]. As a result, child feces management should be promoted among city dwellers in the country.
Conclusion
This study showed that unsafe child stool disposal had spatial variability across survey clusters and regions; it was higher in the northern part of the country. Both the individual-level characteristics (child’s age, wealth index, types of toilet facility) and community-level characteristics (region and community poverty) were statistically significant predictors of unsafe child stool disposal. Hence, the health authorities could tailor effective child stool management programs to mitigate the inequalities identified in this study. It is also better to consider child stool management intervention in existing sanitation activities.
Limitations
Despite we used nationally representative data that can enhance the generalisability of the findings and a multilevel logistic regression model, a model that accounts for the correlated nature of EDHS data. The present study has several limitations. First, due to the secondary nature of the data, the present study was limited by unmeasured confounders such as mother knowledge towards child stool disposal and other community-level factors such as social and cultural norms towards child feces management. Second, self-reported practices can be subject to bias that might underestimate true levels by underreporting socially undesirable behaviors. Additionally, in EDHS self-reported child stool disposal practices have not been validated with objective measurements such as spot check observations. Third, the cross-sectional nature of the survey does not allow the cause-and-effect relationship between independent variables and unsafe child stool disposal.
Data Availability
The data we used which is the 2016 Ethiopian Demographic and Health Survey were obtained from the DHS program (www.dhsprogram.com) but the Dataset Terms of Use do not permit us to distribute this data as per data access instructions (http://dhsprogram.com/data/Access-Instructions.cfm). To get access to the dataset you must first be a registered user of the website (www.dhsprogram.com) and download the 2016 Ethiopian Demographic and Health Survey.
Data availability statement
The data we used which is the ‘2016 Ethiopian Demographic and Health Survey’ were obtained from the DHS program (www.dhsprogram.com)but the ‘Dataset Terms of Use’ do not permit us to distribute this data as per data access instructions (http://dhsprogram.com/data/Access-Instructions.cfm). To get access to the dataset you must first be a registered user of the website (www.dhsprogram.com) and download the 2016 Ethiopian Demographic and Health Survey.
Funding
No fund was received for the present study.
Conflicts of interest
The authors declare that there are no competing interests.
Ethics approval
DHS Programme granted permission to download and use the data for this study after being registered and submitting a request with briefly stated objectives of the study. The Institution Review Board approved procedures for DHS public-use data sets do not in any way allow respondents, households, or sample communities to be identified. There are no names of individuals or household addresses in the data files. The geographic identifiers only go down to the regional level (where regions are typically very large geographical areas encompassing several states/provinces). Each EA (primary sampling unit) has a number in the data file, but their numbers do not have any labels to indicate their names or locations. The detail of the ethical issues has been published in the 2016 EDHS final report, which can be accessed at: http://www.dhsprogram.com/publications [10].
Supporting information
S1_Figure 1: The global spatial autocorrelation based on feature locations and attribute values of unsafe child stool disposal in Ethiopia, EDHS 2016