Determinants of under-five mortality in the high mortality regions of Ethiopia: mixed-effect logistic regression analysis

Abstract

Background

Even though the global under-five mortality rate substantially decreased over time, Sub-Saharan African (SSA) countries including Ethiopia continue to share the huge burden of under-five mortality. Ethiopia showed a substantial reduction in under-five mortality over time but the rate of reduction has varied across regions. Therefore, this study aimed to investigate determinants of under-five mortality in the high mortality regions of Ethiopia.

Methods

A secondary data analysis was done based on the 2016 Ethiopian Demographic and Health Survey (EDHS) data. A total weighted sample of 3446 live births were included for this study. For the determinants of under-five mortality, mixed-effect logistic regression was fitted. The Intra-Class Correlation Coefficient (ICC), and Median Odds Ratio (MOR) were done to assess the presence of a significant clustering effect. The standard binary logistic regression and the mixed-effect logistic regression model were fitted and deviance (−2LL) was used for model comparison as the models were nested models. Variables with a p-value less than 0.2 in the bi-variable mixed-effect binary logistic regression analysis were considered for the multivariable analysis. In the multivariable mixed-effect logistic regression analysis, the Adjusted Odds Ratio (AOR) with the 95% Confidence Interval (CI) were reported to declare the statistical significance and strength of association of under-five mortality and the determinant factors.

Results

Overall, the under-five mortality rate in the high mortality regions of Ethiopia was 74 per 1000 live births and it was highest among twin births (262 per 1000 live births). In the multivariable mixed-effect logistic regression analysis, being having 6 and above births (AOR = 3.66, 95% CI: 1.55, 8.67), preceding birth interval of 2–3 years (AOR = 0.57, 95% CI: 0.41, 0.81) and above 3 years (AOR = 0.35, 95% CI: 0.22, 0.55), being twin (AOR = 5.12, 95% CI: 2.28, 11.46), and being having antenatal care (ANC) visit during pregnancy (AOR = 0.27, 95% CI: 0.16, 0.45) were significant determinants of under-five mortality.

Conclusion

In this study, under-five mortality rate was highest in high mortality regions of Ethiopia. Parity, ANC visit, preceding birth interval, and multiple births were significant predictors of under-five mortality. Therefore, public health interventions that increase maternal health service utilization such as ANC and family planning service utilization to increase birth interval are needed to reduce under-five mortality among these regions of Ethiopia.

Background

Under-five mortality is the death of a child before reaching the age of 5 years [1]. The first 5 years of life are the most crucial years for the physical and intellectual development of children [1].

Even though there is a worldwide decrement in under-five mortality, from 5.9 million deaths in 2015 to 5.3 million in 2018, still there is a high mortality rate in African countries (81 per 1000 live births) including Ethiopia, which is around 7 times higher than in the European region [1, 2]. In Africa, under-five mortality contributes for 14% of the global burden of child mortality and in sub-Saharan Africa, it accounts for nearly 50% of child mortality, while it accounts for just 11% of the world’s population [3]. Sub-Saharan Africa countries had the highest under-five mortality with 1 child in every 13 live birth dying before 5 years of age, which is 15 times higher than in high-income countries and half of all these deaths in 2018 occurred in five countries; India, Nigeria, Pakistan, Ethiopia, and the Democratic Republic of the Congo [3, 4].

Under-five child mortality in Ethiopia is one of the biggest and most difficult problems that should be given emphasis [4]. Millennium development goal (MDG) 4 was aimed at reducing under-five mortality by two-thirds and this was unfinished agenda for most of the middle and lower-income countries [5]. With the current rates of child mortality in Ethiopia achieving the sustainable development goal (SDG), which is reducing under-five mortality to 55 deaths per 1000 live births by 2030 is difficult and studies revealed that there is a significant regional variation in under-five mortality [4].

