Skip to main content

Relative importance of climatic, geographic and socio-economic determinants of malaria in Malawi

Abstract

Background

Malaria transmission is influenced by variations in meteorological conditions, which impact the biology of the parasite and its vector, but also socio-economic conditions, such as levels of urbanization, poverty and education, which impact human vulnerability and vector habitat. The many potential drivers of malaria, both extrinsic, such as climate, and intrinsic, such as population immunity are often difficult to disentangle. This presents a challenge for the modelling of malaria risk in space and time.

Methods

A statistical mixed model framework is proposed to model malaria risk at the district level in Malawi, using an age-stratified spatio-temporal dataset of malaria cases from July 2004 to June 2011. Several climatic, geographic and socio-economic factors thought to influence malaria incidence were tested in an exploratory model. In order to account for the unobserved confounding factors that influence malaria, which are not accounted for using measured covariates, a generalized linear mixed model was adopted, which included structured and unstructured spatial and temporal random effects. A hierarchical Bayesian framework using Markov chain Monte Carlo simulation was used for model fitting and prediction.

Results

Using a stepwise model selection procedure, several explanatory variables were identified to have significant associations with malaria including climatic, cartographic and socio-economic data. Once intervention variations, unobserved confounding factors and spatial correlation were considered in a Bayesian framework, a final model emerged with statistically significant predictor variables limited to average precipitation (quadratic relation) and average temperature during the three months previous to the month of interest.

Conclusions

When modelling malaria risk in Malawi it is important to account for spatial and temporal heterogeneity and correlation between districts. Once observed and unobserved confounding factors are allowed for, precipitation and temperature in the months prior to the malaria season of interest are found to significantly determine spatial and temporal variations of malaria incidence. Climate information was found to improve the estimation of malaria relative risk in 41% of the districts in Malawi, particularly at higher altitudes where transmission is irregular. This highlights the potential value of climate-driven seasonal malaria forecasts.

Background

Malaria is one of the greatest public health problems in Malawi, placing children under the age of five and pregnant women at the highest risk of the disease. It is estimated that about 6 million clinical cases of malaria (for both under and over five years age groups) are reported annually[1], in a population of approximately 15 million[2]. Plasmodium falciparum is the leading malaria-causing parasite in Malawi and accounts for most of the cases reported. The 2010 Malaria indicator survey reported prevalence of 43.3%[3]. In a bid to reduce the spread of the disease, the Malawi government has embarked on different interventions such as vector control through indoor residual spraying and insecticide-treated nets (ITN)[4]. The latter primarily focuses on the high-risk groups of pregnant women and children under the age of five and is the main intervention method in Malawi. The government is now encouraging the use of long-lasting insecticidal nets of which there are mass distributions conducted periodically. Several Non Governmental Organizations (NGOs) and international organizations such as the World Health Organization (WHO), the President’s Malaria Initiative and the Global Fund are involved in malaria intervention activities. However, these NGOs usually work in selected districts and as such, some districts receive no aid. Despite the scaling up of malaria control measures over the last decade, the recent data have not suggested a decrease in the burden of disease[4, 5].

Another challenge that the government faces in the fight against malaria is the detection of the disease at the health facilities, especially those in rural areas without the necessary laboratory facilities for testing. Presumptive malaria diagnosis has been used in the past and this led to over-diagnosis of malaria cases and consequently poor guidance to map regional risk and manage the disease. In 2011 the government of Malawi adopted the WHO recommendation that tests should be carried out in all suspected malaria cases and rapid diagnostic tests (RDT) were subsequently progressively introduced to health facilities across the country[6].

The focus of the government through the National Malaria Control Programme is to achieve universal coverage in the prevention and treatment of malaria. The aim is that by the year 2015, the 2010 levels of malaria morbidity and mortality will be reduced by half[1]. In order to achieve this goal, there is a pressing need for tools that can be implemented rapidly and cost-effectively to mitigate the burden of malaria. The government needs to know which districts are most at risk and when. Defining the spatial distribution of a disease within a country or region allows public health decision makers to identify zones susceptible to epidemics and to target resources toward those areas at greatest risk[7].

Climate variability can affect malaria transmission, both in terms of spatial and seasonal distribution, inter-annual variability and epidemic potential[8]. Rainfall can affect the availability of mosquito breeding, developmental and resting sites[9], while temperature influences the rate of development of immature stages and adult survival rate, biting frequency and the extrinsic incubation periods (the period between infection of the vector and the vector’s ability to infect the next susceptible host) of disease agents[10]. Climate variability and change in the epidemiology of mosquito-borne diseases are complex. While increasing temperature speeds up vector, larvae and parasite biological cycles, high temperatures can increase mosquito and larva mortality[11, 12] and prevent malaria transmission above an upper temperature limit in the range of 33-39°C[13]. Likewise, rainfall can promote transmission by creating ground pools and other breeding sites, but heavy rains can have a flushing effect, cleansing such sites of early stage larvae[14, 15]. The relationships between malaria and climate have been well documented and statistical and dynamical models have been developed to represent them. With such models, climate observations and forecasts could be used to predict epidemics on monthly to seasonal time scales[1623]. In addition to the climatic factors that affect malaria transmission, consideration of the interaction of ecological variables with human behaviour and the urban environment is also important. Predictive models for malaria morbidity risk should also include potentially important non-climate variables that can affect population vulnerability and alter the underlying spatial distribution of infectious disease hazard associated with climate and environment[24].

The socio-economic status of a population can lead to conditions and environments conducive to vector proliferation and enhanced disease transmission. For example, rapid urbanization in Africa has contributed to unprotected water reservoirs, poor housing and lack of sanitation, which have implications for malaria transmission and epidemiology[25]. Poor sanitation in a household can lead to stagnant water holes that can act as breeding sites. Similarly, poor water storage and sanitation can provide breeding sites for mosquitoes around the household. A previous study of malaria risk in Ethiopia found that households with no toilet facilities were more likely to test positive for malaria[26]. That said, it is believed that those living in urban areas in Malawi are at a reduced risk of malaria compared to their rural counterparts. The 2010 Malawi Demographic and Health Survey reported that 30.7% of those surveyed in urban areas had fever in the preceding two weeks before the survey as compared to 35.1% in rural areas[27]. Kelly-Hope and McKenzie[28] report lower transmission intensities in urban African environments relative to rural areas. Urban areas possess qualities that reduce population vulnerability, compared to their rural counterparts. For example, literacy levels are generally higher and access to prevention interventions is greater than in rural settings. A study in Malawi found that the poorest populations in rural areas are not reached by intervention methods and ITN ownership was associated with living in urban areas and higher educational levels[29]. Further explanations posed as to why malaria transmission is lower in urban areas include pollution, which affects larval habitats and their life cycles[25, 30], mosquito avoidance behaviour by urban dwellers[25] and higher population densities resulting in lower biting rates[25, 31]. In rural locations, one room households are likely to be associated with greater risk of malaria compared to houses with more rooms, since such homes are unlikely to have sufficient nets for all members of extended families that frequently occupy such dwellings in rural Malawi. A study in Kenya found that adult mosquito abundance was significantly associated with traditional housing, with more mosquitoes found in grass thatched and mud wall houses[32].

