Skip to main content

Under age five children survival times in Nigeria: a Bayesian spatial modeling approach

Abstract

Background

Nigeria is among the top five countries in the world with the highest under-five mortality rates. In addition to the general leading causes of under-five mortality, studies have shown that disparity in sociocultural values and practices across ethnic groups in Nigeria influence child survival, thus there is a need for scientific validation. This study quantified the survival probabilities and the impact of socioeconomic and demographic factors, proximate and biological determinants, and environmental factors on the risk of under-five mortality in Nigeria.

Methods

The Kaplan-Meier survival curve, Nelson Aalen hazard curve, and components survival probabilities were estimated. The Exponential, Gamma, Log-normal, Weibull, and Cox hazard models in a Bayesian mixed effect hierarchical hazard modeling framework with spatial components were considered, and the Deviance and Watanabe Akaike information criteria were used to select the best model for inference. A \(5\%\) level of significance was assumed throughout this work. The 2018 Nigeria Demographic and Health Survey dataset was used, and the outcome variable was the time between birth and death or birth and the date of interview for children who were alive on the day of the interview.

Results

Findings show that the probability of a child dying within the first two months is 0.04, and the probability of a boy child dying before attaining age five is 0.106, while a girl child is 0.094 probability. Gender, maternal education, household wealth status, source of water and toilet facility, residence, mass media, frequency of antenatal and postnatal visits, marital status, place of delivery, multiple births, who decide healthcare use, use of bednet are significant risk factors of child mortality in Nigeria. The mortality risk is high among the maternal age group below 24 and above 44years, and birth weight below 2.5Kg and above 4.5Kg. The under-five mortality risk is severe in Kebbi, Kaduna, Jigawa, Adamawa, Gombe, Kano, Kogi, Nasarawa, Plateau, and Sokoto states in Nigeria.

Conclusion

This study accentuates the need for special attention for the first two months after childbirth as it is the age group with the highest expected mortality. A practicable way to minimize death in the early life of children is to improve maternal healthcare service, promote maternal education, encourage delivery in healthcare facilities, positive parental attitude to support multiple births, poverty alleviation programs for the less privileged, and a prioritized intervention to Northern Nigeria.

Peer Review reports

Background

The under-five mortality rate is still a persistent problem in the world. According to [1], the global number of under-five deaths was estimated as \(5.2\, (5.0, 5.6)\) million in 2019, which implies that approximately 14, 000 children died every day before attaining their fifth birthday in 2019. Since 1990, the global under-five mortality rate has declined by \(59\%\), dropping from 93 deaths per 1000 live births in 1990 to 38 deaths per 1000 live births in 2019, [2]. Although there is a substantial global reduction in under-five mortality rates since the year 1990, disparities exist across countries. Sub-Saharan Africa region remains the highest-burden region of under-five mortality in the world, with one in thirteen children dying before reaching the fifth birthday thus, placing the region at top risk region of under-five deaths, [2]. According to WHO, half of the global under-five deaths occurred in five countries, including Nigeria.

Nigeria is one of the countries that presented a high under-five mortality rate of 117 per 1, 000 live births in 2019, [3]. The country is listed among the top five countries with the highest under-five mortality rate in 2019. A collective effort of the government, stakeholders, and non-governmental organizations in the fight against under-five mortality to improve the well-being of children has led to a consistent decline in the early child mortality rate in Nigeria, [3]. Despite several interventions to lower the under-five mortality index in Nigeria, the rate is still relatively high [4]. The slow declining pace of the rate may hinder the attainment of the Sustainable Development Goal (SDG) to end preventable under-five death in 2030. The SDG target is to reduce under-five mortality to at least 25 per 1,000 live birth by 2030, [2]. Given the statistics, there must be at least a 50% decline in early childhood mortality in Nigeria to attain the SDG by 2030, [5]. Therefore, the under-five mortality burden in Nigeria calls for rapid attention to scale down and below the global mortality average rate.

Leading causes of under-five death include preterm birth complications, birth trauma, pneumonia, congenital anomalies, diarrhea, and malaria, [2, 6], which are preventable with affordable interventions. However, under-five death in Nigeria has been attributed to additional factors including, disparities in sociocultural values and practices, [7]. Cultural beliefs either positively or negatively influence the health status of children, and consequently introduced heterogeneity in improving the survival and well-being of children. For instance, the large disparity of birth in health facilities, births attended by qualified medical practitioner, number of antenatal and postnatal care utilization across different ethnic groups can be attributed to the heterogeneity in the country, [8,9,10]. Nonetheless, there are enormous geographical variations in the under-five mortality rate in the country. According to [5] report, most under-five deaths in Nigeria occurred in the Northern states. Thus, there is a need for scientific validation of the contribution of socioeconomic and demographic factors increasing the risk of child mortality in Nigeria.

Several types of research have been conducted to contribute to the fight against child mortality in Nigeria. For instance, [11] extensively discussed specific medical solutions to reduce neonatal mortality in Nigeria. The authors suggested the adaptation of imported medical technologies through scientific thinking, as many imported technologies may underperform given the environmental impact. Jones et al. [12] investigated the impact of several interventions on child mortality and estimated the scope of coverage of these interventions. According to the authors, there are sufficient medical knowledge and instruments to reduce child mortality, but the mortality rate is still high because the injection of these interventions into the population has not fully considered the inequalities and heterogeneity among the children. Hence, policies and intervention programs should be aware of the factors responsible for these inequalities. In the attempt to determine these factors, [7] examined the effect of ethnicity on under-five mortality in children and found that the risk of mortality significantly varies between ethnic groups in Nigeria. It was found that the mortality rate is high among Hausa, Fulani, and Kanuri tribes in Nigeria. Angela and Uju [13] investigated the effect of child characteristics on under-five mortality in Nigeria. They found that birth between the maternal ages of 20-24 has a lower risk of mortality. Adeyinka et al. [4] determined the socioeconomic indicators contributing to child mortality in Nigeria, and surprisingly, concluded that attendance of skilled health workers during delivery was associated with an increased risk of neonatal death. Ayoade [14] examined the Spatio-temporal distribution of the under-five mortality rate in Nigeria and found that there is spatial clustering and geographical disparities between states in Nigeria. Recently, [15] examined the effect of spatial distribution and other covariates on under-five mortality deaths in Nigeria.

Many researchers adopted the indirect method of estimating mortality rates in Nigeria, where the survival status of children in a cohort is used to estimate the mortality rate. The implicit assumption of this approach is that the births of a cohort are the number of children born in a period, [16]. This assumption is hardly satisfied as studies have shown a disparity in a maternal age group to the risk of under-five mortality, [17]. Moreover, an approach commonly adopted in these researches in multilevel modeling of under-five mortality rate is to consider the death status of a child as a binary outcome variable, [15]. In this manner, information about the time localization of mortality risk is lost in the process. Thus, give rise to a poor time resolution of under-age five mortality and could lead to biased estimates of the variance, [18, 19].

