Search tips
Search criteria 


Logo of wtpaEurope PMCEurope PMC Funders GroupSubmit a Manuscript
Trop Med Int Health. Author manuscript; available in PMC 2008 January 16.
Published in final edited form as:
PMCID: PMC2202922

Bayesian spatial analysis and disease mapping: tools to enhance planning and implementation of a schistosomiasis control programme in Tanzania


An increasing number of studies have attempted to predict disease risk geographically on the basis of satellite data, although few have been applied to large-scale control programmes. The aim of our study was to predict the spatial distributions of Schistosoma haematobium and Schistosoma mansoni infections to assist planning the implementation of mass distribution of praziquantel as part of an on-going national control programme in Tanzania. Bayesian geostatistical models were developed using parasitological data from 143 schools. In the S. haematobium models, whilst land surface temperature and rainfall were significant predictors of prevalence, they became non-significant when spatial correlation was taken into account. For S. mansoni, distance to water bodies and annual minimum temperature were significant predictors, even in models adjusting for spatial correlation. Spatial correlation occurred over greater distances for S. haematobium than for S. mansoni. Uncertainties in predictions were examined to identify areas requiring further data collection before programme implementation.

Keywords: Spatial distribution, maps, Bayesian analysis, Schistosoma haematobium, Schistosoma mansoni, schistosomiasis, communicable disease control, Tanzania


In implementing any control programme, understanding where at-risk populations live is fundamental for appropriate geographical targeting of resources and cost-effective control. Used appropriately, geographical information systems (GIS), remote-sensing (RS) and spatial analysis have an as-yet unrealised potential as tools for standardized programme surveillance and implementation (Kamel Boulos 2004; Brooker et al. 2005a). Numerous studies have been undertaken using satellite-derived environmental data to predict the distribution, abundance and prevalence of diseases and their vectors, including malaria (Hay et al. 2000; Rogers et al., 2002), leishmaniasis (Elnaiem et al. 2003), filariasis (Lindsay & Thomas 2000), trypanosomiasis (Rogers 2000) and schistosomiasis (Brooker et al. 2001; Malone et al. 2001; Brooker et al. 2002a, 2002b; Moodley et al. 2003; Kabatereine et al. 2004). Although a useful benchmark for future applications of disease mapping in Africa, using accessible and practical methods, these approaches could be improved upon by assessing uncertainty, which is inherent in all aspects of disease mapping, including the data, models, analyses and predictions (Elith et al. 2002) and by incorporating the spatial structure of the data. Maps that include estimates of uncertainty in model outputs can allow more informed and objective decision-making in relation to targeted disease control, as the control programme managers gain greater appreciation of decision risk.

Bayesian methods, which offer a flexible and robust approach, are increasingly being applied to spatial analysis, disease mapping and decision-making (Best et al. 2005). They provide convenient platforms for incorporating spatial correlation and by modelling both the observed data and any unknowns as random variables, they can incorporate uncertainty into the modelling process. Initially derived for use and now frequently employed in small-area analyses of chronic, non-infectious diseases (Best et al. 2005), Bayesian approaches have recently been used to study the geographical distribution of tropical diseases, both at a large scale, including malaria (Gemperli et al. 2004), and at a smaller scale, including filariasis (Alexander et al. 2000), malaria (Diggle et al. 2002), onchocerciasis (Carabin et al. 2003) and Schistosoma mansoni infection (Raso et al. 2005). While the use of these spatial analytical methods is an attractive research objective, they have been little applied in the context of large-scale disease control and there remains a need to document the applicability of GIS, RS and spatial analysis as tools for enhanced planning and implementation of large-scale disease control programmes.

The aim of the present study was to produce accurate, validated prevalence surface maps for both S. haematobium and S. mansoni infections in northwest Tanzania in order to spatially-define an implementation strategy based on mass treatment with praziquantel. An additional aim was to investigate the confidence of the prevalence predictions to inform decisions on whether sufficient data were collected to exclude areas from mass treatment. The applicability of the methods and the implications of the results are discussed in the general context of large-scale disease control programmes.


Context: control programme and study area

The study was conducted in the context of a national schistosomiasis and soil-transmitted helminths (STH) control programme in Tanzania, which was established in 2003 with support from the Schistosomiasis Control Initiative (SCI) and funded by the Bill and Melinda Gates Foundation. Following WHO guidelines, the programme classifies communities on the basis of prevalence of schistosomiasis in school-age children, according to three strategies: 1) schools where prevalence is <10%, praziquantel is to be made available in local health centres, 2) schools where prevalence is 10 - 50%, mass treatment of all school age-children in the community is conducted and 3) schools where prevalence is >50%, mass treatment of all school age children plus other high-risk groups in the community is conducted. Because of the widespread distribution of STH in Tanzania, albendazole is co-administered with praziquantel.

