Correlation between the prevalence of Buruli ulcer and habitat suitability of Belostomatidae (left) and Naucoridae (right) in the wet season in Akonolinga. Colour indicates use of a linear model (Black) or locally weighted scatterplot smoothing (Red), different ways of viewing the correlation. Buruli ulcer was absent from certain villages (grey dots) where habitat suitability for the insects is high. Because this can skew any correlation between habitat suitability and prevalence, we explored the effect of including (thin lines) or excluding (thick lines) these villages. We note that in either case Spearman’s rank correlation coefficient was significant. For Belostomatidae with Buruli ulcer absent villages p = 0.08, without BU absent villages p = 0.03. For Naucoridae with BU absent villages p = 0.04, without BU absent villages p = 0.02.