Skip to main content

Bayesian spatio-temporal modelling and mapping of malaria and anaemia among children between 0 and 59 months in Nigeria



A vital aspect of disease management and policy making lies in the understanding of the universal distribution of diseases. Nevertheless, due to differences all-over host groups and space–time outbreak activities, data are subject to intricacies. Herein, Bayesian spatio-temporal models were proposed to model and map malaria and anaemia risk ratio in space and time as well as to ascertain risk factors related to these diseases and the most endemic states in Nigeria. Parameter estimation was performed by employing the R-integrated nested Laplace approximation (INLA) package and Deviance Information Criteria were applied to select the best model.


In malaria, model 7 which basically suggests that previous trend of an event cannot account for future trend i.e., Interaction with one random time effect (random walk) has the least deviance. On the other hand, model 6 assumes that previous event can be used to predict future event i.e., (Interaction with one random time effect (ar1)) gave the least deviance in anaemia.


For malaria and anaemia, models 7 and 6 were selected to model and map these diseases in Nigeria, because these models have the capacity to receive strength from adjacent states, in a manner that neighbouring states have the same risk. Changes in risk and clustering with a high record of these diseases among states in Nigeria was observed. However, despite these changes, the total risk of malaria and anaemia for 2010 and 2015 was unaffected.


Notwithstanding the methods applied, this study will be valuable to the advancement of a spatio-temporal approach for analyzing malaria and anaemia risk in Nigeria.


Childhood malaria infection has been a major concern, especially in developing countries like Nigeria. The 2021 report from the World Health Organization (WHO) estimated that 241 million malaria cases with 627,000 deaths worldwide [1, 2]. The increase in malaria infection dropped from 81% in 2000 to 59% in 2015 and 56% in 2016 but went up again to 59% in 2020 due to Covid-19 pandemic. Globally, 96% of malaria cases is majorly from 29 countries with Nigeria (27%) topping the list of 6 countries that contribute almost 55% of malaria cases [1, 3]. Children between the age of 0 to 59 months are most vulnerable with an estimate of 213 million to 228 million malaria cases between 2019 and 2020 and a mortality rate of 534,000 to 602,000 in the respective years; 80% of all malaria deaths are among children under age 5 years [1, 4]. Like most other vector borne diseases, malaria is characterized by spatio-temporal variations or changes due to demographic, socio-economic and geographical factors. These covariates can help determine the spatio-temporal patterns of disease and recognize hotspots to aid efficient examination of disease, cost-effective allocation of resources and most importantly, effective disease control [5,6,7,8,9,10].

The spatio-temporal distribution of vector borne diseases is not only determined by environmental factors. Political and state borders are also major determinants because of their involvement in the spatial distribution and enactment of control and prevention programmes [10]. This can be explained by a spatio-temporal study in Northern Thailand, where there was a sharp difference in the malaria prevalence with Myanmar border [7]. Also, only environmental and biological factors cannot justify the differences in local diseases as claimed by Ra et al.[6].

Anaemia is another disease that has become a public health challenge in Nigeria, especially amongst under aged 5 years children. It is a condition that arises as a result of the reduction of hemoglobin in the blood [11, 12]. Globally, more than 273 million children under aged 5 years are affected by anaemia [12]. Sub-Sahara Africa is the most endemic region with about 53.8% childhood anaemia cases [12]. According to the WHO classification, anaemia is considered severe if its prevalence is 40% and above, moderate between 20% and 39.9%, and mild between 5% and 19.9% [12, 13]. Anaemia has become major public health due to its prevalence and effect on child’s health.

Pregnant women and children are most vulnerable to anaemia because of their high requirement of iron. Children between the age of 6–59 months are anaemic if their haemoglobin level is below 11 g/dl. The major causes of anaemia in children are parasitic infection, dietary iron deficiency and inherited disorders but in malaria endemic region, malaria disease is the major cause [14, 15].

