Hierarchical spatial modeling of the presence of Chagas disease insect vectors in Argentina. A comparative approach

We modeled the spatial distribution of the most important Chagas disease vectors in Argentina, in order to obtain a predictive mapping method for the probability of presence of the vector species. We analyzed both the binary variable of presence-absence of Chagas disease and the vector species richness in Argentina, in combination with climatic and topographical covariates associated to the region of interest. We used several statistical techniques to produce distribution maps of presence–absence for the different insect species as well as species richness, using a hierarchical Bayesian framework within the context of multivariate geostatistical modeling. Our results show that the inclusion of covariates improves the quality of the fitted models, and that there is spatial interaction between neighboring cells/pixels, so mapping methods used in the past, which assumed spatial independence, are not adequate as they provide unreliable results.