Progress has been made in mapping and predicting the risk of schistosomiasis using Bayesian geostatistical inference. Applications primarily focused on risk profiling of prevalence rather than infection intensity, although the latter is particularly important for morbidity control. In this review, the underlying assumptions used in a study mapping Schistosoma mansoni infection intensity in East Africa are examined. We argue that the assumption of stationarity needs to be relaxed, and that the negative binomial assumption might result in misleading inference because of a high number of excess zeros (individuals without an infection). We developed a Bayesian geostatistical zero-inflated (ZI) regression model that assumes a non-stationary spatial process. Our model is validated with a high-quality georeferenced database from western Côte d'Ivoire, consisting of demographic, environmental, parasitological and socio-economic data. Nearly 40% of the 3818 participating schoolchildren were infected with S. mansoni, and the mean egg count among infected children was 162 eggs per gram of stool (EPG), ranging between 24 and 6768 EPG. Compared to a negative binomial and ZI Poisson and negative binomial models, the Bayesian non-stationary ZI negative binomial model showed a better fit to the data. We conclude that geostatistical ZI models produce more accurate maps of helminth infection intensity than the spatial negative binomial ones.

The objectives of this study were (1) to examine risk factors for Schistosoma mansoni infection among schoolchildren living in western Côte d'Ivoire, and (2) to carry forward spatial risk prediction and mapping at non-sampled locations. First, demographic and socio-economic data were obtained from 3818 children, aged 6–16 years, from 55 schools. Second, a single stool sample was examined from each child by the Kato-Katz technique to assess infection status of S. mansoni and its intensity. Third, remotely sensed environmental data were derived from satellite imagery and digitized ground maps. With these databases a comprehensive geographical information system was established. Bayesian variogram models were applied for spatial risk modelling and prediction. The infection prevalence of S. mansoni was 38·9%, ranging from 0% to 89·3% among schools. Results showed that age, sex, the richest wealth quintile, elevation and rainfall explained the geographical variation of the school prevalences of S. mansoni infection. The goodness of fit of different spatial models revealed that age, sex and socio-economic status had a stronger influence on infection prevalence than environmental covariates. The generated risk map can be used by decision-makers for the design and implementation of schistosomiasis control in this setting. If successfully validated elsewhere, this approach can guide control programmes quite generally.