disaggregation: An R Package for Bayesian Spatial Disaggregation Modeling

Disaggregation modeling, or downscaling, has become an important discipline in epidemiology. Surveillance data, aggregated over large regions, is becoming more common, leading to an increasing demand for modeling frameworks that can deal with this data to understand spatial patterns. Disaggregation regression models use response data aggregated over large heterogeneous regions to make predictions at fine-scale over the region by using fine-scale covariates to inform the heterogeneity. This paper presents the R package disaggregation, which provides functionality to streamline the process of running a disaggregation model for fine-scale predictions

Encouraging healthier grocery purchases online: a randomised controlled trial and lessons learned

Online supermarket platforms present an opportunity for encouraging healthier consumer purchases. A parallel, double-blind randomised controlled trial tested whether promoting healthier products (e.g. lower fat and lower calorie) on the Sainsbury's online supermarket platform would increase purchases of those products. Participants were Nectar loyalty membership scheme cardholders who shopped online with Sainsbury's between 20th September and 10th October 2017. Intervention arm customers saw advertisement banners and recipe ingredient lists containing healthier versions of the products presented in control arm banners and ingredient lists. The primary outcome measure was purchases of healthier products. Additional outcome measures were banner clicks, purchases of standard products, overall purchases and energy (kcal) purchased. Sample sizes were small due to customers navigating the website differently than expected. The intervention encouraged purchases of some promoted healthier products (spaghetti [B = 2.10, p < 0.001], spaghetti sauce [B = 2.06, p < 0.001], spaghetti cheese [B = 2.45, p = 0.001], sour cream [B = 2.52, p < 0.001], fajita wraps [B = 2.10, p < 0.001], fajita cheese [B = 1.19, p < 0.001], bakery aisle products (B = 3.05, p = 0.003) and cola aisle products [B = 0.97, p < 0.002]) but not others (spaghetti mince, or products in the yogurt and ice cream aisles). There was little evidence of effects on banner clicks and energy purchased. Small sample sizes may affect the robustness of these findings. We discuss the benefits of collaborating to share expertise and implement a trial in a live commercial environment, alongside key learnings for future collaborative research in similar contexts

Mapping malaria by sharing spatial information between incidence and prevalence data sets

As malaria incidence decreases and more countries move towards elimination, maps of malaria risk in low-prevalence areas are increasingly needed. For low-burden areas, disaggregation regression models have been developed to estimate risk at high spatial resolution from routine surveillance reports aggregated by administrative unit polygons. However, in areas with both routine surveillance data and prevalence surveys, models that make use of the spatial information from prevalence point-surveys might make more accurate predictions. Using case studies in Indonesia, Senegal and Madagascar, we compare the out-of-sample mean absolute error for two methods for incorporating point-level, spatial information into disaggregation regression models. The first simply fits a binomial-likelihood, logit-link, Gaussian random field to prevalence point-surveys to create a new covariate. The second is a multi-likelihood model that is fitted jointly to prevalence point-surveys and polygon incidence data. We find that in most cases there is no difference in mean absolute error between models. In only one case, did the new models perform the best. More generally, our results demonstrate that combining these types of data has the potential to reduce absolute error in estimates of malaria incidence but that simpler baseline models should always be fitted as a benchmark

Mapping malaria seasonality in Madagascar using health facility data.

BACKGROUND: Many malaria-endemic areas experience seasonal fluctuations in case incidence as Anopheles mosquito and Plasmodium parasite life cycles respond to changing environmental conditions. Identifying location-specific seasonality characteristics is useful for planning interventions. While most existing maps of malaria seasonality use fixed thresholds of rainfall, temperature, and/or vegetation indices to identify suitable transmission months, we construct a statistical modelling framework for characterising the seasonal patterns derived directly from monthly health facility data. METHODS: With data from 2669 of the 3247 health facilities in Madagascar, a spatiotemporal regression model was used to estimate seasonal patterns across the island. In the absence of catchment population estimates or the ability to aggregate to the district level, this focused on the monthly proportions of total annual cases by health facility level. The model was informed by dynamic environmental covariates known to directly influence seasonal malaria trends. To identify operationally relevant characteristics such as the transmission start months and associated uncertainty measures, an algorithm was developed and applied to model realisations. A seasonality index was used to incorporate burden information from household prevalence surveys and summarise 'how seasonal' locations are relative to their surroundings. RESULTS: Positive associations were detected between monthly case proportions and temporally lagged covariates of rainfall and temperature suitability. Consistent with the existing literature, model estimates indicate that while most parts of Madagascar experience peaks in malaria transmission near March-April, the eastern coast experiences an earlier peak around February. Transmission was estimated to start in southeast districts before southwest districts, suggesting that indoor residual spraying should be completed in the same order. In regions where the data suggested conflicting seasonal signals or two transmission seasons, estimates of seasonal features had larger deviations and therefore less certainty. CONCLUSIONS: Monthly health facility data can be used to establish seasonal patterns in malaria burden and augment the information provided by household prevalence surveys. The proposed modelling framework allows for evidence-based and cohesive inferences on location-specific seasonal characteristics. As health surveillance systems continue to improve, it is hoped that more of such data will be available to improve our understanding and planning of intervention strategies

