Background Ethiopia is a Sub-Saharan country with very high neonatal mortality rates, varying across its regions. The rate of neonatal mortality reduction in Ethiopia is slow, and Ethiopia may not meet the third United Nations sustainable development target by 2030. This study aimed to investigate the spatial variations and contributing factors for neonatal mortality rates in Ethiopia. Methods We analysed data from the 2016 Ethiopian Demographic and Health Survey (EDHS), which used a two-stage cluster sampling technique with a census enumeration area as primary and households as secondary sampling units. A Bayesian spatial logistic regression model using the Stochastic Partial Differential Equation (SPDE) method was fitted accounting for socio-economic, health service-related and geographic factors. Results Higher neonatal mortality rates were observed in eastern, northeastern and southeastern Ethiopia, and the Somali region had higher risks of neonatal mortality. Neonates from frequently drought-affected areas had a higher mortality risk than less drought-affected areas. Application of traditional substances on the cord increased the risk of neonatal mortality (Adjusted Odds Ratio (AOR) = 2.07, 95% Credible Interval (CrI): 1.12 to 4.30) and getting health facility delivery services had a lower odds of neonatal mortality (AOR = 0.60, 95% CrI: 0.37, 0.98). Conclusions Residing in drought-affected areas, applying traditional substances on the umbilical cord and not delivering at health facilities were associated with a higher risk of neonatal mortality. Policy-makers and resource administrators at different administrative levels could leverage the findings to prioritise and target areas identified with higher neonatal mortality rates.
Ethiopia is the second-most populous country in Africa, with a population of more than 112 million and a growth rate of 2.6% in 2019 [36]. The majority (80%) of the people reside in rural areas, with agriculture being the primary income source [37]. The study is based on secondary data analysis of the 2016 Ethiopian Demographic and Health Survey (EDHS) [22]. The EDHS is a cross-sectional survey with a nationally representative sample. The datasets contain socio-economic, neonatal, maternal, geospatial and health service use related variables. The demographic, Global Positioning System (GPS) coordinates and geospatial data were combined using the cluster (enumeration area) code. Permission was granted to access the dataset through the Demographic and Health Surveys Program [38]. The detailed descriptions of the DHS design and sampling procedures can be found elsewhere [39]. The target population of the study were newborns from birth to the 28th day from birth in Ethiopia. The samples were selected using a two-stage cluster sampling technique using a census enumeration area (EA) as the primary sampling unit and households as the secondary sampling unit. An EA or cluster is a geographic area covering, on average, 181 households. In the 2016 EDHS survey, a total of 645 clusters were sampled. For the survey to be cost-effective and produce representative data at a national and sub-national level, the DHS applies an oversampling in regions with a small population and under-sampling in regions with a larger population. Therefore, DHS applies sampling weight to restore the representativeness of the samples and correct the deliberate under-sampling and over-sampling [40]. In the computation of means, totals, and percentages, we applied sample weighting based on the DHS recommendations [41]. The most recent births of mothers were included in the study to identify the factors associated with neonatal mortality because most of the health services use related variables such as antenatal care, place of delivery, and postnatal care were recorded for only the most recent live births. The sampled total number of live births in the past five years of 2016 EDHS was 10,571, and health service use related data were collected from 7,180 most recent live births, of which 7,071 samples were collected from permanently residing respondents. One hundred nine survey participants were visitors and not included in the analysis. After data cleaning, the final analysis included 6,868 samples (see Fig 2). We adopted the Mosley and Chen conceptual framework for child survival [35] to organise variables (see Fig 1). The primary outcome of interest was neonatal mortality, defined as death within 28 days of birth. The explanatory variables at the individual level include multiple births, birth order, child sex, umbilical cord care practice, antenatal care use, place of delivery, duration of pregnancy, counselling about neonatal danger signs and living situation of the mother. Factors such as residing in urban or rural areas, the proportion of postnatal care use per cluster and episodes of drought were considered at the community level (see S1 Table for definitions of variables).