The spread of malaria and anaemia in Nigeria has been a concern to researchers which has led to several studies such as [12, 16,17,18,19,20]. In [21], a quasi-experimental fixed-effect model was used to investigate the effect of malaria on haemoglobin concentration in children under 5 years old. They concluded that there is a strong negative effect on haemoglobin levels among Burkina Faso’s children from malaria infection. Furthermore, [22] stated that anaemia caused by Plasmodium falciparum is a result of excess removal of nonparasitized red blood cells together with the destruction of parasitized red cells immune leading to malfunctioning of the bone marrow. Also, the main cause of mortality and morbidity in children who live in Kenyan malaria hotspots is falciparum malaria [23]. Given these literatures and more, there is a need to monitor the progress of malaria and anaemia in the near future and use the available data to forecast cases in space and for the possible spread of these diseases.

The spatial pattern of diseases and exposures does not explain the temporal variation which is also important and interesting. Besag et al. [24], foremost introduced a spatial pattern which was extended by incorporating a linear time trend for interaction [25]. Knorr‐Held [26], included a non-parametric temporal trend which comprises the time changing effect of predictors. Including disease surveillance studies, spatio-temporal models are mostly used in several fields of science [27]. With the help of Bayesian hierarchical modelling framework, the implementation of these models is made possible. These models accommodate a composite and workable structure in space and time models, with spatio-temporal interaction as the paramount feature. Here, our work is extended from the methods used by Bernardinelli et al. [25], and Knorr‐Held [26] for spatio-temporal framework. By applying the multilevel model analysis [27], we independently investigate the spatio-temporal distribution of malaria and anaemia and their associated risk factors using data from the Nigeria malaria indicator survey (NMIS) for 2010 and 2015.

Data description

This work obtained data from the Nigeria Malaria Indicator Survey (NMIS) which was carried out by the National Malaria Elimination Program (NMEP), National Population Commission (NPopC) and the National Bureau of Statistics (NBS). The data captures the surveys of 2010 and 2015 which were the first and second malaria indicator surveys conducted in Nigeria. The 2015 survey was put into action just a year after the 2010 survey and a year after the development of the new national malaria strategic plan that covers 2014–2020 [28]. The 2 years were used because this survey is usually carried out every 5 years. To find out about the risk of malaria or anaemia disease, two-stage sampling was carried out. Clusters were selected from each urban/rural strata in the first stage and systematic sampling were done for selection of households in the second stage. The data has 12,623 children under age 5 years old in total. 11,172 and 11,072 children were tested for malaria and anaemia out of 12,623 respectively. Finally, the sample size of 9,533 was used for analysis after removing the missing values and this is 75.5% of the original data. Figure 1 shows the map of Nigeria comprising 6 geopolitical zones and their 37 states including the capital territory. It further reveals the location of Nigeria in Africa. While the maps in Figs. 2 and 3 indicate the prevalence of malaria and anaemia in each state.

Fig. 1
figure 1

Location map of Nigeria showing the 6 geopolitical zones and their 37 states including the capital territory

Fig. 2
figure 2

Map of Nigeria showing state rates based on sampling weights of under 5 years old malaria prevalence

Fig. 3
figure 3

Map of Nigeria showing state rates based on sampling weights of under 5 years old anaemia prevalence

In this study, the first dependent variable is the binary response from a child’s RDT outcome while the second dependent variable is the binary response to the anaemic status of a child. For both dependent variables, 1 represents the presence of malaria or anaemia infection and 0 represents no presence of malaria or anaemia infection. The independent variables are the type of place of residence, source of drinking water, type of toilet facility, has electricity, has radio, has television, main floor material, main wall material, main roof material, wealth index, child’s age in months, sex, mother’s highest educational level, state, and region. Table 1 contains the summary of all the variables used in this work.

Table 1 Exploratory data analysis


To begin with, a space–time model was proposed by Bernardinelli et al. [25] using Poisson distribution and the log risk ratio was defined as a linear function of time. Authors expressed the log risk ratio for area \(i\); \(i=1,\dots ,I\) for time \(t\); \(t=1,\dots ,T\) as

$$log\left({\theta }_{it}\right)={\eta }_{it}=\mu +{u}_{i}+{v}_{i}+\left(\beta +{\delta }_{i}\right)\times t.$$

