Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Geographic variation and factors associated with under-five mortality in Ethiopia. A spatial and multilevel analysis of Ethiopian mini demographic and health survey 2019

  • Zemenu Tadesse Tessema ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    zemenut1979@gmail.com

    Affiliations Department of Epidemiology and Preventive Medicine, School of Public Health and Preventive Medicine, Monash University, Melbourne, Australia, Department of Epidemiology and Biostatistics, Institute of Public Health, College of Medicine and Health Sciences, University of Gondar, Gondar, Ethiopia

  • Tsion Mulat Tebeje,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Epidemiology and Biostatistics Unit, School of Public Health, Dilla University, Dilla, Ethiopia

  • Lewi Goytom Gebrehewet

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Epidemiology and Biostatistics, Institute of Public Health, College of Medicine and Health Sciences, University of Gondar, Gondar, Ethiopia

Abstract

Background

The distribution of under-five mortality (U5M) worldwide is uneven and the burden is higher in Sub-Saharan African countries, which account for more than 53% of the global under-five mortality. In Ethiopia, though U5M decreased substantially between 1990 and 2019, it remains excessively high and unevenly distributed. Therefore, this study aimed to assess geographic variation and factors associated with under-five mortality (U5M) in Ethiopia.

Methods

We sourced data from the most recent nationally representative 2019 Ethiopian Mini-Demographic and Health Survey for this study. A sample size of 5,695 total births was considered. Descriptive, analytical analysis and spatial analysis were conducted using STATA version 16. Both multilevel and spatial analyses were employed to ascertain the factors associated with U5M in Ethiopia.

Results

The U5M was 5.9% with a 95% CI 5.4% to 6.6%. Based on the multivariable multilevel logistic regression model results, the following characteristics were associated with under-five mortality: family size (AOR = 0.92, 95% CI: 0.84,0.99), number of under-five children in the family (AOR = 0.17, 95% CI: 0.14, 0.21), multiple birth (AOR = 14.4, 95% CI: 8.5, 24.3), children who were breastfed for less than 6 months (AOR = 5.04, 95% CI: 3.81, 6.67), people whose main roof is palm (AOR = 0.57, 95% CI: 0.34, 0.96), under-five children who are the sixth or more child to be born (AOR = 2.46, 95% CI: 1.49, 4.06), institutional delivery (AOR = 0.57, 95% CI: 0.41, 0.81), resident of Somali and Afar region (AOR = 3.46, 95% CI: 1.58, 7.55) and (AOR = 2.54, 95% CI: 1.10, 5.85), respectively. Spatial analysis revealed that hot spot areas of under-five mortality were located in the Dire Dawa and Somali regions.

Conclusion

Under-five mortality in Ethiopia is high and unacceptable when compared to the 2030 sustainable development target, which aims for 25 per 1000 live births. Breastfeeding for less than 6 months, twin births, institutional delivery and high-risk areas of under-five mortality (Somali and Dire Dawa) are modifiable risk factors. Therefore, maternal and community education on the advantages of breastfeeding and institutional delivery is highly recommended. Women who deliver twins should be given special attention. An effective strategy should be designed for intervention in under-five mortality hot spot areas such as Somali and Dire Dawa.

Background

Under-five mortality is defined as deaths reported at ages 0 to 59 months and includes neonatal, post-neonatal, and child deaths [1]. It is a significant indicator of the socioeconomic, health and environmental conditions, and of national development and health equity and access [2].

Globally, in 1990, the number of under-five deaths was 12,494,000 (93 deaths per 1000 live births). Following an average annual decline of 3.1%, by 2019 it was significantly reduced to 5,189,000 deaths (38 deaths per 1000 live births) representing a 58% reduction [3]. The 2030 sustainable development goals (SDG) aims to end the preventable deaths of newborns and children under 5 years of age, with all countries aiming to reduce under-5 mortality to at least 25 per 1000 live births or lower [4].

The distribution of U5MR worldwide is uneven and the global burden of U5MR is concentrated in two regions, Sub-Saharan Africa and Central and Southern Asia, which accounts for more than 80% of U5MR [5]. The burden is even higher in Sub-Saharan African countries with 74 deaths per 1000 live births, which accounts for more than 53% of the global under-five mortality [6].

