 Research
 Open Access
 Published:
Relative importance of climatic, geographic and socioeconomic determinants of malaria in Malawi
Malaria Journal volume 12, Article number: 416 (2013)
Abstract
Background
Malaria transmission is influenced by variations in meteorological conditions, which impact the biology of the parasite and its vector, but also socioeconomic 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 agestratified spatiotemporal dataset of malaria cases from July 2004 to June 2011. Several climatic, geographic and socioeconomic 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 socioeconomic 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 climatedriven 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 malariacausing 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 insecticidetreated nets (ITN)[4]. The latter primarily focuses on the highrisk 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 longlasting 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 overdiagnosis 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 costeffectively 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, interannual 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 mosquitoborne 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 3339°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[16–23]. 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 nonclimate variables that can affect population vulnerability and alter the underlying spatial distribution of infectious disease hazard associated with climate and environment[24].
The socioeconomic 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]. KellyHope 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 socioeconomic 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 pointreferenced 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 subdistrict level in northern Malawi.
This study extends that of Kazembe[38] by using a spatiotemporal 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 overestimate the risk, provides a general picture of the overall malaria burden and resources required to decrease malaria risk. As well as environmental data, socioeconomic determinants are also considered to better account for spatial variations in malaria risk. The aim of this analysis is to investigate the spatial and interannual variations in malaria morbidity in Malawi and to determine how much, if any, of the interannual variability is due to climate variability relative to other nonclimatic 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[38–44]. Here, the spatiotemporal 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 socioeconomic, geographic and climate effects explain the spatial and interannual 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 km^{2}. 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, lowlying 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 2537°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 1727°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 nonclinically 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 nonclinically 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].
Socioeconomic 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).
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 ERAInterim reanalysis produced by the European Centre for MediumRange 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 locationspecific 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 computingintensive. 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 spatiotemporal variations in malaria risk in Malawi can be accounted for by climate variations is investigated, while accounting for both observed and unobserved nonclimatic 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 agestratified counts of malaria cases per month from July 2004  June 2011. The initial model included the climatic, demographic, socioeconomic 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 nonsignificant interaction terms, quadratic terms and explanatory variables.
Although the GLM accounted for extra variation by the inclusion of climate and nonclimate 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[55–57]. 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 pseudoobservations 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 crossvalidation, 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
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 nonclimate confounding factors, both observed, i.e. cartographic, demographic and socioeconomic 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 loglinear predictor$\mathit{\text{log}}{\rho}_{\mathit{\text{jst}}}=\alpha +{\sum}_{i}{\beta}_{i}{x}_{\mathit{\text{ijst}}}+{\sum}_{i}{\gamma}_{i}{z}_{\mathit{\text{ijst}}}$, where α is the model intercept, β is the parameter associated with climate covariates x_{ ijst } and γ with 'nonclimate’ 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.
Areaspecific 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${\varphi}_{s}\sim N(0,{\sigma}_{\varphi}^{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,${\upsilon}_{s}\sim \text{CAR}({\sigma}_{\upsilon}^{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,${\omega}_{1(t)}=0,\phantom{\rule{2.22144pt}{0ex}}{\omega}_{{t}^{\prime}(t)}\sim \mathrm{N}({\omega}_{{t}^{\prime}(t)1},{\sigma}_{\omega}^{2}),\phantom{\rule{2.22144pt}{0ex}}{t}^{\prime}(t)=2,\dots ,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${\tau}_{1(t)}=0,\phantom{\rule{2.22144pt}{0ex}}{\tau}_{{t}^{\prime}(t)}\sim \mathrm{N}(0,{\sigma}_{\tau}^{2}),\phantom{\rule{2.22144pt}{0ex}}{t}^{\prime}(t)=2,\dots ,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
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 burnin 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 GelmanRubin 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 interannual 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). Figure2cf 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.
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 centralwest 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 13 were found to be statistically significant. Rather than selecting a particular lag, or including all three lags separately, which could result in overfitting, 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 nonlinear effects, along with various interaction terms. Other geographic and socioeconomic 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 overestimating 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.
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].
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 selfdiagnose malaria and selftreat 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 socioeconomic 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 nonsignificant in the mixed effects model. This demonstrates the importance of accounting for spatial heterogeneity and correlation, when analysing geographical data in order to avoid underestimation 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 autocorrelated 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.
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].
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 pseudoobservations 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 interseasonal variability. In order to assess how much additional interseasonal and interannual 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.
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 lakeside communities, prevalence is generally much lower, due to the altitude and lower temperatures[65]. In these circumstances, interannual 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 lowlying southern districts, where climate is conducive to intense transmission, year to year climate perturbations are perhaps not expected to impact morbidity. Instead, changes to socioeconomic 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 socioeconomic 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 districtlevel. 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 subdistrict 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 longterm 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 pointreferenced 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 processbased modelling approaches, spatiotemporal predictions of malaria risk may be possible at a finer spatial scale in areas where data is not recorded. However, there is an ongoing 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 socioeconomic 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 interannual 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 lakeside 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 timespan of health data, this is the first spatiotemporal model for malaria relative risk in Malawi, at the district level. The analysis indicates that a climatebased early warning system could have some value in Malawi’s northern epidemicprone districts and emphasizes the critical requirement for an effective climate monitoring system, and access to high quality climate forecasts. A climatebased 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 20112015: Towards Universal Access: National Malaria Control Programme. 2011, Lilongwe
 2.
Malawi Population Projections: National Statistical Office, Government of Malawi. 2010, Zomba
 3.
Malawi National Malaria Indicator Survey 2010: Ministry of Health. 2010, Lilongwe
 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: 212217. 10.1016/j.actatropica.2011.06.017.
 5.
RocaFeltrer 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: 272278. 10.3201/eid1802.111008.
 6.
Guidelines for use of malaria rapid diagnostic tests (mRDTs) in Malawi: Ministry of Health. 2011, Lilongwe
 7.
KellyHope L, Thomson MC: Climate and infectious diseases. Seasonal Forecasts, Climatic Change and Human Health. 2008, New York: Springer, 3170.
 8.
McMichael AJ, CampbellLendrum 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
 9.
Hunter PR: Climate change and waterborne and vectorborne disease. J Appl Microbiol. 2003, 94: 3746. 10.1046/j.13652672.94.s1.5.x.
 10.
Githeko A, Lindsay S, Confalonieri U, Patz J: Climate change and vectorborne diseases: a regional analysis. Bull World Health Organ. 2000, 78: 11361147.
 11.
Craig MH, Snow RW, le Sueur D: A climatebased distribution model of malaria transmission in subSaharan Africa. Parasitol Today. 1999, 15: 105111. 10.1016/S01694758(99)013964.
 12.
Bayoh MN, Lindsay SW: Temperaturerelated duration of aquatic stages of the Afrotropical malaria vector mosquito Anopheles gambiae in the laboratory. Med Vet Entomol. 2004, 18: 174179. 10.1111/j.0269283X.2004.00495.x.
 13.
Gage K, Burkot T, Eisen R, Hayes E: Climate and vectorborne diseases. Am J Prev Med. 2008, 35: 436450. 10.1016/j.amepre.2008.08.030.
 14.
Reiter P: Climate change and mosquitoborne disease. Environ Health Perspect. 2001, 109: 141161.
 15.
Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007, 2: e114610.1371/journal.pone.0001146.
 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: 5463. 10.1023/A:1011943131643.
 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: 17051706. 10.1016/S01406736(03)133661.
 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: 3710.1186/14752875337.
 19.
Morse AP, DoblasReyes FJ, Hoshen MB, Hagendorn R, Palmer TIMN: A forecast quality assessment of an endtoend probabilistic multimodel seasonal forecast system using a malaria model. Tellus. 2005, 57A: 464475.
 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: 214221.
 21.
Thomson MC, Connor SJ: The development of malaria early warning systems for Africa. Trends Parasitol. 2001, 17: 438445. 10.1016/S14714922(01)020773.
 22.
Thomson MC, DoblasReyes FJ, Mason SJ, Hagedorn R, Connor SJ, Phindela T, Morse AP, Palmer TN: Malaria early warnings based on seasonal climate forecasts from multimodel ensembles. Nature. 2006, 439 (7076): 576579. 10.1038/nature04503.
 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: 6168.
 24.
Lafferty KD: The ecology of climate change and infectious diseases. Ecology. 2009, 90: 888900. 10.1890/080079.1.
 25.
Robert V, Macintyre K, Keating J, Trape JF, Duchemin JB, Warren M, Beier JC: Malaria transmission in urban subSaharan Africa. Am J Trop Med Hyg. 2003, 68: 169176.
 26.
Ayele D, Zewotir T, Mwambi H: Prevalence and risk factors of malaria in Ethiopia. Malar J. 2012, 11: 19510.1186/1475287511195.
 27.
Malawi Demographic and Health Survey 2010: National Statistical Office, Government of Malawi and ICF Macro. 2011, Zomba, Malawi & Calverton, Maryland, USA
 28.
KellyHope L, McKenzie FE: The multiplicity of malaria transmission a review of entomological inoculation rate measurements and methods across subSaharan Africa. Malar J. 2009, 8: 1910.1186/14752875819. doi:10.1186/14752875819
 29.
Mathanga DP, Bowie C: Malaria control in Malawi: are the poor being served?. Int J Equity Health. 2007, 6: 2210.1186/14759276622.
 30.
Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW: Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005, 3: 8190. 10.1038/nrmicro1069.
 31.
Tompkins AM, Ermert V: A regionalscale, high resolution dynamical malaria model that accounts for population density, climate and surface hydrology. Malar J. 2013, 12: 6510.1186/147528751265.
 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: 2935.
 33.
Dunn CE, Le Mare A, Makungu C: Malaria risk behaviours, sociocultural practices and rural livelihoods in southern Tanzania: implications for bednet usage. Soc Sci Med. 2011, 72: 408417. 10.1016/j.socscimed.2010.11.009.
 34.
AppiahDarkwah I, BaduNyarko SK: Knowledge of malaria prevention and control in a suburban community in Accra, Ghana. Int J Trop Med. 2011, 6: 6169.
 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: 109117.
 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: 9610.1186/14752875596.
 37.
Kazembe L, Kleinschmidt I, Holtz T, Sharp B: Spatial analysis and mapping of malaria risk in Malawi using pointreferenced prevalence of infection data. Int J Health Geogr. 2006, 5: 4110.1186/1476072X541.
 38.
Kazembe LN: Spatial modelling and risk factors of malaria incidence in Northern Malawi. Acta Trop. 2007, 102: 126137. 10.1016/j.actatropica.2007.04.012.
 39.
Mabaso M, Vounatsou P, Midzi S, Da Silva J, Smith T: Spatiotemporal analysis of the role of climate in interannual variation of malaria incidence in Zimbabwe. Int J Health Geogr. 2006, 5: 910.1186/1476072X59.
 40.
Craig M, Sharp B, Mabaso M, Kleinschmidt I: Developing a spatialstatistical model and map of historical malaria prevalence in Botswana using a staged variable selection procedure. Int J Health Geogr. 2007, 6: 4410.1186/1476072X644.
 41.
Lowe R, Bailey TC, Stephenson DB, Graham RJ, Coelho CAS, Carvalho MS, Barcellos C: Spatiotemporal modelling of climatesensitive disease risk: Towards an early warning system for dengue in Brazil. Comput Geosci. 2011, 37: 371381. 10.1016/j.cageo.2010.01.008.
 42.
Patil AP, Gething PW, Piel FB, Hay SI: Bayesian geostatistics in health cartography: the perspective of malaria. Trends Parasitol. 2011, 27: 24610.1016/j.pt.2011.01.003.
 43.
Lowe R, Bailey TC, Stephenson DB, Jupp TE, Graham RJ, Barcellos C, Carvalho MS: The development of an early warning system for climatesensitive disease risk with a focus on dengue epidemics in Southeast Brazil. Stat Med. 2013, 32: 864883. 10.1002/sim.5549.
 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 geostatistical models. Malar J. 2012, 11: 16010.1186/1475287511160.
 45.
Climate of Malawi.http://www.metmalawi.com/climate/climate.php,
 46.
Welfare Monitoring Survey 2009: National Statistical Office, Government of Malawi. 2010, Zomba
 47.
Malawi Population Projections 19982023: National Statistical Office, Government of Malawi. 2003, Zomba
 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 ERAInterim reanalysis: Configuration and performance of the data assimilation system. Q J R Meteorol Soc. 2011, 137: 553597. 10.1002/qj.828.
 50.
Farr TG, Kobrick M: Shuttle Radar Topography Mission produces a wealth of data. Eos Trans Amer Geophys Union. 2000, 81: 583585. 10.1029/EO081i048p00583.
 51.
Bivand R, Pebesma E, GómezRubio V: Applied spatial data analysis with R. 2008, New York: Springer
 52.
Diggle P, Tawn J, Moyeed R: Modelbased geostatistics. J R Stat Soc Ser C Appl Stat. 1998, 47: 299350.
 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 0387954570],
 54.
