In many countries health system data remain too weak to accurately enumerate Plasmodium falciparum malaria cases. In response, cartographic approaches have been developed that link maps of infection prevalence with mathematical relationships to predict the incidence rate of clinical malaria. Microsimulation (or 'agent-based') models represent a powerful new paradigm for defining such relationships; however, differences in model structure and calibration data mean that no consensus yet exists on the optimal form for use in disease-burden estimation. Here we develop a Bayesian statistical procedure combining functional regression-based model emulation with Markov Chain Monte Carlo sampling to calibrate three selected microsimulation models against a purpose-built data set of age-structured prevalence and incidence counts. This allows the generation of ensemble forecasts of the prevalence-incidence relationship stratified by age, transmission seasonality, treatment level and exposure history, from which we predict accelerating returns on investments in large-scale intervention campaigns as transmission and prevalence are progressively reduced.