In Ethiopia, U5MR has substantially decreased from 200 per 1000 live births in 1990 to 49 per 1000 live births in 2020 but is not close to Sustainable Development Goal [3]. Additionally, different studies conducted in Ethiopia revealed several socioeconomic, demographic, and geographic or spatial variation. For example, Afar, Somali, and Benishangul Gumuz are the three high-risk regions of Ethiopia for U5MR [7].

Different studies examining U5M previously showed that source of water [8], multiple births [6, 7], region [9, 10], age and sex of household head [1], maternal education status [11, 12], place of delivery [9, 10, 13], maternal age [1416], sex of the child [8, 13, 16], residence [17, 18], family size [9, 15], number of children under five [18, 19], birth order [9, 15], and breastfeeding status [18, 20] have a significant association with under-five mortality.

Although U5M has been declining it remains high and takes the lives of many children, no research on U5M has been made until recently, with the Ethiopian mini–Demographic Health Survey (EMDHS) 2019. Research is required to inform policies and interventions for the prevention of U5M with the efficient allocation of scarce resources to priority areas based. This study aimed to assess the geographic variation and factors associated with under-five mortality in Ethiopia. The results of our study will assist in achieving the SDG, by providing an understanding of the current burden of U5M in Ethiopia, and its determinants, in addition to identifying areas and factors with the highest burden of U5M in Ethiopia.

Methods

Data source

Our source was the 2019 EMDHS, the second mini demographic health survey (DHS) conducted in Ethiopia (a land-locked country located in the Horn of Africa that lies between the 30N and 150N Latitude or 330E and 480E Longitude) [21]. Data collection was conducted from March 21, 2019, to June 28, 2019, the nine regions (Tigray, Afar, Amhara, Oromia, Somali, Benishangul Gumuz, Southern nation nationalities and People region (SNNPR), Harari, and Gambella) and two administrative cities (Addis Ababa and Dire Dawa). The study design was a population-based cross-sectional study. A frame of all census Enumeration areas (EAs) was used as a sampling frame for the 2019 EMDHS. 149,093 EAs were created which cover an average of 131 Households (HHS). A two-stage stratified cluster sampling technique was employed and each region was stratified into urban and rural areas, yielding 21 sampling strata were selected independently in each stratum. In the first stage, 305 clusters (93 urban and 212 rural) were selected with probability proportional to EAs size and with independent selection in each sampling stratum. In the second stage, a fixed number of 30 households per cluster was selected. Finally, women aged 15–49 in 9,150 (2,790 urban and 6,360 rural) households from 305 clusters were selected. The whole procedure of sampling is found in the full 2019 EMDHS report [21].

Study variables

The outcome variable was under-five mortality status, which was categorized as (child alive: Yes = 0 and No = 1). The age was recorded in months. The community-level predictors, were place of residence, region, community place of delivery, community wealth, community media exposure, and community toilet facility. The individual-level predictors were further categorized as socio-demographic and economic factors like the educational level of the mother, sex of the household head, age of household head, number of household members, number of children under the age of five, marital status of the mother, source of water, time to get water, type of toilet facility, household electricity, types of cooking fuel, main floor material, main wall material, and main roof material and maternal and child factors (such as maternal age at first birth, sex of the child), utilization of contraception, order of birth, mode of delivery, duration of breastfeeding and multiple births (Fig 1).

thumbnail
Fig 1. Conceptual framework of under-five mortality in Ethiopia 2019.

https://doi.org/10.1371/journal.pone.0275586.g001

Spatial analysis

The data for spatial analysis was cleaned and merged using STATA version 16 and Microsoft Excel. ArcGIS version 10.8 and saTScan version 9.7 were used for the spatial analysis.

Spatial autocorrelation

Spatial autocorrelation (Global Moran’s I) analysis was conducted to examine whether under-five mortality was dispersed (Moran’s I value closer to -1), clustered (Moran’s I value closer to 1), or randomly distributed (Moran’s I value of 0) in Ethiopia [22].

Spatial interpolation