Breslow N, Clayton D: Approximate inference in generalized linear mixed models. J Am Stat Assoc. 1993, 88 (421): 925.
 55.
Gilks WR, Richardson S, Spiegelhalter DJ: Markov Chain Monte Carlo in Practice. 1996, Boca Raton, Florida, USA: Chapman & Hall/CRC
 56.
Brooks S: Markov chain Monte Carlo method and its application. J R Stat Soc Ser D Statistician. 1998, 47: 69100. 10.1111/14679884.00117.
 57.
Gelman A, Carlin J, Stern H, Rubin D: Bayesian Data Analysis, Second Edition. 2004, Boca Raton, Florida, USA: Chapman & Hall/CRC
 58.
Gelman A, Meng X, Stern H: Posterior predictive assessment of model fitness via realized discrepancies. Stat Sin. 1996, 6: 733759.
 59.
Lowe R: Spatiotemporal modelling of climatesensitive 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: 120. 10.1007/BF00116466.
 61.
Mollie A: Bayesian mapping of disease. Markov Chain Monte Carlo in Practice. 1996, Boca Raton, Florida, USA: Chapman & Hall/CRC, 359379.
 62.
Besag J, Green P, Higdon D, Mengersen K: Bayesian computation and stochastic systems. Stat Sci. 1995, 10: 341. 10.1214/ss/1177010123.
 63.