Following Besag et al. [24] specification, \(\phi ={u}_{i}+{v}_{i}\) are the spatial random effect (convolution model), where \({u}_{i}\)’s are the structured variables and \({v}_{i}\)’s are the unstructured variables. \(\mu\) is the overall mean, \(\beta\) is the universal linear time trend effect and \({\delta }_{i}\) is the random effect for interaction betwixt space and time. For the data to show the time trend, the parameters for the time trend were allocated unclear priors. In terms of the interaction random effect \({\delta }_{i}\), the independent and identically distributed (i.i.d) Gaussian prior was adopted, though alternative prior specification can be given. Based on spatial models, the priors for unstructured and structured spatial effects were described.

On the other hand, Knorr‐Held [26] reformed the earlier method by disabling the parametric limitations. The authors adopted a Binomial distribution for the number of cases in the country \(i\) \(\left(i=1,\dots ,I\right)\) at \(t\) th time \(\left(t=1,\dots ,T\right)\), while the log odds are expressed as:

$$log\left(\frac{{\pi }_{it}}{1-{\pi }_{it}}\right)={\eta }_{it}=\mu +{u}_{i}+{\upsilon }_{i}+{\gamma }_{t}+{\nu }_{t}+{\delta }_{it}$$

where \({\gamma }_{t}\), \({\nu }_{t}\) refers to the temporal random effects that take care of unnamed attributes of year \(t\), and \({\delta }_{it}\); the interaction effects that take care of differences not described by the main effects. Intrinsic conditional autoregression (iCAR) and first order random walk structure were assigned to \({u}_{i}\) and \({\gamma }_{t}\), while independent Gaussian priors were assigned to \({\upsilon }_{i}\) and \({\nu }_{t}\). Based on the temporal effects interaction and spatial effects, the interaction \({\delta }_{it}\) was presumed to have four forms of prior inference.

Bernardinelli et al. [25] and Knorr‐Held [26] performed their parameter estimation under the fully Bayesian approach with the use of Markov chain Monte Carlo (MCMC) through Gibbs sampling techniques. Here, INLA approximation to fully Bayesian estimation was used. Therefore, the method used in this study is discussed as follows.

Let \({y}_{iky}\) be malaria or anaemia status of child \(k\) in state \(i\): \(i=1,\dots ,37\) during year \(y\): \(y=\mathrm{1,2}\). The response outcome variable is a binary response, and it is defined as:

$$Y_{{iky}} = \left\{ {\begin{array}{*{20}c} {1,{\text{Anaemia}}} \\ {0,{\text{No}}\,{\text{anaemia}}} \\ \end{array} } \right.$$

where \({y}_{iky}\) is the binary response outcome and it follows a Bernoulli distribution as.

\({y}_{iky}\sim Bernoulli\left({\theta }_{iky}\right)\), where \({\theta }_{iky}\) are unspecified probabilities associated to the outcome probabilities of the models. The logistic regression model is expressed as;

$$logit\left({\theta }_{iky}\right)={\beta }_{0}+{\eta }_{iky}$$

where \({\beta }_{0}\) is the model intercept and the linear predictor \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta\) with covariate vector \(\mathrm{\rm X}={\left({x}_{iky1},\dots ,{x}_{ikyq}\right)}^{^{\prime}},\) \(\beta =\left({\beta }_{1},\dots ,{\beta }_{q}\right)\) is the vector regression coefficient. We employed the combined formulation of the structured additive regression to permit flexibility where the classical predictor can be expanded to a better flexible additive predictor. Therefore, the structured additive predictor is expanded to spatio-temporal modelling as

$${\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{spat}\left({s}_{i}\right)+{f}_{year}\left(y\right)+{f}_{iy}\left({s}_{i},y\right)$$

where \({f}_{spat}\), \({f}_{year}\), and \({f}_{iy}\) are respectively functions suitable for space, year and space-year interaction. The spatial components \({f}_{spat}\) are disintegrated into two i.e., spatially unstructured \({f}_{unstr}\) and spatially structured \({f}_{str}\) effects. However, \({f}_{year}\) show the random year effects and this is modelled as a first-order random walk or AR(1) according to [29]. While \({f}_{iy}\left({s}_{i},y\right)\) is a space-year interaction.

To ascertain the wellness of the estimators, these seven models were compared as follows.

Model 1: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)\)

Model 2: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)+{\beta }_{y}\)

Model 3: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)+{f}_{year} \left(y\right)\)

Model 4: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)+{f}_{iy}\left({s}_{i},y\right)\)

