Skip to main content


Assessing the effects of mosquito nets on malaria mortality using a space time model: a case study of Rufiji and Ifakara Health and Demographic Surveillance System sites in rural Tanzania

Article metrics



Although malaria decline has been observed in most sub-Saharan African countries, the disease still represents a significant public health burden in Tanzania. There are contradictions on the effect of ownership of at least one mosquito net at household on malaria mortality. This study presents a Bayesian modelling framework for the analysis of the effect of ownership of at least one mosquito net at household on malaria mortality with environmental factors as confounder variables.


The analysis used longitudinal data collected in Rufiji and Ifakara Health Demographic Surveillance System (HDSS) sites for the period of 1999–2011 and 2002–2012, respectively. Bayesian framework modelling approach using integrated nested laplace approximation (INLA) package in R software was used. The space time models were established to assess the effect of ownership of mosquito net on malaria mortality in 58 villages in the study area.


The results show that an increase of 10 % in ownership of mosquito nets at village level had an average of 5.2 % decrease inall age malaria deaths (IRR = 0.948, 95 % CI = 0.917, 0.977) in Rufiji HDSS and 12.1 % decrease in all age malaria deaths (IRR = 0.879, 95 % CI = 0.806, 0.959) in Ifakara HDSS. In children under 5 years, results show an average of 5.4 % decrease of malaria deaths (IRR = 0.946, 95 % CI = 0.909, 0.982) in Rufiji HDSS and 10 % decrease of malaria deaths (IRR = 0.899, 95 % CI = 0.816, 0.995) in Ifakara HDSS. Model comparison show that model with spatial and temporal random effects was the best fitting model compared to other models without spatial and temporal, and with spatial–temporal interaction effects.


This modelling framework is appropriate and provides useful approaches to understanding the effect of mosquito nets for targeting malaria control intervention. Furthermore, ownership of mosquito nets at household showed a significant impact on malaria mortality.


Evidence about the impact of mosquito net ownership on malaria mortality under routine conditions is limited and contradicting. Several studies have been conducted on the effect of this malaria control intervention [15]. Some studies used all cause child mortality as an endpoint [1, 2] while others have focused on indicators of malaria transmission [3, 6] and on changes in treatment outcomes [1, 3]. Other studies have shown a decrease in all-cause mortality associated with insecticide-treated nets (ITN)ownership [1, 2]. Another study found no difference in mortality rates between villages with and without ITNs intervention [5]. However, most of these studies have used a temporal model and only one study identified, on malaria incidence, that considered the spatial–temporal dependency [6]. Therefore, there is a need to assess the effect of mosquito net ownership on malaria mortality in Tanzania.

Insecticide-treated nets (ITNs) have been an effective intervention for preventing malaria transmission by reducing the chances that an individual would be bitten by an infective Anopheles mosquito [7, 8]. Several studies on ITN efficacy or effectiveness conducted have shown continued ability of ITNs protection against malaria transmission [911]. Moreover, studies on pregnancy outcome showed that use of bed nets could reduce miscarriages and stillbirths by 33 % and by 23 % in placental parasitaemia [12]. Overall, the use of ITNs have been effective [1, 7, 13]. Following results from studies on efficacy and effectiveness of ITN use in Tanzania, a market-based national ITN voucher programme to scale up the coverage of bed nets was adopted in 2004 [14]. This programme increased household ownership of at least one ITN from 23 to 92 % between 2004/2005 and 2012 in Tanzania mainland [15].

Long-term evaluation of the impact of malaria prevention and control on malaria mortality is lacking, posing a challenge to malaria endemic countries like Tanzania, where many deaths occur at home [16] where no records are available [17]. Demographic and Health Surveys (DHS) have been used as sources of mortality data but lack the ability to provide causes of deaths. Assumptions about the proportion of malaria deaths to all-cause mortality have been used [4, 18, 19].

Verbal autopsy (VA) is currently an alternative approach that has been used to determine malaria causes of deaths in community-based studies in many of the endemic countries where vital registration data are incomplete or do not exist [20, 21]. VA is an indirect method used in determining the causes of death based on an interview with next of kin or other caregivers. This is done using a standardized questionnaire that obtains information on signs, symptoms, medical history and circumstances preceding death [20]. Studies have evaluated the validity of VA for determining malaria-specific death [22, 23] in HDSS sites. The results concluded that VA methods have an acceptable level of diagnostic accuracy for determining malaria deaths at the community level.

In sub-Saharan African (SSA), VA has been routinely used as a tool to provide information on the burden of the disease using HDSS sites [24, 25]. HDSS has sites established within the INDEPTH network in many SSA countries which continuously collect large amounts of data. The data are collected on spatial and temporal scales and include cause specific deaths using VA, in and out migration, and other demographic information [26]. This study assesses the effect of ownership of mosquito nets on malaria mortality using a Bayesian space time model. Also, the study assessed how environmental factors such as rainfall, temperature and vegetation index affect the impact of mosquito net ownership.


Study area