This study aims to quantify the survival probabilities and the impact of socioeconomic and demographic factors, proximate and biological determinants, and environmental factors on the risk of under-five mortality in Nigeria. A direct estimation method of mortality rates was adopted, where the date of birth, survival status, and ages at death of children was utilized for the analysis. In this manner, the heterogeneity in the mortality risk between ages 0 to 59 months is incorporated for inference. For instance, a heap of under-five deaths a few months after birth can be incorporated into the modeling framework for estimating the risk in these age groups. Moreover, this study accounted for the uncertainty in the geographical heterogeneity of under-five mortality rates between the states in Nigeria. Nonparametric survival models and parametric hazard models within a Bayesian hierarchical modeling framework were used to quantify the survival probabilities and the impact of risk factors and spatial patterns on child survival before the fifth birthday. This would unravel the temporal survival pattern and give insight into the existing survival inequalities among children across the country, which could be valuable for policy-making and strategic

Materials and method

Study area and data collection

Nigeria is a west African country with 36 administrative states and a Federal Capital Territory (FCT). Figure 1 shows the map of Nigeria, which served as the study area. The data used were pooled from the 2018 Nigeria Demographic and Health Survey (NDHS), implemented by the National Population Commission (NPC). The survey took place between 14th August and 29th December 2018. The sampling frame adopted is the population and Housing Census conducted by the NPC in 2006. The Enumeration Areas (EA) in the 2006 Census served as the survey’s primary sampling units. The survey implemented a two-stage stratification sampling technique to include EAs in the selection. In the first stage, samples of 1400 EAs were drawn with probability proportionate to size. The resulting listing of households served as the sampling frame in the second stage. A total of 40,427 women were interviewed successfully, yielding a response rate of 99%. Children born within the last 59 months to the interview day were eligible to be included in the study. The age at death of children and the present age of children alive on the interview day was of prime interest. This analysis included a record of 33, 697 children after data cleaning.

Fig. 1
figure 1

Map of Nigeria

Study variables

The outcome variable is the survival times of under-five children within a five years study frame, which is composed of the time between birth and death, \(t_d\), and the time between birth and the date of interview, \(t_a\), if the child was alive at the interview day. Mathematically, \(t_d={\text {Age at death (in years)}}/{12}\) and \(t_a={[\text {Date of interview}-\text {Date of birth (in years)}]}/{12}\). The independent variables included in the hazard model can be categorized into socioeconomic and demographic factors, proximate and biological factors, and environmental factors. These variables are presented in the Appendix. According to the World Health Organization classification, the environmental factors such as household’s toilet facility and source of water supply, were classified as either improved or unimproved. The improved toilet facility includes flush or pour-flush to a piped sewer system, septic tank or pit latrine, ventilated improved pit latrine, pit latrine with slab, and composting toilet. Similarly, the improved water source includes piped household water connection inside the users dwelling, public taps or standpipes, tube wells or boreholes, protected dug wells, protected springs and rainwater collection, [20].

Data analysis

Nonparametric method

Let \(n_i\) be the number of children included in state i, then, \(n=\sum _{i=1}^{37}n_i=33,697\) is the number of children under study and with \(k\le n\) distinct time of deaths between birth and time of the survey, such that \(t_1 \le t_2 \le ...\le t_k\). Note that there is a possibility of having more than one death at each event time. Denote \(d_i\) as the number of deaths at time \(t_i\), and \(q_i\) as the number of children at risk at time \(t_i\). Thus the Kaplan-Meier estimator of the survival function S(t), and the Nelson-Aalen estimator of the cumulative Hazard function H(t) are given as

$$\begin{aligned} \hat{S(t)}=\prod\limits_{t_i:t_i\le t}^{}\Bigg (1-\frac{d_i}{q_i}\Bigg ),\,\, \hat{H(t)}=\sum\limits_{i:t_i< t}^{}\Bigg (\frac{d_i}{q_i}\Bigg ). \end{aligned}$$
(1)

The \(100(1-\alpha )\%\) confidence interval of S(t) is given as \(\hat{S(t)}\pm z(\alpha /2)\sqrt{V( \hat{S(t)})}\), where

$$\begin{aligned} V(\hat{S(t)})=\left( \hat{S(t)}\right) ^2\sum\limits_{i:t_i< t}^{}\frac{d_i}{q_i(q_i-d_i)}. \end{aligned}$$
(2)

For the description of the computational steps for DHS component survival probabilities see [16, 21].

Parametric method

Let \(C_{ij}\) be a random variable representing the number of days between the birthday of child j in state i and the day of the interview. Let \(t_{ij}\) be the number of days a child lived until the event of interest occurred. Define \(y_{ij}=\min (t_{ij},C_{ij})\), and \(\delta _{ij}\in \{0,1\}\) be a censoring indicator that takes the value 1 if \(t_{ij}\le C_{ij}\) (i.e child j in state i experienced the event of interest (death) and takes the value 0 otherwise). Let Y be a random variable that takes the value \(y_{ij}\) for \(j=1,2,...,n_i\), and \(i=1,...,37\) with density function \(f(y_{ij}|\lambda )\). The survival probability is given as

$$\begin{aligned} S(y_{ij}|\varvec{\lambda }) =\int _{y_{ij}}^{\infty } f(x|\varvec{\lambda })d x, \end{aligned}$$
(3)

which is equivalent to \(S(y_{ij}|\varvec{\lambda })=\frac{f(y_{ij}|\varvec{\lambda })}{h(y_{ij}|\varvec{\lambda })}\), where h is the hazard function of Y. Thus the likelihood function to incorporate the incomplete information is given as

$$\begin{aligned} L(Y|\varvec{\lambda }) = \prod\limits_{i=1}^{37}\prod\limits_{j=1}^{n_i} f(y_{ij}|\varvec{\lambda })^{\delta }S(y_{ij}|\varvec{\lambda })^{1-\delta }=\prod\limits_{i=1}^{37}\prod\limits_{j=1}^{n_i} f(y_{ij}|\varvec{\lambda })^{\delta }\Bigg (\int _{y_{ij}}^{\infty } f(x|\varvec{\lambda })d x\Bigg )^{1-\delta }. \end{aligned}$$
(4)

Let \(\eta\) be a linear predictor, which is a function of demographic and socio-economic variables, and represent \(\eta =g(\lambda )\), where g is a known link function usually determined by the density function of Y. The model is specified through the hazard function of Y. Consequently, the likelihood function is written as \(L(Y|g^{-1}(\eta ))\). In the case of a two parameters model with location and scale parameters, \(\lambda _1\), and \(\lambda _2\) respectively, \(\eta =g(\lambda _1)\) and the likelihood is given as \(L(Y|g^{-1}(\eta ),\lambda _2)\).

In this study, four probability density functions were considered for Y, namely, the Exponential, Gamma, Weibull, and Log-normal probability density functions. In addition, the Cox proportional hazard model was also adopted. The summary of the five competing models and their properties for a single child is presented in Table 1.