Model 5: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)+{\beta }_{y}+{f}_{iy}\left({s}_{i},y\right)\)

Model 6: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)+{f}_{year} \left(y\right)+{f}_{iy}\left({s}_{i},y\right)\)

Model 7: \({\eta }_{iky}={\mathrm{\rm X}}_{iky}^{^{\prime}}\beta +{f}_{str}\left({s}_{i}\right)+{f}_{unstr}\left({s}_{i}\right)+{f}_{1year} \left(y\right)+{f}_{iy}\left({s}_{i},y\right)\)


  • \({\mathrm{\rm X}}_{iky}\) signifies the vector of categorical variables effects for child \(k\) in state \(i\) during year \(y\)

  • \(\beta\) is a vector of regression coefficients

  • \({\beta }_{y}\) means the year-specific fixed effects

  • \({f}_{unstr}\left({s}_{i}\right)\) and \({f}_{str}\left({s}_{i}\right)\) are the unstructured and structured random effects respectively

  • \({f}_{year}\) and \({f}_{1year}\) show the smooth functions of the temporal random effects

  • \({f}_{iy}\left({s}_{i},y\right)\) signifies the spatial-year interaction effect

Model 1 simply takes care of the spatially structured random effects, and this accounts for unobserved significant factors that change spatially transversely over the states and spatially unstructured random effects that caters for undetected variables inside the states. Hereafter, by assuming a categorical variable, it will have a linear effect on malaria and anaemia. The temporal effect is not assumed by this model. Model 2 follow the same pattern but in addition assumes a linear year trend taken by \({\beta }_{y}\). On the other hand, Model 3 contains separable space and year random effect which takes care of the linear effect of categorical variables. Also, Model 4 and Model 1 are parallel to each other but additionally, Model 4 takes care of space and year interaction which captures differences that is not shown by the main effects. Regarding Model 5, the assumption is made on the linear effects of the categorical variables, spatial random effects, linear year trend and space year-interaction. While Model 6 and Model 7 are alike but vary in prior assumptions of the temporal random year effects (\({f}_{year}\), \({f}_{1year}\)). Additionally, both models take on linear effects of categorical variables, spatial random effects of the location, space, and year interaction. In other words, all models take on linear effects of categorical variables through the term \({\mathrm{\rm X}}_{iky}^{^{\prime}}\).

Prior specifications

Here, for the spatio-temporal logistic regression models, the full Bayesian approach was adopted. Diffuse priors were allocated to fixed effects and linear year trend, intrinsic conditional autoregressive (iCAR) was used to model the spatially structured random effects, while the independent and identically distributed (i.i.d) Gaussian prior was assigned spatially unstructured random effects. A first-order random walk was used to model the temporal year random effects \({f}_{1year}\). Nevertheless, it is interesting to note that varied prior specifications for the temporally changeable year effects \({f}_{year}\) were given in the models and penalized spline was given to the spatio-temporal logistic regression model. Also, independent penalized splines for the logistic independent first-order autoregressive model were adopted to model the spatial year-specific effects (interaction).

Parameter estimation

In this research, we discussed the procedure of parameter estimation of the spatio-temporal logistic regression model of malaria and anaemia using the fully Bayesian approach. With due consideration, every unspecified parameter assumed random variables and was given adequate prior distributions. The posterior of the priors is given as:

$$p\left(\varphi ,\psi |y\right)\propto L\left(y|\varphi ,\psi \right)p\left(\varphi ,\psi \right)$$

where \(L\left(y|\varphi ,\psi \right)\) is the likelihood of the penalized spline and \(p\left(\varphi ,\psi \right)\) are the prior distributions of the spatio-temporal logistic regression model. The latent Gaussian field is expressed as \(\varphi =\left\{\left\{\beta \right\},\left\{{\beta }_{y}\right\},\left\{{f}_{str}\left(.\right)\right\},\left\{{f}_{unstr}\left(.\right)\right\},\left\{{f}_{year}\left(.\right)\right\},\left\{{f}_{1year}\left(.\right)\right\},\left\{{f}_{iy}\left(.\right)\right\}\right\}\) while the equivalent hyperparameters are shown as \(\psi =\left\{{\varrho }_{str},{\varrho }_{unstr},{\varrho }_{year},{\varrho }_{1year},{\varrho }_{iy}\right\}\). Conjugate gamma priors \(Gamma\left(1, 0.00005\right)\) were allocated to all hyperparameters while R-integrated nested Laplace approximation (INLA) package was used to estimate the parameters.


