Utilisation de prédictions par géostatistique bayésienne pour estimer les variations locales de la prévalence des infestations par Schistosoma haematobium en Afrique occidentale
Uso de modelos geoestadísticos bayesianos para predecir las variaciones locales de la infección por Schistosoma haematobium en África occidental
Archie CA ClementsI, *; Sonja FirthI; Robert DembeléII; Amadou GarbaIII; Seydou TouréIV; Moussa SackoV; Aly LandouréV; Elisa Bosqué-OlivaVI; Adrian G BarnettVII; Simon BrookerVIII; Alan FenwickVI
IUniversity of Queensland, School of Population Health, Herston Road, Herston, Qld, 4006, Australia
IIProgramme National de Lutte Contre la Schistosomiase, Ministère de la Santé, Bamako, Mali
IIIProgramme National de Lutte Contre la Bilharziose et les Géohelminthes, Ministère de la Santé Publique et de la Lutte Contre les Endémies, Niamey, Niger
IVProgramme National de Lutte Contre la Schistosomiase, Ministère de la Santé, Ouagadougou, Burkina Faso
VInstitut National de Recherche en Santé Publique, Bamako, Mali
VISchistosomiasis Control Initiative, Imperial College, London, England
VIIInstitute for Health and Biomedical Innovation, Queensland University of Technology, Kelvin Grove, Qld, Australia
VIIILondon School of Hygiene and Tropical Medicine, London, England
OBJECTIVE: To predict the subnational spatial variation in the number of people infected with Schistosoma haematobium in Burkina Faso, Mali and the Niger prior to national control programmes.
METHODS: We used field survey data sets covering a contiguous area 2750 × 850 km and including 26 790 school-age children (5-14 years old) in 418 schools. The prevalence of high- and low-intensity infection and associated 95% credible intervals (CrIs) were predicted using Bayesian geostatistical models. The number infected was determined from the predicted prevalence and the number of school-age children in each km².
FINDINGS: The predicted number of school-age children with a low-intensity infection was 433 268 in Burkina Faso, 872 328 in Mali and 580 286 in the Niger. The number with a high-intensity infection was 416 009, 511 845 and 254 150 in each country, respectively. The 95% CrIs were wide: e.g. the mean number of boys aged 10-14 years infected in Mali was 140 200 (95% CrI: 6200-512 100).
CONCLUSION: National aggregate estimates of infection mask important local variations: e.g. most S. haematobium infections in the Niger occur in the Niger River valley. High-intensity infection was strongly clustered in western and central Mali, north-eastern and north-western Burkina Faso and the Niger River valley in the Niger. Populations in these foci will carry the bulk of the urinary schistosomiasis burden and should be prioritized for schistosomiasis control. Uncertainties in the predicted prevalence and the numbers infected should be acknowledged by control programme planners.
OBJECTIF: Prédire les variations spatiales au niveau infranational du nombre de personnes infestées par Schistosoma haematobium au Burkina Faso, au Mali et au Niger, avant la mise en place des programmes nationaux de lutte contre la schistosomiase.
MÉTHODES: Nous avons utilisé un jeu de données d'enquête sur le terrain couvrant une zone contiguë de 2750 x 850 km et 26 790 enfants d'âge scolaire (5-14 ans), répartis dans 418 écoles. La prévalence des schistosomiases de forte et de faible intensité, ainsi que les intervalles de crédibilité à 95 % associés, ont été prédits à l'aide de modèles géostatistiques bayésiens. Le nombre de personnes infestées a été déterminé à partir de la prévalence prédite et du nombre d'enfants d'age scolaire par km².
RÉSULTATS: D'après les prédictions de l'étude, le nombre d'enfants d'âge scolaire atteints d'une schistosomiase de faible intensité serait de 433 268 au Burkina Faso, de 872 328 au Mali et de 580 286 au Niger. S'agissant des enfants fortement infestés, les prédictions donnaient respectivement 416 009 cas pour le Burkina Faso, 511 845 pour le Mali et 254 150 pour le Niger. Les intervalles de crédibilité à 95 % étaient larges : par exemple, le nombre moyen de garçons de 10 à 14 ans infestés au Mali était de 140 200 (ICr à 95 % : 6200-512 100).
CONCLUSION: Les estimations nationales agrégées masquent d'importantes variations locales : par exemple, la plupart des infestations par S. haematobium relevées au Niger étaient apparues dans la Vallée du Niger. Les cas d'infestation lourde étaient très fortement regroupés à l'Ouest et au centre du Mali, au Nord-est et au Nord-ouest du Burkina Faso et dans la Vallée du Niger, au Niger. Les populations de ces foyers supportent la plus grande part de la charge de schistosomiase urinaire et doivent être considérées comme prioritaires dans la lutte contre la schistosomiase. Les planificateurs de programmes de lutte contre cette maladie doivent être conscients des incertitudes qui pèsent sur les prédictions de la prévalence et du nombre de personnes infestées.
OBJETIVO: Predecir la variación territorial subnacional del número de personas infectadas por Schistosoma haematobium en Burkina Faso, Malí y el Níger antes del inicio de los programas de control nacionales.
MÉTODOS: Usamos conjuntos de datos de encuestas sobre el terreno que abarcaron en total a 26 790 niños en edad escolar (5 a 14 años) de 418 escuelas repartidos en una zona de 2750 x 850 km. Mediante modelos geoestadísticos bayesianos se predijeron la prevalencia de infección de alta y baja intensidad y los correspondientes intervalos de credibilidad (ICr) del 95%. El número de personas infectadas se determinó a partir de la prevalencia predicha y del número de niños en edad escolar por km².
RESULTADOS: Las predicciones sobre el número de niños en edad escolar con infección de baja intensidad fueron de 433 268 en Burkina Faso, 872 328 en Malí y 580 286 en el Níger. El número de casos de infección de alta intensidad fue de 416 009, 511 845 y 254 150, respectivamente. Los ICr95% fueron amplios: p.ej., la media de muchachos de 10 a 14 años infectados en Malí fue de 140 200 (ICr95%: 6200-512 100).
CONCLUSIÓN: Las estimaciones totales nacionales de infección ocultan variaciones locales importantes: p.ej., en el Níger la mayoría de las infecciones por S. haematobium se dan en el valle del Río Níger. La infección de alta intensidad se concentra marcadamente en el centro y oeste de Malí, las zonas nororiental y noroccidental de Burkina Faso y el valle del Río Níger en el Níger. Las poblaciones de esos focos soportarán la mayor carga de esquistosomiasis urinaria y deberían ser un objetivo prioritario de la lucha contra la esquistosomiasis. Los planificadores de los programas de control deben tener en cuenta la incertidumbre asociada a la prevalencia predicha y las cifras de personas infectadas.
An accurate estimate of the proportion of a population affected by a disease is important for prioritizing the control of that disease relative to another and for allocating resources to control or prevention programmes. Since high-quality surveillance data for developing countries is frequently lacking, calculations of the burden of a disease are often based on prevalence estimates from cross-sectional surveys, which are rarely randomized or representative of the whole population. For schistosomiasis, caused by trematodes of the genus Schistosoma, a 2003 review reported that 207 million people were infected globally and 779 million people were at risk, the majority in sub-Saharan Africa.1 For most countries in the region, the numbers were derived from prevalence estimates contained in a 1989 report2 that provided national aggregate data even though the prevalence of schistosomiasis is known to be geographically heterogeneous. In addition, previous reports have ignored uncertainties in prevalence estimates and in the size and spatial distribution of the population at risk.
Recently, empirical maps of tropical infectious diseases have been used to improve estimates of the populations infected and at risk at the continental or global level3-5 and, increasingly, to plan and target control programmes. Advances in the production of these maps include geostatistical prediction of the prevalence of infection with Schistosoma haematobium (the aetiological agent of urinary schistosomiasis),6 of other parasitic infections7-10 and of coinfections11 using Bayesian methods.12 The Bayesian approach is advantageous because the effect of covariates and spatial heterogeneity, or clustering, can be modelled simultaneously and the uncertainty of predictions can be assessed.
While the prevalence of schistosome infection has been used for planning large-scale control programmes and for estimating the disease burden, the intensity of the infection is more informative for estimating morbidity, such as urinary tract lesions and anaemia,13-15 and plays a greater role in driving transmission than prevalence. Therefore, maps showing the average infection intensity, or the distributions of low- and high-intensity infection, might provide more effective tools than prevalence maps for estimating the disease burden and developing intervention strategies. An important statistical issue in analysing the intensity of parasitic infections is overdispersion, or aggregation, which occurs when most individuals have few parasites but a few individuals have many.16 Bayesian models have previously been used in the spatial analysis17-19 and prediction20 of the intensity of parasitic infections, due for example to Wuchereria bancrofti and Schistosoma mansoni, with overdispersion in individuals' parasite or egg counts being modelled using the negative binomial distribution.
Burkina Faso, Mali and the Niger, three contiguous countries in the Sahelian zone of western Africa, recently conducted, coordinated, national cross-sectional, school-based, parasitological surveys.21 These surveys were unprecedented in their size and covered approximately 2750 km × 850 km, 26 790 school-age children and 418 schools. We aimed to use survey data to predict subnational spatial distributions of the prevalence of low- and high-intensity S. haematobium infection and to use the prediction maps to calculate the numbers of school-age children infected or at risk. We also aimed to estimate uncertainties in the predicted prevalence and the numbers infected.
Selection of schools and children
The most prevalent parasitic infection in all three study countries is S. haematobium. Programmes supported by the Schistosomiasis Control Initiative were primarily designed to control urinary schistosomiasis21 and our analysis includes only survey data on S. haematobium. Ethical approval for data collection was obtained from St Mary's Hospital Research Ethics Committee in the United Kingdom, the National Public Health Research Institute's (INRSP) scientific committee in Mali, the Ministry of Health's ethics and scientific committees in Burkina Faso and the Ministry of Health's ethical committee in the Niger.
Sample sizes were calculated using historical data from Mali22. It was decided to survey 87 schools in Burkina Faso, 226 in Mali and 215 in the Niger, and to include 60 children at each school. Ultimately, only 418 of the 528 schools were surveyed because remote, sparsely-populated areas had to be excluded for logistical reasons. Geographical coverage was maximized using different spatial stratification methods in the three countries. In Mali and the Niger, sample frames that contained the location of all communities were used. Spatial stratification was performed by overlaying a 1-decimal degree square grid on these countries in the ArcView version 9 (ESRI, Redlands, CA, United States of America) geographical information system. Communities were selected from the cells using simple random selection and, if they had more than one school, the study school was selected using simple random selection when the survey team arrived. In Burkina Faso, lists of schools were available for each province but were not georeferenced. The number of schools selected in each province was weighted according to the area of the province. School sampling was then done in each province using simple random selection.
The surveys were conducted from 2004-2006. The survey team determined the school's coordinates using a global positioning system. If there were fewer than 50 boys or girls in a school, all individuals of that sex were selected because compliance was more difficult if a minority were excluded. If there were more than 50 boys or girls, 30 individuals of that sex were selected using systematic random sampling. Urine and stool samples were collected from each child and processed using standard parasitological methods. The 10-ml urine samples were passed through a filter and examined by microscope in the field. The S. haematobium egg count was recorded and entered into a Microsoft Access database (Microsoft, Redmond, WA, USA).
Predicting the prevalence of infection
Spatial prediction was based on Bayesian geostatistics.23 Rather than modelling egg counts using the negative binomial distribution, we used a multinomial model in which the egg count was categorized as representing: (i) no infection, (ii) low-intensity infection (i.e. 1-50 eggs per 10 ml urine), or (iii) high-intensity infection (i.e. > 50 eggs per 10 ml urine). There were two reasons: (i) expediency, given that in Burkina Faso extremely high egg counts were recorded as > 1000 eggs per 10 ml urine, meaning that the upper tail of the distribution was truncated, and (ii) to facilitate future estimation of the burden of schistosomiasis because existing evidence for related morbidity is based on stratified egg counts, often using WHO definitions of low and high intensity.24
It has been shown that age and sex are associated with the prevalence of urinary schistosomiasis, probably due to physiological differences in susceptibility.25,26 Distance from a perennial inland water body is a plausible risk factor for exposure to schistosomes and subsequent infection because transmission requires contact with the aquatic habitat of intermediate host snails of the Bulinus genus. Distances were derived from an electronic perennial inland water body map obtained from the Food and Agriculture Organization of the United Nations (UN). The effect of temperature and rainfall on the distribution of Bulinus snails is reviewed in Rollinson et al.27 Satellite-derived mean values for the land surface temperature and normalized difference vegetation index (a proxy for rainfall) for 1982-1998 were obtained from the National Oceanographic and Atmospheric Administration's Advanced Very High Radiometer. The initial candidate set of variables included the individual participant variables of sex and age (categorized as 5-9 years and 10-14 years) and the school-level ecological variables of distance from a perennial inland water body, land surface temperature and normalized difference vegetation index. We tested nominal and ordinal multinomial regression models and found that a nominal model provided a better fit. Variables were selected using fixed-effects multinomial regression models in the Stata/SE 10.0 statistical package (StataCorp, College Station, TX, USA). The normalized difference vegetation index was excluded because Wald's P was > 0.2. All remaining variables were selected for inclusion in the spatial model. This model and details of how predictions were made are presented in Box 1. The outputs of Bayesian models, including parameter estimates and spatial predictions, are termed posterior distributions. These distributions fully represent uncertainties associated with estimated values. We summarized the posterior distributions in terms of the posterior mean and 95% credible interval (CrI), within which the true value occurs with a probability of 95%.
Calculating the number infected
An electronic population surface for the study area was obtained from the Global Rural-Urban Mapping Project (GRUMP) alpha version28 and imported into ArcView. In practice, GRUMP is a 30-arc second (1-km²) population raster data set that combines year 2000 census data at a subnational level with an urban-extent mask. In GRUMP, the population is redistributed using an algorithm that assumes a greater proportion is located in urban areas.29
Country-specific population growth rates and the proportions of the population in given sex and age groups (i.e. 5-9 years and 10-14 years) were obtained from the UN Population Division publication World population prospects30 and were used to generate population surfaces for 2005. Surfaces representing the mean and lower and upper 95% CrI limits of the predicted prevalence in each age-sex group were multiplied by the projected 2005 population surfaces using the Spatial Analyst Extension of ArcView, thereby giving the predicted number infected per km² pixel. These numbers were summed for each country.
Calculating the number at risk
By design, the multinomial model gave a predicted prevalence that was non-zero at all locations, even at those where field data indicated the absence of infection. A receiver operating characteristic analysis was conducted to determine the optimal threshold (i.e. where sensitivity = specificity) of the combined predicted prevalence (i.e. low- plus high-intensity infection prevalence) that best discriminated between schools with a zero and non-zero observed prevalence. With this approach, a predicted prevalence threshold of 5.3% gave the best discriminatory performance. A mask, created in the geographical information system to exclude areas with a combined predicted prevalence < 5.3%, was overlaid on the different population surfaces for each country to calculate the numbers and proportions at risk of infection in both the school-age and total population.
Predicted prevalence of infection
The total number of children aged 5-14 years included in the 2004-2006 surveys was 4808 in Burkina Faso, 14 586 in Mali and 7396 in the Niger. The raw prevalence of low-intensity infection was 10.0% (95% confidence interval, CI: 9.2-10.9) in Burkina Faso, 25.2% (95% CI: 24.5-25.9) in Mali and 10.1% (95% CI: 9.4-10.8) in the Niger. For high-intensity infection, the raw prevalence was 8.5% (95% CI: 7.7-9.3) in Burkina Faso, 11.4% (95% CI: 10.9-11.9) in Mali and 3.4% (95% CI: 3.0-3.8) in the Niger. A map of the raw prevalence of S. haematobium infection is presented in Fig. 1.
The spatial model is presented in Table 1. It can be seen from the 95% CrI of the quadratic term for the land surface temperature that the association between this variable and the prevalence of low- or high-intensity infection was not significant. However, the distance from a perennial inland water body was significantly and negatively associated with the prevalence of both low- and high-intensity infection. In addition, the prevalence of low- and high-intensity infection was significantly higher in boys and in children aged 10-14 years than in girls or children aged 5-9 years. The rate of decay of spatial correlation was higher for low-intensity infection than high-intensity infection and the variance of the spatial random effect (i.e. the sill in geostatistical terms) was higher for high-intensity infection than low-intensity infection, which indicates a stronger propensity for spatial clustering for high-intensity infection.
Separate prediction maps were produced for boys and girls and for ages 5-9 years and 10-14 years. The mean predicted prevalences of low- and high-intensity infection in boys aged 10-14 years (the highest prevalence group) are presented in Fig. 2 and Fig. 3, respectively, as illustrative examples. The maps for other age-sex groups (available from the corresponding author) showed the same spatial distribution but a lower predicted prevalence, reflecting the lower odds ratios (ORs) for infection in these groups. There were differences between the spatial distributions for low- and high-intensity infections: low-intensity infection was more widespread but the variation in prevalence was less extreme, with a predicted prevalence between 10% and 50% in most mapped areas of Burkina Faso (excluding the south-west), Mali (excluding the far south) and the Niger (excluding central regions). High-intensity infection had a more restricted spatial distribution, with defined clusters of high prevalence (> 50%) located in a mid-latitude band from western to central Mali, in northern and central Burkina Faso and in the Niger River valley in the Niger. There were large areas of low prevalence (< 5%) in southern, northern and eastern Mali, south-western Burkina Faso and most of the Niger, excluding the Niger River valley. The 95% CrI maps showed wide uncertainty in the predicted prevalence, though the prediction model was clearly able to exclude parts of each country as being at-risk areas for significant transmission of S. haematobium, as indicated by an upper 95% CrI limit for the predicted prevalence of < 5%). Fig. 4 and Fig. 5 (available at: http://www.who.int/bulletin/volumes/87/12/08-058933/en/index.html) illustrate the situation in boys 10-14 years of age; maps for other sex and age groups can be obtained from the corresponding author.
Number infected and at risk
The mean numbers of school-age children with low- and high-intensity infection in the three countries and their associated 95% CrIs are presented in Table 2. The estimated number of school-age children with low-intensity infection was 433 268 in Burkina Faso, 872 328 in Mali and 580 286 in the Niger and the number with high-intensity infection was 416 009 in Burkina Faso, 511 845 in Mali and 254 150 in the Niger. The 95% CrIs for the number infected in each age-sex group were very wide: e.g. the mean number of boys aged 10-14 years infected in Mali was 140 200 (95% CrI: 6200-512 100). Maps of the estimated number of boys aged 10-14 years with low- and high-intensity infection (Fig. 6 and Fig. 7, respectively) showed considerable within-country variation in the burden of schistosomiasis. This was also apparent for boys and girls of both age groups. As was observed with the predicted prevalence, the spatial patterns were identical for each age-sex group, though the proportion infected was uniformly lower for girls and younger boys (maps available from the corresponding author).
We estimated that 3.5 million, 3.4 million and 2.8 million school-age children and a total of 12.5 million (94.5% of the total population), 11.8 million (87.6% of the total) and 9.9 million (70.3% of the total) individuals of all ages were at risk of urinary schistosomiasis in Burkina Faso, Mali and the Niger, respectively. The maps (available from the corresponding author) show that there were parts of each country where people were not at risk and that could be excluded from active surveillance or nationally coordinated schistosomiasis control.
We used robust, contemporary statistical methods in a novel application to predict the spatial distribution of S. haematobium infection. This resulted in estimates of local heterogeneity in high- and low-intensity parasitic infection that could be used in control programme planning. In 2003, Steinmann et al.1 estimated the number of people infected with Schistosoma spp. was 7.8 million in Burkina Faso, 7.8 million in Mali and 3.2 million in the Niger by using prevalence estimates of 60.0%, 60.0% and 26.7% and assuming that the at-risk population was 100% of the estimated population, namely 13.0 million, 13.0 million and 12.0 million in the three countries, respectively. If we assume, as Steinmann et al. did, that the prevalence is the same for all age groups (which overestimates the prevalence since it is usually highest in school-age children),31 our estimates would be 3.0 million (23.0% of the total) in Burkina Faso, 4.8 million (35.4% of the total) in Mali and 3.0 million (21.1% of the total) in the Niger. Clearly, the previously reported number infected was considerably overestimated for Burkina Faso and Mali. In addition, the number at risk was overestimated for all three countries. We have confidence in our conclusions because our data were recent, extensive, randomized and representative. Although our calculations excluded S. mansoni infection, which Steinmann et al. included, its prevalence in our surveys was very low in Burkina Faso and the Niger, at 0.5% and 0.3%, respectively, and moderately low in Mali, at 6.7%. More importantly, Steinmann et al.'s figures overlooked important heterogeneities in the spatial distribution of people infected or at risk of schistosomiasis, which our maps capture.
The uncertainties in the numbers infected with S. haematobium captured by our model were wide despite the large sample size and the geographical coverage of the data. It might be suggested that the large uncertainties shown in our prediction maps limit their utility for decision-making. However, since tools for representing uncertainty in spatial predictions now exist and there is still a need for disease maps in planning control programmes, it is beneficial to acknowledge such uncertainties when interpreting maps for disease control. Moreover, we did not include uncertainties in the projected population or those due to population migration between areas, which are considerable in these three countries. In addition, GRUMP can produce national population totals that do not match UN totals. In 2000, GRUMP population estimates for both Mali and the Niger differed from UN estimates by more than 5%,28 contributing to additional uncertainty in our updated 2005 population estimates. We also assumed that population growth was even across different regions and age groups. The limitations of the ancillary data used in GRUMP calculations have been noted elsewhere,29 but currently it is impossible to quantify the resulting inaccuracies.
Our surveys were spatially stratified, which ensured equal geographical coverage of areas of both high and low population density and made the precision of predicted prevalence estimates more even across the study area. However, less densely populated areas were more strongly represented than they would have been with a nonstratified approach. Because prevalence is likely to differ between low- and high-density areas, sample weighting is necessary to ensure accurate national prevalence estimates; our map calculations are essentially a sample-weighting approach. This is evident from the difference between raw prevalence estimates for low- and high-intensity infections combined, which were 18.5% for Burkina Faso, 36.6% for Mali and 13.5% for the Niger, and map-derived estimates, which were 23.0%, 35.4% and 21.1% in the three countries, respectively. The most striking difference was for the Niger. This arose because the prevalence was highest in the Niger River valley, which has the highest population density, and in the survey this area was underrepresented, from a population perspective, relative to low-density areas.
The different spatial distributions for low- and high-intensity infection are in agreement with Guyatt et al.,16 who demonstrated that the relationship between overall prevalence and the prevalence of high-intensity infection varied geographically across Africa. While low-intensity infection was more widespread, it exhibited less spatial variability than high-intensity infection, which aggregated more into clusters.16 Future control programmes will have the greatest impact if they focus on high-intensity infection clusters.
The maps presented in this report are currently being used by national programme managers as objective decision-support tools for geographically targeting existing resources more efficiently to high-risk communities. They have several other uses. First, national programme managers can use them to argue for resources from governments or international donors, particularly after funding from the Schistosomiasis Control Initiative ends. Second, the maps can be used to formulate and compare different disease control strategies by determining the likely impact on transmission and morbidity. And lastly, they can be used for advocacy and empowerment at the subnational level: local resource needs and priorities for schistosomiasis control, which are often subsumed in aggregate national data, can be presented to national programme managers and government officials. We encourage national programme managers in other countries, and those focussing on other diseases, to conduct spatially stratified disease surveys and undertake mapping of subnational disease distributions to provide evidence for more efficient targeting of resources. ▪
We thank the children, parents, teachers and head teachers who participated in the surveys. We also thank the technicians and support staff who undertook the surveys and the members of national schistosomiasis control programmes who provided administrative and organizational assistance.
Funding: During the course of this study, the Schistosomiasis Control Initiative was supported by the Bill and Melinda Gates Foundation. Simon Brooker is supported by a Career Development Fellowship (081673) from the Wellcome Trust. Archie Clements receives salary support and is a visiting scientist at the Australian Centre for International and Tropical Health, Queensland Institute of Medical Research, Herston, Qld, Australia.
Competing interests: None declared.
1. Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis 2006;6:411-25. PMID:16790382 doi:10.1016/S1473-3099(06)70521-7
2. Utroska JA, Chen MG, Dixon H, Yoon S, Helling-Borda M, Hogerzeil HV, et al. An estimate of global needs for praziquantel within schistosomiasis control programmes. Geneva: World Health Organization; 1989.
3. Brooker S, Clements AC, Bundy DA. Global epidemiology, ecology and control of soil-transmitted helminth infections. Adv Parasitol 2006;62:221-61. PMID:16647972 doi:10.1016/S0065-308X(05)62007-6
4. Guerra CA, Snow RW, Hay SI. Mapping the global extent of malaria in 2005. Trends Parasitol 2006;22:353-8. PMID:16798089 doi:10.1016/j.pt.2006.06.006
5. Hay SI, Guerra CA, Tatem AJ, Noor AM, Snow RW. The global distribution and population at risk of malaria: past, present, and future. Lancet Infect Dis 2004;4:327-36. PMID:15172341 doi:10.1016/S1473-3099(04)01043-6
6. Clements AC, Lwambo NJ, Blair L, Nyandindi U, Kaatano G, Kinung'hi S, et al. Bayesian spatial analysis and disease mapping: tools to enhance planning and implementation of a schistosomiasis control programme in Tanzania. Trop Med Int Health 2006;11:490-503. PMID:16553932 doi:10.1111/j.1365-3156.2006.01594.x
7. Gemperli A, Sogoba N, Fondjo E, Mabaso M, Bagayoko M, Briet OJ, et al. Mapping malaria transmission in West and Central Africa. Trop Med Int Health 2006;11:1032-46. PMID:16827704 doi:10.1111/j.1365-3156.2006.01640.x
8. Gemperli A, Vounatsou P, Kleinschmidt I, Bagayoko M, Lengeler C, Smith T. Spatial patterns of infant mortality in Mali: the effect of malaria endemicity. Am J Epidemiol 2004;159:64-72. PMID:14693661 doi:10.1093/aje/kwh001
9. 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 Cote d'Ivoire. Parasitology 2005;131:97-108. PMID:16038401 doi:10.1017/S0031182005007432
10. Raso G, Vounatsou P, Gosoniu L, Tanner M, N'Goran EK, Utzinger J. Risk factors and spatial patterns of hookworm infection among schoolchildren in a rural area of western Cote d'Ivoire. Int J Parasitol 2006;36:201-10. PMID:16259987 doi:10.1016/j.ijpara.2005.09.003
11. Raso G, Vounatsou P, Singer BH, N'Goran EK, Tanner M, Utzinger J. An integrated approach for risk profiling and spatial prediction of Schistosoma mansoni-hookworm coinfection. Proc Natl Acad Sci USA 2006;103:6934-9. PMID:16632601 doi:10.1073/pnas.0601559103
12. Diggle P, Moyeed R, Rowlingson B, Thompson M. Childhood malaria in the Gambia: a case-study in model-based geostatistics. Appl Stat 2002; 51:493-506.
13. Vester U, Kardorff R, Traore M, Traore HA, Fongoro S, Juchem C, et al. Urinary tract morbidity due to Schistosoma haematobium infection in Mali. Kidney Int 1997;52:478-81. PMID:9264005 doi:10.1038/ki.1997.356
14. Abdel-Wahab MF, Esmat G, Ramzy I, Fouad R, Abdel-Rahman M, Yosery A, et al. Schistosoma haematobium infection in Egyptian schoolchildren: demonstration of both hepatic and urinary tract morbidity by ultrasonography. Trans R Soc Trop Med Hyg 1992;86:406-9. PMID:1440818 doi:10.1016/0035-9203(92)90241-4
15. Warren KS, Mahmoud AA, Muruka JF, Whittaker LR, Ouma JH, Arap Siongok TK. Schistosomiasis haematobia in coast province Kenya. Relationship between egg output and morbidity. Am J Trop Med Hyg 1979;28:864-70. PMID:484768
16. Guyatt HL, Smith T, Gryseels B, Lengeler C, Mshinda H, Siziya S, et al. Aggregation in schistosomiasis: comparison of the relationships between prevalence and intensity in different endemic areas. Parasitology 1994;109:45-55. PMID:8058368 doi:10.1017/S0031182000077751
17. Alexander N, Moyeed R, Stander J. Spatial modelling of individual-level parasite counts using the negative binomial distribution. Biostatistics 2000;1:453-63. PMID:12933567 doi:10.1093/biostatistics/1.4.453
18. Alexander ND, Moyeed RA, Hyun PJ, Dimber ZB, Bockarie MJ, Stander J, et al. Spatial variation of Anopheles-transmitted Wuchereria bancrofti and Plasmodium falciparum infection densities in Papua New Guinea. Filaria J 2003;2:14. PMID:14525619 doi:10.1186/1475-2883-2-14
19. Brooker S, Alexander N, Geiger S, Moyeed RA, Stander J, Fleming F, et al. Contrasting patterns in the small-scale heterogeneity of human helminth infections in urban and rural environments in Brazil. Int J Parasitol 2006;36:1143-51. PMID:16814294 doi:10.1016/j.ijpara.2006.05.009
20. Clements AC, Moyeed R, Brooker S. Bayesian geostatistical prediction of the intensity of infection with Schistosoma mansoni in East Africa. Parasitology 2006;133:711-9. PMID:16953953 doi:10.1017/S0031182006001181
21. Garba A, Toure S, Dembele R, Bosque-Oliva E, Fenwick A. Implementation of national schistosomiasis control programmes in West Africa. Trends Parasitol 2006;22:322-6. PMID:16690357 doi:10.1016/j.pt.2006.04.007
22. Gryseels B. Mission d'Evaluation du report, Lutte contre la schistosomiasis au Mali. Eschborn: Gesellschaft für Technische Zusammenarbeit; 1989.
23. Diggle P, Tawn J, Moyeed R. Model-based geostatistics. Appl Stat 1998; 47:299-350.
24. Coordinated use of anthelminthic drugs in control interventions: a manual for health professionals and programme managers. Geneva: World Health Organization; 2006.
25. Fulford AJ, Webster M, Ouma JH, Kimani G, Dunne DW. Puberty and age-related changes in susceptibility to schistosome infection. Parasitol Today 1998;14:23-6. PMID:17040685 doi:10.1016/S0169-4758(97)01168-X
26. Feldmeier H, Poggensee G, Krantz I. Puberty and age-intensity profiles in schistosome infections: another hypothesis. Parasitol Today 1998;14:435. PMID:17040837 doi:10.1016/S0169-4758(98)01313-1
27. Rollinson D, Stothard JR, Southgate VR. Interactions between intermediate snail hosts of the genus Bulinus and schistosomes of the Schistosoma haematobium group. Parasitology 2001;123 Suppl;S245-60. PMID:11769287 doi:10.1017/S0031182001008046
28. Global Rural-Urban Mapping Project (GRUMP) [Alpha version]. Palisades, NY: Socioeconomic Data and Applications Center, Columbia University. Available from: http://sedac.ciesin.columbia.edu/gpw [accessed on 19 June 2009] .
29. Balk DL, Deichmann U, Yetman G, Pozzi F, Hay SI, Nelson A. Determining global population distribution: methods, applications and data. Adv Parasitol 2006;62:119-56. PMID:16647969 doi:10.1016/S0065-308X(05)62004-0
30. United Nations. World population prospects: the 2006 revision. New York, NY: United Nations, Department of Economic and Social Affairs, Population Division; 2007. Available from: http://www.un.org/esa/population/publications/wpp2006/wpp2006.htm [accessed on 19 June 2009] .
31. Brooker S, Donnelly CA, Guyatt HL. Estimating the number of helminthic infections in the Republic of Cameroon from data on infection prevalence in schoolchildren. Bull World Health Organ 2000;78:1456-65. PMID:11196498Box 1.
(Submitted: 11 September 2008 - Revised version received: 28 January 2009 - Accepted: 4 February 2009 - Published online: 27 July 2009)