StewartIbarra A, Lowe R: Climate and nonclimate drivers of dengue epidemics in southern coastal Ecuador. Am J Trop Med Hyg. 2013, 88: 971981. 10.4269/ajtmh.120478.
 64.
Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences. Stat Sci. 1992, 7: 457472. 10.1214/ss/1177011136.
 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: 37810.1186/1475287510378.
 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: 583639. 10.1111/14679868.00353.
 67.
Larson PS, Mathanga DP, Campbell CH, Wilson ML: Distance to health services influences insecticidetreated net possession and use among six to 59 monthold children in Malawi. Malar J. 2012, 11: 1810.1186/147528751118.
 68.
Amexo M, Tolhurst R, Barnish G, Bates I: Malaria misdiagnosis: effects on the poor and vulnerable. Lancet. 2004, 364: 18961898. 10.1016/S01406736(04)174461.
 69.
Malawi Health Sector Strategic Plan 20112016: Moving towards equity and quality: Ministry of Health. 2011, Lilongwe
 70.
MacNab Y: Hierarchical Bayesian modeling of spatially correlated health service outcome and utilization rates. Biometrics. 2003, 59: 305316. 10.1111/15410420.00037.
 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: 3510.1186/147528751035.
 72.
Hoshen MB, Morse AP: A weatherdriven model of malaria transmission. Malar J. 2004, 3: 3210.1186/14752875332.
 73.
Dzinjalamala F: Epidemiology of malaria in Malawi. Epidemiol of Malawi. 2009, 203: 21
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
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
Below are the links to the 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.
About this article
Cite this article
Lowe, R., Chirombo, J. & Tompkins, A.M. Relative importance of climatic, geographic and socioeconomic determinants of malaria in Malawi. Malar J 12, 416 (2013). https://doi.org/10.1186/1475287512416
Received:
Accepted:
Published:
Keywords
 Malaria
 Climate
 Socioeconomic
 Statistical model
 Confounding factors
 Random effects
 Spatial correlation
Comments
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.