Spatiotemporal mapping of malaria prevalence in Madagascar using routine surveillance and health survey data.

Malaria transmission in Madagascar is highly heterogeneous, exhibiting spatial, seasonal and long-term trends. Previous efforts to map malaria risk in Madagascar used prevalence data from Malaria Indicator Surveys. These cross-sectional surveys, conducted during the high transmission season most recently in 2013 and 2016, provide nationally representative prevalence data but cover relatively short time frames. Conversely, monthly case data are collected at health facilities but suffer from biases, including incomplete reporting and low rates of treatment seeking. We combined survey and case data to make monthly maps of prevalence between 2013 and 2016. Health facility catchment populations were estimated to produce incidence rates from the case data. Smoothed incidence surfaces, environmental and socioeconomic covariates, and survey data informed a Bayesian prevalence model, in which a flexible incidence-to-prevalence relationship was learned. Modelled spatial trends were consistent over time, with highest prevalence in the coastal regions and low prevalence in the highlands and desert south. Prevalence was lowest in 2014 and peaked in 2015 and seasonality was widely observed, including in some lower transmission regions. These trends highlight the utility of monthly prevalence estimates over the four year period. By combining survey and case data using this two-step modelling approach, we were able to take advantage of the relative strengths of each metric while accounting for potential bias in the case data. Similar modelling approaches combining large datasets of different malaria metrics may be applicable across sub-Saharan Africa

Indirect effects of the COVID-19 pandemic on malaria intervention coverage, morbidity, and mortality in Africa: a geospatial modelling analysis.

BACKGROUND: Substantial progress has been made in reducing the burden of malaria in Africa since 2000, but those gains could be jeopardised if the COVID-19 pandemic affects the availability of key malaria control interventions. The aim of this study was to evaluate plausible effects on malaria incidence and mortality under different levels of disruption to malaria control. METHODS: Using an established set of spatiotemporal Bayesian geostatistical models, we generated geospatial estimates across malaria-endemic African countries of the clinical case incidence and mortality of malaria, incorporating an updated database of parasite rate surveys, insecticide-treated net (ITN) coverage, and effective treatment rates. We established a baseline estimate for the anticipated malaria burden in Africa in the absence of COVID-19-related disruptions, and repeated the analysis for nine hypothetical scenarios in which effective treatment with an antimalarial drug and distribution of ITNs (both through routine channels and mass campaigns) were reduced to varying extents. FINDINGS: We estimated 215·2 (95% uncertainty interval 143·7-311·6) million cases and 386·4 (307·8-497·8) thousand deaths across malaria-endemic African countries in 2020 in our baseline scenario of undisrupted intervention coverage. With greater reductions in access to effective antimalarial drug treatment, our model predicted increasing numbers of cases and deaths: 224·1 (148·7-326·8) million cases and 487·9 (385·3-634·6) thousand deaths with a 25% reduction in antimalarial drug coverage; 233·1 (153·7-342·5) million cases and 597·4 (468·0-784·4) thousand deaths with a 50% reduction; and 242·3 (158·7-358·8) million cases and 715·2 (556·4-947·9) thousand deaths with a 75% reduction. Halting planned 2020 ITN mass distribution campaigns and reducing routine ITN distributions by 25%-75% also increased malaria burden to a total of 230·5 (151·6-343·3) million cases and 411·7 (322·8-545·5) thousand deaths with a 25% reduction; 232·8 (152·3-345·9) million cases and 415·5 (324·3-549·4) thousand deaths with a 50% reduction; and 234·0 (152·9-348·4) million cases and 417·6 (325·5-553·1) thousand deaths with a 75% reduction. When ITN coverage and antimalarial drug coverage were synchronously reduced, malaria burden increased to 240·5 (156·5-358·2) million cases and 520·9 (404·1-691·9) thousand deaths with a 25% reduction; 251·0 (162·2-377·0) million cases and 640·2 (492·0-856·7) thousand deaths with a 50% reduction; and 261·6 (167·7-396·8) million cases and 768·6 (586·1-1038·7) thousand deaths with a 75% reduction. INTERPRETATION: Under pessimistic scenarios, COVID-19-related disruption to malaria control in Africa could almost double malaria mortality in 2020, and potentially lead to even greater increases in subsequent years. To avoid a reversal of two decades of progress against malaria, averting this public health disaster must remain an integrated priority alongside the response to COVID-19. FUNDING: Bill and Melinda Gates Foundation; Channel 7 Telethon Trust, Western Australia