Table 2 shows the model fit values for the spatio-temporal logistic regression models of malaria which comprises Deviance Information Criteria (DIC) and the effective number of parameters (DP). Model 7 was chosen as a better model because it gave the least DIC value of 10819.24. Therefore, the presentation of results and interpretations are based on Model 7 which includes both linear and nonlinear effects as well as the spatio-temporal effects.

Table 2 Summary of the model comparisons of malaria

Table 3 provides the adjusted posterior odds ratios estimates (AOR) and 95% confidence interval for the best fitting model mentioned above. Here, the results for significant covariates were discussed alone. Regarding child’s age in months, the results showed an increase in the odds of malaria for all ages. In the same vein, there was a significant increase in the odds of malaria among children whose anaemic status is positive and those who live in rural area. On the other hand, odds of malaria decrease significantly with respect to wealth index and mother’s educational level (Secondary and Higher). Furthermore, household with electricity had significantly lower odds of malaria compared to a household without electricity.

Table 3 Adjusted posterior odd ratios estimates (AOR) of malaria with 95% confidence interval

Figure 4 presents the mapped estimated residual spatial effects of the year 2010 and year 2015. The essence is to study how disease prevalence and risk factors change with time. These maps show unobserved spatial factors that are not captured in the survey or that capture the effects of cultural patterns. From the figure, there was an obvious spatial pattern change over the two years. Though higher concentrations in the two years are scattered, the states inside the Northeast, Northwest and Southwest regions have higher odds of 0.99–2.2 of malaria. In 2010, states in the north-east had higher odds of malaria but had lower odds of malaria later in 2015. .

Fig. 4
figure 4

Maps displaying residual spatial effects of malaria in Nigeria for year 2010 and 2015 obtained from spatio-temporal interaction logistic regression model, i.e., Model 7

Table 4 contains the results of \(DIC\) and the effective number of parameters \(DP\) of the model fit values for the spatio-temporal logistic regression models of anaemia.  From the summary, Model 6 gave the least DIC value of 10330.78. Therefore, the interpretation and presentation of results were based on this model.

Table 4 Summary of the model comparisons of anaemia

Figure 5 shows the posterior relative risk of malaria. Here, there was an obvious change in the relative risk of malaria over the two years. This implies that there has been an increase in the relative risk of malaria from 2010 to 2015.

Fig. 5
figure 5

Depicting estimated posterior relative risk of malaria for the logistic regression best fitting model

Table 5 presents the adjusted posterior odds ratios estimates (AOR) and 95% confidence interval for the best fitting model mentioned above. Discussion on the results was basically on the significant covariates. There was a significant decrease in odds of anaemia for children in age categories 3, 4, 5, and 6 but insignificant for  age category 2. Also, the odds of anaemia decrease significantly with increasing mother’s educational level. Furthermore, sex, household that has radio and television had significantly lower odds of anaemia. On the other hand, child’s age in months, the results showed that there was an increase in the odds of malaria but only significant for children in age categories 2, 3 and 6. While the odds of anaemia increased significantly for malaria rapid test results, source of drinking water, wealth index and type of place of residence.

Table 5 Adjusted posterior odd ratios estimates (AOR) of anaemia with 95% confidence interval

Figure 6 is the graphical representation of child’s age in months and the adjusted odd ratio (AOR) of malaria and anaemia. The relationship between child’s age in months and the adjusted odd ratio of malaria increase significantly across the age group. While anaemia decreases significantly. This might be due to the stimulation of antimalarial immune defenses by malaria antigen in breast milk which reduces malaria risk in infants attributable to breastfeeding [30]. As children are weaned, they are vulnerable to malaria as they have lost maternal immunity and are yet to develop self-immunity against infection [31]. As a result, there will be a decrease in anaemia infection because individuals at some point develop a disease-controlling immunity that makes them asymptomatic [32].

Fig. 6
figure 6

Relationship between child’s age in months and AOR of malaria and anaemia