Table 1 Competing parametric models. \(\Lambda\) is the parameterized Gamma density, and \(\gamma (.,.)\) is the lower incomplete gamma function

In the Cox proportional hazard model [22], the hazard rate is defined as

$$\begin{aligned} h(y_{ij}) = h_0(y_{ij})\exp {(\eta _{ij})}, \end{aligned}$$
(5)

where \(h_0(.)\) is the baseline hazard function. For a given partition of time t of death into K, we assume \(h_0\) to be constant over each time interval. That is

$$\begin{aligned} h_0(t)=\exp (c_k), \,t\in (t_k-1,t_k],\,\,k=1,2,...,K. \end{aligned}$$
(6)

The parameters \(c_k,k=1,2,...,K\) are unknown and estimated from the data. The prior model assigned to \(c_k\) is discussed later in this section.

The performance of each model was evaluated using the Deviance Information Criterion (DIC) and the Watanabe Akaike Information Criterion (WAIC). The lower the criteria values, the better the model, [23]. The best parametric model among the competing models is used to quantify the impact of risk factors and spatial patterns on the hazard rate among under-age five children in Nigeria. Thus, we adopted a Bayesian hierarchical modeling framework.

Given the structural additive linear predictor,

$$\begin{aligned} \eta _{ij}=\mathbf{x}_{ij}^T\boldsymbol{\beta }+\mathbf{z}_{ij}^T\boldsymbol{\psi }+\mathbf{w}_{ij}^T\boldsymbol{\theta }, \end{aligned}$$
(7)

where \(\boldsymbol{\beta }\) is a \(p\times 1\) vector of linear effects, such as the mother’s education level, household wealth index, place of settlement, e.t.c. \(\mathbf{x}\) is a \(p\times 1\) design vector for the linear covariates. \(\boldsymbol{\psi }\) is a \(r\times 1\) vector of non-linear effects to account for variables that varies non-linearly with the time of death, such as the mother’s age. \(\mathbf{z}\) is a \(r\times 1\) design binary vector that links the corresponding effect parameters to the time of death. Finally, \(\boldsymbol{\theta }\) is a \(37 \times 1\) vector of spatial effects to account for heterogeneity between states in Nigeria, and \(\mathbf {w}\) is a \(37\times 1\) vector of covariates to incorporate spatial intercepts.

Prior models

In a Bayesian framework, the parameters in the linear predictor in Equation 7 are considered random, that is, a probability distribution is used to represent the uncertainty in the parameters. Hence, a non-informative multivariate normal distribution was assigned for \(\varvec{\beta }\), that is, \(\varvec{\beta }\sim N_p(\mathbf{0}, 100\mathbf{I})\), where \(\mathbf{I}\) is a \(p\times p\) identity matrix. Since the mothers’ age is equally spaced, a random walk of order two prior distribution was assigned for the non-linear effect, \(\boldsymbol{\psi }\). Let \(\boldsymbol{\psi }=(\psi _1,\psi _2,...,\psi _r)^T\), defined

$$\begin{aligned} \Delta _v=\psi _v-2\psi _{v+1}+\psi _{v+2}\sim N(0,\tau _{\psi }^{-1}),\,\,v=1,2,...,r. \end{aligned}$$
(8)

Thus, the joint prior distribution, given by [24] is written as

$$\begin{aligned} \pi (\boldsymbol{\psi }|\tau _\psi )\propto (\tau _\psi )^{(n-2)/2}\exp \Bigg \{-\frac{\tau _{\psi }}{2}\sum\limits_{v=1}^{r-2}(\Delta _{v})^2\Bigg \}\cdot \end{aligned}$$
(9)

The set of mothers’ ages, R, was ordered at regular intervals, that is, \(R=\{ 15, 16,..., 49\}\). To facilitate estimation, R was re-coded into \(R=\{ 1, 2,..., 35\}\), and \(r=35\). Due to rank deficiency, a sum-to-zero constraint is employed. The random walk two is flexible, invariant under the addition of linear term to \(\boldsymbol{\psi }\), and computationally flexible as it exhibits Markovian property [24].

To account for the spatial effects, the Intrinsic Conditional Autoregressive prior distribution was adopted. Let \(\boldsymbol{\theta }=(\theta _1,\theta _2,...,\theta _{37})\), that is, \(\boldsymbol{\theta }\) accounts for the spatial intercept and heterogeneity between the states. Let \(\boldsymbol{\theta }\) be represented by a Gaussian Markov random field (GMRF) with respect to the Graph \(\mathcal G=(\mathcal V,\mathcal E)\), where \(\mathcal V\) is a set of vertices representing the 37 states and \(\mathcal E\) is the set of edges connecting the states, thus, for \(i\ne i^*\), \(\{i,i^*\}\in \mathcal E\) if and only if state i and \(i^*\) are neighbors. Let W be a \(37\times 37\) adjacency matrix such that \(w_{ii^*}=1\) if \(\{i,i^{\star }\}\in \mathcal E\) and zero otherwise. Let \(\mathbf{D}\) be a diagonal matrix, where each diagonal entry \(d_{ii}\) equals the total numbers of neighbors of state \(i\in \mathcal V\). Therefore, the precision matrix for the GMRF is given as

$$\begin{aligned} Q=\tau _{\theta } \mathbf{D}(\mathbf{I}-\mathbf{W}). \end{aligned}$$
(10)

Notice that by construction, Q is sparse, and hence, \(\boldsymbol{\theta }\) is a Gaussian Markov random field with respect to \(\mathcal G\). That is, \(\theta _{i}\) and \(\theta _{i^*}\) are conditionally independent given \(\boldsymbol{\theta }_{-ii^*}\) if and only if \(\{ i,i^* \}\notin \mathcal E\). Therefore, the prior probability density function is given as

$$\begin{aligned} \pi (\boldsymbol{\theta }|\tau _\theta )\propto \log (|Q|)\Big (-\frac{1}{2}\boldsymbol{\theta }^T Q\boldsymbol{\theta }\Big ). \end{aligned}$$
(11)

For the baseline hazard function Eq. (6), a first-order random walk was assigned for \(c_k,\, k=1,2,...,K\). The definition is analogous to Eq. (8). Hence,

$$\begin{aligned} \Delta _k= & {} c_k-c_{k-1}\sim N(0,\tau ^{-1}_{c}),\nonumber \\ \pi (\boldsymbol{\mathbf c }|\tau _{c})\propto & {} (\tau _c )^{(n-1)/2}\exp \Bigg \{-\frac{\tau _{c}}{2}\sum\limits_{k=1}^{r-1}(\Delta _{k})^2\Bigg \},\\\mathbf c=&(c_1,c_2,...,c_K)^T. \end{aligned}$$
(12)