Education status is a socio-economic factor that may impact malaria prevalence. Education status tends to affect the knowledge about malaria prevention and control among the population. Over recent years there has been emphasis on the idea that improving knowledge about malaria in communities will lead to better use of interventions[33, 34]. A study in India found educated respondents were more knowledgeable about malaria than the illiterate[35]. Literacy has been found to be positively associated with parasitaemia in Kenya[36]. Nevertheless, some of the observed association between disease prevalence and education may also be due to its role as a proxy for poverty.

Concerning interventions, malaria prevalence is expected to reduce as ITN distribution rates increase. The National Malaria Control Programme has been distributing ITNs freely with the purpose of reducing prevalence amongst the most vulnerable groups (children and pregnant women). Prevalence is expected to drop in areas with a high distribution rate, given that this results in greater usage. However, in Malawi, where the use of mosquito nets is the most widely used intervention, owing to its low cost, net usage is still not as high as anticipated, despite several mass distribution campaigns.

Although climate demonstratively defines the temporal limits of the transmission season and spatially demarks endemic, epidemic and malaria free zones, the role that climate variations play in governing year to year variability in malaria morbidity rates is not well understood. Previous modelling efforts in Malawi used spatial analysis to predict and map malaria risk across the country, by modelling point-referenced prevalence of infection data, with topographical and climatic factors as explanatory covariates[37]. Kazembe[38] then went on to profile the spatial variation of malaria risk in under fives, using Bayesian spatial analysis to investigate the possible association of disease risk with environmental factors at the sub-district level in northern Malawi.

This study extends that of Kazembe[38] by using a spatio-temporal dataset stratified by age (under five years and five years and over) for the whole country, spanning several years (July 2004 - June 2011). This improves the ability of the model to inform malaria risk trends in Malawi as a whole. Further, clinical malaria data is used, which, although has the propensity to over-estimate the risk, provides a general picture of the overall malaria burden and resources required to decrease malaria risk. As well as environmental data, socio-economic determinants are also considered to better account for spatial variations in malaria risk. The aim of this analysis is to investigate the spatial and inter-annual variations in malaria morbidity in Malawi and to determine how much, if any, of the inter-annual variability is due to climate variability relative to other non-climatic factors. This is accomplished by using a Bayesian statistical modelling approach, which is increasingly being used for mapping and predicting the risk of infectious diseases and can account for spatial dependence and unknown random effects[3844]. Here, the spatio-temporal model framework developed by Lowe[41, 43] is applied and extended to analyze a newly compiled database of malaria morbidity in Malawi and determine which socio-economic, geographic and climate effects explain the spatial and inter-annual variability.

Methods

Study area

Malawi is a small country in Southern Africa, bordered by Tanzania to the north, Zambia to the west and Mozambique to the south, with an area of about 120,000 km2. Malawi is divided into 28 districts within three administrative regions. The country has a varied topographical landscape with highland areas in the northern districts of Chitipa and Rumphi, where the Nyika Plateau is located. High mountainous areas also include Dedza, the Zomba Plateau and Mulanje Mountain in the south east corner of Malawi. These areas have lower temperatures, which reduce the transmission of the disease. In contrast, low-lying areas are found in the Shire river valley and along Lake Malawi. Districts in these areas experience higher temperatures and generally report higher malaria incidence. Lilongwe and Kasungu plains in the central region are the two biggest plains in the country.

Malawi falls within the tropics between latitudes 9°S and 18°S and longitudes 32°E and 36°E and experiences two distinct rainy and dry seasons. Malaria transmission mainly occurs between November and April during and shortly after the warm and wet months in which 95% of the annual precipitation is recorded. Average temperatures vary between 25-37°C in this period[45], which provides the best environmental conditions for the breeding of the malaria vectors. There is a considerable drop in malaria cases during the dry season between May and October, with mean temperatures varying between 17-27°C.

The majority of the population of Malawi live in rural areas and are involved in small scale subsistence farming. The rural population typically live in traditional dwellings with mud walls and a thatched roof. Overall, illiteracy is still a challenge in Malawi with an adult literacy rate in 2009 of 70%[46].

Data sources

Malaria data

Malaria case information was obtained from the Health Management Information Systems (HMIS) operated by the Ministry of Health (MOH), a database which records routine health data. In this system, both clinically and non-clinically diagnosed malaria cases disaggregated by age are routinely collected from all health facilities across Malawi, aggregated monthly to the district level. The HMIS in its current form was developed in 2002 but only started functioning nationally in 2004. During the first two years of being established, there were many challenges as the system did not function well in most districts. As RDT testing was introduced in 2011, the majority of diagnosed cases are non-clinically diagnosed suspected cases, which is likely to cause a significant overestimation of cases. From the district offices, data is sent to the national central database at the Ministry of Health every quarter. The malaria data used for this study was stratified by age (under five years and five years and over) for the period July 2004 - June 2011, as the reporting year starts in July and ends in June. Population and other demographic indicators were obtained from the National Statistical Office (NSO). The population figures for the districts were obtained from the population projections report by the NSO based on the 1998 population and housing census[47].

Socio-economic data

A database of potential drivers of malaria risk in Malawi was collated for every district for each month between July 2004 - June 2011 (see Table1). Certain variables were calculated using population estimates, for example population density, the proportion of the population living in urban areas and the proportion of health facilities and ITN distribution per inhabitant. Data on urbanization, housing, health care provision, sanitation and literacy levels were obtained from Welfare Monitoring Surveys (WMS)[46] conducted annually by the NSO (see Table1).

Table 1 Source and original resolution of datasets

Climatic and geographic data

Monthly precipitation estimates (units: mm day-1) were derived from the CPC/Famine Early Warning System Daily Rainfall Estimates (RFE 2.0) over Africa, available at a 10km resolution[48]. Temperature estimates (units: °C) were derived from the ERA-Interim reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ECMWF)[49], with a coarser resolution of 80km. Altitude data was obtained from the Shuttle Radar Topography Mission 90m digital elevation dataset[50]. Climatic and topographic data were calculated for each district in Malawi using an interpolation method to relate gridded products to administrative districts[51]. Districts were further grouped into administrative regions (north, central, south) and ecological zones (lakeshore, lowland, highland, highland/lakeshore, highland/lowland).

Statistical analysis

The objective of statistical modelling is to determine a minimal adequate model from the large set of potential models that might be used to describe the given set of data. Selecting few predictors from among a large number of potential candidates is a major challenge and can easily become arbitrary. An explanatory variable should only be included in the model if it significantly improves the fit of the model. A limitation of standard statistical modelling approaches is that they assume independence between survey locations and neglect potential spatial dependency between neighbouring locations due to unobserved common exposures[44]. Estimation of standard errors of explanatory variables is biased if overdispersion or spatial correlation is not taken into account within a model. Geostatistical models take into account spatial correlation by incorporating additional location-specific random effect parameters into a model. Such correlations may arise from factors such as variations in health system performance, intervention coverage or population immunity. Bayesian geostatistical approaches are increasingly used for mapping and predicting the risk of infectious diseases[52]. Ideally every possible combination of variables would be tested and compared in a Bayesian framework. However, this is not presently a viable approach as it is extremely time and computing-intensive. The most practical approach is to reduce the list of potential explanatory variables using general regression selection methods, before moving to a Bayesian context[40].