That is, in some regions of Ethiopia, under-five mortality is still devastating and according to the 2016 Ethiopian Demographic and Health Survey (EDHS) report, under-five mortality per 1000 live birth in afar, Benishangul-Gumuz, and Somali were 125, 98, and 94, respectively [6].

Evidence from different literature showed that low level of maternal education, unsafe drinking water and sanitation, low family income, birth interval, short breastfeeding time, and place of delivery are among the determinants of under-five mortality [4, 7,8,9,10,11,12,13,14,15,16,17,18,19].

Even though Ethiopia is working on the reduction of under-five mortality, through improving coverage, quality and use of skilled care, essential newborn care, and management of preterm and low birth weight, still it is very high in the regions of Afar, Somalia, and Benishangul-Gumuz [5]. Therefore, this study aimed to explore the determinants of under-five mortality in these high-risk regions of Ethiopia. This study might give insight for health professionals, policymakers, and the community in general for appropriate intervention to reduce under-five mortality.

Methods

Data source

Secondary data analysis was conducted based on the EDHS 2016, which is the fourth survey conducted nationally from January 18 to June 27, 2016. The EDHS employed a two-stage stratified sampling technique to select respondents. In the first stage, 645 enumeration areas (EAs) (202 in urban areas and 443 in rural areas) were selected. In the second stage, a fixed number of 28 households per cluster was selected, with an equal probability systematic selection from the newly created household list. The detailed sampling procedure was presented in the full EDHS 2016 report [6].

Afar, Somali, and Benishangul-Gumuz regions were selected since based on the 2016 EDHS report, these regions had the highest rates of under-five mortality. A total weighted sample of 3446 women with a live birth in the 5 years preceding the survey were included for this study. For women who had two or more live births in the 5 years preceding the survey, the last birth was taken (Fig. 1).

Fig. 1

Flow diagram showing the sampling

Study variables

The outcome variable was the death of under-five children, which is a dichotomous variable coded as “0” if the child is alive and “1” if the child is died.

The independent variables were: proximal factors (sex of a child and multiple births), intermediate factors (mother’s level of education, birth order, preceding birth interval, parity, ANC visit, respondent current age, and place of delivery), and distant factors (household wealth index, source of drinking water, sex of household head, place of residency, and presence of toilet facility) (Table 1).

Table 1 Description and measurement of independent variables

Data management and analysis

Stata 14 software was used for extracting, recoding, and analysis of data. The data were weighted using sampling weight, primary sampling unit, and strata before any statistical analysis to restore the representativeness of the data and to get a reliable estimate and standard error. Descriptive statistics were done using frequencies and percentages. The EDHS data has a hierarchical structure; this violates the independent assumptions of a standard logistic regression model. Therefore, to assess the determinant factors, the mixed-effects logistic regression model with the Generalized Linear Mixed Model (GLMM) framework was used.

$$ \log\ \left[\uppi \mathrm{ij}/\left(1-\uppi \mathrm{ij}\right)\right]=\upbeta 0+\upbeta 1\mathrm{X}2\mathrm{ij}+\upbeta 2\mathrm{X}3\mathrm{ij}+\upbeta 3\mathrm{X}2\mathrm{ij}+\cdots +\upbeta \mathrm{nXnij}+\cup 0\mathrm{i}+\upepsilon \mathrm{ij} $$

πij = is probability of under-five mortality

1-πij = is probability of surviving in the first 5 years of birth

β0 = is log odds of the fixed intercept

X1ij… Xnij = are independent variables of individuals and community’s level

β1… βn are effect sizes of individual levels and community-level coefficients

uOj are random errors at cluster level

ϵij random error at individual level estimated by π²/3

Multicollinearity between independent variables was assessed using Variance Inflation Factor (VIF) and each of the independent variables and the mean VIF was less than 5. Intraclass Correlation Coefficient (ICC) and Median Odds Ratio (MOR) were checked to assess whether there was clustering and deviance was used for model comparison. Both bivariable and multivariable mixed-effect logistic regression were done. At the bivariable analysis variables with a p-value < of 0.2 were considered for multivariable analysis and variables with a P-value of < 0.05 in the multivariable analysis were considered as determinants of under-five mortality.

