Spatial and hierarchical Bayesian analysis to identify factors associated with caesarean delivery use in Ethiopia: evidence from national population and health facility data ============================================================================================================================================================================= * Teketo Kassaw Tegegne * Catherine Chojenta * Theodros Getachew * Roger Smith * Deborah Loxton ## Abstract **Background** Caesarean section has a significant role in reducing maternal and neonatal mortality. A linked analysis of population and health facility data is valuable to map and identify caesarean section use and associated factors. This study aimed to identify geographic variation and associated factors of caesarean delivery in Ethiopia. **Method** Linked data analysis of the 2016 Ethiopia Demographic and Health Survey (EDHS) and the 2014 Ethiopian Service Provision Assessment Plus (ESPA+) survey was performed. Spatial analysis was conducted to identify geographic variations and factors associated with caesarean delivery. Hierarchical Bayesian analysis was also performed to identify factors associated with caesarean delivery using the SAS MCMC procedure. **Results** Women’s age and education, household wealth, parity, antenatal care (ANC) visits, and distance to caesarean section facility were associated with caesarean delivery use. Women who had ≥4 ANC visits were 4.67 (95% Credible Interval (CrI): 2.17, 9.43) times more likely to have caesarean delivery compared to those who had no ANC visits. Women who had education and were from rich households were also 2.80 (95% CrI: 1.83, 4.19) and 1.80 (95% CrI: 1.08, 2.84) times more likely to have caesarean deliveries relative to women who had no education and were from poor households, respectively. A one-kilometer increase in distance to a caesarean section facility was associated with an 88% reduction in the odds of caesarean delivery (Adjusted Odds Ratio (AOR) = 0.12, 95% CrI: 0.01, 0.78). Hotspots of high caesarean section rates were observed in Addis Ababa, Dire Dawa, and the Harari region. In addition, women’s age at first childbirth and ≥4 ANC visits showed significant spatially varying relations between caesarean delivery use across Ethiopia. **Conclusion** Caesarean section is a lifesaving procedure, and it is essential to narrow disparities to reduce maternal and neonatal mortality and avoid unnecessary procedures. Keywords * Caesarean delivery * caesarean section * obstetric care * spatial variations * geographic variations * Ethiopia ## Introduction Caesarean section or caesarean delivery is a surgical procedure to deliver a baby, which is often performed when vaginal delivery is not possible or puts the baby and/or the mother at risk (1). A medically indicated caesarean section is effective in saving the lives of mothers and newborns (2). The World Health Organization (WHO) stated that a caesarean section rate above 10% at the population level was not associated with a reduction in maternal and newborn deaths (3, 4). A recent analysis of global trends in caesarean deliveries reported that 21.1% of deliveries were by caesarean section worldwide (5). It is also projected that by 2030, 28.5% of deliveries will be by caesarean section globally (5). Despite an increasing trend in caesarean sections globally, inequalities exist between and within countries (5). Inequalities in caesarean sections could be explained by different demand and supply-side factors. The demand-side factors, including individual, household, and community-level factors, are significantly associated with caesarean delivery use (6-8). The supply-side factors, such as availability and quality of services, facilities’ readiness to provide services, and providers’ competency and caring behaviour, influence women’s use of caesarean sections. A systematic review reported that health facilities in many developing countries are underperforming in these areas (9). In Ethiopia, there was a slight increase in the national rate of caesarean sections from 2000 (0.7%) to 2016 (1.9%) (10). However, there is a gap in the availability and use of caesarean sections depending on geographic locations. It is, therefore, difficult to identify areas with unmet needs and the types of interventions required to meet caesarean section needs. A linked analysis of population and health facility data has a significant role in examining the link between population healthcare needs, service utilization and the health service environment (11, 12). In Ethiopia, the national population and health facility data can be accessed from the Ethiopia Demographic and Health Survey (EDHS) (13) and the Ethiopian Service Provision Assessment Plus (ESPA+) survey (14), respectively. Mapping caesarean delivery use and associated factors using these data would provide valuable information for program managers and decision-makers and help to identify communities requiring special attention. Thus, this study aimed to identify: 1) geographic variations in caesarean delivery use – spatial analysis, and 2) factors associated with caesarean delivery use – multilevel analysis in Ethiopia using the national population and health facility data. ## Methods ### Data sources This study used two data sets: the 2016 EDHS and the 2014 ESPA+ survey. #### Population data The 2016 EDHS was used to obtain national population data. It contained population socio-demographic characteristics, health service use and morbidity/mortality data. The survey used a stratified multistage sampling technique. EDHS clusters, which were census Enumeration Areas in Ethiopia, were used as sampling frames. In addition to population information, the latitude and longitude coordinates of each EDHS clusters was recorded. Data were collected from 645 EDHS clusters (15). In this analysis, we included 6,954 women who gave birth in the five years preceding the survey from 622 EDHS clusters. We excluded 239 women who gave birth during this period from 23 EDHS clusters since they had no geographic coordinates for data linkage. #### Health facility data The 2014 ESPA+ survey provided national health facility data covering different aspects of the availability and readiness of health services. In addition, the data set contained the latitude and longitude coordinates of each health facility. The survey used a combination of census and simple random sampling techniques (16). Due to their importance and limited numbers, the survey included all hospitals. Details of the sampling procedure are discussed elsewhere (12). In this analysis, we included only hospitals proving caesarean sections. Thus, we included only 179 of the 214 hospitals that reported providing caesarean sections. ### Data linking method This study used an Euclidean buffer link to link population and health facility data (12). This is regarded as the method of choice for linking a census of health facilities with population data (11, 12, 17). Details of this method are discussed elsewhere (12). The administrative boundaries of Ethiopia, obtained from Natural Earth (18), were used. ### Health service environment Using the ESPA+ survey, four health service environment variable scores were computed. These were the average distance to the nearest caesarean delivery facilities, comprehensive obstetric care availability score, readiness to provide comprehensive obstetric care score, and a general health facilities readiness score. Using principal component analysis, availability and readiness scores were computed for the nearest hospital providing caesarean delivery. The SCORE procedure in SAS was used to compute availability and readiness scores. The comprehensive obstetric care indices were created using the World Health Organization’s ‘Service Availability and Readiness Indicators’ (19, 20). Details of computing these scores were discussed elsewhere (12). Regarding general service readiness, six service readiness dimensions were used: communication equipment, external supervision, client opinion and feedback, quality assurance, emergency transport, and client latrine. The first two principal components (health facility management system and infrastructure) were used to compute two general service readiness scores (12). For hospitals providing caesarean delivery, indices of comprehensive obstetric care availability and readiness scores were created. Two comprehensive obstetric care availability scores (basic and comprehensive components) were created using seven variables. Similarly, two comprehensive obstetric care readiness scores (equipment and supplies, and skilled personnel) were computed using nine dichotomous variables (12). ### Study variables The outcome variable of this study was caesarean delivery use. A woman was considered to have used caesarean delivery if her most recent birth in the five years preceding the survey was via caesarean section. Sociodemographic and obstetric characteristics and health facility variables were considered independent variables for this study. Since this study is a multilevel analysis, all independent variables included are of three types: level 1, level 2 and level 3. Level 1 or individual-level variables include women’s age, education, and occupation, women’s autonomy in healthcare decision-making, husband’s education and occupation, household wealth, age at first childbirth, parity, antenatal care (ANC) visits, and nature of recent pregnancy (planned or unplanned). Level 2 or EDHS cluster-level variables were general service readiness, comprehensive obstetric care availability and service readiness, and average distance to the nearest caesarean delivery facility. One level 3 or region-level (number of caesarean sections providing hospitals per region) was considered in the analysis. Respondents’ educational status was grouped into two categories: had no education and had an education. Similarly, occupational status was grouped as having work and no work. Respondents who responded not working at the time of the interview or did not work in the last 12 months before the survey were grouped as having no work. The EDHS household wealth quintiles were also grouped into two categories: poor (including lowest, second and middle quintiles) and rich (including fourth and highest quintiles). ANC visits were grouped as no ANC visits, 1 – 3 ANC visits and ≥4 ANC visits. ### Statistical analysis #### Bayesian hierarchical model The EDHS survey employed a multistage cluster sampling technique: women were nested within EDHS clusters and then clusters were nested within regions. Considering this hierarchical nature of the data, a three-stage Bayesian hierarchical model was used. Women’s data within the respective EDHS clusters were also linked with the health facility data. In this model, the probability of caesarean delivery is allowed to vary both between EDHS clusters and between regions. To model the between-EDHS cluster and between-region variability, the logarithms of the odds ratios of caesarean delivery was assumed to have a normal distribution. The means of the normal distribution of log odds ratios across EDHS clusters and regions, therefore, represent the average effect in the EDHS clusters and regions, and the variances represent the variability among EDHS clusters and regions. Bayesian inference allows for the combination of existing knowledge with new information according to established rules of probability. In this study, a non-informative prior that would not have an influence on the posterior distribution was chosen (21, 22). For each fixed effect, a non-informative prior that follows a normal distribution with large variance (*σ*2=1000) and a zero mean (*μ*=0) was used. Similarly, for the random effects, non-informative priors that were set up to follow a normal distribution with zero mean and different variances were used. The variances for the random effects were set up to follow an inverse gamma distribution with a shape parameter of *α*=0.01 and a scale parameter *β*=0.01. The MCMC procedure (PROC MCMC) in SAS was used to estimate the Bayesian hierarchical generalized linear mixed models. The study used a simulation size of 175,000 iterations along with 25,000 burn-in iterations. One of every 25 samples was kept, which gave the Markov chain Monte Carlo (MCMC) sample size of 7,000 for computing the posterior quantities of interest. This analysis procedure enabled the identification of potential factors associated with the rate of caesarean delivery with an equal-tail 95% credible interval. The model building process was started with the unconditional model (a model containing no predictors) and more complex models were gradually built by checking improvements in model fit after each model is estimated (23). A Deviance Information Criterion (DIC) was calculated using the posterior mean estimates of the parameters to select the best fitting model. The unconditional (empty) model was used to calculate the intra-class correlation coefficient (ICC), which estimates how much variation in the use of caesarean delivery exists between EDHS clusters (level-2 units) and between regions (level-3 units) using posterior samples. Three variance estimates were used to calculate EDHS cluster ICC and region ICC. The formula for calculating each of these ICCs is similar to the formula used for two-level ICC calculations. However, the denominator now has three elements (residual variance, cluster variance and region variance) which form the total variance in the model. The two equations for the EDHS cluster ICC (*ICC**C*) and region ICC (*ICC**R*) are below. ![Formula][1] The *ICC**C* tells us how much of the total variation in caesarean delivery exists between EDHS clusters. On the other hand, *ICC**R* indicates the total variation in caesarean delivery that exists between regions (24). In hierarchical generalized linear mixed models, it is assumed that there is no level-1 error variance; however, to calculate the intra-class correlation coefficient, a slight modification is made. The level-1 residual variance (***ε*)** follows a logistic distribution and is standardized with a mean of zero and variance = *π*2/3 (25). Therefore, for three-level random intercept hierarchical generalized linear mixed model with an intercept variance of ![Graphic][2] and ![Graphic][3], the intra-class correlation coefficient (Rho) is given by; ![Formula][4] Using posterior samples from the MCMC simulation, proportional change in variance (PCV), for EDHS clusters and regions, was also estimated to measure changes in area level variance between empty and individual level models and between successive models (26, 27). This was computed using a mathematical equation: ![Graphic][5]; where *V**n*−1 is the neighbourhood variance in the empty model and *V**n*−2 is the neighbourhood variance in the subsequent model. Based on the model fit statistics, allowing EDHS cluster-level variables to vary across regions and adding a region-level variable did not show improvement in model fit statistics. Therefore, a model with six individual-level variables (women’s age, age at first childbirth, education, household wealth, parity, and the number of ANC visits, and one EDHS cluster level variable (average distance to the nearest caesarean delivery facility) was selected. In this model, individual-level variables were allowed to vary across EDHS clusters and regions. #### Spatial analysis ArcGIS 10.6.1 and R version 4.2.0 were used to carry out the spatial analysis. The Ethiopian Polyconic Projected Coordinate System was used to flatten the Ethiopian map (12). Hotspot analysis and spatial regression were performed to identify spatial clusters and factors associated with spatial variations of caesarean delivery use. The unit of spatial analyses were EDHS DHS clusters. Three procedures were followed to carry out the spatial analysis as discussed elsewhere (12). First, the Global Moran’s I statistic, which is a global measure of spatial autocorrelation (28) was carried out. Second, the Incremental Spatial Autocorrelation was used to determine the critical distance at which there is maximum clustering (12). The average distance *(15 kilometres)* at which a feature has at least one neighbour was calculated. This gave 155 kilometres at which clustering of caesarean delivery use peaked. Lastly, the Getis-Ord Gi* statistic was used to identify statistically significant spatial clusters of high caesarean delivery rates (hotspots) and low caesarean delivery rates (cold spots). A False Discovery Rate (FDR) correction method was applied to account for multiple and spatial dependence tests in Local Statistics of Spatial Association (12, 29). Statistical significance was determined based on z-scores and p-values (12). Furthermore, spatial regression was performed to identify key factors behind the observed spatial patterns of caesarean delivery use. Moran eigenvector spatially and non-spatially varying coefficient (M-SNVC) model was used to model spatially and non-spatially varying relations. Eigenvector spatial filtering regression model removes spatially autocorrelated residuals and improves model fit (30, 31). M-SNVC model is also stable and quite robust against spurious correlations compared to the Moran eigenvector spatially varying coefficient model (32). Since the spatial regression was done at the EDHS cluster level, the number of caesarean deliveries per EDHS cluster was first computed. Then we used a transformation-based generalized spatial regression using the *spmoran* package implemented in R (33). Since the count data did not obey the Poisson distribution, the log-Gaussian approximation was applied first to roughly normalize the data, and then the SAL transformation to identify the most likely distribution (33). Finally, the transformed result from the generalized spatial regression model was used in the M-SNVC model to model the spatial relationship between caesarean delivery use and explanatory variables. ### Ethics approval This study used secondary data sets: 2016 EDHS and 2014 ESPA+ collected previously where confidentiality information was maintained. Hence, as we did not collect data directly from participants, recruitment and contacting participants were not required. Ethical approval was obtained from The University of Newcastle, the Ethiopian Public Health Institute and the Measure DHS program to access the datasets. ## Results ### Sociodemographic and obstetric characteristics The sociodemographic and obstetric characteristics of the women were not presented in this paper. These characteristics can be found in our previous publications with different outcomes (34, 35). In this paper, women’s sociodemographic and obstetric characteristics are similar to these papers. ### Hospital characteristics Amongst the 214 hospitals, 200 (93.46%) reported providing normal delivery care. Out of the 200 hospitals, 179 (89.50%) provided caesarean delivery. The national average distance of caesarean delivery providing hospitals to the 2016 EDHS clusters was 33.27 km. The 2016 EDHS sampled clusters in the Somali region had the longest distance to caesarean delivery providing hospitals, 71.36 km. Conversely, on average, EDHS clusters in Addis Ababa were 1.54 km away from caesarean section providing hospitals (Fig 1). ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/10/10/2022.10.07.22280820/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/F1) Figure 1: The average distance of the demographic and health survey clusters to caesarean delivery providing hospitals in Ethiopia, 2016 (n = 179) ### Prevalence of caesarean delivery In Ethiopia, the prevalence of caesarean delivery was found to be 3.65% (urban 12.70%, 1.14% rural). The prevalence of caesarean deliveries were 7.34% in the public and 30.20% in private health facilities. The prevalence of caesarean delivery varied across the different regions and city administrations of the country; the highest caesarean delivery rate was reported in Addis Ababa (21.87%) followed by Harari (11.19%) and Dire Dawa (8.31%). The map (Fig 2) showed that there is regional variation in caesarean delivery rates. ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/10/10/2022.10.07.22280820/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/F2) Figure 2: Caesarean delivery use among pregnant women in Ethiopia, 2016 ### Hospital to the population proportion There is a huge difference in the proportion of hospitals per one-million population in each region. In the Amhara, Oromia and Southern Nations, Nationalities and Peoples (SNNP) regions, approximately one hospital is expected to serve one million people. In the Tigray region, it was estimated that six hospitals serve one million people. In the two city administrations (Addis Ababa and Dire Dawa) and the Harari region, an average of 14 hospitals serve one-million people. The rates of caesarean delivery per EDHS clusters also varied depending on the number of, and distance to, caesarean delivery providing hospitals in each region (Fig 3). ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/10/10/2022.10.07.22280820/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/F3) Figure 3: Distribution of hospitals and caesarean delivery rate in Ethiopia, 2016 ### Determinants of caesarean delivery use The calculated EDHS cluster level intra-class correlation coefficient (*ICC**c*) was 15.95%. This indicated that about 16% of variations in using caesarean delivery exist between EDHS clusters. Similarly, the estimated region level ICC (*ICC**R*) indicated that 33.59% of variations in caesarean section rates exist between regions. This indicated that only 50.46% of variability is accounted for by the differing characteristics of the women, or other unmeasured factors. A practically meaningful proportion of variance in caesarean delivery use exists at the EDHS cluster and region levels, providing support for using three-level analytical model. Women’s age and education, age at first childbirth, household wealth, parity, and ANC visits were significant individual-level caesarean section predictors. A one-year increase in the age of a woman was associated with a 28.13 times increase in the odds of caesarean delivery use. Similarly, a single-year increase in the age of women at first childbirth was associated with a 7.27 times increase in the odds of caesarean delivery. Women who had education were 2.80 times more likely to deliver via caesarean sections compared to those who had no education. A woman from a rich household was 1.80 times more likely to have caesarean delivery relative to a woman from a poor household. Women who had ≥4 ANC visits were 4.67 times more likely to use caesarean delivery compared to those who had no ANC visits. On the other hand, a one-child increase in the number of live births a woman had was associated with a 97% decrease in the odds of caesarean delivery use (Table 1). View this table: [Table 1:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/T1) Table 1: Factors associated with caesarean delivery use among pregnant women in Ethiopia, 2016 At the EDHS cluster level (level 2), only one variable was significantly associated with the use of caesarean delivery. A one-kilometer increase in distance to a caesarean section facility was associated with an 88% reduction in the odds of caesarean delivery (AOR = 0.12, 95% CrI: 0.01, 0.78) (Table 1). Finally, most of the between-EDHS clusters and between-region variances were explained by this model. The proportional change in variance indicated that adding predictors to the empty model explained an increased proportion of variation in caesarean delivery use. The variance estimates between-EDHS clusters decreased from 1.04 in the empty model to 0.05 in the final model. The proportion of the variance of the between-EDHS cluster explained by the final model was 95.19%. Similarly, the variance estimates between-regions decreased from 2.19 in the empty model to 0.14 in the final model. The proportion of the between-region variance explained by the final model was 93.61%. In addition, the empty model showed that 15.95% (*ICC**C*) and 33.5% (*ICC**R*) of the variability in the odds of caesarean delivery use was explained by EDHS-cluster-level and region-level characteristics, respectively. The between-EDHS cluster variability declined over successive models, from 15.95% in the empty model to 1.44% in the final model. Similarly, the between-region variability decreased from 33.59% in the empty model to 4.02% in the final model (Table 2). View this table: [Table 2:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/T2) Table 2: Variations in caesarean delivery use in Ethiopia: random slope and random intercept model ### Hotspots of caesarean delivery use There was evidence of spatial clustering in caesarean delivery use among pregnant women in Ethiopia (Global Moran’s I = 1.15, Z-score = 38.26, P-value < 0.0001). Most of the hotspots were in Addis Ababa, Dire Dawa and Harari, and some parts of the Oromia, Amhara, Somali and SNNP regions. Most of the cold spots were in the Gambela, Benishangul-Gumuz and Tigray regions followed by Amhara and SNNP regions (Fig 4). ![Figure 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/10/10/2022.10.07.22280820/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/F4) Figure 4: Clusters of high and low caesarean delivery rates in Ethiopia, 2016 ### Determinants of spatial variations in caesarean delivery use In the spatial regression analysis, we found that women’s age and education, parity, ANC visits, age at first childbirth, household wealth, and distance to caesarean delivery facilities were associated with spatial variations in caesarean delivery use. Parity was negatively associated with the spatial variations in caesarean delivery use (Table 3). View this table: [Table 3:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/T3) Table 3: Factors associated with the spatial variations of caesarean delivery use in Ethiopia, 2016 In the M-SNVC regression model, we found that the relationship between caesarean delivery use and age at first childbirth and ≥4 ANC visits varied across EDHS clusters (Fig 5 & 6). Caesarean delivery use was positively associated with ≥4 ANC visits across Ethiopia (Fig 5). On the other hand, age at first childbirth showed both negative and positive effects on caesarean delivery use across the country. Age at first childbirth showed a strong negative impact on caesarean delivery use in the Tigray and eastern part of Amhara regions and a strong positive impact in the SNNP region (Fig 6). However, the coefficients of the other variables were estimated constant. ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/10/10/2022.10.07.22280820/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/F5) Figure 5: Spatial variations of caesarean delivery use with spatially varying coefficients of ≥4 ANC visits in Ethiopia, 2016 ![Figure 6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/10/10/2022.10.07.22280820/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2022/10/10/2022.10.07.22280820/F6) Figure 6: Spatial variations of caesarean delivery use with spatially varying coefficients of age at first childbirth in Ethiopia, 2016 ## Discussion This study provides a comprehensive assessment of unmet needs for caesarean section in Ethiopia by location, service type and demographics. Despite an overall increase in caesarean sections in Ethiopia, there were differences by geographic area, and the rate is lower than the WHO recommendation (36). Furthermore, consistent with previous studies (37-42), caesarean sections in private health facilities were higher than in public hospitals. Caesarean sections for medically indicated reasons can effectively prevent maternal and newborn mortality (43, 44). However, overuse of caesarean sections could increase maternal and perinatal morbidity (45). Private health facilities, which are mainly profit-driven, could do unnecessary caesarean sections to protect themselves from unexpected delivery risks (46). Interventions targeting clinical and non-clinical drivers are required to narrow this gap and avoid unnecessary caesarean sections. Hotspots of high caesarean section rates were observed in urban centers, i.e., in Addis Ababa, Dire Dawa and the Harari region. This could be due to differences in the distribution of caesarean delivery facilities. Most emergency obstetric care (EmOC) facilities, both public and private, are concentrated in urban areas. These three cities had the highest proportion of hospitals compared to other administrative regions in Ethiopia. Furthermore, in the spatial regression analysis, the geographic variations in caesarean delivery use were explained by women’s age and education, parity, ANC visits, age at first childbirth, household wealth, and distance to caesarean delivery facilities. Spatially varying relations were also observed between caesarean delivery use and ≥4 ANC visits and age at first childbirth. This finding enables informed decision-making for which communities and health facilities need special attention and where government should spend more money. Women’s use of caesarean sections was influenced by different individual and cluster-level factors. Consistent with previous studies (47-50), a single-year increase in the age of women and age at first childbirth were positively associated with caesarean delivery use. Pregnant women of advanced age are more likely to have pregnancy complications (51-53), comorbidities (54, 55) and higher rates of labour inductions (56). Thus, pregnant women of advanced age should be aware that they might be at increased risk of caesarean section and should prepare in advance. The likelihood of caesarean delivery, consistent with a recent meta-analysis (8), was also positively associated with women’s education. Education could influence women’s decision-making power and delivery mode preference. A systematic review and meta-analysis on women’s caesarean section preference reported that 15.6% of women preferred to give birth via caesarean section (57). It is, therefore, necessary to counsel pregnant women about the potential benefits and risks of caesarean delivery. Consistent with previous studies (7, 8), women from rich households were more likely to have caesarean delivery compared to women from poor households. Women from rich households could afford expenses related to caesarean delivery and transportation. Furthermore, a single unit increase in the number of live births was significantly associated with lower odds of caesarean section. This finding was consistent with previous studies in Ethiopia (58), India (59), Kenya (60) and Rwanda (6). Multiparas could have faster labour progress compared to nulliparous women (61), supported by the effect of parity on uterine activity study (62). This finding indicates that women in their first pregnancy should get more emphasis during their ANC visits and prepare them to plan skilled delivery care. Consistent with studies in Rwanda (6), India (42) and Bangladesh (7), the odds of caesarean sections were higher among women with more ANC visits. As per the WHO recommendations, ANC visits are for risk identification, prevention and management of complications, and health education and promotion (63). Women with pregnancy complications identified during ANC visits could have more ANC visits and increased caesarean section rates. In addition, women could prefer caesarean sections if they are uncertain about their fetal condition (64). Thus, the risks and benefits of caesarean and vaginal deliveries should be discussed during ANC visits to avoid unnecessary caesarean sections and birth complications. Further research is needed to investigate why caesarean sections increase along with ANC visits. Furthermore, a one-kilometer distance increase in caesarean section facility was inversely associated with women’s use of caesarean delivery. Women with medically indicated caesarean sections would not get the service due to limited geographic access to caesarean section facilities (65-68). Thus, limited access to caesarean section facilities could be associated with high maternal and neonatal mortality (9). Despite caesarean service availability and hospital readiness scores did not show significant associations with caesarean delivery, these variables have enormous importance and should not be understated. This non-significant finding could be due to chance or bias, which we could not determine. United Nations recommends that there should be at least five EmOC facilities, including one comprehensive emergency obstetric care facility, per 500,000 population (69), and governments should work to achieve this goal. The findings of this study should be considered with the following limitations in mind. First, removing EDHS clusters without geographic coordinates could under/overestimate caesarean section rate. Second, misclassification errors due to EDHS coordinates displacement (12) could influence the impact of distance on caesarean delivery use. Third, regarding Bayesian analysis, there is no correct way to choose a prior where we decided to use non-informative priors. Finally, Bayesian analysis often comes with high computational costs (memory and time), especially in models with many parameters and multilevel analysis. Despite these limitations, this study has several strengths. This study used national population and health facility data to identify the demand and supply-side determinants of caesarean delivery use, which most studies investigated separately. Bayesian inferences also solve sample size problems, which is the limitation of classical statistics. We also used spatial analysis to identify geographic variations and associated factors in caesarean delivery use. Investigating geographic variations has a significant role in informed decision-making and monitoring and evaluation purposes. ## Conclusion The prevalence of caesarean delivery was low in Ethiopia. There were also significant geographic variations across the country, and different individual and community-level factors were associated with caesarean delivery use. ## Data Availability All data produced in the present work are contained in the manuscript. ## Authors’ contributions TKT, CC, RS, DL conceptualized the design of the analysis. TKT developed and drafted the manuscript. CC, TG, RS and DL participated in critically revising the intellectual contents of the manuscript. All authors read, provided feedback and approved the final manuscript. ## Acknowledgements We thank the University of Newcastle, Australia for offering free access to the digital online library to search the electronic databases that were considered for this analysis. We also thank the Measure DHS Program and the Ethiopian Public Health Institute for providing free access to the data sets used for this analysis. * Received October 7, 2022. * Revision received October 7, 2022. * Accepted October 10, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Karanth L, Abas AB: Maternal and foetal outcomes following natural vaginal versus caesarean section (c-section) delivery in women with bleeding disorders and carriers. Cochrane Database Syst Rev 2021(12). 2. 2.World Health Organization: WHO Statement on Caesarean Section Rates. In. Geneva: WHO/RHR/15.02; 2015. 3. 3.Betran AP, Torloni MR, Zhang J, Ye J, Mikolajczyk R, Deneux-Tharaux C, Oladapo OT, Souza JP, Tunçalp Ö, Vogel JP: What is the optimal rate of caesarean section at population level? A systematic review of ecologic studies. Reprod Health 2015, 12(1):57. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12978-015-0043-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26093498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 4. 4.Ye J, Zhang J, Mikolajczyk R, Torloni M, Gülmezoglu A, Betran A: Association between rates of caesarean section and maternal and neonatal mortality in the 21st century: a worldwide population-based ecological study with longitudinal data. BJOG: An International Journal of Obstetrics & Gynaecology 2016, 123(5):745–753. 5. 5.Betran AP, Ye J, Moller A-B, Souza JP, Zhang J: Trends and projections of caesarean section rates: global and regional estimates. BMJ Global Health 2021, 6(6):e005671. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYm1qZ2giO3M6NToicmVzaWQiO3M6MTE6IjYvNi9lMDA1NjcxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMTAvMTAvMjAyMi4xMC4wNy4yMjI4MDgyMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 6. 6.Kibe PM, Mbuthia GW, Shikuku DN, Akoth C, Oguta JO, Ng’ang’a L, Gatimu SM: Prevalence and factors associated with caesarean section in Rwanda: a trend analysis of Rwanda demographic and health survey 2000 to 2019–20. BMC Pregnancy and Childbirth 2022, 22(1):1–13. 7. 7.Ahmed MS, Islam M, Jahan I, Shaon IF: Multilevel analysis to identify the factors associated with caesarean section in Bangladesh: evidence from a nationally representative survey. Int Health 2022. 8. 8.Islam M, Sathi NJ, Hossain M, Jabbar A, Renzaho A, Islam SMS: Caesarean delivery and its association with educational attainment, wealth index, and place of residence in Sub-Saharan Africa: a meta-analysis. Sci 2022, 12(1):1–14. 9. 9.Paxton A, Bailey P, Lobis S, Fry D: Global patterns in availability of emergency obstetric care. International Journal of Gynecology & Obstetrics 2006, 93(3):300–307. 10. 10.Yisma E, Smithers LG, Lynch JW, Mol BW: Cesarean section in Ethiopia: prevalence and sociodemographic characteristics. The Journal of Maternal-Fetal & Neonatal Medicine 2017:1–6. 11. 11.Burgert-Brucker CR, Prosnitz D: Linking DHS household and SPA facility surveys: Data considerations and Geospatial Methods. DHS Spatial Analysis Reports No. 10: ICF International; 2014. 12. 12.Tegegne TK, Chojenta C, Getachew T, Smith R, Loxton D: Service environment link and false discovery rate correction: Methodological considerations in population and health facility surveys. PLoS ONE 2019, 14(7):e0219860. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 13. 13.Wang W, Winter R, Mallick L, Florey L, Burgert-Brucker C, Carter E: The Relationship between the Health Service Environment and Service Utilization: Linking Population Data to Health Facilities Data in Haiti and Malawi. DHS Analytical Studies No. 51. Rockville, Maryland, USA: ICF International; 2015. 14. 14.Hozumi D, Fronczak N, Minichiello SN, Buckner B, Fapohunda B, Kombe G, Searing H, Ramarao S, Ricca J: Profiles of Health Facility Assessment Methods. Report of the International Health Facility Assessment Network (IHFAN). In.: MEASURE Evaluation, USAID; 2008. 15. 15.Central Statistical Agency (CSA) [Ethiopia], ICF: Ethiopia Demographic and Health Survey 2016. In. Addis Ababa, Ethiopia, and Rockville, Maryland, USA: CSA and ICF; 2016. 16. 16.Ethiopian Public Health Institute, ICF International: Ethiopia Service Provision Assessment Plus (ESPA+) Survey 2014. In. Addis Ababa, Ethiopia and Rockville, Maryland USA: Ethiopian Public Health Institute and ICF International; 2014. 17. 17.Skiles MP, Burgert CR, Curtis SL, Spencer J: Geographically linking population and facility surveys: methodological considerations. Population health metrics 2013, 11(1):14. 18. 18.Free vector and raster map data [[https://www.naturalearthdata.com/downloads/10m-cultural-vectors/](https://www.naturalearthdata.com/downloads/10m-cultural-vectors/)] 19. 19.World Health Organization: Service availability and readiness assessment (SARA): an annual monitoring system for service delivery: reference manual. 2015. 20. 20. World Health Organization: Service Availability and Readiness Assessment (SARA): an annual monitoring system for service delivery. In.: Health Statistics and Information Systems, WHO; 2013. 21. 21.Babapulle MN, Joseph L, Bélisle P, Brophy JM, Eisenberg Mjjtl: A hierarchical Bayesian meta-analysis of randomised clinical trials of drug-eluting stents. 2004, 364(9434):583–591. 22. 22.Ojo OB, Lougue S, Woldegerima WAJPo: Bayesian generalized linear mixed modeling of Tuberculosis using informative priors. 2017, 12(3):e0172580. 23. 23.Ene M, Leighton EA, Blue GL, Bell BA: Multilevel models for categorical data using SAS® PROC GLIMMIX: the basics. In: SAS Global Forum 2015 Proceedings: 2015; 2015: 3430–2015. 24. 24.Hox JJ, Moerbeek M, Van de Schoot R: Multilevel analysis: Techniques and applications: Routledge; 2017. 25. 25.Tom A, Bosker Tasrj, Bosker RJ: Multilevel analysis: an introduction to basic and advanced multilevel modeling: Sage; 1999. 26. 26.Merlo J, Yang M, Chaix B, Lynch J, Råstam L: A brief conceptual tutorial on multilevel analysis in social epidemiology: investigating contextual phenomena in different groups of people. J Epidemiol Community Health 2005, 59(9):729–736. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiamVjaCI7czo1OiJyZXNpZCI7czo4OiI1OS85LzcyOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzEwLzEwLzIwMjIuMTAuMDcuMjIyODA4MjAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 27. 27.Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, Råstam L, Larsen K: A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health 2006, 60(4):290–297. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiamVjaCI7czo1OiJyZXNpZCI7czo4OiI2MC80LzI5MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzEwLzEwLzIwMjIuMTAuMDcuMjIyODA4MjAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 28. 28.Anselin L: Local indicators of spatial association—LISA. Geographical analysis 1995, 27(2):93–115. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=https://doi.org/10.1111/j.1538-4632.1995.tb00338.x&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QU19800001&link_type=ISI) 29. 29.Caldas de Castro M, Singer BH: Controlling the false discovery rate: a new application to account for multiple and dependent tests in local statistics of spatial association. Geographical Analysis 2006, 38(2):180–208. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.0016-7363.2006.00682.x&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236275400005&link_type=ISI) 30. 30.Tiefelsdorf M, Griffith Daje, A P: Semiparametric filtering of spatial autocorrelation: the eigenvector approach. 2007, 39(5):1193–1221. 31. 31.Thayn JB, Simanis JMJAotAoAG: Accounting for spatial autocorrelation in linear regression models using spatial filtering with eigenvectors. 2013, 103(1):47–66. 32. 32.Murakami D: Spatial regression using the spmoran package: Boston housing price data examples. 2020. 33. 33.Murakami D: Transformation-based generalized spatial regression using the spmoran package: Case study examples. arXiv preprint arXiv:210906383 2021. 34. 34.Tegegne TK, Chojenta C, Getachew T, Smith R, Loxton D: Antenatal care use in Ethiopia: a spatial and multilevel analysis. BMC pregnancy and childbirth 2019, 19(1):399. 35. 35.Tegegne TK, Chojenta C, Getachew T, Smith R, Loxton D: Giving birth in Ethiopia: a spatial and multilevel analysis to determine availability and factors associated with healthcare facility births. BJOG: An International Journal of Obstetrics & Gynaecology 2020, 127(12):1537–1546. 36. 36.Betran A, Torloni M, Zhang J, Gülmezoglu A: WHO Statement on caesarean section rates. BJOG: An International Journal of Obstetrics & Gynaecology 2016, 123(5):667–670. 37. 37.Vieira GO, Fernandes LG, de Oliveira NF, Silva LR, de Oliveira Vieira T: Factors associated with cesarean delivery in public and private hospitals in a city of northeastern Brazil: a cross-sectional study. BMC pregnancy and childbirth 2015, 15(1):132. 38. 38. Almeida Sd, Bettiol H, Barbieri MA, Silva AAMd, Ribeiro VS: Significant differences in cesarean section rates between a private and a public hospital in Brazil. Cad Saude Publica 2008, 24:2909–2918. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19082282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262995900020&link_type=ISI) 39. 39.Rebelo F, Da Rocha Cmm, Cortes TR, Dutra CL, Kac G: High cesarean prevalence in a national population-based study in Brazil: the role of private practice. Acta Obstet Gynecol Scand 2010, 89(7):903–908. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/00016349.2010.484044&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20583936&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 40. 40.Mendoza-Sassi RA, Cesar JA, Silva PRd, Denardin G, Rodrigues MM: Risk factors for cesarean section by category of health service. Revista de saúde pública 2010, 44(1):80–89. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20140332&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 41. 41.Beogo I, Rojas BM, Gagnon M-P: Determinants and materno-fetal outcomes related to cesarean section delivery in private and public hospitals in low-and middle-income countries: a systematic review and meta-analysis protocol. Systematic reviews 2017, 6(1):5. 42. 42.Singh P, Hashmi G, Swain PK: High prevalence of cesarean section births in private sector health facilities-analysis of district level household survey-4 (DLHS-4) of India. BMC Public Health 2018, 18(1):613. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 43. 43.Moore B: Appropriate technology for birth. The Lancet 1985, 326(8458):787. 44. 44.Chaillet N, Dubé E, Dugas M, Francoeur D, Dubé J, Gagnon S, Poitras L, Dumont A: Identifying barriers and facilitators towards implementing guidelines to reduce caesarean section rates in Quebec. Bulletin of the World Health Organization 2007, 85:791–797. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2471/BLT.06.039289&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18038061&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250164200015&link_type=ISI) 45. 45.Opiyo N, Kingdon C, Oladapo OT, Souza JP, Vogel JP, Bonet M, Bucagu M, Portela A, McConville F, Downe S: Non-clinical interventions to reduce unnecessary caesarean sections: WHO recommendations. Bulletin of the World Health Organization 2020, 98(1):66. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 46. 46.Han K-T, Kim SJ, Ju YJ, Choi JW, Park E-C: Do hospital characteristics influence Cesarean delivery? Analysis of National Health Insurance claim data. The European Journal of Public Health 2017, 27(5):801–807. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 47. 47.Abebe FE, Gebeyehu AW, Kidane AN, Eyassu GA: Factors leading to cesarean section delivery at Felegehiwot referral hospital, Northwest Ethiopia: a retrospective record review. Reprod Health 2015, 13(1):6. 48. 48.Gutema H, Shimye A: cesarean section and associated factors at mizan aman general hospital, southwest Ethiopia. J Gynecol Obstet 2014, 2(3):37–41. 49. 49.D’orsi E, Chor D, Giffin K, Angulo-Tuesta A, Barbosa GP, Gama AdS, Reis AC: Factors associated with cesarean sections in a public hospital in Rio de Janeiro, Brazil. Cad Saude Publica 2006, 22:2067–2078. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16951878&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 50. 50.Patel RR, Peters TJ, Murphy DJ: Prenatal risk factors for Caesarean section. Analyses of the ALSPAC cohort of 12 944 women in England. Int J Epidemiol 2005, 34(2):353–367. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyh401&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15659468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228978900023&link_type=ISI) 51. 51.Kenny LC, Lavender T, McNamee R, O’Neill SM, Mills T, Khashan AS: Advanced maternal age and adverse pregnancy outcome: evidence from a large contemporary cohort. PLoS ONE 2013, 8(2):e56583. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0056583&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23437176&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 52. 52.Greenberg MB, Cheng YW, Sullivan M, Norton ME, Hopkins LM, Caughey AB: Does length of labor vary by maternal age? American journal of obstetrics and gynecology 2007, 197(4):428. e421-428. e427. 53. 53.Zaki MN, Hibbard JU, Kominiarek MA: Contemporary labor patterns and maternal age. Obstetrics and gynecology 2013, 122(5):1018. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/AOG.0b013e3182a9c92c&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24104787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000330448500013&link_type=ISI) 54. 54.Khan MN, Islam MM, Shariff AA, Alam MM, Rahman MM: Socio-demographic predictors and average annual rates of caesarean section in Bangladesh between 2004 and 2014. PLoS ONE 2017, 12(5):e0177579. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 55. 55.Lin H-C, Sheen T-C, Tang C-H, Kao S: Association between maternal age and the likelihood of a cesarean section: a population-based multivariate logistic regression analysis. Acta Obstet Gynecol Scand 2004, 83(12):1178–1183. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.0001-6349.2004.00506.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15548152&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 56. 56.Ecker JL, Chen KT, Cohen AP, Riley LE, Lieberman ES: Increased risk of cesarean delivery with advancing maternal age: indications and associated factors in nulliparous women. American journal of obstetrics and gynecology 2001, 185(4):883–887. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1067/mob.2001.117364&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11641671&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171734000021&link_type=ISI) 57. 57.Mazzoni A, Althabe F, Liu NH, Bonotti AM, Gibbons L, Sánchez AJ, Belizán JM: Women’s preference for caesarean section: a systematic review and meta-analysis of observational studies. BJOG: an international journal of obstetrics & gynaecology 2011, 118(4):391–399. 58. 58.Yisma E, Smithers LG, Lynch JW, Mol BW: Cesarean section in Ethiopia: prevalence and sociodemographic characteristics. The Journal of Maternal-Fetal & Neonatal Medicine 2019, 32(7):1130–1135. 59. 59. Nupur Shah, Nandita Maitra, Pagi SL: Evaluating role of parity in progress of labour and its outcome using modified WHO partograph. Int J Reprod Contracept Obstet Gynecol 2016, 5(3):860–863. 60. 60.Van Der Spek L, Sanglier S, Mabeya HM, Van Den Akker T, Mertens PL, Houweling TA: Socioeconomic differences in caesarean section–are they explained by medical need? An analysis of patient record data of a large Kenyan hospital. Intern 2020, 19(1):1–14. 61. 61. Anjel Vahratian, Matthew K. Hoffman, James F. Troendle, Zhang J: The impact of parity on course of labor in a contemporary population. Birth 2006, 33(1):12–17. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16499527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 62. 62.Arulkumaran S, Gibb D, Lun K, Heng S, Ratnam S: The effect of parity on uterine activity in labour. BJOG: An International Journal of Obstetrics & Gynaecology 1984, 91(9):843–848. 63. 63.World Health Organization: WHO recommendations on antenatal care for a positive pregnancy experience: WHO; 2016. 64. 64.Gamble JA, Creedy DK: Women’s preference for a cesarean section: incidence and associated factors. Birth 2001, 28(2):101–110. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1523-536X.2001.00101.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11380381&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168872000005&link_type=ISI) 65. 65.Long Q, Kempas T, Madede T, Klemetti R, Hemminki E: Caesarean section rates in Mozambique. BMC pregnancy and childbirth 2015, 15(1):253. 66. 66.Cavallaro FL, Cresswell JA, França GV, Victora CG, Barros AJ, Ronsmans C: Trends in caesarean delivery by country and wealth quintile: cross-sectional surveys in southern Asia and sub-Saharan Africa. Bulletin of the World Health Organization 2013, 91:914–922D. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.2471/BLT.13.117598&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24347730&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F10%2F10%2F2022.10.07.22280820.atom) 67. 67.Hounton SH, Newlands D, Meda N, De Brouwere V: A cost-effectiveness study of caesarean-section deliveries by clinical officers, general practitioners and obstetricians in Burkina Faso. Human Resources for Health 2009, 7(1):34. 68. 68.Fesseha N, Getachew A, Hiluf M, Gebrehiwot Y, Bailey P: A national review of cesarean delivery in Ethiopia. International Journal of Gynecology & Obstetrics 2011, 115(1):106–111. 69. 69.WHO, UNFPA, UNICEF, AMDD: Monitoring emergency obstetric care: a handbook. In. Geneva, Switzerland: WHO; 2009. [1]: /embed/graphic-1.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/inline-graphic-2.gif [4]: /embed/graphic-2.gif [5]: /embed/inline-graphic-3.gif