Fine scale analysis of malaria incidence in under-5: hierarchical Bayesian spatio-temporal modelling of routinely collected malaria data between 2012–2018 in Cameroon

The current study aims to provide a fine-scale spatiotemporal estimate of malaria incidence among Cameroonian under-5, and to determine its associated environmental factors, to set up preventive interventions that are adapted to each health district of Cameroon. Routine data on symptomatic malaria in children under-5 collected in health facilities, between 2012 and 2018 were used. The trend of malaria cases was assessed by the Mann–Kendall (M–K) test. A time series decomposition was applied to malaria incidence to extract the seasonal component. Malaria risk was estimated by the standardised incidence ratio (SIR) and smoothed by a hierarchical Bayesian spatiotemporal model. In total, 4,052,216 cases of malaria were diagnosed between 2012 and 2018. There was a gradual increase per year, from 369,178 in 2012 to 652,661 in 2018. After adjusting the data for completeness, the national incidence ranged from 489‰ in 2012 to 603‰ in 2018, with an upward trend (M–K test p-value < 0.001). At the regional level, an upward trend was observed in Adamaoua, Centre without Yaoundé, East, and South regions. There was a positive spatial autocorrelation of the number of malaria incident-cases per district per year as suggested by the Moran’s I test (statistic range between 0.11 and 0.53). The crude SIR showed a heterogeneous malaria risk with values ranging from 0.00 to 8.90, meaning that some health districts have a risk 8.9 times higher than the national annual level. The incidence and risk of malaria among under-5 in Cameroon are heterogeneous and vary significantly across health districts and seasons. It is crucial to adapt malaria prevention measures to the specificities of each health district, in order to reduce its burden in health districts where the trend is upward.