Bayesian Geostatistical Modeling to Assess Malaria Seasonality and Monthly Incidence Risk in Eswatini

Eswatini is on the brink of malaria elimination and had however, had to shift its target year to eliminate malaria on several occasions since 2015 as the country struggled to achieve its zero malaria goal. We conducted a Bayesian geostatistical modeling study using malaria case data. A Bayesian distributed lags model (DLM) was implemented to assess the effects of seasonality on cases. A second Bayesian model based on polynomial distributed lags was implemented on the dataset to improve understanding of the lag effect of environmental factors on cases. Results showed that malaria increased during the dry season with proportion 0.051 compared to the rainy season with proportion 0.047 while rainfall of the preceding month (Lag2) had negative effect on malaria as it decreased by proportion − 0.25 (BCI: − 0.46, − 0.05). Night temperatures of the preceding first and second month were significantly associated with increased malaria in the following proportions: at Lag1 0.53 (BCI: 0.23, 0.84) and at Lag2 0.26 (BCI: 0.01, 0.51). Seasonality was an important predictor of malaria with proportion 0.72 (BCI: 0.40, 0.98). High malaria rates were identified for the months of July to October, moderate rates in the months of November to February and low rates in the months of March to June. The maps produced support-targeted malaria control interventions. The Bayesian geostatistical models could be extended for short-term and long-term forecasting of malaria supporting-targeted response both in space and time for effective elimination.