Figure 7 displays the mapped estimated adjusted posterior odds of residual spatial effects of the year 2010 and 2015 in Nigeria. In our first work, we focused on the evolution of the geographical variation of anaemia. The map in Fig. 4 represents the estimated residual spatial effects for the two years. The colours for the regions are the same as described above. In both years, the South-South and partially South-West regions showed a higher concentration of anaemia with odds ratio of 0.90–2.2. Notwithstanding, in 2010, some part of North-Central region had a higher concentration of the same odds ratio. Other things seem similar to the previous figure.

Fig. 7
figure 7

Maps displaying residual spatial effects of anaemia in Nigeria for year 2010 and 2015 sprang from the spatio-temporal interaction logistic regression model i.e., Model 6

Figure 8 provides the estimated posterior relative risk of anaemia. There was an increase in the relative risk of anaemia in the two years i.e., the transmission of anaemia infection was high in 2015.

Fig. 8
figure 8

Depicting estimated posterior relative risk of anaemia for the logistic regression best fitting model


This research work applied spatio-temporal models to investigate the relative risks and geographical variation of malaria and anaemia in Nigeria. This research work was carried out with the sole aim of developing and applying the exact statistical models to assess determinant factors and geographical distinctions of malaria and anaemia. In addition, to apply a unified framework of flexible models within Bayesian hierarchical modelling to understand various factors associated with this discreet type of malaria and anaemia prevalence among children from 0 to 59 months in Nigeria. The models considered are an augmentation to classical models which include spatial and spatio-temporal models for identification of geographical variation of year-specific effects. Logistic regression was developed to assess influential factors and state variation of malaria and anaemia prevalence. The structured additive modelling approach gives allowance for different kinds of predictors to be included in classical models in an additive manner by borrowing strength from both parametric and non-parametric models. Integrated nested Laplace approximation was used to investigate the spatio-temporal effect on childhood malaria and anaemia disease with the application of MIS (Malaria indicator survey) datasets in Nigeria. For each model, the Deviance Information Criteria (DIC) were compared, and the best model was used to fit malaria and anaemia data of Nigeria. Among the models considered, the spatio-temporal interaction logistic regression model was chosen as the best model to fit malaria data while for anaemia, model 6 (interaction with one random time effect (autoregressive prior of order 1(AR1))) was chosen. For both diseases, variation can be seen among the Nigeria states and clustering among states with high malaria and anaemia relative risk (RR). Child’s age in month, main wall material, anaemic status, wealth index, mother’s educational level and type of place of residence were significantly related to malaria and anaemia except for the source of drinking water, sex, and household that has radio and television that were significantly related to only anaemia over the two years period, this is in line with [33, 34]. While anaemia is seen as a dominant determinant of malaria, malaria as well as a major determinant of anaemia[35, 36]. In the two years considered, we found out that the states within the northern, southern, and western regions have the higher prevalence of malaria and anaemia. The lowest prevalence of malaria and anaemia was seen in states within the eastern region. Also, we estimated year temporal effects on malaria and anaemia. For both diseases, the plots showed no obvious change in the spread of malaria and anaemia.

There are always limitations in every research. In this study, the main limitation is the number of years available for us to estimate the spatio-temporal trends of malaria or anaemia. This matter poses a hinderance to investigating the trend of malaria or anaemia pandemic during the early years. Also, this study used secondary data from cross-sectional surveys which did not allow the causal relationships to be established. Furthermore, though iron deficiency is one of the major causes of anaemia, there was no information on iron levels in children. Notwithstanding the limitations, the strength of this study is in the use of individual-level malaria RDT results instead of indicators or estimates of malaria or anaemia.


There is room for further investigation on this research work. In this work, FB approach within the Bayesian hierarchical modelling to model malaria and anaemia prevalence in Nigeria. Furthermore, the Bayesian structured additive approach was used to model the determinants of malaria and anaemia. The findings from this work show that there will be likely a reduction in the spread of these diseases if commendations are adequately adhered to. Therefore, government should focus on improving mother’s education and standard of living. Also, pertaining to these diseases, there should be aggressive awareness on social television programs.

Availability of data and materials