A log gamma hyperprior was assigned to \(\log \tau , \log \tau _c,\log \tau _\psi ,\log \tau _\theta\), with hyperparameters \((1,5e^{-5})\). A \(\log\) gamma prior model was assigned to \(\log \phi\) with parameter (1, 0.01), and \(\alpha\) was assigned a PC-prior distribution given as

$$\begin{aligned} \pi (\alpha )=\frac{\lambda }{2}\exp {(-\lambda d(\alpha ))}\left| \frac{\partial d(\alpha )}{\partial \alpha }\right| , \end{aligned}$$
(13)

where \(d(\alpha )=\sqrt{2 KLD (\alpha )}\), and \(KLD(\alpha )=(\Gamma ((1+\alpha )/2)\alpha +\alpha \log (\alpha )-\alpha \gamma +\gamma -\alpha )/\alpha\), \(\gamma\) is the Euler’s constant, and \(\lambda =5\).

Bayesian estimation

To estimate the parameters of interest, the structural additive linear predictor, (Equation 7), is used to formulate a latent variable \(\mathcal X=\{\boldsymbol{\beta },\boldsymbol{\psi },\boldsymbol{\theta }\}\) having a multivariate Gaussian distribution with a sparse precision matrix \(\mathbf{Q}(\boldsymbol{\phi }_1)\). Let \(\boldsymbol{\phi }_2\) be a vector of all the hyperparameters of the likelihood function and the prior distributions, and \(\boldsymbol{\phi }=(\phi _1,\phi _2)\) with joint distribution \(\pi (\boldsymbol{\phi })\). Thus, the joint posterior distribution is given as

$$\begin{aligned} \pi (\mathcal X,\boldsymbol{\phi }|\mathbf{y})\propto \pi (\boldsymbol{\phi })\pi (\mathcal X|\boldsymbol{\phi })L(Y\mid \mathcal {X},\boldsymbol{\phi }), \end{aligned}$$
(14)

where \(L(Y\mid \mathcal {X},\boldsymbol{\phi })\) is the adequate likelihood function given in Equation 4. The interest is to make inferences from the posterior marginal distributions \(\pi (\mathcal X|\mathbf{y})\) and \(\pi (\boldsymbol{\phi }|\mathbf{y})\). The formulation (Equation 14) can be estimated using Integrated Nested Laplace Approximation (INLA) described in [25, 26]. The parameters of interest were estimated using the R-INLA package in R [27]. The R-INLA package has been used as an estimation tool in several analyses across different study fields and has been shown to perform adequately, [28,29,30]. R-INLA estimates the posterior marginal distribution of the parameters of interest for inference. The estimation procedure is fully discussed in [23, 29, 31,32,33]. For the non-linear and spatial effect estimates, the output of the R-INLA provides the posterior marginal distribution for each non-linear variable and each administrative state, which is then summarized and graphically displayed using adequate functions in R. Given that a sum to zero constraints was imposed on the spatial effect, the probability of excess risk was computed from the posterior marginal distribution of the spatial effect, which is given as

$$\begin{aligned} \pi (\theta _i < 0 | \mathbf{y}) = \int _{-\infty }^{0}\pi (\theta _i |\mathbf{y})d \theta _i,\,\,i=1,2,...,37, \end{aligned}$$
(15)

where \(\pi (\theta _i|\mathbf{y})\) is the posterior marginal distribution of the ith administrative state spatial effect. The interpretation of \(\pi (\theta _i|\mathbf{y})\) varies between the models and the link functions adopted in the model specification. For the log-normal and gamma model, a higher probability, \(\pi (\theta _i|\mathbf{y})\), indicates a lower risk of death or a higher probability of a child residing in state i to survive past age five. Whereas, for Cox proportional hazard, exponential, and Weibull model, a lower value of \(\pi (\theta _i|\mathbf{y})\) indicates a lower risk of death or a higher probability a child residing in the state i to survive past age five. Moreover, the integral in Equation 15 was approximated using Monte Carlo samples from the marginal posterior distribution. That is

$$\begin{aligned} \pi (\theta _i<0 | \mathbf{y}) \approx \frac{\sum\nolimits_{l=1}^{N}I(\theta _{il}< 0)}{N}, \end{aligned}$$
(16)

where \(\theta _{il}\) is the lth sample from \(\pi (\theta _i |\mathbf{y})\), I(.) is an indicator function and \(N=10,000\).

Results

Among the 33, 697 children included in the analysis, the estimated censoring rate is \(90\%\). It indicates that for every 1000 live birth, 100 children never attained their fifth birthday. The descriptive statistics of the data used are shown in the Appendix.

Nonparametric method

Figure 2a presents the overall Kaplan-Meier survival probability estimates without taking any covariate effects into account. The middle black line indicates the mean survival probability, while the red lines indicate the 95% confidence interval of the survival probability. The curve shows a sharp drop in the survival probabilities among children within the first two months, and further consistently drops at a constant rate until the age of 25th months. It slowly descends further until a child attains 40 months, and then levels off at a probability of 0.88. The probability of a child surviving past the first two months declines significantly, and consequently, the probability of dying within this period increases with a magnitude of at least \(1-0.96=0.04\). In Figure 2b, the gender effect was taken into account. Though the survival curves for both males and females relatively have a similar pattern, the survival probability of male children is lower compared with that of female children. While the survival curve for females levels off at 0.906, the curve for male children leveled off at 0.894 before the fifth birthday. This implies that male children are more likely to die, with a probability of \(1-0.894=0.106\) compared with female children with a probability \(1-0.906=0.094\). Figure 3a presents the cumulative hazard rate computed using the Nelson Aalen estimator of cumulative hazard curve, and Figure 3b presents the cumulative component probabilities of a child dying before reaching 59 months at each age group. In Figure 3a, there is a sharp increase in the cumulative hazard rates within the first two months after birth. The curve then increased at a constant pace until 25 months and then gradually levels off until 59 months. In Figure 3b, the highest component probability occurred within the first months after birth. It is followed by the probability associated with a child crossing from the age group \(3-5\) to \(12-23\) months and \(12-23\) to \(24-26\) months.

Fig. 2
figure 2

a Kaplan-Meier survival probability estimate, b Kaplan-Meier survival probability estimate accounting for gender

Fig. 3
figure 3

a Nelson Aalen cumulative hazard curve, and b cumulative component mortality probability estimate. The highest hazard rate and highest component mortality probability are within the first two months after birth

Parametric method

Table 2 shows the DIC and WAIC information criteria. Based on these measures, the log-normal model is the most appropriate model for the analysis, and thus, only the results of the log-normal model are presented. However, the presented result is consistent for all the fitted models. It is important to mention that the interpretation of the models’ effects varies between the models considered and depends on the link function adopted (see Table 1). For the log-normal model and the link function considered in this study, an increased value of model effects implies decreased hazard rate and increased survival probability. Thus, factors with lower effects (coefficients) are considered risk factors for under-five mortality.

Table 2 Model adequacy. The Log-normal model has the least DIC and WAIC, and thus, it is the best model among the competing models