The study was carried out using data from Rufiji and Ifakara Rural Health and Demographic Surveillance System (HDSS) sites [27, 28]. The sites are primarily rural with majority of the population relying on subsistence farming or fishing. The Rufiji HDSS covers 1813 km2 that comprises 33 villages of Rufiji District, Coast region. Ifakara HDSS covers an area of 2400 km2 across which covers 25 villages (13 in Kilombero and 12 in Ulanga districts) in Morogoro region south-eastern part of the country. These two HDSS sites were selected because they are among the HDSS sites, which continuously collect large amounts of data longitudinally on defined geographical areas for malaria specific cause of death in Tanzania. Furthermore, these HDSS sites continue to have high malaria prevalence. For example, microscopy estimates for health facilities in 2012 in Rufiji were 19.2 % in high season and 7.2 % in low season, and, in Ifakara they were 9.4 % in high season and 4.2 % in low season 4.2 % [29]. In these two HDSS, malaria is the major contributor among deaths in children below 5 years of age [30, 31].

Study design and data collection

The analysis used longitudinal data collected from Rufiji and Ifakara HDSS sites. Yearly malaria deaths were extracted from the Rufiji and Ifakara HDSS database covering a period of 1999–2011 and 2002–2012 respectively. The two HDSS sites have consistently been recording pregnancies, pregnancy outcomes, deaths and migrations by visiting households once every 4 months since 1997 in Ifakara HDSS and 1998 in Rufiji HDSS. Also, the HDSS sites collect information on socioeconomic status including household ownership of mosquito nets once per year during visit round of update. The mosquito nets in this study include untreated and treated nets (ITNs). The study period for this particular study covers the social marketing and national scale up of ITNs in the country [32]. The mosquito net ownership in HDSS sites has been increased over the years. In Rufiji HDSS the ownership of at least one mosquito net increased from 25 % of households in 1999 to 95 % in 2011 while in Ifakara HDSS the ownership of at least one mosquito net has been high coverage since 2001 (Fig. 1). The details on how data collection, management and credibility in HDSS are discussed in detail elsewhere [27, 28].

Fig. 1

Malaria mortality rates and percentage of households with at least one mosquito nets in Rufiji and Ifakara HDSS. a For Rufiji HDSS; b for Ifakara HDSS

Verbal autopsy procedure

The WHO and INDEPTH network [33] standardized VA questionnaire was adapted and used for data collection on causes of death. In the HDSS, deaths were captured during sequential rounds of data collection updates. Then the trained HDSS field interviewers visited the deceased’s home after a grieving period to administer a verbal autopsy questionnaire. An interview was administered to relatives or caregivers who were closely associated with the deceased during the period leading to his or her death. The questionnaire assessed the identity of the deceased and established the sequence of events leading to death, including symptoms and signs of the illness before death. Verbal autopsy was carried out since 1998 in Rufiji HDSS and 2002 in Ifakara HDSS. The verbal autopsy forms are independently reviewed by two physicians according to a list of causes of death based on the 10th revision of the international classification of diseases (ICD-10). A third physician is asked to code the cause of death in the case of discordant results. If there is disagreement among the three physicians, the death is coded as “undetermined” cause [34]. Causes of death (main, immediate, and/or contributing) are coded to be consistent with the ICD-10 [35]. Malaria deaths are coded as direct where malaria is the underlying cause of death or indirect where malaria is one of several diseases leading to death but where the death is attributed by a different cause [13]. In this analysis direct and indirect malaria deaths were included.

Confounding variables

The analysis used mean rainfall, mean temperature, mean normalized difference vegetation index (NDVI) and age as confounding variables in assessing the effect of mosquito net ownership on malaria mortality.

Climate data includes monthly mean rainfall and temperature (maximum and minimum) with these measures obtained from Tanzania Meteorological Authority (TMA) head office in Dar es Salaam. TMA provides meteorological services, weather forecasts, climate services, and warnings including daily forecast information for each region in Tanzania [36]. The climate data are collected through the different gauges located in different stations in each district. In Rufiji District, Utete, Ikwiriri and Kibiti weather station records climate data. Kibiti weather station data were used for Rufiji HDSS site and the Kilombero Agricultural Training and Research Institute (KATRIN) weather station data were used for villages in Ifakara HDSS for Kilombero district. Other villages used Mahenge weather station data for villages located in Ulanga district.

Remote sensing data,which included normalized difference vegetation index (NDVI)was extracted from moderate resolution imaging spectroradiometer (MODIS) on board NASA’s terra satellite. The vegetation index was processed from MODIS (MOD13A3) using monthly composite images at a 1 × 1 km resolution. Vegetation indices are used for global monitoring of vegetation conditions and are used in products displaying land cover and land cover changes. ERDAS Imagine software version 10.1 was used for processing the satellite images, and ArcGIS version 10 was used for spatial analyses. For each village within the study area, mean NDVI were calculated each year to link with malaria mortality for each village. The village boundaries were overlaid with the raster environmental data (i.e., MODIS) and calculated the weighted mean value for each village based on the proportion of the village area.

Data processing and analysis

