Background: Neonatal mortality remains a serious public health concern in developing countries including Ethiopia. Ethiopia is one of the countries with the highest neonatal mortality in Africa. However, there is limited evidence on the incidence and predictors of neonatal mortality at the national level. Therefore, this study aimed to investigate the incidence of neonatal mortality and its predictors among live births in Ethiopia. Investigating the incidence and predictors of neonatal mortality is essential to design targeted public health interventions to reduce neonatal mortality. Methods: A secondary data analysis was conducted based on the 2016 Ethiopian Demographic and Health Survey (EDHS) data. A total weighted sample of 11,022 live births was included in the analysis. The shared frailty model was applied since the EDHS data has hierarchical nature, and neonates are nested within-cluster, and this could violate the independent and equal variance assumption. For checking the proportional hazard assumption, Schoenfeld residual test was applied. Akakie Information Criteria (AIC), Cox-Snell residual test, and deviance were used for checking model adequacy and for model comparison. Gompertz gamma shared frailty model was the best-fitted model for this data since it had the lowest deviance, AIC value, and the Cox-Snell residual graph closet to the bisector. Variables with a p-value of less than 0.2 were considered for the multivariable Gompertz gamma shared frailty model. In the multivariable Gompertez gamma shared frailty model, the Adjusted Hazard Ratio (AHR) with a 95% confidence interval (CI) was reported to identify significant predictors of neonatal mortality. Results: Overall, the neonatal mortality rate in Ethiopia was 29.1 (95% CI: 26.1, 32.4) per 1000 live births. In the multivariable Gompertz gamma shared frailty model; male sex (AHR = 1.92, 95% CI: 1.52, 2.43), twin birth (AHR = 5.22, 95% CI: 3.62, 7.53), preceding birth interval less than 18 months (AHR = 2.07, 95% CI: 1.51, 2.85), small size at birth (AHR = 1.64, 95% CI: 1.24, 2.16), large size at birth (AHR = 1.53, 95% CI: 1.16, 2.01) and did not have Antenatal Care (ANC) visit (AHR = 2.10, 95% CI: 1.44, 3.06) were the significant predictors of neonatal mortality. Conclusion: Our study found that neonatal mortality remains a public health problem in Ethiopia. Shorter birth interval, small and large size at birth, ANC visits, male sex, and twin births were significant predictors of neonatal mortality. These results suggest that public health programs that increase antenatal care service utilization should be designed to reduce neonatal mortality and special attention should be given for twin births, large and low birth weight babies. Besides, providing family planning services for mothers to increase birth intervals and improving accessibility and utilization of maternal health care services such as ANC is crucial to improve neonatal survival.
A secondary data analysis was conducted based on the EDHS 2016 data. The 2016 EDHS survey was the fourth survey conducted in Ethiopia, situated in the Horn of Africa. Ethiopia has 9 regional states (Afar, Amhara, Benishangul-Gumuz, Gambela, Harari, Oromia, Somali, Southern Nations, Nationalities, and People’s Region (SNNPR) and Tigray) and two Administrative Cities (Addis Ababa and Dire-Dawa). The EDHS used a stratified two-stage cluster sampling technique selected in two stages using the 2007 Population and Housing Census (PHC) as a sampling frame. Stratification was achieved by separating each region into urban and rural areas. In total, 21 sampling strata have been created. In the first stage, 645 Enumeration Areas (EAs) (202 in the urban area) were selected with probability selection proportional to the EA size and independent selection in each sampling stratum. In the second stage, on average, 28 households were systematically selected. A total weighted sample of 11,022 live births within 5 years preceding the survey were included. The detailed sampling procedure was presented in the full EDHS 2016 report [32]. The outcome variable for this study was neonatal survival status categorized as being alive (coded as 0) or died (coded as 1). Neonatal mortality is defined as the death of live birth within 28 days of life. Age at death was recorded in days if the child died within 28 days of delivery. The independent variables considered for this study were categorized as socio-demographic and economic variables (residence, region, religion, maternal education, husband education, maternal occupation, sex of household head, distance to the health facility and wealth status), child-related factors (sex of neonate, type of birth, preceding birth interval and birth size), and maternal healthcare services related factors (ANC visit, place of delivery, mode of delivery, parity, birth order, wanted pregnancy, and health insurance coverage). The data were weighted using sampling weight, primary sampling unit, and strata before any statistical analysis to restore the representativeness of the survey and take into account the sampling design to get reliable statistical estimates. The sampling statisticians determine how many samples are needed in each stratum to get reliable estimates, in EDHS, some regions were oversampled, and some regions were under-sampled. So, to get statistics that are representative of Ethiopia, the distribution of neonate in the sample need to be weighted (mathematically adjusted) such that it resembles the true distribution in Ethiopia by using sampling weight (v005), primary sampling unit (v021) and strata (v022). The descriptive and summary statistics were conducted using STATA version 14 software. The EDHS data has a hierarchical structure, and therefore neonates are nested within a cluster/EAs. This violates the traditional regression model assumption, which is the independence of observations and equal variance across clusters. We have checked whether there was clustering or not by running the frailty model (random effect survival model). EA was used as a random effect (clustering variable). The theta parameter was used to assess whether there is significant clustering or not. The EDHS data were collected at two-level at individual and at the community level. Therefore, neonates in the same cluster are more of share similar characteristics than neonates in another cluster. The theta (frailty parameter) was significant at the null model (θ = 0.45, 95% CI: 0.22, 0.83). It showed that there was unobserved heterogeneity or shared frailty that needs to be taken into account to get a reliable estimate. Schoenfeld residual test was applied to check the Proportional Hazard (PH) assumptions, and the PH assumption was fulfilled (p-value> 0.05). For model selection, log-likelihood ratio test, deviance (−2LL), Akaike Information Criteria (AIC), and Cox-Snell residual plot were applied. Cox-Snell Residual test is the difference between an observed data point and a predicted or fitted value. A Cox-Snell residual considers the distribution and estimated parameters from the lifetime regression model. A model with the highest values of log-likelihood and the lowest value of AIC was the best-fitted model. Nested parametric models in generalized gamma (Exponential, Weibull, lognormal) were compared based on deviance, and non-nested models (Gompertz and log-logistic) were compared using AIC. Deviance, AIC, and Cox-Snell residual graph showed that the Gompertz gamma shared frailty model had the lowest value and the closet graph to the bisector. Therefore, the Gompertz gamma shared frailty model was the best-fitted model for the data. Variable with a p-value less 0.20 in the uni-variable gamma shared frailty analysis were included in the multivariable analysis. We estimate the hazard ratio and 95% confidence interval. In the multivariable analysis, the Adjusted Hazard Ratio (AHR) with 95% Confidence Interval (CI) was used to declare significant predictors of neonatal mortality. Permission for data access was obtained from 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 identifier.
N/A