The under-five mortality was known in enumerated areas, while in areas that were not selected, the under-five mortality rates were predicted. Spatial interpolation was applied using the geostatistical ordinary Kriging spatial interpolation technique to predict under-five mortality from existing sample data points to un-sampled areas [23].

Spatial scan statistics

The scan analysis was performed using SaTscan, based on the Bernoulli test for cases (child is not alive) and controls (child is alive). The upper limit used was the default maximum spatial cluster size of less than 50% of the population, allowing both small and large clusters to be detected, while clusters that contained more than the maximum limit with the circular shape of the window were avoided [24]. Most likely clusters were identified using p-values and likelihood ratio tests, which is the ratio of the likelihood of the alternative hypothesis (higher activity level inside the window) over the likelihood of the null hypothesis (same activity level inside and outside).

Data management and analysis

We used STATA version 16 and R statistical software version 4.0.5 to analyze the data. A total of 31 variables were retained for the analysis. Residents who were not De jure were excluded which affected under-five mortality, as they could not respond to most of the socio-demographic and economic characteristics even though they could answer the maternal and child characteristics. This exclusion changed our sample size from 5,753 to 5,695.

The outcome variable was re-coded to (child alive: Yes = 1 and No = 0). Four community variables were generated by taking the individual variable, calculating their proportion, and dichotomizing them based on their mean or median according to their distribution.

In the end, we had 29 predictor variables, of which 6 were community-level predictor variables and 23 were individual-level predictor variables.

Based on EMDHS, respondents in the same cluster showed similar outcomes or functions at the same level and the data has a hierarchical structure. This made binary logistic regression not the most appropriate as it violates the assumption of independence of the residuals. Instead, a model that considers clustering effect should be used [7, 10]. Multilevel logistic modeling separates the within-cluster effects from the between-cluster effects [25]. Therefore, to assess the predictors associated with U5M, a non-weighted multilevel logistic regression model was used.

Bivariable multilevel logistic regression was used to screen each predictor variable for a p-value less than 0.2. Significant variables were included in multivariable multilevel logistic models. Twenty predictor variables (3 of which were community-level predictor variables) were included in the multivariable analysis. In the multivariable analysis, a p-value less than 0.05 was considered a factor associated with U5M.

The first model fitted was the null model (intercept model), which contained the outcome variable only (under-five mortality status) with the cluster number. The intra-cluster correlation (ICC) was used to assess whether there was a random effect. An ICC of 0.130 which meant there was a minimum of 13% under-five mortality was explained by between-cluster differences. We found that 87% of under-five mortality was explained by within-cluster differences, which was not negligible. The second model was fitted using the outcome, the cluster and the individual-level predictor variables only. The probability of U5M was predicted as a function of individual-level predictors.

For the third model the outcome variable, the cluster number, and the community-level predictor variables were accounted for. Then the final model was fitted by taking both the individual-level and the community-level predictor variables into account.

The models were compared by using a log-likelihood statistic, where the best model was selected based on smallest deviance.

Ethical consideration

Permission for data access was obtained from a major demographic and health survey through an online request from http://www.dhsprogram.com. The data used for this study were publicly available with no personal identifiers.

Results

Socio-demographic background and maternal and child health factor

A total of 5,695 live births were included in this study, of which 2,942 (51.66%) were male, 5,528 (97.07%) were single children, 2,937 (51.57%) were delivered at home and 6.04% of the children were born via cesarean section (Table 1).

thumbnail
Table 1. Sociodemographic, maternal, and child health factors.

https://doi.org/10.1371/journal.pone.0275586.t001

Out of the total children’s mothers, most were married (5,346, 93.87%) and 3,135 (55.05%) had no education, only 291 (5.11%) attained higher education, and 51.62% were of poor economic status. 64.27% had access to improved water; 30.99% travel more than 30 minutes, and 50.85% travel less than 30 min to get water because they do not have water on-premises. The majority of household heads were male, 4,561 (80.09%), and the most common age ranges of the household head was 35–50 years (Table 1).

At the community level, the majority were rural residents 4,389 (77.07%), and 3,457 (60.70%) in the community practiced home delivery. The number of communities that used unimproved toilets was 3,175 (55.75%). And 3,391 residents of the population (59.54%) were not exposed to the media (Table 2).