In the first year of implementation (2005), five regions1 in coastal Tanzania and six regions in northwest Tanzania will be targeted for mass drug administration. Our study was conducted as part of the intervention in northwest Tanzania. For logistical reasons, data collection was limited to a 670 × 530 km area incorporating Kagera, Mwanza, Shinyanga and Tabora regions (see inset, Figure 1), with the aim of making spatial predictions in all six regions (including Mara and Kigoma).

Figure 1
Geographical distribution of schools surveyed in northwest Tanzania and prevalence of (a) Schistosoma haematobium and (b) S. mansoni infections.

Data collection

Parasitological data were collected from at least 60 children (range 60-65) in each of 143 primary schools in northwest Tanzania. Sample sizes (including numbers of schools and children within schools to be surveyed) were determined using simulation studies, based on existing data (Lwambo et al. 1999). Lists of schools from the regional basic education statistics reports were used as the sampling frame. Coordinates of the schools were not available and stratification was conducted according to districts. No two schools were selected from the same ward to ensure good geographical coverage of the district. Sampling started with grade 6 and proceeded down through the grades until 30 boys and 30 girls were selected. If surplus students were present, systematic random sampling was employed. Overall, students were aged 6-22 years, with a mean, median and mode age of 14 and 95% of students were aged between 11 and 17. Urine and stool samples were collected from each student and processed. Egg counts were conducted according to standard parasitological procedures on a single slide from each urine sample and two Kato-Katz slides from each stool sample. Ethical approval was obtained from the Ethical Review Board of National Institute of Medical Research (NIMR), Tanzania and the National Health Service Local Research Ethics Committee (NHS-LREC) of St Mary’s Hospital, London, UK (EC Number: 03 36).

Variable selection