Linear effect

The linear effect estimates are presented in Table 3, showing the posterior mean, standard deviation (sd), the lower (2.5%), and upper (97.5%) quantile of the 95% credible interval. The reference category for the categorical variables is indicated as “Reference”. The non-inclusion of zero in the credible interval was used to determine significance.

Table 3 Posterior estimates of the fixed effects

The result revealed that female children are significantly less likely to die before their fifth birthday compared with male children, which is the counterpart category. This result is consistent with the Kaplan-Meier estimates. Children with mothers who achieved at least a primary level of education are more likely to attain age five compared with children whose mothers did not have any formal education. Moreover, as mothers’ education increases, the more likely a child will celebrate the fifth birthday. A similar pattern is found for the household wealth quantile. Children whose household wealth quantile is at least poorer are significantly more likely to attain the fifth birthday compared with the poorest category. As the wealth quantile increases, the more likely a child would survive the fifth birthday.

Children living in a household with improved water and toilet facility are significantly more likely to celebrate their fifth birthday compared with their counterpart, which is unimproved water and toilet facility. The result shows that children living in rural areas are significantly more likely to attain the fifth birthday compared with their counterparts living in urban settlements. Children whose parents pay attention to at least one mass media means are significantly more likely to attain age five. Children whose parents visit antenatal care more than five times before birth is not significantly different from those whose parents had 2-4 visits before childbirth. Surprisingly, children whose parents do not visit antenatal care are more likely to attain age five compared with the base category. However, for postnatal care utilization, this is not the case. Children whose parent frequently visits postnatal care are significantly more likely to attain age five.

Children not being breastfed are significantly more likely to die before their fifth birthday. Unexpectedly, the result indicated that home-delivered children are less likely to die before age five compared with their counterparts who were delivered in a health facility. Children whose mothers have been working (employed) in the last twelve months are more likely to attain age five compared with children whose mothers have not been working in the last twelve months before the survey. For the type of religion, children whose religion is Catholic, Islam, Traditional, and Others are more likely to attain their fifth birthday compared with the base category, which is other Christians. Children whose birth was multiples are more likely to die before age five. However, the birth position is not significant.

The result indicates that the mode of delivery does not significantly contribute to the survival of a child. As expected, children who took tetanus injections are less likely to die before age five compared with children who do not. Surprisingly, children whose household has a bednet are more likely to die before the fifth birthday compared with children whose household does not have a bednet. However, children who sleep under bednets are less likely to die before age five. Children whose both parents or spouses decide the health care use are more likely to attain the fifth birthday compared with the children whose health care use is only decided by the respondent alone. Moreover, children whose health care use is decided by someone else are not significantly different from the base category. Access to a healthcare facility does not significantly improve the likelihood that a child would survive age five.

Spatial effect

Figure 4 presents the posterior estimate of the spatial effect of a child dying before the fifth birthday. In Figure 4a the reddish States represent the states with a higher likelihood of under-five mortality, while the greenish States are regions with a lower likelihood relative to the overall country’s average. Figure 4b presents the probability of elevated risk of a child dying before attaining age five. In general, the northern region of Nigeria is at a higher risk of under-five mortality. In other words, the likelihood of child mortality is more severe in the North. Among the Northern states, Kebbi, Kaduna, and Jigawa states have above 0.70 probability of elevated risk of under-five mortality compared with the overall average risk. That is, the likelihood of child mortality in these states is above \(70\%\) higher compared with the country’s average. Invariably, Kebbi, Kaduna, and Jigawa can be tagged as child mortality hotspots in Nigeria. These are followed by Adamawa, Gombe, Kano, Kogi, Nasarawa, Plateau, and Sokoto states. However, turning attention to the southern region of Nigeria, Ekiti, Imo, and Osun states have below 0.25 probability of having an elevated risk of a child dying before attaining age five compared with the overall average risk. It indicates that the likelihood of under-five mortality in these states compared with the country’s average is below \(25\%\).

Fig. 4
figure 4

a The spatial effect of state base survival risk, and b is the probability of elevated risk of mortality before attaining the fifth birthday. The red states indicate regions with a higher mortality risk than the overall average. States with higher probability indicate regions with a higher probability of under-five death before attaining the fifth birthday

Non-linear effect

Figure 5 presents the non-linear effects of maternal age and a child’s birth weight on child mortality before attaining the fifth birthday. In general, the shape of the effect of maternal age and birth weight on child mortality is an “n” shape. In Figure 5a the likelihood that a child will attain the fifth birthday increases gradually as the maternal age increases from 15 years and attains a peak at the age of 27 years, and maintains the peak for further 6 years, and then gradually decreases through age 35 to age 49 years. This indicates that there is a higher likelihood of child mortality among young women and older women compared with women within the age range of 24 to 44 years. Turning attention to Figure 5b, the likelihood of a child to attain the fifth birthday gradually increases from a birth weight of 0.5 kg until it attains a peak at 3kg and maintains the same effect for an additional 0.5kg. The likelihood gradually decreases until a birth weight of 6kg. In other words, the likelihood of child mortality is higher among birth weights less than 2.5kg and those higher than 4.5kg.

Fig. 5
figure 5

a The effect of maternal age (years), and b child’s birth weight (Kg) to under-five mortality before attaining the fifth birthday. The risk is higher at the beginning and end of maternal age, and birth weight below 2.5kg and above 4.5kg

Discussion

This study considered the time-to-death outcomes of under-five children in Nigeria using both parametric and nonparametric methods. The results of the adopted nonparametric methods were consistent, and it indicates a higher risk of under-five mortality within the first two months after birth and thus, should be considered as a hot period in time localization. Improving the well-being of children at this stage could significantly scale down the child mortality index in Nigeria. Moreover, among all the parametric models considered, the log-normal model stood out as the most suitable model to explain the variations in the time-to-death outcomes variable. Findings show that gender, maternal education, household wealth status, source of water and toilet facility, residence, mass media, frequency of antenatal and postnatal visits, marital status, place of delivery, multiple births, who decide healthcare use, use of bednet are significant risk factors of child mortality. Findings also indicate a disparity in geographical regions, across mothers’ ages, and birth weights.

Findings show that the survival probability of under-five children has the highest drop within the first two months after birth. That is, children within the first two months of birth have the highest expected risk of death. This finding corroborates the result conducted in rural Burkina Faso, [34], and Ethiopia, [35]. Hence, this study suggests that intervention programs set to mitigate mortality should give priority to children within the first two months as this could effectively scale down the overall under-five mortality index in Nigeria. Findings indicate that male children have an elevated risk of dying before their fifth birthday compared with female children in the same age group. This finding is consistent with biological results that male children are biologically weaker and more susceptible to diseases compared with female children, [36, 37]. Also, it may suggest that male children are more vulnerable to health inequalities which is more noticeable in the early stages after birth, [38]. Findings showed that attainment of higher maternal education increases the survival chances of children. This finding is expected as educated parents are more likely to have a positive attitude towards parenting, often seek expert advice, and involve in a child’s physical and mental growth, [39, 40]. Thus, this finding suggests that maternal education is among the key important markers to mitigating the high mortality rate in Nigeria.

