the objective of this study was to produce maps that could be integrated into the sci - supported national intervention strategies and that explicitly represent uncertainties in spatial predictions so that national control managers could judge the quality of the evidence upon which the strategies will be based .
the sci - supported programs involve mass distribution of praziquantel ( for urinary and intestinal schistosomiasis ) and albendazole ( for soil - transmitted helminths ) .
the parasitic infection with the highest prevalence is urinary schistosomiasis , caused by flukes ( schistosoma hematobium ) , and the programs are planned to control this disease ( 2 ) .
parasitologic data were collected in coordinated school - based field surveys in burkina faso , mali , and niger ( figure 1 ) during 20042006 ( preintervention ) by using standardized protocols ( available on request ) .
the collated dataset covered a spatially contiguous area , 2,750 km  850 km , and included the infection status of 27,939 school - age children in 418 randomly selected locations .
infection status was defined according to egg count determined by microscopic examination of urine samples ; > 1 s. hematobium eggs indicated infection .
prevalence of infection with schistosoma hematobium at 418 survey locations in burkina faso , mali , and niger , 20042006 .
spatial prediction was based on a logistic regression model ( table ) , constructed by using the software winbugs , version 1.4.2 ( mrc biostatistics unit , cambridge and imperial college , london , uk ) .
the model had infection status as the binary outcome variable , age and sex of the survey participants as individual - level fixed effects , and distance from perennial water body ( derived from electronic maps obtained from the food and agriculture organization ) and land surface temperature ( lst ; with a quadratic term ; see hay et al .   for details on how these data were derived ) as survey location  level fixed effects .
the model also included a geostatistical random effect for residual spatial clustering of infection prevalence ( 11 ) .
values for the fixed effects are odds ratios ; note the odds ratios for the climate variables are on a common scale , where the variables were standardized to have a mean = 0 and sd = 1 .
the reference group for sex was boys and for age was 68 y. the number of children found to be infected with s. haematobium was modeled by using a binomial distribution described by the proportion infected and the total number sampled in each survey location .
the proportion infected was modeled by using logistic regression with an intercept , covariates ( sex , age , distance to perennial water body , land surface temperature , and a quadratic term for land surface temperature ) , and a random effect that described spatial correlation ( i.e. , clustering ) .
model outputs were distributions ( termed posterior distributions ) that can be summarized by using the mean , sd , and 95% cri ( representing the range of values that contains the true value with a probability of 95% ) .
a prevalence map for the study area was constructed , using the model , by predicting infection prevalence at the centroids of cells of a 0.15  0.15 decimal degree ( 18 km  18 km ) grid .
this model was implemented with the spatial.unipred command of winbugs ( details are provided in the technical appendix ) .
estimates from bayesian models are distributions ( termed posterior distributions ) that represent the probability of each of a range of plausible values being true for the variable being modeled . to quantify the uncertainties surrounding the model predictions
, we plotted the probability of each prediction location having a prevalence > 50% , rather than mean predicted prevalence at each location .
the probabilities were calculated from the posterior distributions of predicted prevalence at each location ( i.e. , if 95% of the posterior distribution of predicted prevalence was > 0.5 , the probability of prevalence > 50% at that location was 95% ) .
cross - validation was done by randomly allocating survey locations to 3 groups and undertaking 3 separate runs of the model ; 1 of the 3 groups was sequentially omitted , and predicted prevalence at the omitted locations was determined by using the model .
predicted prevalence was compared with observed prevalence , dichotomized , according to a 50% observed prevalence threshold .
the comparison statistic was the area under the curve ( auc ) of the receiver operating characteristic , and a value of > 0.7 was considered to indicate acceptable predictive ability .
an average auc was calculated across the 3 model runs . in the final model ( table )
, statistically significant correlations suggested that infection prevalence was higher in older boys and increased with proximity to perennial bodies of water , but no association was found between prevalence and lst .
the range over which spatial correlation was > 5% ( chosen to indicate statistically important spatial correlation ) was 177 km , indicating the approximate radius of clusters .
results of the validation analysis showed an average auc of 0.86 , indicating that the model had an acceptable predictive performance .
bayesian probability maps were produced for each sex and age group , but for illustrative purposes we present predicted probability of prevalence > 50% in boys ages 1316 years ( the group with the highest infection prevalence ; figure 2 ) .
large clusters of prediction locations with a high probability ( i.e. , > 50% ; indicative of low uncertainty ) of prevalence being > 50% intervention threshold were located in a mid - latitudinal band across mali , running from western to central regions , and in the niger river valley region of niger .
smaller clusters were located in various parts of southern and eastern mali , northwestern and northeastern burkina faso , and south - central niger . predicted probability of prevalence of infection with schistosoma hematobium being > 50% in burkina faso , mali , and niger in boys ages 1316 years ; results are based on a bayesian geostatistical model .
the red areas had a low degree of uncertainty that predicted prevalence was > 50% , and the blue areas had a high degree of uncertainty that predicted prevalence was > 50% .
future schistosomiasis control plans should acknowledge uncertainties such as those presented in figure 2 . a possible approach would be to introduce a second threshold for the level of uncertainty that a location is above the intervention prevalence threshold ; if the uncertainty is greater than this second threshold , then the location is excluded until new evidence is obtained that confirms prevalence is above or below the intervention prevalence threshold .
this second uncertainty threshold should be determined by the quantity of resources available for disease control and the level of decision risk deemed appropriate .
in addition to providing an evidence base for distributing resources in 3 west african countries as part of the sci - supported national control programs , the maps presented here have a potential role in maintaining sustainability of schistosomiasis control after sci support ends ( sci is funded through 2009 ) . they can be used as advocacy tools for channeling funds to high - risk populations in the affected countries and , in the likely event that money for schistosomiasis control in these countries becomes more limited after sci support ends , they can be used to ensure that scarce governmental resources are distributed as efficiently as possible .
national coordinators who might face accountability for targeted ( i.e. , unequal ) distribution of resources will benefit from the defendable , scientifically sound methods presented in this article . by focusing on uncertainty in spatial predictions , more flexible tools for disease control can be developed that allow the geographic dimensions of control programs to be scaled and modified according to available resources and acceptable levels of decision risk .
 we aimed to map the probability of schistosoma haematobium infection being > 50% , a threshold for annual mass praziquantel distribution . 
 parasitologic surveys were conducted in burkina faso , mali , and niger , 20042006 , and predictions were made by using bayesian geostatistical models . 
 clusters with > 50% probability of having > 50% prevalence were delineated in each country . 