Environmental data included satellite-derived mean land surface temperature (LST) and normalized difference vegetation index (NDVI) for 1982-1998, obtained from the National Oceanographic and Atmospheric Administration’s (NOAA) Advanced Very High Radiometer (AVHRR), elevation, obtained from an interpolated digital elevation model from the Global Land Information System (GLIS) of the United States Geological Survey (, annual rainfall and the location of inland water-bodies. These data were imported into the GIS ArcMap version 8.1 (ESRI, Redlands, CA) and linked according to location to the parasitological data.

Although intensity of schistosome infection has greater relevance to transmission dynamics and morbidity (Anderson & May 1991), prevalence of infection, based on microscopic examination of schistosome eggs in stool or urine samples, remains the most widely used indicator of infection status and is employed by WHO in making recommendations for control. Binomial logistic regression models were constructed in Stata/SE version 8.1 (Stata Corporation, College Station, Texas). Univariable logistic regression was conducted and variables with Wald’s P >0.2 were excluded from further analyses. Collinearity was investigated among all possible pairs of potential predictor variables and if any pair had a correlation coefficient >|0.9|, the member of the pair that was thought less likely to be biologically important was excluded. With the remaining variables, backwards-stepwise binary logistic regression was conducted using Wald’s P >0.1 as the exit criterion and P ≤0.05 as the entry criterion. Non-linear relationships were examined using scatter-plots and entry of categorised predictor variables into the models. In the final model for S. haematobium, minimum LST (categorised into <35°C, 35-39°C and >39°C) and annual rainfall (categorised into <1050 mm and ≥1050 mm) were included as predictors. In the final model for S. mansoni, distance to perennial water bodies (categorised into <0.04, 0.04 - 0.1, 0.1 - 0.4 and >0.4 decimal degrees) and annual minimum temperature (as a continuous variable) were included.

Spatial and non-spatial Bayesian modelling

Spatially-explicit binomial logistic regression was undertaken in WinBUGS version 14 (MRC Biostatistics Unit, Cambridge, UK) with the predictor variables identified above entered as fixed-effects and the spatial structure of the residuals modelled as a Gaussian stochastic process using a geostatistical design (Diggle et al. 1998). Separate models were developed for S. haematobium and S. mansoni, including a Bayesian geostatistical model with covariates and a Bayesian geostatistical model without covariates for each of the two parasitic infections. The reason for including/excluding the covariates was to determine if this had an affect on the fit of the model and to determine whether the covariates accounted for part or all of the spatial correlation in the parasitological data. Additionally, binomial logistic regression models were constructed with the selected predictor variables included as fixed-effects, but without explicitly considering the spatial dependence structure of the data. Again, this was to determine the impact of including the spatial random effect on the fit of the model and also to determine the affect of accounting for spatial correlation on the credible intervals of the coefficients of the selected covariates. Detailed descriptions of the structure of the Bayesian geostatistical models and the process of model assessment are described in Appendix 1.

The deviance information criterion (DIC) statistic was calculated for the Bayesian geostatistical models with and without covariates and the fixed-effects models to determine statistically if addition of the geostatistical component and/or covariates improved model fit (models with a lower DIC statistic are considered to demonstrate a better fit). Maps of the summary statistics of the posterior distributions of predicted prevalence and the spatial random effects were then constructed using ArcMap version 8.1. More details on how spatial predictions were made at non-sampled locations are provided in Appendix 2.

Model validation

Validation was undertaken using independently-collected parasitological data from 54 schools in Magu district (Lwambo et al. 1999). Predicted prevalence at validation locations was compared to observed prevalence (the gold standard) using receiver operating characteristics (ROC) analysis, a method widely applied to diagnostic test evaluation and recently applied to validation of regression models (Brooker et al. 2002b). The gold standard threshold used for the S. haematobium model was ≥0.5, which enabled an assessment of the accuracy of the selected S. haematobium model in delineating high priority areas where all at risk groups will receive mass treatment from lower priority areas where only school-age children will receive mass treatment. The gold standard threshold used for the S. mansoni model was ≥0.1, which enabled an assessment of the accuracy of the selected S. mansoni model in delineating areas receiving mass treatment from areas not receiving mass treatment. Neither the 0.1 threshold for the S. haematobium model nor the 0.5 threshold for the S. mansoni model could be tested because the validation dataset contained no schools with a S. haematobium prevalence less than 0.1 and only two schools with a S. mansoni prevalence greater than 0.5. The statistic used for the comparison was the area under the curve (AUC), which relates to the ROC curve, a plot of sensitivity versus one minus specificity, where values greater than 0.9 indicate an extremely well-fitting model, values greater than 0.7 indicate a moderately well-fitting model and values approaching 0.5 indicating a model that is no improvement on random allocation of test status.

Application of spatial predictions to an intervention plan

The SCI treatment programme does not discriminate between S. haematobium and S. mansoni infections as praziquantel is effective against both species. The prevalence surfaces for S. haematobium and S. mansoni were, therefore, combined to produce a single intervention map, where contours equating to median predicted prevalence of 0.1 and 0.5, were placed in order to delineate areas with different intervention strategies. The uncertainty in the model predictions was also taken into consideration in planning the intervention. It was decided not to exclude areas where the 95% Bayesian credible intervals suggested that, despite a low median predicted prevalence, high prevalence of either infection may occur.


The observed spatial distributions of infection prevalence based on the 143 schools surveyed are presented in Figure 1. The schools with the highest prevalence of S. haematobium infection were located south of Lake Victoria in Mwanza and Shinyanga regions whereas, for S. mansoni, they were located close to the shores of Lake Victoria, particularly in the south-western part of the lake.

Model selection and validation

The Bayesian spatial and non-spatial models for S. haematobium and S. mansoni are presented in Tables Tables11--2.2. Using the DIC, the S. haematobium model containing the geostatistical component but no covariates was the best-fitting, most parsimonious model. Although the covariates were significant when included in the model without the geostatistical component, they became non-significant when it was included. The value of the decay parameter for spatial correlation (ϕ) was 0.2, indicating that a high degree of spatial correlation of S. haematobium prevalence was evident between locations with relatively large distances separating them (Figure 2).

Figure 2
Distance decay of spatial correlation for Schistosoma haematobium in northwest Tanzania. Note: at the equator, 1 decimal degree equates to approximately 120 kilometres.
Table 1
Bayesian models for prevalence of Schistosoma haematobium in northwest Tanzania based on 143 school, 2004. LST refers to land surface temperature and DIC refers to deviance information criterion.
Table 2
Bayesian models for prevalence of Schistosoma mansoni in northwest Tanzania, based on 143 schools, 2004. DIC refers to deviance information criterion.

The S. mansoni model with both environmental covariates and the geostatistical component was the best-fitting, most parsimonious model according to the value of the DIC statistic. The odds ratios indicated that locations adjacent to perennial water-bodies had a much higher prevalence than locations at intermediate distances, which in turn had a much higher prevalence than the reference category (locations at greater than 0.4 decimal degrees from the nearest perennial water body). There was also a positive association between S. mansoni prevalence and annual minimum temperature. In contrast to the S. haematobium model, the value of the decay parameter for spatial correlation (ϕ) was 2.8, indicating that spatial correlation of S. mansoni prevalence occurred over much shorter distances (Figure 3).

Figure 3
Distance decay of spatial correlation for Schistosoma mansoni in northwest Tanzania. Note: at the equator, 1 decimal degree equates to approximately 120 kilometres.

Validation of the selected S. haematobium model using an observed prevalence threshold of ≥0.5 gave an AUC of 0.76 indicating a moderately good predictive performance of the selected model in the district from which the validation data were obtained. The AUC for the selected S. mansoni model using an observed prevalence threshold of ≥0.1 was 0.95, indicating an extremely good predictive performance.

Spatial predictions: Schistosoma haematobium

The outputs of the model are distributions indicating the likely range of values of prevalence at each prediction location and as such, give an indication of the uncertainty surrounding the predictions. Maps of the median and 95% Bayesian credible intervals for the posterior distributions of predicted S. haematobium prevalence are presented in Figure 4. An area of high median predicted prevalence (>50%) was apparent directly south of Lake Victoria (but not adjacent to it). Most of the western part of the study area was predicted to have low median prevalence (<10%). On examination of the 95% Bayesian credible intervals, it was clear that S. haematobium is not likely to occur at levels above 10% in the northwest of the study area (Kagera region) and Urambo district. However, with the wide credible intervals in the south-western part of the intervention area (Kigoma region), predictions could not be made with a sufficient degree of certainty in that area.

Figure 4
Prevalence predictions for Schistosoma haematobium in northwest Tanzania: a) posterior median; b) lower 95% Bayesian credible interval limit; c) upper 95% Bayesian credible interval limit.