This study identified a discrepancy in the risk of mortality among households’ wealth status. That is, the risk of mortality is most severe among the poorest wealth quantile, and the risk consistently decreases with an increase in wealth quantile. The most prominent increase in under-five survival was found for a one-step wealth quantile increase from the rich to the richest. However, for intervention, improving household wealth quantile from rich to richest seems impractical, rather programs could focus on the improvement from poorest to poorer quantile by providing basic living conditions, such as access to good shelter, roads, and electricity being that the improvement from poorest to poorer wealth quantile is the second most prominent increase in under-five survival probability. This finding is consistent with existing studies, [41, 42]. Wealth inequalities, non-diversification of the economy, unemployment, and low education financing in Nigeria are the major key factors contributing to increasing poverty, [43], and consequently, increase the risk of child mortality in Nigeria. This finding identifies poverty as a leading cause of mortality. Poverty alleviation programs to improve the livelihoods of households are one of the main remedies to child mortality in Nigeria.

Findings also indicated that access to an improved water source and toilet facility lowers the risk of mortality. However, the improvement is more prominent for the toilet facilities. Thus, policies to safely managed improved household water sources and toilet facilities are required to mitigate child mortality, [44]. Surprisingly, the frequency of antenatal care visits does not improve the survival probabilities of children. However, the survival probability increases for postnatal visits. It could be associated with late first-time antenatal utilization or irregular visitation. Women tend to wait towards the delivery period or complications during gestation before utilizing antenatal care services. Thus, educating mothers on the importance of early antenatal utilization could scale down the under-five mortality rate. Similarly, home delivery in Nigeria significantly improves the survival probabilities of a child to celebrate the fifth birthday. This finding is contrary to the common rationale that babies delivered in a health facility are more likely to attain their fifth birthday since medical experts are readily available to give the first medical support and advice. According to [8], the prevalence of non-utilization of health facilities for childbirth is about 62% of the total birth in Nigeria. It could be associated with poor quality of health facility, caring, cost, and cultural beliefs that lack the acceptability of health care services, and the belief that healthcare facility is only necessary during birth complications, [45]. Thus, intervention programs that give financial support and prenatal education are needed to encourage the frequent use of health facilities but not only when there are complications during childbirth.

Findings indicate that maternal occupation since a year before the survey significantly lowers the risk of child mortality. This finding corroborates the previous study, [46]. However, the reverse may be the case if only children under breastfeeding were considered, [47]. Findings indicate that married couples increase the survival probability of children attaining their fifth birthday. Nonetheless, divorced, deceased spouses or separated parents also increase survival compared to never-married parents. It could be linked to the elevated likelihood of children of unmarried parents suffering the consequence of limited resources, poverty, and poor schooling, [48, 49]. This study also found that multiple births significantly increase the risk of under-five mortality compared with a singleton. It could be associated with the consequence of shared limited family resources among the children. This finding suggests that multiple births jeopardize the survival of under-five children. Hence, children of multiple births should be given special treatment in addition to positive involvement in parental support. Findings also indicate that immunization with Tetanus strongly increases the survival of children. This finding is in line with the result of a previous study, [50]. Thus, immunization with Tetanus and its benefits on the under-five survival should be well publicized and target parents with low maternal education attainment. This study found that household with bednets does not lower the risk of child mortality. However, the usage of insecticide-treated bednets among children significantly lowers the risk of mortality in children. The magnitude of the effect of bednet use in the improvement of child survival is equivalent to the magnitude of immunization with Tetanus. This finding suggests that the non-usage of insecticide-treated bednet increases the risk of mortality, and it should be considered as a significant risk factor. Thus, intervention programs that donate and distribute insecticide-treated bednets to households in Nigeria should take it as an additional duty to educate parents on the importance of the use of bednets on the survival of children. This finding is consistent with a previous study in Ghana, [51]. This study found that both parent and husband’s decision on the mother’s health use significantly improves child survival compared to mother decision alone. Moreover, a someone-else decision on the mother’s use of health services increases the risk of mortality. This finding suggests that decisions on mother’s health use outside the parents could increase the risk of mortality. Hence, at least spouses (fathers) should be involved in the decision-making on health use. The role of the father to a child’s health plays a crucial part in child survival, as studies have shown that the father’s involvement is related to positive child health outcomes, [52].

This study also substantiates the effect of maternal age and birth weight on child mortality. Findings indicate that the risk of mortality among under-age five children increases for maternal age below 24 and above 44 years. This finding suggests that maternity age in the mentioned disjoint age intervals is a risk factor for child mortality. For younger women, this could be attributed to a lack of experience with a child’s nutritional needs or could be attributed to the biological and social mechanism, [53]. Whereas for older women, it could be attributed to the prevalence of complications during childbirth and preterm delivery, [54]. Similarly, the likelihood of under-age five mortality is high at the lowest birth weights, improves within the most frequent average birth weights, and then picks again at extremely higher birth weights. Specifically, it was found that birth weight less than 2.5kg and above 4.5kg elevates the risk of underage five mortality. It could be attributed to the prevalence of subnormal growth and neurodevelopmental complications among children with low birth weight, [55]. These findings are consistent with previous studies, [56, 57].

This study identified significant spatial heterogeneity in the risk of under-five mortality before attaining the fifth birthday. Findings show a clear divide between North and South Nigeria, placing the Northern states at higher risk regions compared with the South. Specifically, among the Northern states in Nigeria, this study identified Kebbi, Kaduna, and Jigawa as the hotspot of child mortality in Nigeria. Additionally, Adamawa, Gombe, Kano, Kogi, Nasarawa, Plateau, and Sokoto states were also identified as high-risk regions. This finding corroborates the work of [15].

This study has some limitations. Firstly, regional factors that directly influence the survival probabilities of children were not available in the data. These factors include variables related to food security and child nutritional intake, which can play key roles in a child’s survival. Moreover, there is a nonavailability of data on the cause of child death to account for the risk of unavoidable deaths and its impact on the research results. Secondly, the adopted model could suffer from the omission of diseased children from the survey and the missed information on the exact day of death, which consequently may affect the point survival probability estimates. Lastly, the multi-stage design of the survey was not incorporated in the estimation stage. Future work can consider such implementation to verify whether there is a sufficient improvement in the results relative to the complexity it introduces into the model.

Conclusion