Prevalence of under-five mortality

The under-five mortality was 5.9% (95% CI 5.4%, 6.6%) and it varied across regions ranging from 2.13% in Addis Ababa to 9.82% in the Somali region. Out of the total, there were 918 (15.79%) children who died having been breastfed for less than 6 months (Table 3).

Spatial autocorrelation

The global Moran’s index was 0.09, z-score 3.5, and p-value <0.001, suggesting that the spatial distribution under-five mortality was not at random in the 2019 EMDHS (Fig 2).

Hotspot analysis

A high proportion of U5M was observed around the regions of Somali and Dire Dawa (Fig 3) [26].

Spatial interpolation

The spatial interpolation or prediction analysis result showed that the under-five mortality ranged from 0% to 50%. The areas predicted to have the highest risk were located in Somali, Dire Dawa, Benishangul Gumuz, Gambella, and SNNPR (Fig 4) [26].

thumbnail
Fig 4. Empirical bayesian kriging prediction map of under-five mortality in Ethiopia 2019.

(Source: https://open.africa/dataset/africa-shapefiles).

https://doi.org/10.1371/journal.pone.0275586.g004

Spatial scan statistics analysis

As observed in Fig 5, significant clusters were identified in the Somali region. Located within 7.384788 N, 45.908936 E with 612.75 km radius. Children living in the primary cluster were 62% more likely to die before they celebrate their fifth birthday when compared to those who lived outside the window (relative risk (1.62), likelihood ratio (7.36), and p-value (0.03) (Fig 5) [26].

Factors associated with under-five mortality

After filtering, all variables in bivariable multilevel logistic regression, 4 models were fitted (null model, individual-level model, community-level model, and both individual- and community-level model). The ICC was 0.13 with 95% CI: 0.08 to 0.2, which is not negligible. The fourth model was the best-fitting model with the lowest deviance. When the number of household members increased by one, the odds of under-five mortality decreased by 8.5% (AOR = 0.92, 95% CI: 0.84, 0.99). The odds of under-five death decreased by 82.8% (AOR = 0.17, 95% CI: 0.14, 0.21) as the number of under-five children in the family increased by one. Children who were part of multiple births had 14 times higher odds of under-five mortality (AOR = 14.4, 95% CI: 8.5, 24.3) than singletons. Children who were breastfed for less than 6 months had 5 times higher odds of experiencing under-five mortality (AOR = 5.04, 95% CI: 3.81, 6.67) than those who were breastfed. People whose main roof was palm/bamboo had 43% (AOR = 0.57, 95% CI: 0.34, 0.96) decreased odds of having under-five deaths in the family than those with corrugated iron. Under-five children who were the sixth born or more had 2.5 (AOR = 2.46, 95% CI: 1.49, 4.06) times higher odds of under-five mortality compared to the second and third born. Place of delivery was found to have a significant association with U5M; those who were delivered in institutions had 43% less chance of experiencing U5M (AOR = 0.57, 95% CI: 0.41, 0.81) compared to those delivered at home. Being a resident of the Somali and Afar region increased the odds of under-five death by 3.5 (AOR = 3.46, 95% CI: 1.58, 7.55) and 2.5 (AOR = 2.54, 95% CI: 1.10, 5.85) times, respectively, compared to the Tigray region (Table 4).

thumbnail
Table 4. Multivariable multilevel logistic regression analysis of both individual and community-level factors associated with under-five mortality in Ethiopia, EMDHS 2019.

https://doi.org/10.1371/journal.pone.0275586.t004

The residual plot shows that the model performed adequately, with the Kolmogorov smirnov (KS) test being insignificant (p-value = 0.29), and also the deviation was not significant (Fig 6).

thumbnail
Fig 6. residual diagnostic plot for the final model of under-five mortality in Ethiopia 2019.

https://doi.org/10.1371/journal.pone.0275586.g006

Discussion

The study mainly focused on the spatial distribution of under-five mortality and on factors that influence in Ethiopia. The findings of this study revealed that the prevalence of under-five mortality was 5.9% (95% CI: 5.4, 6.6). The multivariable multilevel logistic regression model result showed that family size, number of under-five children, birth order, breastfeeding, multiple birth, region, and institutional delivery were associated with under-five mortality in Ethiopia. The spatial analysis result revealed that Somali, Dire Dawa and Eastern Oromia were high-risk areas for under five mortality in Ethiopia.

Increased number of family size resulted in decreased odds of under-five death. This finding is consistent with previous studies conducted in Ethiopia [16], Eastern Nigeria [27], and Ghana [28]. This may be due to the mother of the children being more experienced in taking good and adequate care of the children as the family size increases, in addition to there being more family members that can care for under-five children. However, this is in contrast to another study conducted in Ethiopia, which demonstrated that the odds of under-five mortality increased as the number of under-five children increased [18]. This may indicate that, as the number of under-five children increases, it may become difficult to fulfill their basic needs, and family exhausted to look after the children.

Birth order had significant effect on under-five mortality. This study indicates that, as order of birth increased, the likelihood of child death also increased. This finding is consistent with a study conducted in rural parts of Ethiopia [9]. However, it does not support another study conducted in Ethiopia which found that being firstborn increased the odds of under-five mortality twofold when compared to those with five and above birth orders [18]. This may be because firstborns are associated with younger and less experienced mothers, which may be related to increased mortality. In addition, firstborns are more likely to be hospitalized as a result of congenital malformations and perinatal conditions in their early childhood [29]. This changes for firstborn children as they get older, while younger brothers and sisters are more likely to be hospitalized for injuries and avoidable conditions, indicating less parental care [30].

Being breastfed for less than six months increased the odds of under-five mortality compared with those who were breastfed for six and more months. This is consistent with a finding in Ethiopia [18]. Breastfeeding is known as a preventive method for reducing child mortality and as well as a prevention for delayed growth [31, 32].

Being multiple birth children increased the odds of under-five mortality when compared to singleton children. This finding is in line with several studies conducted in Ethiopia at different times, and study conducted in rural parts of Ethiopia [6, 7, 9, 16, 18]. This may be due to multiple births imposing short- and long-term medical complications as well as adding increased burden on the mother and the family [33].

Facility or institutional delivery reduced the odds of under-five mortality by 43% as compared to home delivery. This finding is consistent with studies conducted in Ethiopia [9, 10] and the Democratic Republic of Congo [13]. This is because home delivery increases the death of under-five children due to a higher risk of intrapartum and postpartum complications. These complications include exposure to common vaccine-preventable diseases if there is absence of vaccination, and a lower likelihood of attending prenatal and postpartum visits [10, 34].

The odds of under-five mortality were higher for residents of Afar and Somali regions compared to the Tigray region. This coincides among other things, with Somali and Afar regions having the first and second lowest coverage of children with all basic as well as second and first prevalence respectively [35, 36]. Additionally, the Somali region was one of the regions that were in the primary cluster from our spatial scan statistics analysis where under-five mortality is more likely to occur. This is not supported by studies conducted in Ethiopia based on the 2011 and 2016 EDHS [9, 10], where other regions were found to have a significant association with under-five mortality. This discrepancy may be a result of the different number of enumerated areas (clusters) used in the different studies, which may have led to different samples, or different statistical models used for analysis.

Conclusion

Under-five mortality in Ethiopia was 5.9% of livebirth. Being breastfed for more than 6 months, twin birth, institutional delivery and high-risk areas of U5M (Somali and dire Dawa) were modifiable risk factors. Therefore, maternal and community education on the advantage of breastfeeding and institutional delivery is highly recommended. Women who deliver multiple births should be given special attention. Further, effective strategies should be designed for interventions on hot spot areas of U5M, namely Somali and Dire Dawa.

Strengths and limitations of the study

We used EMDHS 2019 data which is more inclusive and generalizable as Enumeration Areas (EAs) are taken from most parts of the country. The larger sample size increased the power of the study and helped to make geographical comparisons. As this was cross-sectional, factors associated are not necessarily of causal effects, since results do not show a temporal relationship. Additionally, potentially important variables such as antenatal care visits, postnatal care, and preceding birth interval were excluded from the study due to their high rate of missing values, hence the study may missing important predictive factors associated with U5M.

Acknowledgments

We acknowledge the DHS program for providing us with the dataset. Elizabeth Korevaar reviewed the English language of the manuscript.

References

  1. 1. Croft T, Marshall A, Allen C. Guide to DHS statistics. Rockville: ICF; 2018. 2018.
  2. 2. Woldeamanuel BT, Aga MA. Count models analysis of factors associated with under-five mortality in Ethiopia. Global Pediatric Health. 2021;8:2333794X21989538. pmid:33623812
  3. 3. Ezbakhe F, Pérez-Foguet A. Child mortality levels and trends. Demographic Research. 2020;43:1263–96.
  4. 4. Johnston R. Arsenic and the 2030 Agenda for sustainable development. Arsen Res Glob Sustain-Proc 6th Int Congr Arsen Environ AS. 2016;2016:12–4.
  5. 5. Li Z, Hsiao Y, Godwin J, Martin BD, Wakefield J, Clark SJ, et al. Changes in the spatial distribution of the under-five mortality rate: Small-area analysis of 122 DHS surveys in 262 subregions of 35 countries in Africa. PloS one. 2019;14(1):e0210645. pmid:30668609
  6. 6. Fenta SM, Fenta HM, Ayenew GM. The best statistical model to estimate predictors of under-five mortality in Ethiopia. Journal of Big Data. 2020;7(1):1–14.
  7. 7. Worku MG, Teshale AB, Tesema GA. Determinants of under-five mortality in the high mortality regions of Ethiopia: mixed-effect logistic regression analysis. Archives of Public Health. 2021;79(1):1–9.
  8. 8. Ayele DG, Zewotir TT, Mwambi HG. Structured additive regression models with spatial correlation to estimate under-five mortality risk factors in Ethiopia. BMC public health. 2015;15(1):1–12. pmid:25884813
  9. 9. Gebremichael SG, Fenta SM. Under-Five Mortality and Associated Risk Factors in Rural Settings of Ethiopia: Evidences from 2016 Ethiopian Demographic and Health Survey. Advances in Public Health. 2020;2020.
  10. 10. Geremew BM, Gelaye KA, Melesse AW, Akalu TY, Baraki AG. Factors Affecting Under-Five Mortality in Ethiopia: A Multilevel Negative Binomial Model. Pediatric Health, Medicine and Therapeutics. 2020;11:525. pmid:33408551
  11. 11. Ghilagaber G, Antai D, Kandala N-B. Modeling spatial effects on childhood mortality via geo-additive Bayesian discrete-time survival model: A case study from Nigeria. Advanced techniques for modelling maternal and child health in Africa: Springer; 2014. p. 29–48.
  12. 12. Singh A, Masquelier B. Continuities and changes in spatial patterns of under-five mortality at the district level in India (1991–2011). International journal of health geographics. 2018;17(1):1–18.
  13. 13. Kandala N-B, Mandungu TP, Mbela K, Nzita KP, Kalambayi BB, Kayembe KP, et al. Child mortality in the Democratic Republic of Congo: cross-sectional evidence of the effect of geographic location and prolonged conflict from a national household survey. BMC public health. 2014;14(1):1–9. pmid:24649944
  14. 14. Anafcheh T, Yaghoubi Doust M, Mojadam M, Mirkazemi R, Khafaie MA. Temporal and spatial distribution of under-five mortality and factors associated with multiple cases of under-five deaths within a family in the rural area of Khuzestan, Southern Iran. Scientific reports. 2018;8(1):1–8.
  15. 15. Ayele DG, Zewotir TT. Childhood mortality spatial distribution in Ethiopia. Journal of Applied Statistics. 2016;43(15):2813–28.
  16. 16. Dejene T, Girma E. Social determinants of under-five mortality in Ethiopia: Event history analysis using evidence from Ethiopian Demographic and Health Survey (EDHS). 2013.
  17. 17. Fikru C, Getnet M, Shaweno T. Proximate determinants of under-five mortality in Ethiopia: using 2016 nationwide survey data. Pediatric health, medicine and therapeutics. 2019;10:169. pmid:31908566
  18. 18. Zewudie AT, Gelagay AA, Enyew EF. Determinants of under-five child mortality in Ethiopia: analysis using Ethiopian demographic health survey, 2016. International Journal of Pediatrics. 2020;2020. pmid:33029153
  19. 19. Abir T, Agho KE, Page AN, Milton AH, Dibley MJ. Risk factors for under-5 mortality: evidence from Bangladesh Demographic and Health Survey, 2004–2011. BMJ open. 2015;5(8):e006722. pmid:26297357
  20. 20. Pretorius CE, Asare H, Genuneit J, Kruger HS, Ricci C. Impact of breastfeeding on mortality in sub-Saharan Africa: a systematic review, meta-analysis, and cost-evaluation. European Journal of Pediatrics. 2020;179(8):1213–25. pmid:32592027
  21. 21. Ethiopian Public Health Institute—EPHI, Federal Ministry of Health—FMoH, ICF. Ethiopia Mini Demographic and Health Survey 2019. Addis Ababa, Ethiopia: EPHI/FMoH/ICF; 2021.
  22. 22. Fotheringham ASRP. The SAGE handbook of spatial analysis. Los Angeles; London: SAGE Publications; 2009.
  23. 23. Bhunia G, Shit P, Maiti R. Comparison of GIS-based interpolation methods for spatial distribution of soil organic carbon (SOC). J Saudi Soc Agric Sci 17: 114–126. 2018.
  24. 24. Kulldorff M. A spatial scan statistic. Communications in Statistics-Theory and methods. 1997;26(6):1481–96.
  25. 25. Sommet N, Morselli D. Keep calm and learn multilevel logistic modeling: A simplified three-step procedure using stata, R, Mplus, and SPSS. International Review of Social Psychology. 2017;30:203–18.
  26. 26. GeoPortal I. OpenAfrica shapfile 2019 [https://open.africa/dataset/africa-shapefiles.
  27. 27. Azuike E, Onyemachi P, Amah C, Okafor K, Anene J. Determinants of under-five mortality in South-Eastern Nigeria. J Community Med Public Health Care. 2019;6:049.
  28. 28. Aheto JMK. Predictive model and determinants of under-five child mortality: evidence from the 2014 Ghana demographic and health survey. BMC public health. 2019;19(1):1–10.
  29. 29. Kurinczuk J, Hollowell J, Boyd P, Oakley L, Brocklehurst P, Gray R. The contribution of congenital anomalies to infant mortality. Inequalities in Infant Mortality Project Briefing Paper 4. Oxford: National Perinatal Epidemiology Unit, University of Oxford; 2010.
  30. 30. Julihn A, Soares FC, Hammarfjord U, Hjern A, Dahllöf G. Birth order is associated with caries development in young children: a register-based cohort study. BMC Public Health. 2020;20(1):1–8.
  31. 31. Ahmed AH. The International Code of Marketing of Breast-milk Substitutes: Update on the Global Implementation. Journal of Human Lactation. 2020;36(4):803–7. pmid:32853085
  32. 32. WHO U. Global nutrition targets 2025: breastfeeding policy brief (WHO/NMH/NHD14. 7). Geneva: World Health Organization. 2014.
  33. 33. Elster N. Less is more: the risks of multiple births. Fertility and Sterility. 2000;74(4):617–23.
  34. 34. Goyal S, Kortsmit K, Cox S, D’Angelo DV, Romero L, Henderson ZT, et al. Prevalence of Home Births and Associated Risk Profile and Maternal Characteristics, 2016–2018. Obstetrics & Gynecology. 2020;136(6):1195–203.
  35. 35. Dessalegn M, Ayele M, Hailu Y, Addisu G, Abebe S, Solomon H, et al. Gender inequality and the sexual and reproductive health status of young and older women in the Afar region of Ethiopia. International Journal of Environmental Research and Public Health. 2020;17(12):4592. pmid:32604777
  36. 36. Oladeji O, Oladeji B, Omer MD, Farah AE, Ameda IM, Gera R, et al. Exploring opportunities to enhance effectiveness of mobile health and nutrition strategy for providing health and nutrition services amongst pastoralists in Somali region, Ethiopia. African Journal of Primary Health Care & Family Medicine. 2021;13(1):1–7. pmid:33881333