Results

Sociodemographic characteristics of study participants

A total weighted sample of 3446 was included in this study. A majority (80.9%) of the mothers were without formal education. Approximately, 87% of study participants were from rural areas and 43.67% were from the Somali region. Most (85.95%) of respondents were followers of the Muslim religion. About 19.38% of respondents were from households with safe drinking water supplies and the majority (58.5%) of participants were from households with no toilet facility. Concerning wealth status, 76% of respondents were from households with a poor wealth status. Most (82.3%) of mothers gave their birth at home and 70.3% of mothers did not have ANC visit during their last pregnancy. Regarding birth order, 30% of mothers had a birth order greater than six and 41% of mothers gave birth with a preceding birth interval less than 2 years. Moreover, 52.3% of children were males and the great majority (98.8%) were singletons (Table 2).

Table 2 Descriptive characteristics of the respondents in Ethiopia, 2016 (N = 3446)

Prevalence of under-five mortality in the high-risks region of Ethiopia

Two hundred fifty-seven children (7.46, 95%CI; 6.6, 8.3%) died before 5 years of age from total live birth with a mortality rate of 74 per 1000 live births. Among mothers with no ANC visit during their last pregnancy, there is higher under-five mortality (95 deaths per 1000 live births) relative to mothers with ANC visits (26 death per 1000 live birth). In multiple births, the rate of under-five mortality was higher (27 deaths per 1000 live births) than in singletons (8 deaths per 1000 live births). The under-five mortality rate is also differed with wealth status, falling from 80 deaths to 53 deaths per 1000 live births in poor vs rich households. In rural and urban areas, the under-five mortality rate was 78 and 47 deaths per 1000 live births, respectively. In households with untreated water sources, the under-five mortality rate was higher (78 deaths per 1000 live births) than in those with a safe water source (60 deaths per 1000). Among households with no standard toilet facility, under-five mortality was higher (79 deaths per 1000 live births) compared to those with a standard toilet facility (67 deaths per 1000 live births) (Table 3).

Table 3 The prevalence of under-five mortality in Ethiopia, 2016 (N = 3446)

Determinants of under-five mortality in the high-risks region of Ethiopia

Model comparison/random effect analysis

Log-likelihood Ratio (LLR) test and deviance were checked, and the mixed effect logistic regression model was chosen since it had the smallest value of deviance. Furthermore, the ICC value was 0.14 and this revealed about 14% of the variability of under-five mortality was attributed due to the difference between communities/clusters. The median odds ratio was 1.49, indicating that if we randomly select one child from a higher-risk cluster (cluster with higher under-five mortality) the odds of under-five mortality is 1.49 times compared to those from a lower-risk cluster. Moreover, the Log-likelihood ratio test was significant (X² value = 4.63, P value = 0.0157). These all informed us to choose a mixed-effect logistic regression model (GLMM) over the basic model (Table 4).

Table 4 Model comparison and fitness for the assessment of determinants of under-five mortality

Fixed effects analysis

Determinants of under-five mortality in the high-risk regions of Ethiopia were analyzed using a mixed-effect logistic regression model. Variables with a p-value < 0.2 in the bivariable analysis were eligible for the multivariable analysis. In the multivariable analysis, parity, preceding birth interval, birth type, and ANC visit were significant determinants of under-five mortality. The odds of under-five mortality was 3.66 (AOR = 3.66: 95%CI; 1.55, 8.68) times higher among mothers with parity of greater than six compared with those mothers having children of three and less in the family. Mothers with preceding birth intervals of two to 3 years and above 3 years had 43% (AOR = 0.57: 95%CI; 0.41, 0.81) and 65% (AOR = 0.35: 95% CI; 0.22, 0.55) lower odds of under-five mortality as compared to mothers with preceding birth interval of fewer than 2 years. Being multiple births had 5.1 (AOR = 5.1: 95%CI; 2.28, 11.46) times higher odds of under-five mortality compared to a single birth. Moreover, there was a 63% (AOR = 0.27 95%CI; 0.16, 0.45) lower odds of under-five mortality among mothers with ANC visit during their last pregnancy compared to their counterparts (Table 5).