In this paper, the extent to which spatio-temporal variations in malaria risk in Malawi can be accounted for by climate variations is investigated, while accounting for both observed and unobserved non-climatic confounding factors, spatial heterogeneity and correlation. First, a maximal 'fixed effects’ model was fitted within a negative binomial generalized linear model (GLM) framework[53], to assess the relation between potential predictors and age-stratified counts of malaria cases per month from July 2004 - June 2011. The initial model included the climatic, demographic, socio-economic and cartographic variables described above, with relevant lags and polynomial terms. Categorical variables, to account for the annual cycle, administrative regions and ecological zones, with associated interaction terms, were also tested. With the assistance of a stepwise model selection procedure based on the Akaike Information Criterion (AIC), the model was simplified by removing non-significant interaction terms, quadratic terms and explanatory variables.

Although the GLM accounted for extra variation by the inclusion of climate and non-climate variables and factors, such as the annual cycle and ecological zones, there was still a large proportion of the variance that was unexplained. Consequently, a generalized linear mixed model (GLMM)[54] was adopted. The GLMM is an extension of the GLM that allows for additional variation in the response arising from unobservable random effects. The inclusion of random effects introduces an extra source of variability (a latent effect) into the model to capture the impact of unknown/unobserved confounding factors, such as variations in health care provision, unequally distributed aid or variations in population immunity. Spatially unstructured random effects can assist in modelling overdispersion, previously allowed for solely via the single scale parameter in the negative binomial GLM, while spatially structured random effects allow for correlated heterogeneity between districts. Parameters in a GLMM can be estimated using a Bayesian framework, where parameter uncertainty is accounted for by assigning prior distributions to the parameters. Hierarchical models can be created by parameterizing prior distributions with unknown 'hyperparameters’ which have their own 'hyperprior’ distribution. Markov Chain Monte Carlo (MCMC methods) make estimation of parameters in Bayesian models a practical feasibility[5557]. This is because associated MCMC sampling yields samples from full posterior predictive distributions, which automatically incorporate all components of variance at the different levels in the model and therefore, provide a full assessment of prediction uncertainty.

When assessing complex Bayesian models, it can be useful to use posterior predictive distributions as reference distributions for comparison to observed data[58]. The posterior predictive distribution of the response is obtained by simulating new pseudo-observations using samples from the posterior distribution of the parameters in the model. The distribution of estimated values can then be compared to observed values. This approach is an alternative to cross-validation, where the model is fitted with part of the data and the remaining observations are compared to the posterior predictive distribution calculated from the sample used for fitting.

Model formulation

Let y jst be counts of malaria cases in each age group (j = 1,2, where age group 1 represents five years and over and age group 2 represents under five years), district (s = 1,…,27) and month (t = 1,…,84). After conducting preliminary tests to assess the presence of overdisperion in the count data[59], it was assumed that y jst arises from a negative binomial distribution

y jst | μ jst NegBin ( μ jst , κ ) log μ jst = log e jst + log ρ jst ρ jst = i = 1 p exp ( θ i x ijst ) ,

where μ jst is the corresponding distribution mean, which is equal to the expected number of cases e jst multiplied by the unknown relative malaria risk ρ jst for a given age group j, district s and time t. κ is the scale (or overdispersion) parameter and θ represents the parameters associated with fixed and random effects included in the model parametrization. Note that population effects are accounted for by including the expected number of cases e jst (i.e. the population within each district, multiplied by the overall malaria risk) as an offset (see[41] for more details). The model equation can then be rearranged such that the relative risk ρ jst is equivalent to the standardized morbidity ratio (SMR), where SMR = y jst /e jst , i.e. the ratio of observed to expected cases within a district at a given time. Then, the most suitable estimate of the relative risk ρ jst (or SMR) is sought via a linear combination of climate covariates (temperature and precipitation) and non-climate confounding factors, both observed, i.e. cartographic, demographic and socio-economic covariates, or unobserved (using random effects) that might explain variations in malaria risk.

Initially for the model selection stage 'fixed effects’, both continuous and categorical, were included in the log-linear predictorlog ρ jst =α+ i β i x ijst + i γ i z ijst , where α is the model intercept, β is the parameter associated with climate covariates x ijst and γ with 'non-climate’ covariates z ijst . Next, random effects were included (hence mixed effects model). As the model parameters were estimated in a Bayesian model framework, prior distributions were specified for all parameters.

Area-specific random effects that are divided into spatially unstructured ϕ s and structured υ s components are often termed the 'convolution prior’, ϕ s  + υ s [60, 61]. Spatial heterogeneity was introduced by assigning exchangeable location specific random effects using a Gaussian distribution with zero mean and large variance for the unstructured prior ϕ s N(0, σ ϕ 2 ). Spatial clustering and correlation were accounted for by assigning a conditional intrinsic Gaussian autoregressive model (CAR) to the spatially structured prior, which takes the neighbourhood structure of the districts into account, υ s CAR( σ υ 2 )[62]. Note that Likoma island, located in Lake Malawi, was excluded from the analysis to facilitate the creation of the neighbourhood structure (hence 27 districts were modelled rather than 28).

Autocorrelated random effects for each calendar month were included to account for the annual cycle of malaria. Since only part of the malaria annual cycle may be attributable to climatic conditions, the inclusion of this effect allowed the model to account for other potential seasonal confounding variables, such as seasonal population movements[63]. Thus, climate variables are retained in the GLMM only if they add additional information that improves the fit of the model. This seasonal term is included as a structured first order autoregressive month effect to account for temporal serial correlation in malaria transmission (e.g., malaria relative risk in one month may depend on the risk in the previous month). The month effect was assigned a random walk or first difference prior distribution, in which each effect is derived from the immediately preceding effect, ω 1 ( t ) =0, ω t ( t ) N( ω t ( t ) - 1 , σ ω 2 ), t (t)=2,,12[55]. To account for the apparent trend in the data and unobserved confounding factors, an exchangeable unstructured prior was assigned to the year effect with year 1 (July 2004 - June 2005) set to zero and subsequent years assigned a Gaussian distribution with zero mean and large variance τ 1 ( t ) =0, τ t ( t ) N(0, σ τ 2 ), t (t)=2,,7[61].

Diffuse gamma hyperpriors were assigned to the precisions (1/σ2) for the spatial and temporal random effects. The specification for the relative risk is then

log ρ jst = α + i β i x ijst + i γ i z ijst + ϕ s + υ s + ω t ( t ) + τ t ( t ) .

MCMC simulation was used to produce samples of model parameter values from their joint posterior distribution. Two parallel MCMC chains were generated, each of length 25,000 with a burn-in of 20,000 and thinning of 10 to obtain 1000 samples from the joint posterior distribution (see[43]). Convergence was assessed by inspecting plots of traces of simulations for individual parameters and monitoring the Gelman-Rubin diagnostic[64]. Finally, posterior predictive distributions were generated to compare model predictions to observations.

Results & discussion

Potential drivers

