Background: Over the past few decades, maternal and child mortality had drawn the attention of governments and policymakers. Institutional delivery has been among the implementations needed to reduce maternal and child mortality. The fact that the problem was persisted intensified studies to research for more factors. Thus, the current study was intended for further analyses of EMDHS to identify the magnitude, spatial patterns, and predictors of institutional delivery. Methods: A cross-sectional survey data from EMDHS 2019 was analyzed involving 5488 reproductive-age women regarding institutional deliveries. We presented descriptive statistics using mean, standard deviations, and proportions. To check the nature of the distribution of institutional delivery, we applied the global Moran’s I statistics. Getis-Ord Gi statistics was applied to detect spatial locations, and we applied spatial interpolation to predict unknown locations of institutional delivery using the Ordinary Kriging method. Kulldorff’s SatScan was also applied to identify the specific local clustering nature of institutional delivery using the Bernoulli method. We applied multilevel binary logistic regression for the scrutiny of individual and community-level factors. We applied P < 0.25 to include variables in the model and P < 0.05 to declare associations. AOR with 95% CI was used to describe variables. Results: The prevalence of institution/facility delivery was 2666.45(48.58%) in the survey. The average number of children was 4.03 ± 2.47, and most women in this survey were in the age range of the 25-29 years (31.84%) and 30–34 years (21.61%). Women who learned primary education (AOR = 1.52; 95% CI 1.20–1.95), secondary education (AOR = 1.77; 95% CI 1.03–3.07), and higher education (AOR = 5.41; 95% 1.91–15.25), while those who can read and write sentences (AOR = 1.94; 95% 1.28–2.94), Rich (AOR = 2.40 95% CI 1.82–3.16), and those followed 1–2 ANC (AOR = 2.08; 95% CI 1.57–2.76), 3 ANCs (AOR = 3.24; 95% CI 2.51–418), and ≥ 4 ANCs (AOR = 4.91; 95% CI 3.93–6.15) had higher odds of delivering at health institutions. Conclusion: The institutional delivery was unsatisfactory in Ethiopia, and there were various factors associated differently across the different regions. Pastoralist regions showed high home delivery than institutions which invites further interventions specific to those regions. Factors like age, highest education level achieved, preceding birth interval, literacy status, wealth status, birth order, regions, and rural residences were all affected institutional delivery so that interventions considering awareness, access, and availability of the services are vital.
Ethiopia is the country located at (3o-14oN, 33o – 48°E). The country had undertaken four standard Demographic Health Surveys (EDHS). The country started EDHS in the year 2000 and conducted every five years since then. There were also two Ethiopian Mini Demographic Health Surveys (EMDHS) conducted in 2014 and 2019. EMDHS usually conduct between the standard EDHS (two to three years) after the EDHS conducted. The 2019 EMDHS is the second nationwide mini survey conducted in the country. In Ethiopian DHS, data has been collected using a two-level multistage stratified cluster sampling to pick eligible respondents from rural and urban areas. For the current analysis, we used Ethiopia Mini Demographic Health Survey (EMDHS) 2019 data. All nine regions and two city administrations were involved in the data collection. The regions were further categorized as agrarian (Benishangul-Gumuz Amhara, Southern Nations, Nationalities, and People Gambela, Oromia, Harari, Region (SNNPR), and Tigray), pastoralists (Afar and Somali), and city administrations (Addis Ababa and Dire-Dawa) contextually. We retrieved the data from the DHS website: (www.dhsprogram.com) after the measure program allowed us to download the datasets. The weighted sample became 5488 women who had live births in the last five years before the survey. They conducted the interview on the permanent residents and visitors who stayed the day before the survey in the residences, and it was a face-to-face manner [30]. The outcome variable for this study was the health institutions/facilities delivery, which was coded as “0” if the women gave birth at home and “1” if the women gave birth at a health facility. Institutions/facilities delivery was stated as the births at health institution/facility within five years afore the survey. Maternal education, maternal age, religion, ANC follow-up, sex of household head, literacy, the total number of children, birth order, preceding birth interval, the timing of 1st ANC visit, wealth index, and marital status were the variables. Region and place of residence. Before conducting the descriptive data analysis, we weighted the data to adjust the non-proportional allocation of samples to strata and regions. Then, descriptive statistics were presented using weighted and unweighted frequencies, mean ± (standard deviations), and percentage, while all analyses were performed using STATA version 15 (STATA Corporation. IC., TX, USA). The mean-variance inflation factor also was checked to be 3.53, which was in the acceptable range. For spatial analysis, we used ArcGIS 10.7 that determined the clustering, dispersion, and random distribution nature of the institutional delivery. Moran’s I output lies between (− 1 to + 1). The values close to − 1 indicated dispersed institutional delivery, and those closes to + 1 indicated clustering distribution. After discovering significant global autocorrelation, we tested the local Getis Ord statistics to identify the areas with high and low institutional deliveries [31]. For statistical optimization of the weight, the Ordinary Kriging spatial interpolation method was applied, and enabled us to make the prediction of institutional delivery for un-sampled areas of the country. SaTScan Version 9.6 software was used for the local cluster detection. A circular window that moves systematically throughout the study area was used to identify a significant clustering of institutional delivery. We presented the results of primary and secondary observed clusters using log-likelihood (LL) and p-value < 0.05. Since the data from country representative surveys are usually clustered or have hierarchical structure, we applied multilevel analysis. We went through four consecutive models building strategies to identify felicitous predictors of institutional delivery in the country. Model 0 is an empty/null (the intercept only model) existed before addition of the predictors. Model 1 (fixed effect model) included all individual-level variables that were initially significant at p-value of < 0.25 to determine the level of variance explained by the model. Model 2 (random effect model) included cluster-level (community -level) variables and model 3 (the mixed effect model) was the final model in which both the individual and community level variables introduced to judge final model performance. The log of the probability of the institutional delivery was modeled using multilevel binary logistic regression as: logπij1-πij=β0+β0Xij+β2Zij+uij; where, i and j are the level 1 (individual) and level 2 (community) units; X and Z refer to individual and community-level variables, in sequence. Πij is the probability of the institutional delivery for the ith mother in the jth community. We resolute random effect using Intra-community Correlation (ICC), ICC=σ2aσ2a+σ2b; where, σ2a is the community level variance and σ2b indicates individual level variance. The variance (σ2b) is equal to π23 which is the fixed value. Likelihood Ratio (LR) test for model comparison and deviance (−2LL) for the goodness of fit check were calculated, while Median Odds Ratio (MOR) and Proportional Change in Variance (PCV) were also estimated [32]. Finally, the mixed effect model, which included both fixed and random effect variables were fitted. To include variable in the model p-value < 0.25 and to declare association p-value< 0.05 were used. AOR with 95% CI was also used to articulate the results.
N/A