The number of malaria deaths in each village were extracted from the HDSS database and aggregated yearly in order to link with mosquito net ownership which is also surveyed yearly in the HDSS sites. The percentages of household’s mosquito net ownership at village level were calculated each year. Time at risk (person-years) contributed by each person was calculated from 1st January for each year until exit. Exit from the study was due to migration (outside the HDSS area), death or end of the study. In a case where a person migrated to a different household location within the study area, time at risk was computed separately for new location and added to the total time at risk. The outcome of interest is total yearly death due to malaria at village level for specific age groups. The age was categorized into under five and five and above. The outcome variable was aggregated at village level to capture heterogeneity for mortality in the study area [37, 38]. The malaria mortality rates were calculated by dividing the number of deaths by the person-years of observation and were expressed per 1000 person-years (py). Figures 2, 3, 4, and 5 display malaria mortality by village in each HDSS site.

Fig. 2

Annual malaria mortality rate trends for under-five by village in Rufiji HDSS

Monthly environmental and climate data were aggregated and averaged to obtain mean values to link with the mosquito net ownership and malaria death data. Mean annual rainfall (millimetre), mean annual temperature (°C), and mean annual NDVI were used as confounder variables because they are associated with malaria transmission [39, 40]. The percentage of households ownership at least one mosquito netat household was used also as main explanatory variable.

In order to highlight spatio-temporal trends of malaria mortality in the study areas, direct estimates of malaria mortality rates were calculated. The malaria mortality rates were calculated by dividing the number of deaths by the person-years of observation and were expressed per 1000 person-years (py).

Space time hierarchical model

This study propose a spatial and temporal hierarchical model that considers the effects of predictors variables as described in detail [41]. These type of models usually have to model data, processes, and parameters [42]. The study used space time models that includes spatial and time random effects as described [43]. Covariate informations considered in this study for village i at year j were mosquito net ownership as exposure variable and confounder variables (rainfall, temperature and NDVI). For this study i the number of village in the study, I = 33 for Rufiji and 25 for Ifakara HDSS and j for years, J = 13 for Rufiji HDSS and J = 11 for Ifakara HDSS.

The processes are modelled by common specifications where both spatial component and temporal random time effects as first autoregression AR (1) are considered (see Additional file 1). In a broad view, random effects specifications for the spatial random terms include the unstructured noise term that follows a normal distribution and structure term which modelled with the set of villages that adjacent to other villages by weight of neighboring village. Weight is considered to be 1 when villages i and j share the same boundary and 0 otherwise. The prior distribution of the random effect term is a conditional autoregressive prior used to model the spatial dependence. The temporal random effect determines the prior distribution consider a first order auto-regression AR (1) given by the partial autocorrelation function. The distribution of the interaction term is characterized by a precision matrix obtained as the Kronecker product of the precisions of unstructured and structured noise terms [44].

Ownership of at least one mosquito net was considered as the major malaria control strategy in this study. Additional confounding variables incorporated in the space time model included annual mean rainfall, and average annual temperature, and mean annual NDVI. Models considered included a non-spatial and temporal effect model (M 1 ), a model with spatial and temporal random terms (M 2 ) and a model with spatial and temporal random terms, and interaction (M 3 ) (see Additional file 1). Model (M 2 ) and (M 3 ) included spatial random effects which arising from a Gaussian stationary process with covariance matrix capturing correlation between any pair of villages as a function of their distance. Additional models were considered to check for multi-colinearity using variance inflation factors for environmental variables. The model was checked for statistical interactions between mosquito net ownership and confounder variables.

In this analysis, R programming language (R Development Core Team, 2009) using the INLA package was used to fit three models defined M 1 to M 3. Latent models available in INLA package such as Besag model [45], an independent random noise (iid) model, and first order auto-regression AR (1) [46] were used to fit models to data. Besag model is frequently used approach to spatial disease modelling in small areas [45]. INLA package is a computational approach used to approximate Bayesian inference based on an efficient combination of Laplace approximations and numerical integration. Unlike Markov Chain Monte Carlo (MCMC), the INLA method does not sample from the posterior distribution [47]. Posterior distributions are approximated with a closed form expressions, reducing convergence and mixing problems. INLA is suitable for Bayesian hierarchical models for large number of unknown parameters that follow a Gaussian Markov random field and a small number of hyper-parameters.

Model selection

Considering the dynamic nature of malaria, an important issue was whether ownership of mosquito nets has an impact on malaria death. Similar to Arab et al. [30], deviance information criterion (DIC) was used to assess model goodness of fit and selection for best model that produced lower values of the DIC. DIC is a generalization of the Akaike and Bayesian information criterion [48]. Model estimates were then exponentiated to obtain the incidence rate ratios (IRR) for each predictor.

Ethical approval

The Ethical clearance was granted by the Ifakara Health Institute (IHI)’s Institutional Review Board (IRB), Tanzania and Medical Research Coordinating Committee (MRCC) of the National Institute for Medical Research (NIMR) for the establishment of Rufiji and Ifakara HDSS. For each household visit, verbal consent was sought from the respondent.


Table 1 shows the annual malaria mortality rate in Rufiji and Ifakara HDSS. In Rufiji HDSS, malaria mortality rate decreased from 11.2 per 1000 person-years in 1999 to 5.3 per 1000 person-years in 2011 for under five children age, and from 3.5 per person-years in 1999 to 2.2 per 1000 person years in 2011 for all ages. In Ifakara HDSS, results show the decrease in malaria mortality rate from 6.9 per 1000 person-years in 2009 to 4.0 per 1000 person-years in 2012 for under five children age. Furthermore, the malaria mortality rate for all age decreased from 2.3 per person-years in 2008 to 1.2 per 1000 person years in 2012. However, there was a delay in the decrease in mortality due to malaria rate for both under five children and all age in Ifakara HDSS as shown in Table 1.