This study quantified the risk of under-five mortality in the form of survival probabilities before a child attains the fifth birthday, and also identified the risk factors elevating the under-five mortality rate in Nigeria. These findings call for action on the identified risk factors to mitigate the under-five mortality rate in Nigeria. This study accentuates the need for special attention for the first two months after childbirth as it is the age group with the highest expected mortality. A practicable way to minimize death in the early life of children is to improve maternal healthcare service, promote maternal education, encourage delivery in healthcare facilities, positive parental attitude to support multiple births, poverty alleviation programs for the less privileged, and a prioritized intervention to Northern Nigeria. The identified hotspot could be used for policy-making and as a guide towards the distribution of intervention resources. Additionally, it could serve for the construction of a localized framework for intervention program evaluation. The maps produced in this study could be placed beside another map, such as poverty, unemployment rate, or insurgence maps to identify possible correlation across the country.

Availability of data and materials

The Nigeria Demographic and health survey data used in this study are available and accessible on the DHS program website (http://dhsprogram.com).

Abbreviations

DIC:

Deviance Information criterion

EA:

Enumeration Areas

NDHS:

Nigeria Demographic and Health Survey

NPC:

National Population Commission

R-INLA:

R-Integrated Nested Laplace Approximation

UNICEF:

United Nations Children’s Fund

WAIC:

Watanabe-Akaike Information Criterion

WHO:

World Health Organization

References

  1. UNICEF. Levels and trends in child mortality. Technical report, eSocialSciences. 2020.

  2. World Health Organization. (2021). Children: improving survival and well-being. 2020. https://www.who.int/news-room/factsheets/detail/children-reducing-mortality. Accessed 1 Nov 2021.

  3. Unicef. UNICEF Data: Monitoring the situation of children and women. 2020. https://data.unicef.org/country/nga/. Accessed 1 Nov 2021.

  4. Adeyinka DA, Muhajarine N, Petrucka P, Isaac EW. Inequities in child survival in Nigerian communities during the sustainable development goal era: insights from analysis of 2016/2017 multiple indicator cluster survey. BMC Public Health. 2020;20(1):1–18.

    Article  Google Scholar 

  5. Unicef, et al. Nigeria multiple indicator cluster survey 2016-2017. Technical report, NBS, UNICEF. 2018.

  6. Abraha HE, Belachew AB, Ebrahim MM, Tequare MH, Adhana MT, Assefa NE. Magnitude, trend, and causes of under-five mortality from kilite-awlaelo health demographic surveillance database, northern Ethiopia, 2009–2017. BMC Public Health. 2020;20(1):1–8.

    Article  Google Scholar 

  7. Adedini SA, Odimegwu C, Imasiku EN, Ononokpono DN. Ethnic differentials in under-five mortality in Nigeria. Ethn Health. 2015;20(2):145–62.

    Article  PubMed  Google Scholar 

  8. Adedokun ST, Uthman OA. Women who have not utilized health service for delivery in Nigeria: who are they and where do they live? BMC Pregnancy Childbirth. 2019;19(1):1–14.

    Article  Google Scholar 

  9. PMNCH. Women’s children’s and adolescents’ health. 2021.

  10. Olopha PO, Fasoranbaku AO, Gayawan E. Spatial pattern and determinants of sufficient knowledge of mother to child transmission of hiv and its prevention among Nigerian women. PLoS ONE. 2021;16(6):0253705.

    Article  Google Scholar 

  11. Amadi HO, Kawuwa MB. Reducing early neonatal mortality in Nigeria-the solution. Sel Top Neonatal Care. 2018;1:221.

    Google Scholar 

  12. Jones G, Steketee RW, Black RE, Bhutta ZA, Morris SS, Group BCSS, et al. How many child deaths can we prevent this year? Lancet. 2003;362(9377):65–71.

    Article  PubMed  Google Scholar 

  13. Angela C, Uju O. Determinants of under-five mortality in Nigeria: an application of cox proportional hazard and cox frailty models. IOSR J Math Ver I [Internet]. 2015;11(4):2278–5728.

    Google Scholar 

  14. Ayoade MA. Spatio-temporal patterns of under 5 mortality in Nigeria. SSM-Popul Health. 2018;6:116.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Fagbamigbe AF, Nnanatu CC. Modelling the spatial distribution and the factors associated with under-five mortality in Nigeria. Spat Demogr. 2021;10(2):1–28.

    Google Scholar 

  16. DHS. Guide to DHS statistics DHS-7. 2021. https://dhsprogram.com/data/Guide-to-DHS-Statistics/Early_Childhood_Mortality.htm. Accessed 1 Nov 2021.

  17. Sinha S, Aggarwal AR, Osmond C, Fall CH, Bhargava SK, Sachdev HS. Maternal age at childbirth and perinatal and under-five mortality in a prospective birth cohort from delhi. Indian Pediatr. 2016;53(10):871–7.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Royston P, Altman DG, Sauerbrei W. Dichotomizing continuous predictors in multiple regression: a bad idea. Stat Med. 2006;25(1):127–41.

    Article  PubMed  Google Scholar 

  19. MacCallum RC, Zhang S, Preacher KJ, Rucker DD. On the practice of dichotomization of quantitative variables. Psychol Methods. 2002;7(1):19.

    Article  PubMed  Google Scholar 

  20. WHO/UNICEF. Joint Monitoring Programme for Water Supply and Sanitation: Progress on sanitation and drinking-water: 2010 update. Geneva: World Health Organization; 2010.

  21. United Nations. Mortality estimates from major sample surveys: towards the design of a database for the monitoring of mortality levels and trends. New York: United Nations, Department of Economic and Social Affairs; 2011.

  22. Cox DR. Regression models and life-tables. J R Stat Soc Ser B (Methodol). 1972;34(2):187–220.

    Google Scholar 

  23. Egbon OA, Somo-Aina O, Gayawan E. Spatial weighted analysis of malnutrition among children in Nigeria: A Bayesian approach. Stat Biosci. 2021;13(3):1–29.

    Article  Google Scholar 

  24. Lindgren F, Rue H. On the second-order random walk model for irregular locations. Scand J Stat. 2008;35(4):691–700.

    Article  Google Scholar 

  25. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent gaussian models by using integrated nested laplace approximations. J R Stat Soc Ser B (Stat Methodol). 2009;71(2):319–92.

    Article  Google Scholar 

  26. Martino S, Riebler A. Integrated nested laplace approximations (inla). Wiley StatsRef Stat Ref Online. 2014;1–19.

  27. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2019. https://www.R-project.org/.

  28. Lindgren F, Rue H, Lindström J. An explicit link between gaussian fields and gaussian markov random fields: the stochastic partial differential equation approach. J R Stat Soc Ser B (Stat Methodol). 2011;73(4):423–98.

    Article  Google Scholar 

  29. Martins TG, Simpson D, Lindgren F, Rue H. Bayesian computing with inla: new features. Comput Stat Data Anal. 2013;67:68–83.

    Article  Google Scholar 

  30. Egbon OA, Belachew AM, Bogoni MA. Risk factors of concurrent malnutrition among children in ethiopia: a bivariate spatial modeling approach. All Life. 2022;15(1):512–36.

    Article  Google Scholar 

  31. Ruiz-Cárdenas R, Krainski ET, Rue H. Direct fitting of dynamic models using integrated nested laplace approximations-inla. Computational Statistics & Data Analysis. 2012;56(6):1808–28.

    Article  Google Scholar 

  32. Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with r-inla. Spat Spatio-Temporal Epidemiol. 2013;4:33–49.

    Article  Google Scholar 

  33. Palmi-Perales F, Gomez-Rubio V, Martinez-Beneito MA. Bayesian multivariate spatial models for lattice data with inla. 2019. arXiv preprint arXiv:1909.10804.

  34. Schoeps A, Gabrysch S, Niamba L, Sié A, Becher H. The effect of distance to health-care facilities on childhood mortality in rural burkina faso. Am J Epidemiol. 2011;173(5):492–8.

    Article  PubMed  Google Scholar 

  35. Ayele DG, Zewotir TT, Mwambi H. Survival analysis of under-five mortality using cox and frailty models in Ethiopia. J Health Popul Nutr. 2017;36(1):1–9.

    Article  Google Scholar 

  36. Costa JC, da Silva ICM, Victora CG. Gender bias in under-five mortality in low/middle-income countries. BMJ Glob Health. 2017;2(2):000350.

    Article  Google Scholar 

  37. Pongou R. Why is infant mortality higher in boys than in girls? a new hypothesis based on preconception environment and evidence from a large sample of twins. Demography. 2013;50(2):421–44.

    Article  PubMed  Google Scholar 

  38. Egbon OA, Belachew AM, Bogoni MA. Modeling spatial pattern of anemia and malnutrition co-occurrence among under-five children in Ethiopia: A Bayesian geostatistical approach. Spat Spatio-Temporal Epidemiol. 2022;43:100533.

    Article  Google Scholar 

  39. Andriano L, Monden CW. The causal effect of maternal education on child mortality: Evidence from a quasi-experiment in malawi and uganda. Demography. 2019;56(5):1765–90.

    Article  PubMed  Google Scholar 

  40. Cleland JG, Van Ginneken JK. Maternal education and child survival in developing countries: the search for pathways of influence. Soc Sci Med. 1988;27(12):1357–68.

    Article  PubMed  CAS  Google Scholar 

  41. Chao F, You D, Pedersen J, Hug L, Alkema L. National and regional under-5 mortality rate by economic status for low-income and middle-income countries: a systematic assessment. Lancet Glob Health. 2018;6(5):535–47.

    Article  Google Scholar 

  42. Cha S, Jin Y. Have inequalities in all-cause and cause-specific child mortality between countries declined across the world? Int J Equity Health. 2020;19(1):1–13.

    Article  Google Scholar 

  43. Ucha C. Poverty in Nigeria: Some dimensions and contributing factors. Glob Majority E-J. 2010;1(1):46–56.

    Google Scholar 

  44. Ezeh OK, Agho KE, Dibley MJ, Hall J, Page AN. The impact of water and sanitation on childhood mortality in Nigeria: evidence from demographic and health surveys, 2003–2013. Int J Environ Res Public Health. 2014;11(9):9256–72.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Konje ET, Hatfield J, Kuhn S, Sauve RS, Magoma M, Dewey D. Is it home delivery or health facility? community perceptions on place of childbirth in rural northwest tanzania using a qualitative approach. BMC Pregnancy Childbirth. 2020;20(1):1–11.

    Article  Google Scholar 

  46. Akinyemi JO, Solanke BL, Odimegwu CO. Maternal employment and child survival during the era of sustainable development goals: Insights from proportional hazards modelling of Nigeria birth history data. Ann Glob Health. 2018;84(1):15.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Bankole A. Maternal employment and infant mortality: an examination of the role of breastfeeding as an intermediate factor. Biol Soc J Eugenics Soc. 1989;6(1):19–26.

    CAS  Google Scholar 

  48. Clark C, Hamplová D. The impact of mother’s marital status on child mortality in Sub-Saharan Africa: An analysis of birth and marital histories. Canada: Canada Research Chair in Youth, Gender and Global Health; 2011.

    Google Scholar 

  49. McLanahan S, Sandefur G. Growing up with a Single Parent. What Hurts, What Helps. USA: Harvard University Press; 1994.

  50. Blencowe H, Lawn J, Vandelaer J, Roper M, Cousens S. Tetanus toxoid immunization to reduce mortality from neonatal tetanus. Int J Epidemiol. 2010;39(suppl-1):102–9.

    Article  Google Scholar 

  51. Afoakwah C, Nunoo J, Andoh FK. Effect of insecticide-treated bed net usage on under-five mortality in northern ghana. Malar J. 2015;14(1):1–6.

    Article  Google Scholar 

  52. Nair S, Chandramohan S, Sundaravathanam N, Rajasekaran AB, Sekhar R. Father involvement in early childhood care: insights from a mel system in a behavior change intervention among rural indian parents. Front Public Health. 2020;8:516.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Finlay JE, Özaltin E, Canning D. The association of maternal age with infant mortality, child anthropometric failure, diarrhoea and anaemia for first births: evidence from 55 low-and middle-income countries. BMJ Open. 2011;1(2):000226.

    Article  Google Scholar 

  54. Cavazos-Rehg PA, Krauss MJ, Spitznagel EL, Bommarito K, Madden T, Olsen MA, Subramaniam H, Peipert JF, Bierut LJ. Maternal age and risk of labor and delivery complications. Matern Child Health J. 2015;19(6):1202–11.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Hack M, Klein NK, Taylor HG. Long-term developmental outcomes of low birth weight infants. Future Childr. 1995;176–196.

  56. Kim Y-N, Choi D-W, Kim DS, Park E-C, Kwon J-Y. Maternal age and risk of early neonatal mortality: a national cohort study. Sci Rep. 2021;11(1):1–9.

    Google Scholar 

  57. Wilcox AJ, Russell IT. Birthweight and perinatal mortality: II. On weight-specific mortality. Int J Epidemiol. 1983;12(3):319–25.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

Osafu Augustine Egbon and Mariella Ananias Bogoni acknowledge the support from the Brazilian institution CAPES. Francisco Louzada acknowledges the support from the Brazilian institution CNPq and FAPESP. The authors would like to acknowledge the Demographic and Health Survey program for granting access to the data used in the study.

Funding

No funding

Author information

Authors and Affiliations

Authors

Contributions

OAE conceived the study idea, BTB and MAB sourced the data, OAE and MAB performed the computation, FL and BTB validated the model. All authors participated in the manuscript writing, read and accepted the final manuscript revision.

Corresponding author

Correspondence to Bayowa Teniola Babalola.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in this study involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Egbon, O.A., Bogoni, M.A., Babalola, B.T. et al. Under age five children survival times in Nigeria: a Bayesian spatial modeling approach. BMC Public Health 22, 2207 (2022). https://doi.org/10.1186/s12889-022-14660-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-022-14660-1

Keywords