A precursory view of the temporal and spatial variation of the malaria SMR in Malawi is given in Figure1 and2, along with potential driver variables. Figure1a shows the temporal series of malaria SMR from July 2004 - June 2011 for the under five year and five year and over age categories. A strong annual cycle is apparent, with the peak in the early months of the year. Figure1b shows the corresponding monthly average precipitation and temperature. The known lag between the malaria transmission season and the rains is clearly apparent. The inter-annual variability in the peak SMR is superimposed on an upward trend over the period. While changes in climate and environmental conditions cannot be ruled out, it is far more likely that this trend is a result of the improved levels of reporting that resulted as districts moved to and became familiar with the electronic based reported system that was introduced in 2004. Figure2a and2b show the overall malaria SMR (for under fives and five years and over respectively) in each district over the whole time period (84 months). Figure2c-f shows the ecological zones, mean altitude, population density, proportion of households with one room for sleeping, the mean ITN distribution rate over the seven year period and the number of health facilities per 1000 inhabitants, respectively.

Figure 1
figure 1

Malaria SMR and average climate in Malawi for the period July 2004 - June 2011. (a) Malaria standardised morbidity ratios (SMR) for the under five (dashed curve) and five years and over (solid curve) age categories and (b) average precipitation (solid bars) and average temperature (dashed curve) in Malawi for the period July 2004 - June 2011.

Figure 2
figure 2

Spatial distribution of malaria SMR, geographic and socio-economic indicators across Malawi for the period July 2004 - June 2011. Map of (a) malaria SMR for under fives, (b) malaria SMR for five years and over, (c) ecological zones, (d) mean altitude, (e) population density, (f) proportion of households with only one room for sleeping, (g) mean ITN distribution rate and (h) the number of health facilities per 1000 inhabitants in each district over the period July 2004 - June 2011.

It is interesting to note that the spatial distribution of SMR for the under fives category broadly reflects the map of prevalence produced by the Malaria Atlas Project (MAP) Bayesian analysis of survey data[65], with higher SMR rates along the western shoreline of lake Malawi and central-west lowlands of the southern part of the country. The adult distribution of SMR reflects the same pattern to a certain extent, with the exception of the high SMR in Mzimba district which was due to an outbreak in April 2006. The greatest contrast to the MAP data is that significant cases are still reported in the northern districts of Chitipa, Rumphi and Karonga, for which the MAP analysis reports very low prevalence rates, although it should be recalled that the MAP survey data is relatively sparse in Chitipa and Karonga districts. In the southern most district of Malawi, Nsanje, a relatively high SMR is observed along with a high proportion of household with only one room for sleeping and also low altitude. Note that the health facility rate is also higher than in surrounding districts.

Fixed effects model

Using the negative binomial GLM framework specified above (i.e. with only fixed effects), exploratory analyses were conducted to find the best time lags between climate variables and malaria. At the 0.05 level of significance, precipitation and temperature covariates lag 1-3 were found to be statistically significant. Rather than selecting a particular lag, or including all three lags separately, which could result in over-fitting, these variables were combined into three month average precipitation and temperature variables, lagged two months previous to the malaria month of interest. Quadratic terms related to these climate covariates were also tested in order to capture possible non-linear effects, along with various interaction terms. Other geographic and socio-economic variables (listed in Table1) were tested.

As there may be other spatially varying determinants of disease transmission risk, such as soil type, land cover and land use, for which data was not available, a low order varying proxy location parameter was incorporated in the GLM analysis by including polynomial functions of longitude and latitude, thereby treating location as a continuous predictor. Several categorical variables were also tested including age group (under five years and five years and over) to account for different vulnerability for children under the age of five; calendar month, to model the annual cycle and avoid over-estimating relationships between malaria and climatic variables due to seasonality; year, to account for the trend in the data; region (north, central and south) to account for possible inequalities between the management of malaria between these administratively defined regions; and ecological zone (lakeshore, lowland, highland, highland/lakeshore, highland/lowland). This preliminary analysis was assisted by use of a model selection algorithm based on the AIC stepwise regression[53]. In order to compare the fixed effects model with subsequent models including random effects, model parameters were estimated within a Bayesian framework using MCMC and were considered to be statistically significant if their 95% credible interval did not contain zero. Note that all continuous variables were first standardized to zero mean and unit variance to aid MCMC convergence. Table2 shows the parameter estimates and 95% credible interval for the continuous explanatory variables that were retained in the model as statistically significant, and resulted in the lowest deviance information criterion (DIC)[66] (indicating goodness of fit). Certain variables, such as the proportion of households with no toilet facilities, the proportion of the adult population who never attended school and literacy rates were not found to be statistically significant and did not improve the fit of the model. Note that the categorical variables age, month, year, region and ecological zone were also found to be statistically significant.

Table 2 Parameter estimates for statistically significant continuous explanatory variables for selected fixed effects model (GLM)

Mixed effects model

As standard statistical regression models assume independent (or at least uncorrelated) observations, a mixed effects model (GLMM) was fitted in a Bayesian hierarchical modelling framework to account for spatial and temporal dependence. The continuous and categorical variables selected in the fixed effects model were included. However, a random walk was introduced in the annual cycle, to account for dependencies between one month to the next and exchangeable random year effects were assigned to account for the trend and unobserved confounding factors. Spatial heterogeneity and correlation were accounted for using the convolution prior, described in the Model formulation section. Interestingly, all the previously selected variables, except average temperature, rainfall (quadratic relation), age group, the ITN distribution rate and proportion of health facilities per inhabitant, ceased to be statistically significant (i.e. the credible interval contained zero).

Figure3 shows the kernel density estimates for the marginal posterior distributions for the statistically significant parameters associated with the variables rainfall and rainfall squared, temperature and temperature squared, the proportion of health facilities per inhabitant and the ITN distribution rate. As in the fixed effects model, a statistically significant quadratic relation between average rainfall during the proceeding three months and malaria risk was found (see Figures3a and3b). Although temperature was statistically significant, once confounding factors were accounted for, a quadratic relationship between temperature and malaria risk in Malawi was not found to be statistically significant (see Figure3d). This is likely due to the monthly average temperature range (15.8 - 28.9°C) in Malawi over the time period not exceeding values at which mosquito activity is suppressed[13].

Figure 3
figure 3

Kernel density estimates for significant explanatory variables. Kernel density estimates for the marginal posterior distributions for the parameters associated with (a) average precipitation, (b) precipitation squared, (c) average temperature (d) temperature squared, (e) health facilities per inhabitant and (f) ITN distribution rate.

As in the fixed effects model, the number of health facilities per inhabitant was positively associated to malaria relative risk, as was the ITN distribution rate (see Figures3e and3f). At first this appears contrary to expectations, since improved access to medical treatment and preventative measures is expected to reduce the parasite burden in the population and hence, the transmission intensity, resulting in an inverse relationship. The positive relationship observed most likely reflects more frequent reporting and greater distribution of ITNs where more health facilities are present, along with possible targeting of ITN distribution by donors in highly burdened districts. This result potentially highlights the gap between ITN possession and proper use[67]. According to Amexo et al.[68] 70% of people in Africa self-diagnose malaria and self-treat at home. It is likely that the proportion of the population reporting promptly increases with proximity to health centres. The construction of new health centres is determined by the distance that members of a community have to walk to an already existing facility. The Government of Malawi recommends that the population should live within an 8km radius of a health facility. New health centres are constructed in order to attain this target. In a study conducted by the Ministry of Health, some districts were found to be better served than others. For example, 51% of the population in sparsely populated Chitipa district lives more 8km from a health facility. In Blantyre, Chiradzulu, Mulanje and Zomba less than 5% of the population reside more than 8km from a health facility[69].