Spatial predictions: Schistosoma mansoni

Maps of the median and 95% Bayesian credible intervals for the posterior distributions of predicted S. mansoni prevalence are presented in Figure 5. According to the median prediction, all high-prevalence areas (>50%) were located adjacent to the shores of Lake Victoria, particularly along the southern shore and the islands in the southern part of the lake. Intermediate prevalence was predicted along the shores of other perennial water bodies (e.g. Lake Tanganyika) and the eastern and western shores of Lake Victoria. Examination of the 95% Bayesian credible intervals demonstrated a high level of certainty that S. mansoni occurs at a high prevalence along the southwest shore of Lake Victoria and the islands in the southwest part of the lake, but that high prevalence may also occur adjacent to other parts of the lake and other water-bodies, despite a low to intermediate median prevalence prediction. It also demonstrated with a high level of certainty that a low prevalence of S. mansoni occurs at large distances from perennial water bodies (e.g. in Tabora and Shinyanga regions). Examination of an interpolated surface of the median posterior values of θi, the geostatistical component of the model, may give an indication of the factors that have an important role in determining the distribution of S. mansoni that were not included in the model (Figure 6). Interestingly, the area of low residual S. mansoni prevalence approximates the area of high S. haematobium prevalence.

Figure 5
Prevalence predictions for Schistosoma mansoni in northwest Tanzania: a) posterior median; b) lower 95% Bayesian credible interval limit; c) upper 95% Bayesian credible interval limit.
Figure 6
Geostatistical random-effect component (θi) of the spatially-explicit logistic regression model for Schistosoma mansoni in northwest Tanzania.

Intervention plan

The intervention map, combining both predicted S. haematobium and S. mansoni risk, is presented in Figure 7. Wide credible intervals in the spatial predictions were mainly apparent in the districts in Mara and Kigoma regions from which no schools were sampled. As all of Mara region was placed in the intervention zone, no further evidence was considered necessary, but as all of Kigoma region fell outside the intervention zone, it was decided not to exclude Kigoma without conducting further data collection and model validation in this region. Additionally, it is believed that no transmission of S. mansoni occurs in Lake Tanganyika and model validation in Kigoma region will include an assessment of S. mansoni prevalence along the lake.

