Introduction
Despite the progress made in reducing the global under-five mortality rate by 59% from 93 deaths per 1000 live births in 1990 to 38 deaths per 1000 live births in 2019, millions of children continue to die yearly in resource-limited countries (UNICEF, 2021). Among the 5.2 million under-five deaths recorded in 2019, 2.4 million deaths (47%) occurred during the neonatal period (first 28 days of life), making the neonatal period the time of highest mortality risk for children (UNICEF, 2021). With 27 deaths per 1000 live births, neonatal mortality rates (NMR) remain highest in the region of sub-Saharan Africa (UNICEF, 2021). Moreover, the inter-country disparity in neonatal mortality is more profound for Nigeria, as it ranks second after India in the global neonatal mortality burden league table (UNICEF, 2021). In Nigeria, 1 in 28 newborns did not survive the first 28 days of life, resulting in an estimated 270,000 neonatal deaths, which accounted for 11% of the global burden of neonatal deaths in 2019 (UNICEF, 2021). The risk of neonatal death in Nigeria (36 deaths per 1000 live births) is twice the global NMR of 17 deaths per 1000 live births (UNICEF, 2021).
Although recent reviews of health indicators in Nigeria suggest moderate improvement in the survival of under-five children, NMR has stagnated in recent years (Akinyemi, Bamgboye, and Ayeni, Reference Akinyemi, Bamgboye and Ayeni2015; Morakinyo and Fagbamigbe, Reference Morakinyo and Fagbamigbe2017; Ayoade, Reference Ayoade2018). A further concern is the masking of health inequities at the subnational level by aggregating the country’s performance at the national level (Simpson’s paradox), when tracking the child health-related Sustainable Development Goals (SDG). Understanding the geographical heterogeneity of mortality during the neonatal period – the most crucial period of early child’s development (Ezeh et al., Reference Ezeh, Agho, Dibley, Hall and Page2015; Naline and Viswanathan, Reference Naline and Viswanathan2017), within Nigeria, can provide additional information needed for local-level planning and allocation of resources to areas where they are needed most (underserved population).
Previous studies have addressed social determinants of neonatal mortality (Akinyemi, Bamgboye, and Ayeni, Reference Akinyemi, Bamgboye and Ayeni2015; Kayode et al., Reference Kayode, Grobbee, Amoakoh-Coleman, Ansah, Uthman and Klipstein-Grobusch2017; Morakinyo and Fagbamigbe, Reference Morakinyo and Fagbamigbe2017; Neal, Channon, and Chintsanya, Reference Neal, Channon and Chintsanya2018), but more remains to be known about the spatial variations of neonatal mortality in Nigeria. Also, there is a dearth of information on gender differences in NMR across urban–rural areas and geographical regions. According to Akinyemi et al. (Reference Akinyemi, Bamgboye and Ayeni2015), NMR varied by regions and was highest in the northern region of Nigeria from 1990 to 2013. The regional variations have been traced to the differences in socioeconomic, cultural, and environmental factors (Akinyemi, Bamgboye, and Ayeni, Reference Akinyemi, Bamgboye and Ayeni2015). Drawing on the World Health Organization (WHO) Commission of Social Determinants of Health (SDH), gender and spatial dimensions are important determinants of population health (World Health Organization, 2008).
With respect to the roles of gender (a social and cultural construct) and sex (biological identity) on child survival, epidemiological studies have provided contradictory findings. While some studies have linked worse survival outcomes among male children to biological disadvantages (Boco, Reference Boco2014; Gebretsadik and Gabreyohannes, Reference Gebretsadik and Gabreyohannes2016; Morakinyo and Fagbamigbe, Reference Morakinyo and Fagbamigbe2017), others have reported excess of girl-child mortality due to gender discrimination, especially in terms of selective termination of female fetuses and newborns, and neglect of nutrition and health care for the girl-child (Costa, da Silva, and Victora, Reference Costa, da Silva and Victora2017). Also, evidence suggests that the population residing in socioeconomically disadvantaged areas such as rural residence experience worse health outcomes because of high levels of poverty, inaccessibility to quality health care, and inadequate social infrastructure (McMichael, Reference McMichael2000; Morakinyo and Fagbamigbe, Reference Morakinyo and Fagbamigbe2017). In contrast, some studies have noted urban area disadvantage for under-five mortality due to air pollution, overpopulation, and waste disposal crisis (Van de Poel, O’Donnell, and Van Doorslaer, Reference Van de Poel, O’Donnell and Van Doorslaer2007; Antai and Moradi, Reference Antai and Moradi2010; Kimani-Murage et al., Reference Kimani-Murage, Fotso, Egondi, Abuya, Elungata, Ziraba, Kabiru and Madise2014).
It is fundamental that policymakers should address the gender bias and rural–urban disparity in child survival and social determinants of neonatal mortality across the states and regions in Nigeria. In this study, the primary objective was to determine the patterns and determinants of geographical clustering of neonatal mortality in the state and geopolitical zones in Nigeria. The secondary objective was to assess gender inequity in neonatal mortality between urban and rural communities across the zones in Nigeria.
Methods
Study area
Nigeria, the most populated country in sub-Saharan Africa, is located in West Africa. It comprises six geopolitical zones (i.e. North-West [NW], North-East [NE], North-Central [NC], South-West [SW], South-East [SE], and South-South [SS]), which are further divided into 36 states and Federal Capital Territory (FCT) (Figure 1). There are more than 250 ethnic groups, which are divided into three major ethnic groups. Predominantly, there are Yoruba in the SW, Igbo in the SE, and Hausa in the Northern Nigeria.
States and Zones in Nigeria (National Bureau of Statistics, 2013).
Study design and data sources
This is a cross-sectional study that utilised full birth history, along with maternal and household data files, obtained from the 2016/2017 Nigeria Multiple Indicator Cluster Survey (MICS) (UNICEF MICS, 2018). The MICS is a national population-based survey conducted by trained interviewers in Nigeria, with support from UNICEF Headquarters, New York, to provide estimates of maternal and child health indices for the country (UNICEF and NBS, 2017). A detailed explanation of the methodology used for MICS has been described in the full report (UNICEF and NBS, 2017). With a complex, multi-stage stratified cluster sampling technique, data were collected from 36,176 women aged 15–49 years between September 2016 and January 2017 in the 36 states and FCT of Nigeria (UNICEF and NBS, 2017). The 37 states (including FCT) in the six geopolitical zones corresponded to the sampling strata. The strata were then sorted into rural and urban areas. Overall, 33,901 households from 2239 enumeration areas, otherwise referred to as primary sampling units (PSU), were covered during fieldwork. The PSU, which was used as a proxy for the community, was defined as an administrative/enumeration area with homogeneous population characteristics. Of the 36,176 women selected for interview, the response rate was 95% (34,376 women). To minimise recall bias and accurately estimate neonatal mortality, this study considered children born alive within the last five years preceding the survey commencement (i.e. September 2011–September 2016) and defined neonatal death as a death that occurred within 28 days of birth. The cohort was selected to ensure that the analysis exclusively included neonates delivered in recent years. After removing the data of children whose survival outcomes, dates of birth, and deaths were not documented, an analysis was conducted on 29,786 neonates (corresponding to a weighted sample of 30,924 neonates) delivered to 18,497 women.
Theoretical background
This study utilised the Mosley–Chen framework (Mosley and Chen, Reference Mosley and Chen1984), programmatic experience of the authors, and evidence from the literature to identify the relevant social determinants of neonatal deaths in Nigeria. The framework underscores that childhood mortality in resource-limited countries results from the complex interrelationships of multiple biological and social factors at the child, maternal, household, and community levels. In this study, the hypothesis is that substantial variations in neonatal mortality exist across rural/urban communities, states, and geopolitical zones of Nigeria. Also, the geographical and gender inequities in neonatal survival are expected to vary based on the impact of other social determinants of health on the child–mother dyads. In this study, the gender roles depicted by boys and girls are used to show how people are viewed and expected to act based on societal norms, not solely due to biology (i.e. males and females) but also because of environmental factors and upbringing.
Variables
Dependent variable
The outcome variable – neonatal survival status – was generated from information on child survival outcome, age at death, and current age of living children and divided into two categories: alive (coded as 0) and dead (coded as 1).
Independent variables
The selected independent variables were informed by the Mosley–Chen framework, programmatic experience of authors, evidence from literature, and availability of variables in the MICS dataset. The variables were layered across child, maternal, household, and community levels (see Table 1 for details). From the variables collected in the MICS dataset, the housing condition index was generated by applying a principal components analysis (PCA) to reduce variables on the quality of the roof, exterior wall, and floor; overall, Kaiser–Meyer–Olkin measure of adequacy was 0.7, and p-value (Bartlett’s test of sphericity) <0.001. The first component was selected based on an eigenvalue of 2.03, explaining 67.7% of the total variance. With the median value as the cut-off point, the housing condition index variable was coded as adequate and inadequate. Also, maternal media exposure, access to drinking water, sanitation, and indoor pollution were generated. Exposure to media was defined as the frequency that mothers were exposed to at least a source of mass media – newspaper/magazine, radio, and television (1: almost every day [high exposure]; 2: at least once a week [moderate exposure]; 3: less than once a week [low exposure]; 4: not at all). Household sanitation, sources of drinking water, and cooking fuel (proxy for indoor pollution) were re-categorised into improved and unimproved, as defined by the WHO/United Nations Children’s Fund (UNICEF) Joint Monitoring Programme for Water Supply, Sanitation and Hygiene (WHO, UNICEF, and JMP, 2018). Furthermore, community-contextual variables were derived from maternal-level variables (community level of maternal education) and household-level variables (community infrastructural development) (see Table 1 for details).
Sociodemographic Characteristics of Study Participants, 2016/2017 MICS, Nigeria
In union refers to a boarder range of relationships including both married and unmarried partnerships (e.g. cohabitation).
a Data available for only women with a live birth in the two years prior to the survey. ANC: antenatal care.
Statistical analyses
This study employed artificial intelligence technique – backpropagation feedforward multilayer perceptron (MLP) neural network and geospatial analyses. Descriptive statistics were initially generated for all variables using StataTM version 15.1 software (College Station, Texas) (Stata version 15.1, 2017). The bivariate associations between the outcome variable (neonatal deaths) and independent variables were assessed with the Chi-square test. The significance level was set at two-tailed α=5%. The complex survey design commands in StataTM software and sampling weights were applied to account for the hierarchical sampling and unequal selection probabilities of samples. Early NMR (risk of death from birth to 6 days of life) and late NMR (risk of death from 7 days to 27 days after birth) were also computed. Furthermore, strip plots were used to visualise the state distributions of early and late mortality rates based on their zones. The gender inequity gaps in NMR, disaggregated by urban–rural residence across the six geographical zones, were visualised with equiplots.
Geospatial analysis
Initially, the spatial dependencies of NMR were assessed across the states in Nigeria by using geometric centroids and a spatial arc distance of 407 km in a Nigeria shapefile that was obtained from the United Nations Office (United Nations Office for the Coordination of Humanitarian Affairs, 2017). A distance-based weight matrix was generated to ensure that all the states were interconnected – a condition for spatial analysis (Anselin, Reference Anselin2005). The symmetry of connectivity histogram and connectivity map were used to assess the suitability of the spatial weights. The units of spatial analysis were states. As proposed by Anselin (Reference Anselin2005), univariate global spatial autocorrelation was performed, and the degree of similarity (i.e. spatial clustering) of NMR across the states was further assessed. A global Moran’s I index of a positive value indicates spatial clustering, while a negative value indicates spatial dispersion (Anselin, Reference Anselin2005). Also, univariate local indicator spatial autocorrelation (LISA) cluster and significance maps were generated to identify the states with high NMR as denoted by hot spots, and the states with low NMR (cold spots). The statistically significant hot spots were determined by the grouping of states with high NMRs and vice versa for cold spots.
In the second step, the artificial intelligence technique, specifically the MLP neural network was employed to identify the important predictors of neonatal mortality, stratified by zones. Predictors with ≥50% normalised importance were deemed major contributors to NMR. For details of MLP, see ‘Artificial neural network’, below).
The final step involved determining spatial dependencies between the key predictors identified in step 2 and NMR by generating bivariate LISA cluster maps. To visualise the combined effects of the variables that were spatially autocorrelated with NMR from the bivariate LISA maps (i.e. proportion of mothers who had children with previous birth interval <2 years, birth order >3, young maternal age at birth [<20 years] and no maternal exposure to mass media), PCA was used to reduce these variables to a composite variable referred to as socio-behavioural index. High values of the composite index indicate poor social behaviours (i.e. higher position in a family birth order, births closer together, young maternal age at birth, no maternal access to mass media). With the singular value decomposition (SVD) method and z-score transformation (i.e. mean of zero and variance of one), the first component explained 74.9% of the total variance, and its eigenvalue was 3.0. The most important advantage of SVD is its robustness against outliers (Anselin, Reference Anselin2020). As a result, the combined spatial pattern of the variables that constituted the first principal component of NMR was visualised with a cluster map. The statistical significance of spatial autocorrelations was tested by running 999 Monte Carlo simulations with a p-value<0.05. The spatial analysis was implemented in GeoDa software version 1.14 (GeoDa on Github, no date).
Artificial neural network
With a view of identifying the key predictors of neonatal mortality at zonal and national levels, backpropagation feedforward MLP neural networks were implemented in the IBM SPSS neural networks software version 21.0 (IBM, 2012). The unit of analysis for the neural network was at the level of children, stratified by zones. Figure 2 shows the flow chart of data preprocessing for the MLP neural network (see Appendix for details). With the dataset randomly partitioned into training (70%), testing (20%), and holdout (10%) sets, the MLP learning algorithm (gradient descent algorithm) used the training set to learn the inherent pattern of the dataset. The testing and holdout sets were used for model validation. The MLP neural networks were validated based on values of the error function (cross-entropy), area under receiver operating characteristic curve (AUROC), and accuracy rate from holdout samples (IBM Corporation, 2012). The missing data were represented by a specific number (99) and labelled as missing before running the analyses. This approach was undertaken to ensure that all available information was considered. The missing categories were not reported in the results section to avoid ambiguity.
Flow Chart of Multilayer Perceptron Neural Network Modelling.
Ethical considerations
Ethical clearances were obtained earlier by the UNICEF MICS team from the National Health Research Ethics Committee, Nigeria, before the survey commencement. In addition, this present study was exempted from ethical review by the University of Saskatchewan Behavioural Ethics Committee (ID no. 904) as datasets were de-identified of the respondents’ personal information. The participants’ anonymity and confidentiality are assured.
Results
Table 1 shows the sociodemographic characteristics of the study participants. Of the weighted sample of 30,924 live births included in the study, 50.8% were boys, 96.0% were single births, and 60.9% were children who had preceding birth interval ≥2 years.
The NMR was 37.9 deaths per 1000 live births. Most of the neonatal deaths (85.9%) occurred within the first 7 days of life, translating to 32.3 deaths per 1000 live births (early NMR), and declined to 5.7 deaths per 1000 live births during the late neonatal period (7–27 days). Although not statistically significant, the NMR was slightly higher in rural areas (39.3 deaths per 1000 live births) than in urban areas (34.7 deaths per 1000 live births); p=0.262. Higher NMR was observed for boys (43.2 deaths per 1000 live births), compared to girls (32.5 deaths per 1000 live births); p<0.001.
Geographical variations of neonatal mortality rates
Figures 3 and 4 indicate that there were wide variations in NMR across the states and zones in Nigeria. NMR was highest in the NW zone (44.1 deaths per 1000 live births) and lowest in the SS zone (22.1 deaths per 1000 live births). Specifically, Kano State had the highest NMR (67.4 deaths per 1000 live births), followed by Niger State (58.2 deaths per 1000 live births). The lowest NMR was observed in Edo State (7.9 deaths per 1000 live births). The distribution of early and late NMR for states based on their regions is presented in Figure 4a and b, respectively.
A. Spatial Distribution of NMR, 2016/2017 Nigeria MICS B. Univariate LISA cluster map for NMR, 2016/2017 Nigeria MICS.
State Level A. Early Neonatal Mortality Rates, 2016/2017 Nigeria MICS B. Late Neonatal Mortality Rates, 2016/2017 Nigeria MICS.
As shown in Figure 3b, significant clustering of NMR was observed across the states in Nigeria (global Moran’s I index=0.1, p=0.02). The univariate LISA cluster map for NMR further revealed that 13.5% of the states were high-high clusters (hot spots) and located majorly in the NW zone (Jigawa, Katsina, Kano, and Zamfara States) and to a lesser extent in the NC zone (Niger State). The high-high clusters depict the states with statistically significantly higher NMR than the national average. Contrarily, low-low clusters (cold spots) were the states with statistically significantly lower NMR than the national average. The cold spots (24.3%) were concentrated in the SS zone (Akwa Ibom, Bayelsa, Cross River, Delta, and Rivers States), and SE zone (Abia, Anambra, Ebonyi, and Enugu States). However, there were some outliers – high-low and low-high clusters. The high-low clusters were formed by significant clustering of two states (5.4%), where NMRs were high and adjacent to states with low NMRs. The high-low clusters were formed by Benue State (NC zone) and Imo State (SE zone). Two states (5.4%) were identified as low-high clusters in the NW zone (Kaduna and Sokoto States) – states with low NMRs and were neighbours to states with high NMRs. The spatial correlogram that shows the changes in spatial autocorrelation of NMR with distance is shown in Figure A1 in the Appendix. Also, the spatial patterns of early and late NMRs are shown in Figures A2 and A3 in the Appendix).
Magnitude of gender inequity in neonatal mortality rate across urban–rural residence
As shown in Figure 5, mortality rates were highest among boys residing in the rural NW zone (54.2 deaths per 1000 live births) and lowest among boys in the urban SS zone (10.8 deaths per 1000 live births). For girls, the highest NMR was observed in urban NC (39.9 deaths per 1000 live births) and lowest in rural SS (18 deaths per 1000 live births).
Absolute Gender Inequity in Neonatal Mortality Rate Across the Urban–Rural Residence of Geographical Regions, 2016/2017 Nigeria MICS.
Figure 5 shows the absolute inequity (i.e. risk difference) between boys and girls disaggregated by urban–rural residence across the geographical zones. The absolute difference was largest in urban SS and urban SW zones (20 deaths per 1000 live births) and lowest in urban NC (0.4 deaths per 1000 live births). Overall, the gender differences in NMR tended to be larger in the rural North; however, the urban South was observed to have large differences. Except for NW and SW zones, NMR among girls was higher in urban areas. However, NMR among boys was generally higher in rural areas in all zones.
Determinants of neonatal mortality across the zones
From the predictive MLP neural nets, there were zone-specific determinants of neonatal mortality in Nigeria (Table 2). Overall, based on the normalised importance values, multiple births (100%), previous birth interval (53.9%), and birth order (51.3%) were identified as the major contributors to neonatal mortality in Nigeria. The zonal-level analysis also found similar evidence that these three factors were common contributors across all the zones. Except for the NW zone, maternal age at birth appeared consistently across the zones (Table 2). Also, maternal mass media exposure was observed as a top contributor to neonatal mortality in the southern part of the country (Table 2).
Zonal Comparison of NMR, Progress Towards SDG Targets, Top Contributors to Neonatal Mortality and Model Performance
a Current deficit in SDG 3 target (%): difference, in percentage, the SDG target, and the estimated rate for NMR as of 2016/2017.
b Normalised importance: equivalent to regression coefficient. This is the neural network classification of independent variables based on their strength of association with the outcome variable.
c Model accuracy: the predictive performance of the trained neural network on previously unseen datasets (i.e. validation data sets).
To establish spatial relationships of the most important determinants identified from MLP – that is, multiple births, previous birth interval, birth order, maternal age at birth, and maternal access to mass media – bivariate LISA cluster maps were generated (Figure 6 and Appendix Figure A4). There was no statistically significant global spatial association between multiple births and spatial lag of NMR in Nigeria (Moran’s I index=-0.1, p=0.05). The local spatial association between multiple births and NMR is presented in Appendix Figure A4. However, births closer together (less than two years gap) and increasing birth order (>3) were significantly associated with the spatial clustering of NMR: Moran’s I index=0.1, p=0.01 and Moran’s I index=0.2, p=0.003, respectively. Also, Moran’s I index for NMR and deliveries by adolescent mothers indicates significant clustering (I=0.21, p=0.001) and no maternal exposure to mass media (I=0.1, p=0.008).
Bivariate LISA Cluster Maps A. NMR and Previous Birth Interval <2 Years B. NMR and Birth Order >3. C. NMR and Deliveries by Adolescent Mothers D. NMR and No Maternal Exposure to Mass Media.
Figure 6a suggests that there were 5(13.5%) states that formed high-high clusters, implying the clustering of states with high NMR and high percentage of children who were born less than two years apart. The high-high clusters were in the NE zone (Bauchi and Yobe States) and NW zone (Jigawa, Kano, and Katsina States). Also, SE zone (Abia and Anambra States), SS zone (Bayelsa, Delta, and Edo States), NC zone (Kogi and Kwara States), and SW zone (Lagos, Ogun, and Ondo States) clearly indicate clustering of states with significantly low NMRs and low percentage of children who were born less than two years apart (low-low clusters). Figure 6b–d displays bivariate LISA cluster maps for the association between NMR and increasing birth order (>3), deliveries by adolescent mothers, and no maternal exposure to mass media.
The multivariate cluster map also shows evidence of significant positive spatial autocorrelation (I=0.2, p=0.002) between increasing birth order, births closer together, young maternal age at birth, no maternal access to mass media, and NMR in Nigeria (Figure 7). The multivariate map shows that high NMR was spatially correlated with the high values of the derived socio-behavioural index (hot spots) in NE (Bauchi and Yobe), NW (Jigawa, Kano, Katsina, Kebbi, and Zamfara), and NC (Plateau). The SE (Abia, Anambra, Ebonyi, Enugu), SS (Akwa Ibom, Bayelsa, Cross River, Delta, Edo, Rivers), NC (Kogi, Kwara), and SW (Lagos, Ogun, Ondo) were identified as the cold spots – low NMR was spatially correlated with low socio-behavioural index. The high-low clusters were found in NC (Benue), SW (Ekiti, Osun, Oyo), and SE (Imo). The high-low clusters were the states with significant spatial correlation between high NMR and low socio-behavioural index. However, the low-high clusters indicated states with significantly low NMR and high socio-behavioural index. The low-high clusters were formed by NE (Adamawa, Borno, Gombe) and NW (Kaduna and Sokoto).
Multivariate (PCA) LISA Cluster Maps of Neonatal Mortality Rate and Socio-Behavioural Index (Previous Birth Interval <2 Years, Birth Order >3, Deliveries by Adolescent Mothers, No Maternal Exposure to Mass Media).
NMR: Neonatal mortality rate; PCA: Principal components analysis.
Model assessment
The neural network models produced AUROC ranging from 0.69 to 0.95 (Table 2). The AUROC for the overall (national) model was 0.71, translating to an accuracy rate of 96.0%. With similar cross-entropy errors and accuracy rates across the training, testing, and holdout samples, there is evidence to suggest that the neural networks were not overtrained.
Discussion
In this study, NMR was 37.9 deaths per 1000 live births in Nigeria, and most of the deaths (85%) occurred within the first week of life. This is an indication of a missed opportunity for improving child survival in the early phase of a child’s development in Nigeria. Among the zones and states, NC (Benue, Nasarawa, and Niger), NE (Bauchi and Yobe), NW (Kano, Kebbi, and Zamfara), and SW (Ekiti, Osun, and Oyo) were observed to have higher NMRs than the national rate. More importantly, Kano State contributed 15.5% to all neonatal deaths in the country. Also, the findings show that children delivered as part of multiple births, and those with high socio-behavioural index (i.e. delivered later in birth order (>3), born within two years after preceding births, delivered to adolescent mothers, and those without access to mass media) had a higher risk of neonatal death in Nigeria. However, there was no clear evidence of spatial dependence of multiple births on the geographic pattern of NMR. Also, there were interregional disparities in the absolute gender-specific NMR between urban and rural areas. This study clearly shows the need to implement targeted interventions to reduce the gender gap between urban and rural residences across the geographical zones in Nigeria. Gender inequity was worse in the rural areas of northern Nigeria, while it was worse in the urban areas of southern Nigeria. NMR was disproportionately higher among girls in urban areas (except NW and SW zones). Conversely, boys had higher mortality risks in the rural areas for all the zones. The zones with the most equitable NMR were the urban NC region (0.4), followed by urban NW (3.3) and rural SE (3.4).
The geospatial analyses suggest that there was a huge disparity in neonatal mortality within Nigeria. The states in the NW and NC zones had higher NMR and clustered to form hot spots of neonatal mortality. However, a distinct pattern of lower NMR was observed towards the southern part (i.e. SS and SE) – cold spots. The finding that NW and NC had higher NMR is consistent with the pattern reported in the 2016/2017 Nigeria MICS report (UNICEF MICS, 2018). Previous surveys conducted before the Boko Haram insurgency started in Nigeria (i.e. 2009) showed that NMR was worst in the NE zone (National Population Commission (NPC) [Nigeria] and ORC, 2004; National Population Commission (NPC) [Nigeria] and ORC Macro, 2010). It is important to reinforce the consequences of conflicts on maternal and child health in the affected region (Omole, Welye, and Abimbola, Reference Omole, Welye and Abimbola2015; Howell et al., Reference Howell, Waidmann, Birdsall, Holla and Jiang2020). In line with the ongoing insurgencies in the northern region (Brechenmacher, Reference Brechenmacher2019), hot spots formed in the NW and NC zones might be in part due to neglect of the internally displaced people and a decline in the quality of health services in that area (especially for vulnerable women). According to a report by the Displacement Tracking Matrix, NW and NC zones have also been affected by the protracted humanitarian crisis, made worse due to inter-community violence and banditry (Global Data Institute Displacement Tracking Matrix, 2020). The recommendation is for the Government of Nigeria to prioritise maternal and child health services among vulnerable populations and to develop innovations aimed at improving child health outcomes in the northern region.
In the southern part of Nigeria, maternal mass media exposure was observed to be a social determinant of neonatal survival. These findings at least hint that mass media exposure among women contributed to a reduction of NMR in southern Nigeria. This implies that more policies should be shifted towards implementing context-specific strategies in the states and zones. Culturally appropriate reproductive, maternal, and child health messages targeted to women in the hot spots (northern Nigeria) via mass media campaigns are expected to increase maternal demands for quality preventive and curative services for children. In the same manner, the states that reported better exposure to mass media by women were found to have lower percentages of children who were delivered within two years before the previous birth. Notably, previous birth interval and birth order were observed to be consistent across the zones. Although recent evidence indicates that short previous birth interval is a major driver of childhood mortality, the conclusions have been mixed (Zhu et al., Reference Zhu, Rolfs, Nangle and Horan1999; Conde-Agudelo and Belizán, Reference Conde-Agudelo and Belizán2000; Kwarteng Acheampong and Eyram Avorgbedor, Reference Kwarteng Acheampong and Eyram Avorgbedor2017). However, most studies emphasise ‘maternal depletion hypothesis’ for the elevated childhood mortality risks arising from short preceding birth intervals (Conde-Agudelo and Belizán, Reference Conde-Agudelo and Belizán2000; Davanzo et al., Reference Davanzo, Hale, Razzaque and Rahman2008). In line with the idea of strong competition among siblings, this study shows that children born later in birth order were likely to experience neonatal mortality.
The impact of young maternal age at birth on neonatal mortality differed across the zones. In 21.6% of the states, mortalities were markedly elevated among states with a high percentage of neonates delivered by adolescent mothers. These states were located in northern Nigeria – Bauchi, Jigawa, Kano, Katsina, Kebbi, Plateau, Yobe, and Zamfara. Further analysis revealed that in these eight states that formed clusters of high adolescent mothers and high NMR, neonatal deaths were highest (83.9%) among women who were first married or in union during adolescence stage. Contrarily, NMR was low in 37.8% of the states with a low percentage of neonates delivered by adolescent mothers – located in southern Nigeria – Abia, Akwa Ibom, Anambra, Bayelsa, Cross River, Delta, Ebonyi, Edo, Enugu, Kogi, Lagos, Ogun, Ondo, and Rivers. Child marriage is a serious violation of human rights that has been closely linked to maternal and child mortality (UNICEF, 2021). In light of the reported findings of UNICEF (2021), child marriage affects about 650 million girls and women globally. According to UNICEF, child marriage is prevalent (44%) in the NW and NC regions of Nigeria – ranking 11th globally (Girls Not Brides, 2018; UNICEF, 2018a). Evidence suggests that the key drivers of child marriages in Nigeria include inadequate girl-child education, political/economic ties, gender norms, poverty, and violence against girls (Girls Not Brides, 2018). These findings highlight the need for stakeholders to sensitise the communities in Nigeria on the 2003 Child’s Rights Act (Nigeria: Act No. 26 of 2003, Child’s Rights Act, 2003, 2003) which prohibits forced and child marriages.
What is more striking from this study is that the determinants of neonatal mortality were not uniform across the six geographical zones. Across the zones, the most important contributors to NMR were previous birth interval (NC), birth order (NE), multiple births (NW and SE), polygamy family (SS), and maternal exposure to mass media (SW). This finding indicates that both broad and targeted strategies may be necessary to alleviate the NMR in Nigeria. Nigeria requires impactful policy actions to address the social determinants of neonatal mortality because of the gender, urban–rural, state, and zonal differences in the patterns of neonatal mortalities. To achieve this overarching goal will require the engagement of community members, decision and policymakers, and research institutions. More prominently, the states and zones are at different levels of progress towards achieving the SDG targets.
To the best of the authors’ knowledge, this is the first known published literature that utilised spatial analysis and artificial neural networks to cast new light on the urban–rural, state, and zonal variations of social determinants of neonatal mortality in Nigeria. This study is in line with the aspirations of SDG 3 (good health and well-being of children), SDG 5 (promote gender equality), SDG 10 (reduce inter[intra]-country inequality), and SDG 17 (increase the availability of high-quality, timely, and reliable data disaggregated by social determinants of health). Despite the strengths, some limitations exist. Owing to the cross-sectional design of this study, causal arguments should not be made. Rather, the findings should be considered within the context of associations. Another major source of limitation is recall bias/poor memory which could be due to self-reporting by women. This is not likely to markedly affect the study findings because data were limited to live births five years prior to the survey commencement.
Overall, this study found a considerable spatial clustering of NMR in Nigeria (majorly driven by young maternal age at birth, short birth interval, low maternal exposure to mass media, and increasing birth order). This study also found that, with the exception of NW and SW zones, NMR was higher for girls than boys in urban areas. However, NMR was higher for boys than girls in rural areas across all zones. This highlights the need for the country to develop and implement state and region-specific child survival initiatives to address the high rates of neonatal mortality, especially in northern Nigeria (i.e. hot-spot zones). Due to the secondary and quantitative nature of this study, it is challenging to offer explanations for the observed patterns. Interpreting the variation in the impact of social determinants of health across different regions in Nigeria requires a nuanced understanding of the local context, such as cultural, economic, and healthcare factors. Further qualitative studies are therefore needed to explain why some states were in clusters of outliers (high-low and low-high) and reasons for varying patterns of gender inequity across the rural and urban areas.
Data availability statement
Data may be found on the UNICEF MICS website (UNICEF MICS 2018).
Acknowledgements
This study was made possible by doctoral support for DAA through the College of Medicine Graduate Award from the University of Saskatchewan, Canada. The content is solely the responsibility of the authors and does not necessarily reflect the official views of the College of Medicine, University of Saskatchewan, Canada.
Author contribution
DAA conceived the study, analysed and interpreted the data, and wrote the first draft of the paper. NM assisted in the design and data interpretation and critically reviewed the manuscript. NM supervised this study. All authors read and approved the final manuscript.
Funding statement
The authors received no specific funding for this work.
Competing interests
No potential conflict of interest was reported by the authors.
APPENDIX
Artificial neural network
With a view of identifying the key predictors of neonatal mortality at regional and national levels, backpropagation feedforward MLP neural networks were implemented in the IBM SPSS neural networks software version 21.0 (IBM, 2012). Much like the human brain, the architecture of MLP is a collection of several artificial neurones that are connected by their weights and assembled into an input layer, at least one hidden layer, and an output layer (IBM Corporation, 2012). The input layer receives the inputs and performs the calculations via the neurones before onward transmission to the hidden layer. The hidden layer is the ‘black box’ connecting the input and output layers. The output layer receives information from the hidden layer to produce the final results (IBM Corporation, 2012). In addition to the linear activation function of MLP, it can also account for the nonlinear relationship between the inputs and outputs, hence producing more accurate results, compared to the traditional statistical methods (IBM Corporation, 2012).
Data preprocessing
With the dataset randomly partitioned into training (70%), testing (20%), and holdout (10%) sets, the MLP learning algorithm (gradient descent algorithm) used the training set to learn the inherent pattern of the dataset. The testing and holdout sets were used for model validation. The model-building process in this study used two hidden layers. The input layer factors were the independent variables (stratified by zones). The output layer consisted of one factor – neonatal survival status (dependent variable). After eliminating the bias units, the input layer had 74 units, while the output layer had 2 units. The number of units in the hidden layers was automatically computed. The number of units in the hidden layer 1 and hidden layer 2 were 13 and 10, respectively. The initial process generated a weighted sum and bias of input and hidden layers. The weighted sum and bias were activated through the hyperbolic tangent (tanh) activation function and softmax activation function (IBM Corporation, 2012). The computational efficiency of the models was optimised at a learning rate (η) of 0.4 and momentum (μ) of 0.8. The gradient descent allows for weight correction (i.e. tuning) through the backpropagation of errors across the networks using an iterative process (IBM Corporation, 2012). The number of epochs (ρ) that allowed convergence of input data was 100. Each epoch involved complete iteration over a batch of training set to adjust the weights and minimise errors of the network with an error rate (δ) of 0.00001.
Spatial Correlogram of Neonatal Mortality Rate, 2016/2017 Nigeria MICS
The spatial autocorrelation of NMR dropped to 0.12 between 172.5 and 345 km. Beyond 345 km, NMRs were spatially dispersed, ranging from -0.18 to -0.01. The spatial autocorrelation was zero at 409.8km (i.e. absence of systematic spatial variation in NMR).
Univariate Local Indicator Spatial Autocorrelation (LISA) Cluster Map of Early Neonatal Mortality Rate, 2016/2017 Nigeria MICS.
Global Moran’s I index=0.02, p-value=0.173
There was no significant spatial dependence of early NMR
Univariate Local Indicator Spatial Autocorrelation (LISA) Cluster Map of Late Neonatal Mortality Rate, 2016/2017 Nigeria MICS.
Global Moran’s I index=0.3, p-value=0.001
A statistically significant moderate global spatial clustering of late NMR was observed
Bivariate Local Indicator Spatial Autocorrelation (LISA) Cluster Cap of Neonatal Mortality Rate and Multiple Births, 2016/2017 Nigeria MICS.
Global Moran’s I index=-0.1, p-value=0.05
Open access