Figure4 shows the multiplicative effect of the two components of the convolution prior to the model. The key feature of the convolution prior is that it allows the assessment of relative contributions of unstructured heterogeneity and spatial clustering to the overall variation of the area effects[70]. From Figure4 it is evident that spatial heterogeneity is the dominant cause of overdispersion in Malawi. The spatially unstructured random effect ϕ s accounts for residual overdispersion in districts that is not attributable to spatial correlation between districts. Here, the spatial correlation component has a minimal yet significant contribution to the convolution prior. Although other geographic and socio-economic covariates such as altitude, longitude, latitude, ecological zone, region, and proportion of the population in each district residing in traditional housing were significant in the fixed effects model, they became non-significant in the mixed effects model. This demonstrates the importance of accounting for spatial heterogeneity and correlation, when analysing geographical data in order to avoid under-estimation of the credible intervals of model covariates. The structure of the random spatial component of the model provides a combined measure of the various potential risk factors that might contribute to the underlying spatial variation in malaria risk. The advantage is that only two hyperpriors are estimated for the precisions of spatial random effects, rather than numerous parameters for each different fixed effect. This results in a more parsimonious model, containing few strong predictors that are more easily interpretable. Figure5 shows the contribution of the auto-correlated annual cycle and random yearly effects to the malaria relative risk, stratified by age group (under five years and five years and over) over the period July 2004 to June 2011. These effects help account for the annual cycle in malaria, that could be attributed to climate and/or seasonal population movements, and the overall upward trend that could be the result of improved reporting over the years as the health facilities became accustomed to the newly established HMIS.

Figure 4
figure 4

Multiplicative contribution of spatially unstructured and structured random effects to malaria relative risk. Spatial distribution of multiplicative contribution of posterior mean spatially (a) unstructured ϕ ̂ s and (b) structured υ ̂ s random effects.

Figure 5
figure 5

Multiplicative contribution of temporally unstructured and structured random effects to malaria relative risk. Temporal distribution of auto-correlated random month effects ω ̂ t ( t ) (accounting for annual cycle) and random year effects τt(t) (to account for unexplained trend) for under five and five years and over age categories.

The role of climate in estimating malaria relative risk

Figure6 shows a surface of the multiplicative contribution of climate variables to malaria relative risk in Malawi. Given varying average precipitation and temperature values, the maximum relative risk is found at the maximum temperature of 28°C and a precipitation rate of 6.24 mm day-1. This result is supported by other studies, for example, a quadratic relationship between malaria incidence and rainfall was found in Botswana[20]. This likely relates to the wash out of first stage larvae from breeding sites by intense rainfall. These effects have been included in some dynamical models of malaria transmission either implicitly[71] or explicitly[31].

Figure 6
figure 6

Multiplicative contribution of climate variables to malaria relative risk. Surface of malaria relative risk given varying average precipitation and temperature values. Note that the maximum relative risk is found at the maximum temperature of 28°C and a precipitation threshold of 6.24 mm day-1.

To assess the predictive ability of the mixed effects model, posterior predictive distributions of malaria relative risk were obtained for each district and month. New pseudo-observations were simulated by drawing random values from a negative binomial distribution with mean and scale parameter estimated using samples from the posterior distribution of the parameters in the model. To summarize this information, the observed and posterior predictive mean malaria risk (SMR) estimates were aggregated across space. Figure7 shows scatter plots and time series of observed versus predicted malaria SMR for the 84 month period for age groups five years and over (upper panel) and under five years (lower panel). To evaluate the model across space, the root mean squared error (RMSE), a measure of the difference between model predicted and observed values, was calculated over the 84 month period for each district. For both age groups there is an overall positive agreement between predicted and observed space aggregated malaria risk (Figures7a and7e). Both time series plots for the two age groups averaged over the whole of Malawi show that the model is able to capture the inter-seasonal variability. In order to assess how much additional inter-seasonal and inter-annual variability is explained by the climate covariates, climate was removed from the model and a prediction without climate was superimposed on the plot (Figures7b and7f). In general, when averaging across the country, little or no improvement in malaria relative risk estimation is achieved by the addition of climate covariates in the model. Figure7 (c and g) show maps comparing observed to predicted malaria risk in each district. Relatively low values of RMSE are found, particularly in the southern districts of Malawi. In general, crude maps of SMRs are subject to considerable random error, particularly if the population count within a district is low. Therefore, visual attention is drawn to areas where rates are based on the least stable estimates, for example, the large but sparsely population Mzimba district (north region). The inclusion of random effects in the model framework makes it possible to 'borrow strength’ from neighbouring districts, resulting in spatial smoothing of risk surfaces. Therefore, the greater difference between observed and predicted malaria risk for Mzimba may in fact highlight the poor reliability of the data in this area.

Figure 7
figure 7

Observed versus predicted malaria SMR in space and time. Scatter plot (a, e) and time series (b, f) of space aggregated observed versus posterior predictive mean malaria SMR for the 84 month time period. Root mean squared error (RMSE) of observed and posterior mean malaria SMR for the 27 districts of Malawi for the period July 2004 - June 2011 (c, g) for the five years and over (upper panel) and under five years (lower panel) age groups. The lower the RMSE, the better the model fit. Difference between RMSE for the model including climate information and RMSE for a model fit without climate information (d, h). Districts with negative values of RMSE clim - RMSE noclim (white) suggest that climate information improves the model in these areas. Districts with positive values of RMSE clim - RMSE noclim (grey) suggest that climate information does not improve the model.

To assess whether the inclusion of climate information in the model could improve model estimation of malaria relative risk at the district level from year to year, the RMSE of the model excluding climate covariates, RMSE noclim , was subtracted from the RMSE of the model including climate covariates, RMSE clim . Areas where RMSE clim  - RMSE noclim  < 0, highlighted in grey in Figure7 (d and h), indicate that climate information improves the estimation of malaria relative risk, as the inclusion of these covariates results in a smaller difference between the model predicted values and the observations. Therefore, according to this model, climate information could help estimate malaria relative risk in under fives and five years and over several months ahead in the districts: Blantyre, Chikhwawa, Chitipa, Machinga, Mzimba, Nkhata Bay, Nkhotakota, Nsanje, Ntchisi, Rumphi and Salima. This represents 41% of the districts in Malawi.

It is interesting to note that many of the districts, where climate is found to improve the model, are located in the north of the country. In this region, away from the lake-side communities, prevalence is generally much lower, due to the altitude and lower temperatures[65]. In these circumstances, inter-annual variability in climate can intermittently lead to years with wetter or warmer conditions, which can result in more intense malaria transmission. In contrast, in the low-lying southern districts, where climate is conducive to intense transmission, year to year climate perturbations are perhaps not expected to impact morbidity. Instead, changes to socio-economic conditions and interventions may dominate. In addition, it should be recalled that the spatial resolution of the district scale data may average out climatic effects if districts include widely varying terrain.