Table 1 Malaria mortality trend and environmental factors in Rufiji and Ifakara HDSS

Figure 1 indicates the percentages of households that owned at least one mosquito net and malaria mortality rate in the two HDSS sites. The percentages have been increased while the malaria mortality rate decreased in Rufiji HDSS. Ifakara HDSS has high percentage of households with at least one mosquito net ownerships since 2002.

Malaria mortality rate by village in the selected villages in HDSS sites are shown in Figs. 2, 3, 4 and 5. The results show that in most of the villages, malaria mortality rates have been decreasing and few villages show high malaria mortality for exception years. Malaria mortality rate for under-five children were high in Mangwi village in 2005 and 2011, Nyamwimbe in 2007 and 2011 compared to other years and villages (Fig. 2). In 2005 and 2007 there was high malaria mortality in Mangwi and Nyamwimbe in 2011 for all age (Fig. 3). In Ifakara, HDSS, malaria mortality declined from 2008 to 2012 for both all age and under-five children. High malaria mortality for under five was observed in Nakafuru in 2008 and Igota village in 2011 (Fig. 4). In addition, higher malaria mortality for all age was observed in 2004 in Miwangani village and Mpofu village in 2005 (Fig. 5).

Fig. 3

Annual malaria mortality rate trends for all age by village in Rufiji HDSS

Fig. 4

Annual malaria mortality rate trends for under-five by village in Ifakara HDSS

Fig. 5

Annual malaria mortality rate trends for all age by village in Ifakara HDSS

Three space time models were implemented for malaria mortality and ownership of mosquito nets and environmental factors. Table 2 shows the deviance information criterion for fitted three models in Rufiji HDSS. Based on DIC values, the model with spatial and temporal random effect terms performed best among the models considered (DIC = 1880. 076), followed by the model with spatial and temporal interaction with DIC = 1902.217. The model without spatial and temporal effect had the worst performance (DIC = 1933.614). Also these models were fitted to malaria mortality for under five children, model with spatial and temporal random effect performed better compared to other two models (DIC = 1483.748). Based on these results, it is clearly observed that including spatial and temporal random effect terms is necessary.

Table 2 Adjusted estimated effect of ownership of mosquito net on malaria mortality in three models for Rufiji HDSS

Table 2 also shows that mosquito net ownership at household was associated with malaria mortality in Rufiji HDSS. The results in Table 2 show that an increase of 10 % ownership of at least one mosquito net at households in village is associated with (IRR = 0.948, 95 % CI: 0.917–0.977) 5 % decrease of malaria risk for all age per year and (IRR = 0.946, 95 % CI: 0.909–0.982) 5 % for under five children per year. The results show an increase of 100 mm for rainfall corresponds approximately to 6.6 % (IRR = 1.066, 95 % CI: 1.031–1.107) increase of malaria risk for all age per year. For under-five children, increase of 100 mm for rainfall corresponds to increase of (IRR = 1.056, 95 % CI: 1.031–1.107) 5. 6 % in malaria risk for under-five children per year.

Table 3 shows incidence rate ratio for the regression coefficients in three models for Ifakara HDSS with DIC. The model with spatial and temporal random effect terms performed best among the models considered (DIC = 1226.043) and (DIC = 1047.046) for all age and under five rates, respectively. Model comparison showed that model with spatial and temporal random effect had a small DIC value and therefore was the best fitted model.

Table 3 Adjusted estimated effect of ownership of mosquito net on malaria mortality in three models for Ifakara HDSS

The results show that ownership of at least one mosquito net at household was associated with malaria mortality. In Table 3, the results show that an increase of 10 % ownership of at least one mosquito net at household in village is associated with (IRR = 0.879, 95 % CI: 0.806–0.959) 12 % and (IRR = 0.899, 95 % CI: 0.816–0.995) 10 % decrease of malaria risk for all age and under-five children per year respectively. Also for an increase of 100 millimetres for rainfall there is approximately (IRR = 1.008, 95 % CI: 0.986–1.030) 1 % and (IRR = 1.018, 95 % CI: 0.988–1.039) 1.3 % increase of malaria risk for all age and under-five children per year, respectively. All models (Table 3) indicated mean annual temperature, and mean NDVI were positive predictors of increased annual malaria mortality in contrast to ownership of mosquito nets which had a reducing effect. The model without spatial and temporal effect underestimated the malaria mortality risk. The effect of mosquito net ownership was not affected by the effect of environmental factors (rainfall, temperature and mean NDVI). There was no statistical significance for interaction of mosquito net ownership and environmental factors in the study areas.

The posterior mean risks, represent smoothed values of the raw standardized malaria mortality ratios and, therefore, give a smoothed picture of what is going on with the distribution of malaria mortality in the study area. Model two with spatial and temporal random effect used to predict year 2011 for Rufiji HDSS and 2012 for Ifakara HDSS. The observed malaria mortality and predicted malaria mortality for the 2011 in Rufiji HDSS and 2012 in Ifakara HDSS respectively using model with spatial and temporal random terms (see Additional files 2, 3). The results from each of the Additional files is that the villages with high observed malaria mortality tend to have a high predicted malaria mortality risk; this is consistent with the results in Figs. 2, 3, 4, and 5 for malaria mortality rate.