Maps and metrics of insecticide-treated net access, use, and nets-per-capita in Africa from 2000-2020.

Insecticide-treated nets (ITNs) are one of the most widespread and impactful malaria interventions in Africa, yet a spatially-resolved time series of ITN coverage has never been published. Using data from multiple sources, we generate high-resolution maps of ITN access, use, and nets-per-capita annually from 2000 to 2020 across the 40 highest-burden African countries. Our findings support several existing hypotheses: that use is high among those with access, that nets are discarded more quickly than official policy presumes, and that effectively distributing nets grows more difficult as coverage increases. The primary driving factors behind these findings are most likely strong cultural and social messaging around the importance of net use, low physical net durability, and a mixture of inherent commodity distribution challenges and less-than-optimal net allocation policies, respectively. These results can inform both policy decisions and downstream malaria analyses

Global estimation of anti-malarial drug effectiveness for the treatment of uncomplicated Plasmodium falciparum malaria 1991-2019.

BACKGROUND: Anti-malarial drugs play a critical role in reducing malaria morbidity and mortality, but their role is mediated by their effectiveness. Effectiveness is defined as the probability that an anti-malarial drug will successfully treat an individual infected with malaria parasites under routine health care delivery system. Anti-malarial drug effectiveness (AmE) is influenced by drug resistance, drug quality, health system quality, and patient adherence to drug use; its influence on malaria burden varies through space and time. METHODS: This study uses data from 232 efficacy trials comprised of 86,776 infected individuals to estimate the artemisinin-based and non-artemisinin-based AmE for treating falciparum malaria between 1991 and 2019. Bayesian spatiotemporal models were fitted and used to predict effectiveness at the pixel-level (5 km × 5 km). The median and interquartile ranges (IQR) of AmE are presented for all malaria-endemic countries. RESULTS: The global effectiveness of artemisinin-based drugs was 67.4% (IQR: 33.3-75.8), 70.1% (43.6-76.0) and 71.8% (46.9-76.4) for the 1991-2000, 2006-2010, and 2016-2019 periods, respectively. Countries in central Africa, a few in South America, and in the Asian region faced the challenge of lower effectiveness of artemisinin-based anti-malarials. However, improvements were seen after 2016, leaving only a few hotspots in Southeast Asia where resistance to artemisinin and partner drugs is currently problematic and in the central Africa where socio-demographic challenges limit effectiveness. The use of artemisinin-based combination therapy (ACT) with a competent partner drug and having multiple ACT as first-line treatment choice sustained high levels of effectiveness. High levels of access to healthcare, human resource capacity, education, and proximity to cities were associated with increased effectiveness. Effectiveness of non-artemisinin-based drugs was much lower than that of artemisinin-based with no improvement over time: 52.3% (17.9-74.9) for 1991-2000 and 55.5% (27.1-73.4) for 2011-2015. Overall, AmE for artemisinin-based and non-artemisinin-based drugs were, respectively, 29.6 and 36% below clinical efficacy as measured in anti-malarial drug trials. CONCLUSIONS: This study provides evidence that health system performance, drug quality and patient adherence influence the effectiveness of anti-malarials used in treating uncomplicated falciparum malaria. These results provide guidance to countries' treatment practises and are critical inputs for malaria prevalence and incidence models used to estimate national level malaria burden