Although the model is able to identify the relative importance of climatic, geographic and socio-economic determinants of malaria in Malawi, this study has several limitations. Firstly, malaria data was only available for the whole of Malawi at the relatively coarse spatial resolution of the district-level. One advantage of using aggregated data is to alleviate problems of misreporting due to variations in the diagnostic capabilities and reporting practices between individual health facilities. However, given this limitation and the coarse resolution of the climate data, the model formulated is unable to capture sub-district variations in malaria, which are likely influenced by localized meteorological and social conditions.

Secondly, the statistical analysis is limited by the short time period for which malaria data is available. Nevertheless, the information gained from the spatial component allows some inference to be drawn as the role of climate and other factors in the transmission of malaria in Malawi which can also be used for improvement and evaluation of dynamical transmission model[31, 72]. As time goes by and the established monitoring system provides a clearer long-term picture of malaria transmission, more detailed temporal and spatial information can be included in the model.

Thirdly, malaria data was obtained from case records at clinics, health facilities and hospitals. Therefore, malaria cases treated at home are missed. As malaria is endemic in much of the country, data incurs the omission of asymptomatic cases. Both home treated and asymptomatic cases result in underestimation. However, overestimation can occur when health facilities report suspected cases of malaria that are not clinically confirmed. In fact, it is estimated that only 37% of fever cases in children in Malawi are actually due to malaria[73]. Although there is no specific measure included in the model to account for misreporting over time, the health care provision indicator gives a broad indication of the accuracy of recording, while the spatial and temporal random effects account for potential unmeasured variability.

This work compliments the spatial risk maps of malaria that have been produced using finer scale point-referenced prevalence of infection data for Malawi[37]. An advantage of using point referenced data obtained from surveys such as the malaria indicator survey (MIS) is that malaria in children under five years of age is clinically diagnosed using RDTs, thus preventing over estimation of cases. The MIS for Malawi was first conducted in 2010 and will be repeated every two years, thus allowing the incorporation of temporal effects in such models. By combining geostatistical and process-based modelling approaches, spatio-temporal predictions of malaria risk may be possible at a finer spatial scale in areas where data is not recorded. However, there is an on-going need for continued collaboration between statisticians, mathematical modellers, meteorologists, public health decision makers and stakeholders in Malawi to construct models and interpret model results. Despite the limitations of this case study, due to the relatively coarse spatial resolution, short time series and data quality issues, with careful model selection this sophisticated modelling framework could serve as a useful tool to understand the relationship between climate, geographic and socio-economic conditions and malaria burden in other countries.

Conclusion

The main contribution of this paper is the collation of a unique dataset of potential spatial and temporal drivers of malaria in Malawi and the use of a sophisticated modelling procedure to determine the most important of these drivers. An initial model was selected that contained statistically significant fixed effects. After accounting for spatial heterogeneity and correlation, the mixed effects model was reduced to contain a few predictors that are easily interpretable, including average temperature and rainfall. Including climate information improves the estimation of inter-annual variations in malaria relative risk in 41% of the districts in Malawi, some of which are located in the north highland regions that are subject to lower and intermittent malaria transmission intensity (with the exception of lake-side communities). In the southern region, where malaria transmission is more intense, climate improved the model’s capability to represent year to year variations in malaria relative risk in only a few districts.

While this analysis has the common caveats associated with reliability and limited time-span of health data, this is the first spatio-temporal model for malaria relative risk in Malawi, at the district level. The analysis indicates that a climate-based early warning system could have some value in Malawi’s northern epidemic-prone districts and emphasizes the critical requirement for an effective climate monitoring system, and access to high quality climate forecasts. A climate-based malaria decision support system could be invaluable for the Malawi National Malaria Control Programme to be able to annually plan their locally targeted control interventions and manage scarce health resources.

Abbreviations

AIC:

Akaike Information Criterion

CAR:

Conditional autoregressive (model)

DIC:

Deviance information criterion

ECMWF:

European Centre for Medium Range Weather Forecasts

GLM:

Generalized linear model

GLMM:

Generalized linear mixed model

HMIS:

Health Management Information Systems

ITN:

Insecticide treated net

MAP:

Malaria Atlas Project

MCMC:

Markov Chain Monte Carlo

MIS:

Malaria indicator survey

MOH:

Ministry of Health

NGOs:

Non Governmental Organizations

NSO:

National Statistical Office

RDT:

Rapid diagnostic tests

RSME:

Root mean square error

SMR:

Standardized morbidity ratio

WHO:

World Health Organization

WMS:

Welfare Monitoring Survey.