Data is available in Demographic Health Survey. Public access to databases is closed (requires a password) ( We received permission from DHS to use the data via a password supplied by DHS Admin. (Authorization letter is attached as supplementary material).


  1. WHO. World malaria report 2021. Geneva: World Health Organization; 2021.

    Google Scholar 

  2. Steketee RW, Choi M, Linn A, Florey L, Murphy M, Panjabi R. World Malaria Day 2021: commemorating 15 years of contribution by the United States President’s Malaria Initiative. Am J Trop Med Hyg. 2021;104:1955–9.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Rahi M, Das P, Sharma A. Malaria elimination in India requires additional surveillance mechanisms. J Public Health (Oxf). 2022;44:527–31.

    Article  Google Scholar 

  4. Douglas NM, Burkot TR, Price RN. Malaria eradication revisited. Int J Epidemiol. 2022;51:382–92.

    Article  PubMed  Google Scholar 

  5. Ahmad F, Goparaju L, Qayum A. Studying malaria epidemic for vulnerability zones: multi-criteria approach of geospatial tools. J Geosci Environ Protect. 2017;5:30–53.

    Article  Google Scholar 

  6. Ra PK, Nathawat MS, Onagh M. Application of multiple linear regression model through GIS and remote sensing for malaria mapping in Varanasi District, India. Health Sci J. 2012;6:731.

    Google Scholar 

  7. Childs D, Cattadori I, Suwonkerd W, Prajakwong S, Boots M. Spatiotemporal patterns of malaria incidence in northern Thailand. Trans R Soc Trop Med Hyg. 2006;100:623–31.

    Article  PubMed  CAS  Google Scholar 

  8. Mabaso ML, Vounatsou P, Midzi S, Da Silva J, Smith T. Spatio-temporal analysis of the role of climate in inter-annual variation of malaria incidence in Zimbabwe. Int J Health Geogr. 2006;5:20.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Clements AC, Barnett AG, Cheng ZW, Snow RW, Zhou HN. Space-time variation of malaria incidence in Yunnan province, China. Malar J. 2009;8:180.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Jana S, Fu SH, Gelband H, Brown P. Jha P99. Malar J. 2022;21:90.

    Article  PubMed  PubMed Central  Google Scholar 

  11. WHO. Nutrional anaemias: tools for effective prevention and control. Geneva, World Health Organization, 2017. Accessed 24 Feb 2022.

  12. Hailegebreal S, Nigatu AM, Mekonnen ZA, Endehabtu BF. Spatio-temporal distribution and associated factors of anaemia among children aged 6–59 months in Ethiopia: a spatial and multilevel analysis based on the EDHS 2005–2016. BMJ Open. 2021;11: e045544.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Kassebaum NJ, Collaborators GA. The global burden of anemia. Hematol Oncol Clin North Am. 2016;30:247–308.

    Article  PubMed  Google Scholar 

  14. Leal LP, Batista Filho M, de Lira PIC, Figueiroa JN, Osório MM. Temporal trends and anaemia-associated factors in 6-to 59-month-old children in Northeast Brazil. Public Health Nutr. 2012;15:1645–52.

    Article  PubMed  Google Scholar 

  15. Sirait M, Flora R, Anwar C, Zulkarnain M, Fajar NA, Faisya AF. Risk factors for the incidence of anemia in elementary-school children living in malaria-endemic regions. 2022;25:235-244.

  16. Okunlola OA, Oyeyemi OT. Spatio-temporal analysis of association between incidence of malaria and environmental predictors of malaria transmission in Nigeria. Sci Rep. 2019;9:17500.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Okunlola OA, Oyeyemi OT, Lukman AF. Modelling the relationship between malaria prevalence and insecticide-treated bed net coverage in Nigeria using a Bayesian spatial generalized linear mixed model with a Leroux prior. Epidemiol Health. 2021;43: e2021041.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Akinyemi YC. Spatiotemporal patterns and determinants of reproductive health care utilization in Nigeria: 2008–2018. Papers Applied Geogr. 2021;7:453–81.

    Article  Google Scholar 

  19. Jaja C. Spatiotemporal patterns and determinants of reproductive health care utilization in Nigeria: 2008–2018. Papers Applied Geogr. 2021;7:1–29.

    Google Scholar 

  20. Borderon M, Ludwig A, Brazeau S. Vulnerable populations. In: Brazeau S, Ogden NH, editors. Earth observation, Public Health and One Health. Wallingford: Cabi Digital Library; 2022.

    Google Scholar 

  21. Starck T, Bulstra CA, Tinto H, Rouamba T, Sie A, Jaenisch T, Bärnighausen T. The effect of malaria on haemoglobin concentrations: a nationally representative household fixed-effects study of 17,599 children under 5 years of age in Burkina Faso. Malar J. 2021;20:416.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Ekvall H. Malaria and anemia. Curr Opin Hematol. 2003;10:108–14.

    Article  PubMed  Google Scholar 

  23. Newton C, Warn P, Winstanley P, Peshu N, Snow R, Pasvol G, et al. Severe anaemia in children living in a malaria endemic area of Kenya. Trop Med Int Health. 1997;2:165–78.

    Article  PubMed  CAS  Google Scholar 

  24. Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Statit Math. 1991;43:1–20.

    Article  Google Scholar 

  25. Bernardinelli L, Clayton D, Pascutto C, Montomoli C, Ghislandi M, Songini M. Bayesian analysis of space—time variation in disease risk. Stat Med. 1995;14:2433–43.

    Article  PubMed  CAS  Google Scholar 

  26. Knorr-Held L. Bayesian modelling of inseparable space-time variation in disease risk. Stat Med. 2000;19:2555–67.

    Article  PubMed  CAS  Google Scholar 

  27. Mdakane BP. Spatial and Spatio-temporal modelling and mapping of self-reported health among individuals between the ages of 15-49 years in South Africa. Thesis, University of Kwazulu-Natal, 2019

  28. NMCP. Nigeria MIS final report. Federal Republic of Nigeria. Abuja; 2015.

  29. DiMaggio C. Spatial epidemiology notes. Applications and Vignettes in R. 2014. Accessed 24 Feb 2022.

  30. Van Den Elsen LW, Verhasselt V, Egwang T. Malaria antigen shedding in the breast milk of mothers from a region with endemic malaria. JAMA Pediatr. 2020;174:297–8.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Schumacher R-F, Spinelli E. Malaria in children. Mediterr J Hematol Infect Dis. 2012;4: e2012073.

    Article  PubMed  PubMed Central  Google Scholar 

  32. White NJ. Anaemia and malaria. Malar J. 2018;17:371.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Ugwu CLJ, Zewotir T. Spatial heterogeneity and determinants of childhood anaemia in Nigeria. Res Sq. 2020.

    Article  Google Scholar 

  34. Roberts DJ, Matthews G, Snow RW, Zewotir T, Sartorius B. Investigating the spatial variation and risk factors of childhood anaemia in four sub-Saharan African countries. BMC Public Health. 2020;20:126.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Weze K, Abioye AI, Obiajunwa C, Omotayo M. Spatio-temporal trends in anaemia among pregnant women, adolescents and preschool children in sub-Saharan Africa. Public Health Nutr. 2021;24:3648–61.

    Article  PubMed  Google Scholar 

  36. Rumisha SF, Shayo EH, Mboera LE. Spatio-temporal prevalence of malaria and anaemia in relation to agro-ecosystems in Mvomero district, Tanzania. Malar J. 2019;18:228.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


The authors are grateful to the University of KwaZulu-Natal, Durban for operational and infrastructural support and DHS for granting the author access to NDHS 2010 and 2015 data.


No funding.

Author information

Authors and Affiliations



JUI: conceptualization, methodology, data curation, writing—original draft preparation; HM: conceptualization and supervision; A-KI: reviewing and editing, supervision. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jecinta U. Ibeji or Henry Mwambi.

Ethics declarations

Ethics approval and consent to participate

Secondary data was used for this study. The data is available for research purposes from the Demography and Health Survey website. Therefore, formal ethical approval is not applicable for this study. All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors do declare that there is no competing interest.

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Ibeji, J.U., Mwambi, H. & Iddrisu, AK. Bayesian spatio-temporal modelling and mapping of malaria and anaemia among children between 0 and 59 months in Nigeria. Malar J 21, 311 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Spatio-temporal
  • Heterogeneity
  • Bayesian Hierarchical
  • Deviance Information Criteria
  • Risk Ratio and R-integrated nested Laplace approximation (INLA)