The malaria mortality rate in the study area decreased by 52.7 and 37.1 % for under-five and all age respectively between 1999 and 2011 in Rufiji HDSS. Also, malaria mortality rate decreased from 2009 in Ifakara HDSS. The space time model suggests that mosquito net ownership is a protective risk for malaria mortality in the study area. The findings show that a 10 % increase in the ownership of at least one mosquito net in the study areas was associated with 6 and 10 % reduction of malaria mortality for under-five children and 5 and 12 % for all ages in Rufiji and Ifakara HDSS, respectively. This result is in line with other previous studies [1, 49] in the same areas which investigated the impact of mosquito net ownership to all cause child mortality. This analysis adds to the existing literature by providing evidence of the effect of mosquito net ownership on malaria mortality in small scale geographic areas; previous studies were predominantly in all cause mortality. The effect of mosquito net ownership in this study was higher in Ifakara HDSS (about 12 % decrease) compared to Rufiji HDSS (5 %). The possible explanation for this is that a number of key malaria control interventions have been implemented in Ifakara HDSS to increase ITN coverage. For example, the KINET project was piloted to promote and subsidize ITNs between 1997 and 1999 [50]. This was also followed subsequently by scale-up nationwide within the frame of the Tanzania National Voucher Scheme (TNVS) [51].

The findings from this study support the continued scale-up of mosquito nets to increase coverage within households in Tanzania. This is an important measure to prevent malaria and is one of the core strategies of World Health Organization’s Roll Back Malaria programme [52]. It also emphasizes the importance of ongoing and future efforts to maintain mosquito net coverage to vulnerable groups and the population as a whole. Furthermore, it also suggests that the massive effort to scale up mosquito net coverage over the past decade has paid off with measurable impacts and that it is possible for health systems to increase coverage of interventions and affect health outcomes over a relatively short period of time. Continued coordinated efforts between local and national governments, international organizations, funding agencies, and researchers are needed to make sure that mosquito nets are reaching all populations at risk for malaria. With the relatively large impact of mosquito nets on child mortality, The findings of this study also support the continued emphasis on malaria control, including the push towards malaria elimination, as a way of improving reduction of malaria deaths at the household and country level.

The findings show an association between rainfall and malaria mortality in under-five age and all ages. The findings of this study are similar to previous assessments of weather related mortality [1, 53, 54]. Rainfall is widely considered to be a major driver of inter-annual variability of malaria incidence in Africa. About 90 % of the deaths that occurred in sub-Saharan Africa are believed to be due to malaria [55]. The possible explanation for this is that rainfall provides breeding sites for mosquitoes and increases the humidity which enhances their survival and therefore increases the spread of the disease [56]. In addition, excess rainfall is known to create considerable flooding which limits access to care and increases the risk of other diseases such as diarrhoea, which in turn increases physical vulnerability and risk of dying from malaria infection, especially in pregnant women and in children under five [57].

Findings of this study indicate that models with spatial and temporal random effect terms performed better for goodness of fit and the effect of mosquito net ownership was higher compared to other models. This shows the necessity for inclusion of spatial and temporal terms in the Bayesian model framework; which is supported by other studies [5860] that also described the potential use of spatial and temporal terms in the model. The possible explanation is because of the complex dependence patterns over space and over time of the occurrence of malaria deaths, and the inherent large stochastic variability due to rare events. Estimating time trends in each area separately is not efficient enough because it will be difficult to establish a baseline pattern separately for each area.

Strengths and limitations

This study used datasets from Health and Demographic Surveillance System sites which continuously registered vital demographic events in a geographical defined area. The study findings from health and demographic surveillance systems data provides information to policy makers and program managers which can be translated into policy and practice for targeting malaria control interventions. Although, this study used VA data collected in HDSS sites, few studies have used VA data for investigating malaria cause specific mortality [20, 61, 62], verbal autopsy has great potential for countries like Tanzania where more people die outside of health facility care where no records are available. Also, VA has been shown to provide the best results to obtain the specific causes of death in most of SSA [63] and widely used [22].

This study has some limitations that need to be considered in interpreting the findings. First, presence of at least one mosquito net was considered as a proxy for use of mosquito net in the household. Information about exact use of the mosquito net was not collected during the study because significantly more investment is required for capturing this information. Secondly, there is a potential risk of misclassification about the causes of death where the sensitivity and specificity of the VA technique is relatively low for assessing cause of death and this may lead to underestimate or overestimate for malaria death. Thirdly, there are other possible factors associated with malaria deaths that were not used as confounders in this study and were not available in HDSS database, such as the availability of anti-malarial treatment.


This study used Bayesian framework modelling to assess the effect of mosquito net ownership on malaria mortality in rural Tanzania. The proposed hierarchical model takes into account spatial and temporal dependencies effect and environmental factors as confounders. It provides a useful tool for health environmental investigations and other factors. Results demonstrate that the proposed modelling approach is robust and can be useful in understanding the effect of ownership of mosquito net on the malaria mortality. Additionally, the model can be applied to analyze spread of other infectious diseases and can be extended to account for other factors such socio-economic factors as well as other important factors such as access to health, information on drug availability and drug resistance.