References

  1. Malaria Strategic Plan 2011-2015: Towards Universal Access: National Malaria Control Programme. 2011, Lilongwe

    Google Scholar 

  2. Malawi Population Projections: National Statistical Office, Government of Malawi. 2010, Zomba

    Google Scholar 

  3. Malawi National Malaria Indicator Survey 2010: Ministry of Health. 2010, Lilongwe

    Google Scholar 

  4. Mathanga DP, Walker ED, Wilson ML, Ali D, Taylor TE, Laufer MK: Malaria control in Malawi: current status and directions for the future. Acta Trop. 2012, 121: 212-217. 10.1016/j.actatropica.2011.06.017.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Roca-Feltrer A, Kwizombe CJ, Sanjoaquin MA, Sesay SS, Faragher B, Harrison J, Geukers K, Kabuluzi S, Mathanga DP, Molyneux E, Chagomera M, Taylor T, Molyneux M, Heyderman RS: Lack of decline in childhood malaria, Malawi, 2001–2010. Emerg Infect Dis. 2012, 18: 272-278. 10.3201/eid1802.111008.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Guidelines for use of malaria rapid diagnostic tests (mRDTs) in Malawi: Ministry of Health. 2011, Lilongwe

    Google Scholar 

  7. Kelly-Hope L, Thomson MC: Climate and infectious diseases. Seasonal Forecasts, Climatic Change and Human Health. 2008, New York: Springer, 31-70.

    Chapter  Google Scholar 

  8. McMichael AJ, Campbell-Lendrum DH, Corvalán CF, Ebi KL, Githeko AK, Scheraga JD, Woodward A: Climate change and human health: risks and responses. 2003, Geneva: World Health Organ, 322pp

    Google Scholar 

  9. Hunter PR: Climate change and waterborne and vector-borne disease. J Appl Microbiol. 2003, 94: 37-46. 10.1046/j.1365-2672.94.s1.5.x.

    Article  Google Scholar 

  10. Githeko A, Lindsay S, Confalonieri U, Patz J: Climate change and vector-borne diseases: a regional analysis. Bull World Health Organ. 2000, 78: 1136-1147.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.

    Article  CAS  PubMed  Google Scholar 

  12. Bayoh MN, Lindsay SW: Temperature-related duration of aquatic stages of the Afrotropical malaria vector mosquito Anopheles gambiae in the laboratory. Med Vet Entomol. 2004, 18: 174-179. 10.1111/j.0269-283X.2004.00495.x.

    Article  CAS  PubMed  Google Scholar 

  13. Gage K, Burkot T, Eisen R, Hayes E: Climate and vectorborne diseases. Am J Prev Med. 2008, 35: 436-450. 10.1016/j.amepre.2008.08.030.

    Article  PubMed  Google Scholar 

  14. Reiter P: Climate change and mosquito-borne disease. Environ Health Perspect. 2001, 109: 141-161.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Githeko AK, Ndegwa W: Predicting malaria epidemics in the Kenyan highlands using climate data: a tool for decision makers. Glob Change & Hum Health. 2001, 2: 54-63. 10.1023/A:1011943131643.

    Article  Google Scholar 

  17. Hay SI, Were EC, Renshaw M, Noor AM, Ochola SA, Olusanmi I, Alipui N, Snow RW: Forecasting, warning, and detection of malaria epidemics: a case study. Lancet. 2003, 361: 1705-1706. 10.1016/S0140-6736(03)13366-1.

    Article  PubMed Central  PubMed  Google Scholar 

  18. DaSilva J, Garanganga B, Teveredzi V, Marx SM, Mason SJ, Connor SJ: Improving epidemic malaria planning, preparedness and response in Southern Africa. Malar J. 2004, 3: 37-10.1186/1475-2875-3-37.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Morse AP, Doblas-Reyes FJ, Hoshen MB, Hagendorn R, Palmer TIMN: A forecast quality assessment of an end-to-end probabilistic multi-model seasonal forecast system using a malaria model. Tellus. 2005, 57A: 464-475.

    Article  Google Scholar 

  20. Thomson MC, Mason SJ, Phindela T, Connor SJ: Use of rainfall and sea surface temperature monitoring for malaria early warning in Botswana. Am J Trop Med Hyg. 2005, 73: 214-221.

    PubMed  Google Scholar 

  21. Thomson MC, Connor SJ: The development of malaria early warning systems for Africa. Trends Parasitol. 2001, 17: 438-445. 10.1016/S1471-4922(01)02077-3.

    Article  CAS  PubMed  Google Scholar 

  22. Thomson MC, Doblas-Reyes FJ, Mason SJ, Hagedorn R, Connor SJ, Phindela T, Morse AP, Palmer TN: Malaria early warnings based on seasonal climate forecasts from multi-model ensembles. Nature. 2006, 439 (7076): 576-579. 10.1038/nature04503.

    Article  CAS  PubMed  Google Scholar 

  23. Ceccato P, Ghebremeskel T, Jaiteh M, Graves P, Levy M, Ghebreselassie S, Ogbamariam A, Barnston A, Bell M, del Corral J, Connor SJ, Fesseha I, Brantly EP, Thomson MC: Malaria stratification, climate, and epidemic early warning in Eritrea. Am J Trop Med Hyg. 2007, 77: 61-68.

    PubMed  Google Scholar 

  24. Lafferty KD: The ecology of climate change and infectious diseases. Ecology. 2009, 90: 888-900. 10.1890/08-0079.1.

    Article  PubMed  Google Scholar 

  25. Robert V, Macintyre K, Keating J, Trape JF, Duchemin JB, Warren M, Beier JC: Malaria transmission in urban sub-Saharan Africa. Am J Trop Med Hyg. 2003, 68: 169-176.

    PubMed  Google Scholar 

  26. Ayele D, Zewotir T, Mwambi H: Prevalence and risk factors of malaria in Ethiopia. Malar J. 2012, 11: 195-10.1186/1475-2875-11-195.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Malawi Demographic and Health Survey 2010: National Statistical Office, Government of Malawi and ICF Macro. 2011, Zomba, Malawi & Calverton, Maryland, USA

    Google Scholar 

  28. Kelly-Hope L, McKenzie FE: The multiplicity of malaria transmission a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009, 8: 19-10.1186/1475-2875-8-19. doi:10.1186/1475-2875-8-19

    Article  PubMed Central  PubMed  Google Scholar 

  29. Mathanga DP, Bowie C: Malaria control in Malawi: are the poor being served?. Int J Equity Health. 2007, 6: 22-10.1186/1475-9276-6-22.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW: Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005, 3: 81-90. 10.1038/nrmicro1069.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Tompkins AM, Ermert V: A regional-scale, high resolution dynamical malaria model that accounts for population density, climate and surface hydrology. Malar J. 2013, 12: 65-10.1186/1475-2875-12-65.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Zhou G, Munga S, Minakawa N, Githeko AK, Yan G: Spatial relationship between adult malaria vector abundance and environmental factors in western Kenya highlands. Am J Trop Med Hyg. 2007, 77: 29-35.

    PubMed  Google Scholar 

  33. Dunn CE, Le Mare A, Makungu C: Malaria risk behaviours, socio-cultural practices and rural livelihoods in southern Tanzania: implications for bednet usage. Soc Sci Med. 2011, 72: 408-417. 10.1016/j.socscimed.2010.11.009.

    Article  PubMed  Google Scholar 

  34. Appiah-Darkwah I, Badu-Nyarko SK: Knowledge of malaria prevention and control in a sub-urban community in Accra, Ghana. Int J Trop Med. 2011, 6: 61-69.

    Article  Google Scholar 

  35. Sharma AK, Aggarwal OP, Chaturvedi S, Bhasin SK: Is education a determinant of knowledge about malaria among Indian tribal population?. J Commun Dis. 2003, 35: 109-117.

    CAS  PubMed  Google Scholar 

  36. Coldren RL, Prosser T, Ogolla F, Ofula VO, Adungo N: Literacy and recent history of diarrhoea are predictive of Plasmodium falciparum parasitaemia in Kenyan adults. Malar J. 2006, 5: 96-10.1186/1475-2875-5-96.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Kazembe L, Kleinschmidt I, Holtz T, Sharp B: Spatial analysis and mapping of malaria risk in Malawi using point-referenced prevalence of infection data. Int J Health Geogr. 2006, 5: 41-10.1186/1476-072X-5-41.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Kazembe LN: Spatial modelling and risk factors of malaria incidence in Northern Malawi. Acta Trop. 2007, 102: 126-137. 10.1016/j.actatropica.2007.04.012.

    Article  PubMed  Google Scholar 

  39. Mabaso M, 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: 9-10.1186/1476-072X-5-9.

    Article  Google Scholar 

  40. Craig M, Sharp B, Mabaso M, Kleinschmidt I: Developing a spatial-statistical model and map of historical malaria prevalence in Botswana using a staged variable selection procedure. Int J Health Geogr. 2007, 6: 44-10.1186/1476-072X-6-44.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Lowe R, Bailey TC, Stephenson DB, Graham RJ, Coelho CAS, Carvalho MS, Barcellos C: Spatio-temporal modelling of climate-sensitive disease risk: Towards an early warning system for dengue in Brazil. Comput Geosci. 2011, 37: 371-381. 10.1016/j.cageo.2010.01.008.

    Article  Google Scholar 

  42. Patil AP, Gething PW, Piel FB, Hay SI: Bayesian geostatistics in health cartography: the perspective of malaria. Trends Parasitol. 2011, 27: 246-10.1016/j.pt.2011.01.003.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Lowe R, Bailey TC, Stephenson DB, Jupp TE, Graham RJ, Barcellos C, Carvalho MS: The development of an early warning system for climate-sensitive disease risk with a focus on dengue epidemics in Southeast Brazil. Stat Med. 2013, 32: 864-883. 10.1002/sim.5549.

    Article  PubMed  Google Scholar 

  44. Raso G, Schur N, Utzinger J, Koudou B, Tchicaya E, Rohner F, N’Goran E, Silué K, Matthys B, Assi S: Mapping malaria risk among children in Cote d’Ivoire using Bayesian geo-statistical models. Malar J. 2012, 11: 160-10.1186/1475-2875-11-160.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Climate of Malawi.http://www.metmalawi.com/climate/climate.php,

  46. Welfare Monitoring Survey 2009: National Statistical Office, Government of Malawi. 2010, Zomba

    Google Scholar 

  47. Malawi Population Projections 1998-2023: National Statistical Office, Government of Malawi. 2003, Zomba

    Google Scholar 

  48. Love T: The climate prediction center rainfall algorithm version 2. Tech. rep., NOAA Climate Prediction Center, available athttp://www.cpc.noaa.gov/products/fews/ 2002.,

  49. Dee D, Uppala S, Simmons A, Berrisford P, Poli P, Kobayashi S, Andrae U, Balmaseda M, Balsamo G, Bauer P: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Q J R Meteorol Soc. 2011, 137: 553-597. 10.1002/qj.828.

    Article  Google Scholar 

  50. Farr TG, Kobrick M: Shuttle Radar Topography Mission produces a wealth of data. Eos Trans Amer Geophys Union. 2000, 81: 583-585. 10.1029/EO081i048p00583.

    Article  Google Scholar 

  51. Bivand R, Pebesma E, Gómez-Rubio V: Applied spatial data analysis with R. 2008, New York: Springer

    Google Scholar 

  52. Diggle P, Tawn J, Moyeed R: Model-based geostatistics. J R Stat Soc Ser C Appl Stat. 1998, 47: 299-350.

    Article  Google Scholar 

  53. Venables WN, Ripley BD: Modern Applied Statistics with S. 2002, New York: Springer, fourth edition, [http://www.stats.ox.ac.uk/pub/MASS4] [ISBN 0-387-95457-0],

    Book  Google Scholar 

  54. Breslow N, Clayton D: Approximate inference in generalized linear mixed models. J Am Stat Assoc. 1993, 88 (421): 9-25.

    Google Scholar 

  55. Gilks WR, Richardson S, Spiegelhalter DJ: Markov Chain Monte Carlo in Practice. 1996, Boca Raton, Florida, USA: Chapman & Hall/CRC

    Google Scholar 

  56. Brooks S: Markov chain Monte Carlo method and its application. J R Stat Soc Ser D Statistician. 1998, 47: 69-100. 10.1111/1467-9884.00117.

    Article  Google Scholar 

  57. Gelman A, Carlin J, Stern H, Rubin D: Bayesian Data Analysis, Second Edition. 2004, Boca Raton, Florida, USA: Chapman & Hall/CRC

    Google Scholar 

  58. Gelman A, Meng X, Stern H: Posterior predictive assessment of model fitness via realized discrepancies. Stat Sin. 1996, 6: 733-759.

    Google Scholar 

  59. Lowe R: Spatio-temporal modelling of climate-sensitive disease risk: towards an early warning system for dengue in Brazil. PhD thesis,. University of Exeter 2010

  60. Besag J, York J, Mollié A: Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math. 1991, 43: 1-20. 10.1007/BF00116466.

    Article  Google Scholar 

  61. Mollie A: Bayesian mapping of disease. Markov Chain Monte Carlo in Practice. 1996, Boca Raton, Florida, USA: Chapman & Hall/CRC, 359-379.

    Google Scholar 

  62. Besag J, Green P, Higdon D, Mengersen K: Bayesian computation and stochastic systems. Stat Sci. 1995, 10: 3-41. 10.1214/ss/1177010123.

    Article  Google Scholar 

  63. Stewart-Ibarra A, Lowe R: Climate and non-climate drivers of dengue epidemics in southern coastal Ecuador. Am J Trop Med Hyg. 2013, 88: 971-981. 10.4269/ajtmh.12-0478.

    Article  PubMed Central  PubMed  Google Scholar 

  64. Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences. Stat Sci. 1992, 7: 457-472. 10.1214/ss/1177011136.

    Article  Google Scholar 

  65. Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, Tatem AJ, Hay SI: A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011, 10: 378-10.1186/1475-2875-10-378.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Spiegelhalter D, Best N, Carlin B, van der Linde A: Bayesian measures of model complexity and fit. J R Stat Soc Series B Stat Methodol. 2002, 64: 583-639. 10.1111/1467-9868.00353.

    Article  Google Scholar 

  67. Larson PS, Mathanga DP, Campbell CH, Wilson ML: Distance to health services influences insecticide-treated net possession and use among six to 59 month-old children in Malawi. Malar J. 2012, 11: 18-10.1186/1475-2875-11-18.

    Article  PubMed Central  PubMed  Google Scholar 

  68. Amexo M, Tolhurst R, Barnish G, Bates I: Malaria misdiagnosis: effects on the poor and vulnerable. Lancet. 2004, 364: 1896-1898. 10.1016/S0140-6736(04)17446-1.

    Article  PubMed  Google Scholar 

  69. Malawi Health Sector Strategic Plan 2011-2016: Moving towards equity and quality: Ministry of Health. 2011, Lilongwe

    Google Scholar 

  70. MacNab Y: Hierarchical Bayesian modeling of spatially correlated health service outcome and utilization rates. Biometrics. 2003, 59: 305-316. 10.1111/1541-0420.00037.

    Article  PubMed  Google Scholar 

  71. Ermert V, Fink AH, Jones AE, Morse AP: Development of a new version of the Liverpool Malaria Model. I. Refining the parameter settings and mathematical formulation of basic processes based on a literature review. Malar J. 2011, 10: 35-10.1186/1475-2875-10-35.

    Article  PubMed Central  PubMed  Google Scholar 

  72. Hoshen MB, Morse AP: A weather-driven model of malaria transmission. Malar J. 2004, 3: 32-10.1186/1475-2875-3-32.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Dzinjalamala F: Epidemiology of malaria in Malawi. Epidemiol of Malawi. 2009, 203: 21-

    Google Scholar 

Download references

Acknowledgements

This study was funded by the EU project QWeCI (Quantifying Weather and Climate Impacts on health in developing countries; funded by the European Commissions Seventh Framework Research Programme under the grant agreement 243964). The authors are grateful to Harry Gombachika and Chris Moyo for helping to establish the collaboration between QWeCI and the Malawi Ministry of Health.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Rachel Lowe.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RL collated the data, devised the model framework and wrote the manuscript. JC provided data and local expert advice, helped collate the data and contributed to the manuscript. AMT helped compile the climatic data, design the analysis and contributed to the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License(http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lowe, R., Chirombo, J. & Tompkins, A.M. Relative importance of climatic, geographic and socio-economic determinants of malaria in Malawi. Malar J 12, 416 (2013). https://0-doi-org.brum.beds.ac.uk/10.1186/1475-2875-12-416

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1475-2875-12-416

Keywords