Table 5 Multivariable mixed effect analysis for assessing the determinants of under-five mortality in high-risk regions of Ethiopia, 2016

Discussion

The SDG is aimed at reducing preventable deaths of newborns and children under the age of five [9, 16]. In this study, 7.46% (95%CI; 6.6, 8.3%) of under-five children had died before celebrating their fifth birthday with a mortality rate of 74 per 1000 live birth. This finding was relatively higher than another study done in Ethiopia [10] and lower than other studies conducted in Ethiopia [2, 11]. The disparity may be attributed to variations in sample size and study setting, as the sample size in our study was relatively small compared to others. Besides, the discrepancy or the higher prevalence found in this study might be due to the poor socioeconomic status of the participants in these high-risk regions of Ethiopia [10].

The multivariable mixed-effects analysis showed that type of birth, ANC visit, preceding birth interval, and parity had significant associations with under-five mortality.

The odds of under-five mortality among multiple births were higher as compared to singleton births and this is consistent with different studies conducted in Ethiopia [2, 3, 11, 12]. This might be because multiple births can lead to growth retardation and prematurity, which are the main risk factor for under-five mortality [13]. In addition, being twins or triplets might result in under-nutrition because of inadequate breast milk as well as infections due to inappropriate formula feeding and feeding of cow’s milk.

In the present study, under-five mortality increases as parity increases. This is consistent with a study carried out in Ethiopia [2] but in contrast to other studies [7, 14]. This might be because mothers with higher parity might be busy caring for many children in the family and may give limited attention to the health of their last child. Besides, mothers with higher parity had enough number of children so they might be careless in seeking appropriate care from the health facility for this child [17]. Additionally, in this study, there were higher odds of under-five mortality among children born from mothers with a previous birth interval of fewer than 2 years. This finding is in line with other studies [2,3,4]. This might be as a result of short birth interval might lead to premature birth and limited family resources because of the increased number of children in the family [11, 16, 19]. Also, mothers with short preceding birth interval might be exposed to malnutrition and this, in turn, affect her child in getting appropriate nutrients from the mother and might expose the child to malnutrition and infections, which finally ends up with death.

Moreover, children of women who had ANC visits during pregnancy had a lower risk of mortality relative to those of mothers with no ANC visits. This finding is consistent with the study conducted in Ethiopia [4] but in contrast with other studies [3, 15]. This might be because mothers with appropriate utilization of maternal health services such as ANC visits might induce appropriate information for the mother regarding the health of the newborn and appropriate feeding practices to prevent infections [18]. In addition, mothers with ANC visits might have a probability of taking postnatal care that is crucial for searching for childhood health problems and taking appropriate intervention timely.

Strength and limitation of the study

This study was based on nationally representative data with a large sample size. Besides, it was based on an appropriate statistical approach (mixed-effect analysis) to accommodate the community or cluster level variability of under-five mortality. Moreover, since it is based on the national survey data the study has the potential to give insight for policy-makers and program planners to design appropriate intervention strategies both at national and regional levels. However, this study had limitations in that the EDHS is mostly based on respondents’ self-report and might have the possibility of recall bias. Besides, since the study was based on the information available in the EDHS data, there might be residual confounding factors.

Conclusion

In the high-risk region of Ethiopia the prevalence of under-five mortality is high. Under-five mortality in these regions was significantly associated with preceding birth interval, type of birth, parity, and ANC visits. Thus, special emphasis should be given to those mothers with no ANC visit, preceding birth interval of fewer than 2 years, mothers with higher parity, and multiple births.