Reduction in mortality due to malaria requires more attention given to climate change such as rainfall and other non- climatic factors that affect malaria deaths such as ownership of mosquito nets. Ownership and use of mosquito nets should be a continuous strategy in the study areas and Tanzania as a whole for better prevention of malaria transmission and death.


  1. 1.

    Alba S, Nathan R, Schulze A, Mshinda H, Lengeler C. Child mortality patterns in rural Tanzania: an observational study on the impact of malaria control interventions. Int J Epidemiol. 2013;43:204–15.

  2. 2.

    Atieli HE, Zhou G, Afrane Y, Lee M-C, Mwanzo I, Githeko AK, Yan G. Insecticide-treated net (ITN) ownership, usage, and malaria transmission in the highlands of western Kenya. Parasit Vectors. 2011;4:113.

  3. 3.

    Mugittu K, Ndejembi M, Malisa A, Lemnge M, Premji Z, Mwita A, et al. Therapeutic efficacy of sulfadoxine-pyrimethamine and prevalence of resistance markers in Tanzania prior to revision of malaria treatment policy: Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthase mutations in monitoring in vivo resistance. Am J Trop Med Hyg. 2004;71:696–702.

  4. 4.

    Eisele TP, Larsen DA, Walker N, Cibulskis RE, Yukich JO, Zikusooka CM, et al. Estimates of child deaths prevented from malaria prevention scale-up in Africa 2001–2010. Malar J. 2012;11:93.

  5. 5.

    Lindblade KA, Eisele TP, Gimnig JE, Alaii JA, Odhiambo F, ter Kuile FO, et al. Sustainability of reductions in malaria transmission and infant mortality in western Kenya with use of insecticide-treated bednets: 4–6 years of follow-up. JAMA. 2004;291:2571–80.

  6. 6.

    Villalta D, Guenni L, Rubio-Palis Y, Arbel´aez RR. Bayesian space-time modelling of malaria incidence in Sucre state-Venezuela. AstA Adv Stat Anal. 2013;97:151–71.

  7. 7.

    Lengeler C. Insecticide-treated nets for malaria control: real gains. Bull World Health Organ. 2004;82:84.

  8. 8.

    Phillips-Howard PA, Nahlen BL, Kolczak MS, Hightower AW, ter Kuile FO, Alaii JA, et al. Efficacy of permethrin-treated bed nets in the prevention of mortality in young children in an area of high perennial malaria transmission in western Kenya. Am J Trop Med Hyg. 2003;68:23–9.

  9. 9.

    Okoyo C, Mwandawiro C, Kihara J, Simiyu E, Gitonga CW, Noor AM, et al. Comparing insecticide-treated bed net use to Plasmodium falciparum infection among schoolchildren living near Lake Victoria, Kenya. Malar J. 2015;14:515.

  10. 10.

    Smithuis FM, Kyaw MK, Phe UO, van der Broek I, Katterman N, Rogers C, et al. The effect of insecticide-treated bed nets on the incidence and prevalence of malaria in children in an area of unstable seasonal transmission in western Myanmar. Malar J. 2013;12:363.

  11. 11.

    Lindblade KA, Mwandama D, Mzilahowa T, Steinhardt L, Gimnig J, Shah M, et al. A cohort study of the effectiveness of insecticide-treated bed nets to prevent malaria in an area of moderate pyrethroid resistance, Malawi. Malar J. 2015;14:31.

  12. 12.

    Gamble C, Ekwaru PJ, Garner P, Kuile FO. Insecticide-treated nets for the prevention of malaria in pregnancy: a systematic review of randomised controlled trials. PLoS Med. 2007;4:e107.

  13. 13.

    Rowe A, Richard WS. Predictions of the impact of malaria control efforts on all cause child mortality in Sub-Saharan Africa. Am J Trop Med Hyg. 2007;77(suppl 6):48–55.

  14. 14.

    Hanson K, Nathan R, Marchant T, Mponda H, Jones C, Bruce J, et al. Vouchers for scaling up insecticide-treated nets in Tanzania: methods for monitoring and evaluation of a national health system intervention. BMC Public Health. 2008;8:205.

  15. 15.

    Tanzania Commission for AIDS (TACAIDS), Zanzibar AIDS Commission (ZAC), National Bureau of Statistics (NBS), Office of the Chief Government Statistician (OCGS) and Macro. International Inc. Tanzania HIV/AIDS and Malaria Indicator Survey 2011–2012. Dar es Salaam; 2013.

  16. 16.

    De Savigny D, Mayombana C, Mwagen E, Masanja H, Minhaj A, Mkilindi Y, et al. Care-seeking patterns for fatal malaria in Tanzania. Malar J. 2004;3:27.

  17. 17.

    Adult Morbidity and Mortality Project. Adult Morbidity and Mortality Project Report. Tanzania: End of Phase I. Dar es Saaam; 1997.

  18. 18.

    Gemperli A. Development of spatial statistical methods for modelling point-referenced spatial data in malaria epidemiology. Basel: Basel University; 2003.

  19. 19.

    Kazembe LN, Chirwa TF, Simbeye JS, Namangale JJ. Applications of Bayesian approach in modelling risk of malaria-related hospital mortality. BMC Med Res Methodol. 2008;8:6.

  20. 20.

    Deressa W, Fantahun M, Ali A. Malaria-related mortality based on verbal autopsy in an area of low endemicity in a predominantly rural population in Ethiopia. Malar J. 2007;6:128.

  21. 21.

    Chandramohan D, Quigley M. Effect of misclassification of causes of death in verbal autopsy: can it be adjusted? Int J Epidemiol. 2001;30:509–14.

  22. 22.

    Mpimbaza A, Filler S, Katureebe A, Kinara SO, Nzabandora E, Quick L, et al. Validity of verbal autopsy procedures for determining malaria deaths in different epidemiological settings in uganda. PLoS One. 2011;6:e26892.

  23. 23.

    Setel PW, Whiting DR, Hemed Y, Chandramohan D, Wolfson LJ, Alberti KG, et al. Validity of verbal autopsy procedures for determining cause of death in Tanzania. Trop Med Int Health. 2006;11:681–96.

  24. 24.

    Mirza NM, Macharia WM, Wafula EM, Agwanda RO, et al. Verbal autopsy: a tool for determining cause of death in a community. East Africa Med J. 1990;67:693–8.

  25. 25.

    Soleman N, Chandramohan D, Shibuya K. Verbal autopsy: current practices and challenges. Bull World Health Organ. 2006;84(3):239–45.

  26. 26.

    INDEPTH Network. Population and health in developing countries population, health and survival at INDEPTH sites. vol 1. Ottawa: International Development Research Centre; 2002.

  27. 27.

    Geubbels E, Amri S, Levira F, Schellenberg J, Masanja H, Nathan R. Health and Demographic Surveillance System profile. The Ifakara rural and urban Health and Demographic Surveillance System (Ifakara HDSS). Int J Epidemiol. 2015;44:848–61.

  28. 28.

    Mrema S, Kante AM, Levira F, Mono A, Irema K, de Savigny D. Health and Demographic Surveillance System profile. The Rufiji health and demographic surveillance system (Rufiji HDSS). Int J Epidemiol. 2015;44:472–83.

  29. 29.

    Masanja IM, Selemani M, Amuri B, Kajungu D, Khatib R, Kachur S, et al. Increased use of malaria rapid diagnostic tests improves targeting of anti-malarial treatment in rural Tanzania: implications for nationwide rollout of malaria rapid diagnostic tests. Malar J. 2012;11:221.

  30. 30.

    Kanté AM, Nathan R, Helleringer S, Sigilbert M, Levira F, Masanja H, et al. The contribution of reduction in malaria as a cause of rapid decline of under-five mortality: evidence from the Rufiji Health and Demographic Surveillance System (HDSS) in rural Tanzania. Malar J. 2014;13:180.

  31. 31.

    Kama MK. Clustering of mortality among children under five Years due to malaria at the Ifakara demographic surveillance site in Tanzania. Johannesburg: University of the Witwatersrand; 2008.

  32. 32.

    Hanson K, Marchant T, Nathan R, Mponda H, Jones C, Bruce J, et al. Household ownership and use of insecticide-treated nets among target groups after implementation of a national voucher programme in the United Republic of Tanzania: plausibility study using three annual cross sectional household surveys. BMJ. 2009;339:b2434.

  33. 33.

    INDEPTH Network. INDEPTH standardized verbal autopsy questionnaire. Accessed 28 Sep 2015.

  34. 34.

    Mwageni E, Momburi D, Juma Z, Irema M, Masanja H, the TEHIP, et al. Rufiji DSS. In: INDEPTH monograph series: population and health in developing countries. Ottawa; 2002.

  35. 35.

    WHO. International statistical classification of diseases and related health problems (ICD-10). Geneva: WHO; 1999.

  36. 36.

    Tanzania meteorological authority website. Accessed 30 Sep 2015.

  37. 37.

    Lutambi AM, Alexander M, Charles J, Mahutanga C, Nathan R. Under-five mortality: spatial-temporal clusters in Ifakara HDSS in South-eastern Tanzania. Glob Health Action. 2010;3:32–41.

  38. 38.

    Shabani J, Lutambi AM, Mwakalinga V, Masanja H. Rufiji Health and Demographic Surveillance System in rural Tanzania. Glob Health Action. 2010;3(10):3402.

  39. 39.

    Dhiman RC, Pahwa S, Dash A. Climate change and malaria in India interplay between temperature and mosquitoes. Regional Health Forum. 2008;12:27–31.

  40. 40.

    Mabaso MLH, 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.

  41. 41.

    Arab A, Jackson MC, Kongoli C. Modelling the effects of weather and climate on malaria distributions in West Africa. Malar J. 2014;13:126.

  42. 42.

    Berliner L. Hierarchical bayesian time series models. In: maximum entropy and bayesian methods. Netherlands: Springer; 1996.

  43. 43.

    Blangiardo M, Cameletti M. Spatial and spatio-temporal bayesian models with R-INLA. New Jersey: Wiley; 2015.

  44. 44.

    Knorr-Held L, Rasser G. Bayesian detection of clusters and discontinuities in disease maps. Biometrics. 2000;56:13–21.

  45. 45.

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

  46. 46.

    Pearson K. The problem of random walk. Nature. 1905;72:318.

  47. 47.

    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:319–92.

  48. 48.

    Zhu L, Carlin BP. Comparing hierarchical models for spatio-temporally misaligned data using deviance information creterion. Stat Med. 2000;19:2265–78.

  49. 49.

    Rumisha SF. Modelling the seasonal and spatial variation of malaria transmission in relation to mortality in Africa. Basle: Basle University; 2013.

  50. 50.

    Schellenberg JA, Abdulla S, Nathan R, Mukasa O, Marchant TJ, Kikumbih N, et al. Effect of large-scale social marketing of Insecticide treated nets on child survival in rural Tanzania. Lancet. 2001;357:1241–7.

  51. 51.

    Magesa SM, Lengeler C, de Savigny D, Miller JE, Njau RJA, Kramer K, et al. Creating an “enabling environment” for taking insecticide treated nets to national scale: the Tanzanian experience. Malar J. 2005;4:34.

  52. 52.

    WHO. World malaria report. Geneva: World Health Organization; 2008 (2009).

  53. 53.

    Connor S, Thomson M, Molyneux D. Forecasting and prevention of epidemic malaria: new perspectives on an old problem. Parassitologia. 1999;41:439–48.

  54. 54.

    Mrema S, Shamte A, Selemani M, Masanja H. The influence of weather on mortality in rural Tanzania: a time-series analysis 1999–2010. Glob Health Action. 2012;5:33–43.

  55. 55.

    WHO. World malaria report. Geneva: World Health Organization; 2013.

  56. 56.

    Reiter P. Climate change and mosquito borne diseases. Environ Health Perspect. 2001;109:223–34.

  57. 57.

    Wang L, Kanji S, Bandyopadhyay S. The health impact of extreme weather events in sub-Saharan Africa. New York: World Bank; 2009.

  58. 58.

    Branscum AJ, Perez AM, Johnson WO, Thurmond MC. Bayesian spatiotemporal analysis of foot-and-mouth disease data from the Republic of Turkey. Epidemiol Infect. 2008;136:833–42.

  59. 59.

    Abellan JJ, Richardson S, Best N. Use of space time models to investigate the stability of patterns of disease. Environ Health Perspect. 2008;116:1111–9.

  60. 60.

    Popoff E. An approximate spatio-temporal bayesian model for Alberta wheat yield. Vancouver: The University of British Columbia; 2014.

  61. 61.

    Sarkar J, Murhekar MV, Shah NK, van Hutin Y. Risk factors for malaria deaths in Jalpaiguri district, West Bengal, India: evidence for further action. Malar J. 2009;8:133.

  62. 62.

    Kaatano GM, Mashauri FM, Kinung’hi SM, Mwanga JR, Malima RC, Kishamawe C, et al. Patterns of malaria related mortality based on verbal autopsy in Muleba District, north-western Tanzania. Tanzan J Heal Res. 2009;11:210–8.

  63. 63.

    Kaufman JS, Asuzu MC, Rotimi CN, Johnson OO, Owoaje EE, Cooper RS. The absence of adult mortality data for sub-Saharan Africa: a practical solution. Bull World Health Organ. 1997;75:389–95.

