Background Globally, about 810 women die every day due to pregnancy and its related complications. Although the death of women during pregnancy or childbirth has declined from 342 deaths to 211 deaths per 100, 000 live births between 2000 and 2017, maternal mortality is still higher, particularly in sub-Saharan Africa and South Asia, where 86% of all deaths occur. Methods A secondary analysis was carried out using the 2014 Ghana Demographic and Health Survey. A sample total of 4, 290 women who had a live birth in the 5 years preceding the survey was included in the analysis. GIS software was used to explore the spatial distribution of unskilled birth attendance in Ghana. The Geographic Weighted Regression (GWR) was employed to model the spatial relationship of some predictor of unskilled birth attendance. Moreover, a multilevel binary logistic regression model was fitted to identify factors associated with unskilled birth attendance. Results In this study, unskilled birth attendance had spatial variations across the country. The hotspot, cluster and outlier analysis identified the concerned districts in the north-eastern part of Ghana. The GWR analysis identified different predictors of unskilled birth attendance across districts of Ghana. In the multilevel analysis, mothers with no education, no health insurance coverage, and mothers from households with lower wealth status had higher odds of unskilled birth attendance. Being multi and grand multiparous, perception of distance from the health facility as not a big problem, urban residence, women residing in communities with medium and higher poverty level had lower odds of unskilled birth attendance. Conclusion Unskilled birth attendance had spatial variations across the country. Areas with high levels of unskilled birth attendance had mothers who had no formal education, not health insured, mothers from poor households and communities, primiparous women, mothers from remote and border districts could get special attention in terms of allocation of resources including skilled human power, and improved access to health facilities.
The Ghana Demographic and Health Survey (GDHS) used a standard Demographic and Health Survey (DHS) model questionnaire developed by the Measure DHS programme. This study used the most recent DHS data and a cross-sectional study design. The DHS are national surveys carried out every five years in low-and middle-income countries globally. The surveys concentrate on maternal and child health, physical activity, sexually transmitted infections, fertility, health insurance, tobacco use, and alcohol consumption. They provide data to monitor the demographic and health profiles of the respective countries. For the study, women with birth history who had given birth up to five years before the survey were included. Only the last birth of the women aged 15–49 years preceding the survey was included in the study. A sample of 4,290 women with complete data required for our analysis participated in this study. Permission to use the data set was given by the MEASURE DHS following the assessment of our concept note. This study is a secondary analysis of the de-identified 2014 Ghana Demographic and Health Survey (GDHS), a publicly available dataset. Therefore, ethics approval or consent to participate is not applicable. The datasets are freely available to the public at www.measuredhs.com. Outcome variable. The primary outcome variable was unskilled birth attendance. The outcome variable was derived from the response to the question “who assisted with the delivery?” Responses were categorized under Health Personnel and Other Person. Health personnel included doctor, nurse, nurse/midwife, auxiliary midwife, and other people consisting of a traditional birth attendant, traditional health volunteer, community/village health volunteer, neighbours/ friends/relatives, etc. For this study, unskilled birth attendance referred to births assisted by a traditional birth attendant, traditional health volunteer, community/village health volunteer, neighbours/ friends/relatives, other [11]. Explanatory variables. Eleven explanatory variables were used. These were grouped into individual and community level variables. The individual characteristics consist of level of education, age, parity, health insurance coverage, wealth status, distance to a health facility, and media exposure. The community-level characteristics comprised the type of residence, community socioeconomic status, community literacy, and region of residence. We employed both descriptive and inferential analytical approaches. First, we computed the proportion of women who utilized the service of unskilled birth attendants during delivery. This ensued with bivariate analysis between individual characteristics (level of education, age, parity, health insurance coverage, wealth status, distance to a health facility, and media exposure), community characteristics (Type of residence, Community socioeconomic status, Community literacy, Region of residence) and utilization of unskilled birth attendants (see Table 1). Following the hierarchical nature of the data set, the Multilevel Logistic Regression Model (MLRM) was employed. This comprises fixed effects and random effects [23]. The fixed effects of the model were gauged with binary logistic regression, which resulted in odds ratios (ORs) and adjusted odds ratios (aORs) (see Table 2). Model 1 was an empty table, where model 2 looked at the relationship between the individual variable and the outcome variable. Model 3 looked at the relationship between community variables and the outcome variable. Model 4 was the complete model that looked at the relationship with both the individual and community variables and the outcome variable. The random effects, on the other hand, were assessed with Intra-Cluster Correlation (ICC) [23] (Table 3). The sample weight (v005/1,000,000) was applied in all the analyses to control for over and under-sampling. All the analyses were carried out using STATA version 14. Source: GDHS, 2014 Source: GDHS, 2014 *p<0.05 **p<0.01 *** p<0.001 Source: GDHS, 2014 We assessed the fitness of all the models with the Likelihood Ratio (LR) test. The presence of multicollinearity between the independent variables was checked before fitting the models. The variance inflation factor (VIF) test revealed the absence of high multicollinearity between the variables (Mean VIF = 2.28). In the conduct of the survey, instead of mapping outhouses in which the data were collected, clusters were mapped to protect the actual identity and location of respondents [24]. These clusters are developed to suit the district-level data, making it easy to merge the household records with spatial data. During the data collection period, there were 216 administrative districts in Ghana; however, not all districts had respondents drawn from for the survey. This aided in the merger of the data gathered with the district shapefiles obtained from the Department of Geography and Regional Planning, University of Cape Coast, Ghana. This was done to permit the analysis to be made at a district level. The data is best analysed at the district level since the information is more representative at the cluster level [24]. This study extracted the required variables from the 2014 GDHS. The extracted data maintained the mapped clusters information. This mapped cluster information was used to help join the extracted non-spatial data to the coordinates gathered for the clusters. All the data required (GDHS data and 216 district boundary) were projected into the projected coordinate system of Ghana Meter Grid to aid in the spatial analysis. The extracted GDHS data were merged with coordinate, and a spatial join was undertaken to transfer the cluster point to the 216-district boundary (polygon) layer using ArcMap version 10.5. This activity enabled us to easily identify and trace where each case is located within a district. It was identified that some of the district boundaries had more than one cluster. In such cases, the data from the clusters were aggregated, and their means were computed to represent the respective district they fell within [24]. With regards to the geospatial analyses, four spatial statistical tools were applied to analyse the data. These tools were spatial autocorrelation (Global Moran’s I), hot spot analysis (Getis-Ord G), outlier and cluster analysis, and Geographically Weighted Regression. The spatial autocorrelation was used to assess whether unskilled birth attendance in Ghana had a clustering or dispersion pattern at the district level. This study hypothesized that unskilled birth attendance is randomly distributed across various districts in the country. The null hypothesis is rejected if a calculated p-value is small (95% confidence interval), which implies an unlikely situation that the observed spatial pattern results from random processes [24]. Further, hot spot analysis (Getis-Ord G) was used to ascertain statistically significant spatial variations in unskilled birth attendance [24, 25]. This analysis was conducted to determine districts with high prevalence against areas of the low prevalence of unskilled birth attendance. In addition, an outlier and cluster analysis was conducted to identify districts that appeared as outliers. Outlier districts could either be a hot spot district that is surrounded by cold spot districts and vice-versa. The geographically weighted regression (GWR) modelling was conducted after ascertaining the hot spot and cluster and outlier analysis of unskilled birth attendance, the geographically weighted regression (GWR) modelling was conducted. This spatial regression modelling was performed to identify which explanatory variables best account for the observed spatial patterns of unskilled birth attendance [25]. To be specific, the GWR uses the OLS coefficient from the clusters concerning its nearest neighbours in modelling the predictability of the explanatory variable. The output shows how the strength of each explanatory variable changed across space. Therefore, maps of the statistically significant coefficients were generated.