Figure 7
Intervention map for northwest Tanzania with prevalence contours defining areas to be excluded from mass treatment (prevalence of Schistosoma mansoni or S. haematobium <0.1), areas to receive mass treatment of school-age children (prevalence of ...


As new initiatives for the control of helminths and other parasites are underway there is a need to geographically stratify areas by risk to guide treatment interventions. Fulfilment of this requirement has been facilitated by developments in GIS/RS and spatial analytical methods. Our study employed statistically robust Bayesian geostatistical models to predict the spatial distributions of S. haematobium and S. mansoni infections in northwest Tanzania, enabling the targeted implementation of an on-going national schistosomiasis control programme. To the authors’ knowledge, this is the first report where spatially-explicit analytical methods have been applied to every stage of the planning and implementation of a large-scale infectious disease control programme.

The application of Bayesian modelling in the current study conferred considerable advantages over traditional, frequentist modelling approaches, in that the spatial dependence structure of the observed data was able to be incorporated in the modelling process in an intuitive and interpretable way and the uncertainties surrounding the predicted prevalence estimates could be fully appreciated, enabling objective decisions to be taken as to the need for further data collection. Failure to account for spatial correlation in the data can lead to inflated significance of regression coefficients and overestimation of the confidence of predictions (Thomson et al. 1999), whereas uncertainty in geographical data may be associated with adverse consequences arising from erroneous decisions taken using the data. Agumya and Hunter (2002) discussed a number of ways in which the decision risks associated with geographical data uncertainty can be managed, including avoidance, reduction, retention and transfer of risk. We took a risk-avoidance approach by deciding not to rely on model predictions in areas where the Bayesian credible intervals suggested that schistosomiasis could be a problem despite low posterior median estimates of prevalence.

In the non-spatial model for S. haematobium and the spatial and non-spatial models for S. mansoni, temperature, rainfall and distance to perennial water bodies were important predictors of either or both parasites. These statistical relationships are consistent and interpretable with the known biology of freshwater snails, the intermediate hosts for schistosomes: S. haematobium is transmitted by snails from the genus Bulinus and S. mansoni is transmitted by the genus Biomphalaria (Brown 1994; Sturrock 1993). Bi. choanomphala in Lake Victoria and Bi. sudanica, which is found adjacent to the lake, are the main intermediate snail hosts of S. mansoni in the study area (Lwambo et al. 1999). This distribution of Biomphalaria spp. hosts explains, at least in part, the high prevalence of S. mansoni in schools at the lakeshore and the low prevalence away from the lakeshore. In contrast, S. haematobium has three snail intermediate hosts in the area: Bulinus nasutus, Bu. africanus and Bu. globosus, which occupy a mosaic of habitats throughout the area such that transmission of the parasite is widespread (Webbe 1962). Malacological studies have demonstrated that the population dynamics of both Bulinus spp. and Biomphalaria spp. are sensitive to temperature and rainfall (Sturrock 1993; Appleton 1978; Brooker & Michael 2000).

It is well recognised that schistosomiasis, like many parasitic diseases, has a focal distribution that is associated with spatial correlation in observed parasitological data. The importance of accounting for spatial correlation was highlighted by the fact that the coefficients for the predictors in the S. haematobium model became non-significant and although remaining significant, the credible intervals surrounding the coefficients in the S. mansoni model were much wider after accounting for spatial correlation. Incorporation of the spatial dependence structure of the data made it apparent that, despite the known biological importance of the environmental covariates, the statistical relationships in the non-spatial models were not supported by the data and spurious significant relationships between the covariates and prevalence would have been accepted had spatial correlation not been considered. The different range of spatial correlation for the two species of schistosome at least partly reflects differing ecologies of the two parasites, most probably due to the different environmental requirements of their respective intermediate hosts.

The statistical technique used in the current study required that a number of assumptions be made - the most significant of which is that the spatial correlation structure does not vary across the study area (the assumption of stationarity), although methods have been developed recently to deal with non-stationarity in a Bayesian geostatistical framework (Gemperli 2003). Additionally, the method used has not been extended to incorporate potential anisotropy (i.e. where spatial correlation is stronger in one or more directions). The main practical drawback of the approach used in the current study in comparison to previous (non-Bayesian) efforts remains the relative complexity and computational difficulty of Bayesian modelling, despite advances such as Gibbs sampling, improved computer hardware with ever-faster processors and the development of the freely-available WinBUGS programme.

An important issue is data quality and the relationship it has with the accuracy of predictions and, ultimately, the population-level impact of the intervention on the outcomes of interest (in our case, reducing morbidity, mortality and the socioeconomic consequences of schistosomiasis). Often, spatial analyses are conducted retrospectively using data from a range of different sources, collected during surveys that have divergent goals (e.g. the MARA/ARMA database (MARA/ARMA 1998)) and the data are therefore not optimal for the purpose of disease mapping. In the current study, spatial analysis and disease mapping were considered as integral components of the planning process and the study and sampling designs were specifically aimed at obtaining statistically efficient, high-quality data for the production of accurate prevalence maps. One potential outcome of demonstrating the practical benefits of high-quality maps to a large scale control programme such as SCI is that future programmes will pay more attention to the need for collecting optimal data for subsequent spatial analysis.

As well as employing a more robust analytical approach to predict disease prevalence, an important aim of disease-mapping in the current study was to define target populations for a national schistosomiasis control programme. However, disease mapping is one of a number of alternative approaches that may be employed to achieve this aim. These alternatives include targeting all individuals without any attempt to define high and low prevalence areas, lot quality assurance sampling (where the minimum number of individuals required to give an assessment of whether prevalence is above a pre-determined “trigger” threshold, with a given confidence, is sampled in every location) (Brooker et al. 2005b), targeting all communities within a specified distance of a known infection source (e.g. Lake Victoria for S. mansoni) (Brooker 2002) and the use of morbidity questionnaires (for self-reported urinary schistosomiasis) (Lengeler et al. 2002). An important area of further research to support (or refute) the credibility of spatial analysis and disease-mapping as tools for planning large-scale interventions is to compare them to the alternative methods described above using economic criteria, such as those used in cost-effectiveness analyses. With each alternative, the relative costs (e.g. monetary costs, human resources and time) and the relative efficacy (in terms of the proportion of locations correctly classified, a function of the sensitivity and specificity of each method) need to be balanced in the context of the wider logistic and operational issues. It will also be important to determine how different targeting approaches alter the allocation of financial resources at local levels and whether this substantially affects the health impact of the programme.

In conclusion, geographic targeting has tremendous potential to enhance the cost-effectiveness of national disease control programmes. Bayesian geostatistical analysis has proved a powerful and statistically robust tool for identifying high prevalence areas in a heterogeneous and imperfectly-known environment. The maps we derived using these methods are now being employed at the district level in Tanzania to deliver the intervention where it is needed most.


We sincerely extend thanks to the children and school teachers who took part in the parasitological surveys, to staff at the National Institute for Medical Research (Petro Mnyeshi, Silinus Nyanda, John Igogote, Philbert Kashangaki and Sahani Kaboya) for conducting the parasitological examinations, to District and Regional Health and Education Officers in northwest Tanzania for assisting survey teams, to staff of the school health programme located in the central offices of the Tanzanian Ministry of Health (Mohamed Mmole, Ali Ngomelo, Isaak Njau, Edith Ufwinki and Augusta Assey) and Ministry of Education (Mrs Urassa), for their administrative assistance, to Dr Ali Mzige, head of Preventive Services, Tanzanian Ministry of Health, for his support of the SCI programme, to Prof. Charles Kihamia of Muhimbili University College of Health Sciences for providing scientific and administrative support and to Dr John Changalucha, director of the National Institute for Medical Research, Mwanza for supporting the surveys. We also thank Dr Nicky Best of Imperial College London for giving advice on Bayesian statistical modelling, Sarah Whawell, Profs. Christl Donnelly and Roy Anderson of Imperial College London and Profs. Mike Kenward and Simon Cousens of London School of Hygiene and Tropical Medicine for statistical and epidemiological advice on planning the data collection and Dr Simon Hay of Oxford University for providing access to the satellite data. Funding for the activities of the Schistosomiasis Control Initiative (SCI) in Tanzania is generously provided by the Bill and Melinda Gates Foundation. Simon Brooker is supported by a Wellcome Trust Advanced Training Fellowship (073656).

Appendix 1: Model structure and iterative model assessment

The Bayesian geostatistical models were of the form:



where Yi is the observed number positive at location i, ni is the number tested at location i, pi is predicted prevalence at location i, α is the intercept, NβN×xN,i is a vector of N predictor variables measured at each location i multiplied by their coefficients and θi, the residual spatial component, is defined by a powered exponential spatial correlation function:


where dij are the distances between pairs of points i and j, ϕ is the rate of decline of spatial correlation with distance and κ is the degree of spatial smoothing (1). Non-informative priors were specified for the intercept (uniform prior with bounds - ∞ and ∞) and the coefficients (normal prior with mean = 0 and precision = 1 × 10-4). The prior distribution of κ was uniform with upper and lower bounds set at 0.5 and 1.50 respectively. The prior distribution of ϕ was also uniform with upper and lower bounds set at 0.1 and 6.0, which gave possible values for spatial correlation of 0.99 - 0.55 with a separating distance of 0.1 decimal degrees (the minimum distance between observed data points) and possible values for spatial correlation of 0.00 - 0.55 with a separating distance of 6.0 decimal degrees (the maximum distance between observed data points), assuming κ = 1.0.

Three chains of the models were run consecutively. A burn-in of 1,000 iterations was allowed, followed by 10,000 iterations where values for the intercept and coefficients were stored. Diagnostic tests for convergence of the stored variables were undertaken, including visual examination of history and density plots of the three chains and visual analysis of the Brooks, Gelman and Rubin statistic (2). Convergence was successfully achieved after 10,000 iterations for the S. mansoni models and 40,000 iterations for the S. haematobium models. The chains were also examined for autocorrelation by visual assessment of the in-built autocorrelation function of WinBUGS. Considerable autocorrelation was apparent for the intercept and some coefficients and it was decided to reduce autocorrelation by thinning subsequent sampling by only storing every third iteration. A total of 10,000 iterations were then run, giving 10,000 stored values from the posterior distribution of each variable (3 chains × 10,000 iterations/3). Descriptive statistics for the posterior distributions of each variable were calculated and analysed. In order to determine whether sufficient iterations had been conducted to fully describe the posterior distributions, Monte Carlo error (MCE)/SD was calculated for each variable and if MCE/SD was less than 0.05, it was decided that sufficient iterations had been conducted. With all stored variables, MCE/SD was less than 0.05 after 10,000 values were generated and stored.

  1. Spiegelhalter D, Thomas A, Best N & Lunn D (2003) WinBUGS user manual. Available at
  2. Brooks SP & Gelman A (1998) Alternative methods for monitoring convergence of iterative simulations. Journal of Computational Graphical Statistics 7, 434-455.

Appendix 2. Spatial predictions

A set of prediction locations was specified according to a 0.05×0.05 decimal degree square grid where the raw data suggested that most of the spatial variation of each infection occurred and a 0.1×0.1 decimal degree square grid covering the rest of the study area, where little or no infection was apparent. Grid dimensions were calculated according to computational limits and to give a meaningful prediction density for the intervention. The values of the predictor variables were then determined for each of the prediction locations using ArcMap version 8.1.

Predicted prevalence was calculated for each prediction location using the Bayesian model that gave the lowest DIC statistic and the “Spatial.unipred” function of WinBUGS (1), which solves the model equation at each prediction location given the values of each covariate at the prediction location and the distance between prediction locations and observed data locations.

Thomas A, Best N, Lunn D, Arnold R & Spiegelhalter D (2004) Geobugs user manual. Available from


1Region is the first administrative division of Tanzania. Each region is divided into districts (the second administrative division) and each district into wards (the third administrative division).


  • Agumya A, Hunter GJ. Responding to the consequences of uncertainty in geographical data. International Journal of Geographical Information Science. 2002;16:405–417.
  • Alexander N, Moyeed R, Stander J. Spatial modelling of individual-level parasite counts using the negative binomial distribution. Biostatistics. 2000;1:453–463. [PubMed]
  • Anderson RM, May RM. Infectious Diseases of Humans: dynamics and control. Oxford: Oxford University Press; 1991.
  • Appleton CC. Review of Literature on abiotic factors influencing the distribution and life-cycles of Bilharziasis intermediate host snails. Malacological Review. 1978;11:1–25.
  • Best N, Richardson S, Thomson A. A comparison of Bayesian spatial models for disease mapping. Statistical Methods in Medical Research. 2005;14:35–59. [PubMed]
  • Brooker S, Michael E. The potential of geographical information systems and remote sensing in the epidemiology and control of human helminth infections. Advances in Parasitology. 2000;47:246–288. [PubMed]
  • Brooker S, Hay SI, Issae W, et al. Predicting the distribution of urinary schistosomiasis in Tanzania using satellite sensor data. Tropical Medicine and International Health. 2001;6:998–1007. [PMC free article] [PubMed]
  • Brooker S. Schistosomes, snails and satellites. Acta Tropica. 2002;82:209–216. [PubMed]
  • Brooker S, Beasley M, Ndinaromtan M, et al. Use of remote sensing and a geographical information system in a national helminth control programme in Chad. Bulletin of World Health Organization. 2002a;80:783–789. [PubMed]
  • Brooker S, Hay SI, Bundy DAP. Tools from ecology: useful for evaluating infection risk models. Trends in Parasitology. 2002b;18:70–74. [PMC free article] [PubMed]
  • Brooker S, Clements ACA, Bundy DAP. Global mapping and control of soil-transmitted helminthiasis. In: Hay SI, Graham AJ, Rogers DJ, editors. Global mapping of infectious diseases: methods, examples and emerging applications. London: Academic Press; 2005a. (in press)
  • Brooker S, Kabatereine NB, Myatt M, Stothard JR, Fenwick A. Rapid assessment of Schistosoma mansoni: the validity, applicability and cost-effectiveness of the Lot Quality Assurance Sampling method in Uganda. Tropical Medicine and International Health. 2005b;10:647–658. [PMC free article] [PubMed]
  • Brown DS. Freshwater Snails of Africa and their Importance. London: Taylor and Francis; 1994.
  • Carabin H, Escalona M, Marshall C, Vivas-Martinez S, Botto C, Joseph L, et al. Prediction of community prevalence of human onchocerciasis in the Amazonian onchocerciasis focus: Bayesian approach. Bulletin of World Health Organization. 2003;81:482–490. [PubMed]
  • Diggle PJ, Tawn JA, Moyeed RA. Model-based geostatistics. Applied Statistics. 1998;47:299–350.
  • Diggle P, Moyeed R, Rowlingson B, Thompson M. Childhood malaria in the Gambia: a case-study in model-based geostatistics. Applied Statistics. 2002;51:493–506.
  • Elith J, Burgman MA, Regan HM. Mapping epistemic uncertainties and vague concepts in predictions of species distributions. Ecological Modelling. 2002;157:313–329.
  • Elnaiem DE, Schorscher J, Bendall A, et al. Risk mapping of visceral leishmaniasis: the role of local variation in rainfall and altitude on the presence and incidence of kala-azar in eastern Sudan. American Journal of Tropical Medicine and Hygiene. 2003;68:10–17. [PubMed]
  • Gemperli A. Development of spatial statistical methods for modelling point-referenced spatial data in malaria epidemiology. University of Basel; 2003. pp. 111–134. PhD thesis.
  • Gemperli A, Vounatsou P, Kleinschmidt I, Bagayoko M, Lengeler C, Smith T. Spatial patterns of infant mortality in Mali: the effect of malaria endemicity. American Journal of Epidemiology. 2004;159:64–72. [PubMed]
  • Hay SI, Omumbo JA, Craig MH, Snow RW. Earth observation, geographical information systems and Plasmodium falciparum malaria in sub-Saharan Africa. Advances in Parasitology. 2000;47:174–215. [PMC free article] [PubMed]
  • Kabatereine NB, Brooker S, Tukahebwa EM, Kazibwe F, Onapa A. Epidemiology and geography of Schistosoma mansoni in Uganda: implications for planning control. Tropical Medicine and International Health. 2004;9:372–380. [PubMed]
  • Kamel Boulos MN. Towards evidence-based, GIS-driven national spatial health information infrastructure and surveillance services in the United Kingdom. International Journal of Health Geographics. 2004;3 Available at [PMC free article] [PubMed]
  • Lengeler C, Utzinger J, Tanner M. Questionnaires for rapid screening of schistosomiasis in sub-Saharan Africa. Bulletin of World Health Organization. 2002;80:235–242. [PubMed]
  • Lindsay SW, Thomas CJ. Mapping and estimating the population at risk from lymphatic filariasis in Africa. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2000;94:37–45. [PubMed]
  • Lwambo NJS, Siza JE, Brooker S, Bundy DAP, Guyatt H. Patterns of concurrent infection with hookworm and schistosomiasis in school children in Tanzania. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1999;93:497–502. [PubMed]
  • Malone JB, Yilma JM, McCarroll JC, Erko B, Mukaratirwa S, Xinyu Zhou. Satellite climatology and the environmental risk of Schistosoma mansoni in Ethiopia and East Africa. Acta Tropica. 2001;79:59–72. [PubMed]
  • MARA/ARMA Towards and atlas of malaria risk in Africa, first technical report of the MARA/ARMA collaboration. 1998. Available at http//
  • Moodley I, Kleinschmidt I, Sharp B, Craig M, Appleton C. Temperature-suitability maps for schistosomiasis in South Africa. Annals of Tropical Medicine and Parasitology. 2003;97:617–627. [PubMed]
  • Raso G, Matthys B, N’Goran EK, Tanner M, Vounatsou P, Utzinger J. Spatial risk prediction and mapping of Schistosoma mansoni infections among schoolchildren living in western Côte d’Ivoire. Parasitology. 2005;131:1–12. [PubMed]
  • Rogers DJ. Satellites, space, time and the African Trypanosomiases. Advances in Parasitology. 2000;47:129–171. [PubMed]
  • Rogers DJ, Randolph SE, Snow RW, Hay SI. Satellite imagery in the study and forecast of malaria. Nature. 2002;415:710–715. [PMC free article] [PubMed]
  • Sturrock RF. The Intermediate Hosts and Host-Parasite Relationships. In: Jordan P, Webbe G, Sturrock RF, editors. Human schistosomiasis. Wallingford: CAB International; 1993. pp. 33–85.
  • Thomson MC, Connor SJ, D’Alessandro U, et al. Predicting malaria infection in Gambian children from satellite data and bed net use surveys: the importance of spatial correlation in the interpretation of results. American Journal of Tropical Medicine and Hygiene. 1999;61:2–8. [PubMed]
  • Webbe G. The transmission of Schistosoma haematobium in an area of Lake Province, Tanganyika. Bulletin of World Health Organization. 1962;27:59–85. [PubMed]