Download references

Authors’ contributions

MS conceived the study idea, designed the study, carried out the statistical analyses, interpretation of the results, and drafted the manuscript. AML, SM, AS, KY, and ASM participated in the design of the study and reviewed the manuscript. AML, MJM and MCYM provided technical support for interpreting results and reviewed it for intellectual content. All authors read and approved the final manuscript.


We are thankful to the members of staff at the Rufiji and Ifakara HDSS and to all field enumerators who collect data for the HDSS. We are grateful to the Tanzania Commission of Science and Technology (COSTECH) and Grand Challenges Canada HASA Project through MEDA for funding this study as well as INDEPTH Network for its continued support to the Rufiji and Ifakara HDSS.

Competing interests

The authors declare that they have no competing interests.

Author information

Correspondence to Majige Selemani.

Additional files


Additional file 1. Appendix.


Additional file 2. Observed and predicted malaria mortality rate using model 2 for 2011in Rufiji HDSS. A1 Observed malaria mortality rate for under-five, A2 predicted malaria mortality rate for under-five, B1 Observed malaria mortality rate for all age and B2 Predicted malaria mortality for all age.


Additional file 3. Observed and predicted malaria mortality rate using model 2 for 2012in Ifakara HDSS. A1 Observed malaria mortality rate for under-five, A2 predicted malaria mortality rate for under-five, B1 Observed malaria mortality rate for all age and B2 Predicted malaria mortality for all age.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Selemani, M., Msengwa, A.S., Mrema, S. et al. Assessing the effects of mosquito nets on malaria mortality using a space time model: a case study of Rufiji and Ifakara Health and Demographic Surveillance System sites in rural Tanzania. Malar J 15, 257 (2016) doi:10.1186/s12936-016-1311-9

Download citation


  • Space time model
  • Malaria mortality


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.