Skip to main content


Agent-based models of malaria transmission: a systematic review



Much of the extensive research regarding transmission of malaria is underpinned by mathematical modelling. Compartmental models, which focus on interactions and transitions between population strata, have been a mainstay of such modelling for more than a century. However, modellers are increasingly adopting agent-based approaches, which model hosts, vectors and/or their interactions on an individual level. One reason for the increasing popularity of such models is their potential to provide enhanced realism by allowing system-level behaviours to emerge as a consequence of accumulated individual-level interactions, as occurs in real populations.


A systematic review of 90 articles published between 1998 and May 2018 was performed, characterizing agent-based models (ABMs) relevant to malaria transmission. The review provides an overview of approaches used to date, determines the advantages of these approaches, and proposes ideas for progressing the field.


The rationale for ABM use over other modelling approaches centres around three points: the need to accurately represent increased stochasticity in low-transmission settings; the benefits of high-resolution spatial simulations; and heterogeneities in drug and vaccine efficacies due to individual patient characteristics. The success of these approaches provides avenues for further exploration of agent-based techniques for modelling malaria transmission. Potential extensions include varying elimination strategies across spatial landscapes, extending the size of spatial models, incorporating human movement dynamics, and developing increasingly comprehensive parameter estimation and optimization techniques.


Collectively, the literature covers an extensive array of topics, including the full spectrum of transmission and intervention regimes. Bringing these elements together under a common framework may enhance knowledge of, and guide policies towards, malaria elimination. However, because of the diversity of available models, endorsing a standardized approach to ABM implementation may not be possible. Instead it is recommended that model frameworks be contextually appropriate and sufficiently described. One key recommendation is to develop enhanced parameter estimation and optimization techniques. Extensions of current techniques will provide the robust results required to enhance current elimination efforts.


Malaria, alongside HIV and tuberculosis, is considered one of the “big three” infectious diseases of humans. The global response to malaria transmission has been significant, with Plasmodium falciparum malaria eliminated from 79 countries from 1979 to 2010 [1]. Modelling suggests that 70% of the reduction in malaria cases in sub-Saharan Africa (SSA) between 2000 and 2015 was attributable to the implementation of intervention strategies [2]. Key interventions included insecticide-treated bed nets (ITNs), artemisinin-based combination therapy (ACT) and indoor residual spraying (IRS). Often, field data are used as evidence for the efficacy and cost-effectiveness of selected interventions; however, these methods can be resource intensive, or have prohibitive ethical barriers. In such situations, mathematical simulation is increasingly used to provide further insights.

Infectious disease modelling of malaria has existed for over a century [3], with the dominant paradigm being the Ross–Macdonald models used by the Global Malaria Eradication Programme (GMEP) from 1955 to 1969 [4, 5]. These are examples of compartmental transmission dynamic models, in which the simulated human population consists of groups of individuals in disease states such as “susceptible”, “exposed”, “infectious” and “recovered”. More recent compartmental models of malaria provide insights into risk-stratification of populations, multiple mosquito populations, and waning immunity; however, the majority of models still closely resemble the Ross–Macdonald configuration [6]. Comparing a variety of modelling approaches can provide robustness of results, and highlight areas for development of modelling techniques [7]. As such, comparing alternative model frameworks may accelerate learning about disease transmission and control.

One approach that has significant potential is the use of agent-based models (ABMs). There appears to be no universal definition of ABMs; this review includes any model that explicitly models individual actions and responses, and associates with each individual respective state variables and parameters. In the following review, all published ABMs are stochastic in nature.

As ABMs focus on the individual, they afford flexibility in modelling factors, such as spatial heterogeneity (e.g. host movement, heterogeneous implementation of interventions) and stochasticity (e.g. inter-patient variability in time of infection, time to recovery, and location of infection). Compartmental models of malaria transmission do exist that incorporate either stochasticity of individual infections [8] or spatial heterogeneity [9]. However, in areas of increasing spatial variation, compartmental models may face convergence issues, or provide no more insight than alternative model structures. In low-transmission environments, where patient variability is more pronounced, ABMs can better represent the stochasticity of disease progression and transmission than compartmental models, where (to some resolution) people are grouped together and treated as interchangeable, so individual agent behaviour cannot be determined. Models that accommodate patient individuality and spatial variation can help fill knowledge gaps [7] about transmission heterogeneities important in malaria elimination strategies.

The flexibility of agent-based approaches also allows models to be constructed to address practical questions relating to malaria control and elimination in specific local contexts [7, 10]. This is advantageous because identifying optimal local intervention strategies can provide a strong evidence-based framework for National Malaria Control Programmes (NMCPs). ABMs can be constructed to resemble such specific settings closely, due to their flexibility in altering model attributes to reflect local individual characteristics and geographical factors.

As more is learned about malaria transmission, the complexity of the questions asked increases, which in turn calls for more nuanced models. The role of mathematical models continues to grow as both technical expertise and computing power increase. With the increasing capacity for modelling to assist in malaria elimination programmes, a review of the published literature for ABMs of malaria transmission was performed. Analysis included characterization of the structure of existing models, the factors influencing malaria transmission modelled, and the methods of data use and output analysis. The approaches used were highlighted and ideas for advancing the field proposed.


Search strategy and selection criteria

A systematic literature review was performed, consistent with the PRISMA statement [11]. Database searches of OVID Medline, CINAHL Plus and OVID Embase were performed as outlined in Additional file 1.

The search strategy aimed to return publications referring to each of the following three concepts in their subject heading, keywords list, title or abstract:

  1. 1.

    Terms relating to malaria, malaria vaccines, anti-malarials (including individual medications), and relevant Plasmodium species

  2. 2.

    Terms relating to epidemiology, demography, infectious disease outbreaks or transmission, epidemics and key outbreak model parameters (e.g. basic reproduction number)

  3. 3.

    Terms relating to agent-based, individual-based or microsimulation models.

No limits were placed on publication type, language, location, dates or publication status. Further to database searching, additional articles were found through expert knowledge, and backward and forward citation searching (the latter using Google Scholar) of included articles.

The abstracts of all returned studies were assessed for suitability, with those that mentioned an agent-based, individual-based, or microsimulation model of malaria transmission selected for full-text review. Full-text articles were excluded if no agent-based model of hosts or vectors was described or used, or if no components of malaria transmission (such as human disease, vector biting, or interventions) were explicitly modelled. Models of vector life cycles or ecology without malaria-specific elements were excluded. Articles that described comparison or ensemble modelling of pre-existing models were considered separately, and were not “included studies” for the purposes of this review.

Data collection and analysis

Due to the anticipated variation within the literature, a data extraction form was not developed prospectively. Instead, for included publications, key characteristics of each model were identified, with common themes and unique properties of models characterized. After model themes were identified, a classification for key components of models was created (see Additional file 2). The properties of each model were collated as discussed in the “Results” section.


The search yielded 406 abstracts potentially meeting inclusion criteria; 137 were selected for full-text assessment, with 90 papers included in the final analysis. Search results are presented diagrammatically using the PRISMA template (available at: in Additional file 3. The table in Additional file 2 provides detailed information about each model.

Model construction

Development of models and model families

Figure 1 outlines the timeline of models in the included papers. 27 original model frameworks have been published, with 15 models reported once to date, and 12 ‘baseline’ models collectively resulting in 63 further publications. Of the 27 original ABMs, six models were developed by extending compartmental or within-host models to an individualized framework [12,13,14,15,16,17]. Model progression was largely similar across the largest model frameworks, including OpenMalaria [18], EMOD [19] HYDREMATS [20], and work from Imperial College [14], such that the topic of each paper progressed logically from previous works. Extensions included simulating interventions [21], spatial mapping [22], or embedding “sub-modules” of key disease concepts into pre-existing works, such as host infectivity [23], mortality [24], and potential impacts of climate change [25]. When interventions were adequately simulated, cost-effectiveness analysis (CEA) was sometimes performed [26,27,28].

Fig. 1

Overview of included studies. The most recent studies are towards the bottom of the diagram. Solid arrows indicate papers directly linked to one another. Dotted arrows indicate models by the same lead author that are not related in methodology. Black solid or dotted boxes indicate models that did not meet inclusion criteria or are not agent-based, respectively, that guided the creation of a later ABM. Papers named in bold are spatially explicit. Boxes coloured orange, green and purple include human agents, mosquito agents or both, respectively

Model structure and frameworks

Methods of modelling at the individual level are influenced by factors including the question at hand, characteristics of the agent, and interactions of interest. Certain ABM frameworks naturally arise from the above concepts; in malaria modelling, key considerations include the choice of agent and whether to focus on disease states or transmission. While these choices are not mutually exclusive, three broad methods of individual simulation were common in the literature, each with differing degrees of agent autonomy. First, models developed from a compartmental structure typically used probabilities in place of flow rates to determine whether an individual transitioned to a new disease state at a given time step (e.g. Gurarie and McKenzie [13], and McKenzie et al. [29]). Using this method, the success or failure of a Bernoulli trial generally dictated disease progression. Second, and particularly in models focused on host parasite densities (see the OpenMalaria simulations [18]), distributions of key variables were sampled to generate differences between agents, with temporal disease state changes governed by a set of equations. In the third method, the specific actions of individuals, for example blood meal searching, resting and oviposition, were simulated according to a process represented by a flow chart (e.g. Pizzitutti et al. [30], Zhu et al. [31] or the EMOD framework [19]). Figure 2 presents a hypothetical flow chart of vector actions, to illustrate a typical set of transitions available to a mosquito.

Fig. 2

A hypothetical ‘decision tree’ approach to modelling mosquito agents. At each time step t, mosquitoes will check their status, determining their subsequent action(s) i with probability pi of success. Individuals may have identical or differing probabilities of success for each task

The diversity of methods led to a broad range of agent constructions, affecting the number and type of variables, the variation between individuals, and the type of agents required. Because methods to construct ABMs are so broad, two examples are provided below for comparison. These two models are examples of ABMs being tailored to the question at hand and highlights the variation in model construction.

Consider an early OpenMalaria simulation of the relationship between the entomological inoculation rate (EIR) and force of infection [18]. Human agents had individual ages, leading to an age-adjusted EIR. Individual mosquitoes were not explicitly modelled; instead, EIR was used to guide biting rates [32]. Two individuals of the same age may have a different number of infections on a particular day, but these numbers are drawn from the same distribution [33]. As OpenMalaria grew to answer new questions, these individuals also had different gametocyte densities [34], pyrogenic thresholds [35], and morbidity and mortality [24].

In contrast, Zhu et al. [21, 31, 36] developed a spatial ABM to investigate the impact of the location of food, hosts, and resting and breeding sites on malaria transmission. This required simulation of a physical landscape, with both mosquito and human locations guided by movement patterns. Mosquito agents had a sequence of actions to be undertaken, and therefore required variables capturing their current hunger, breeding state, and flight distance since last meal. The number of bites for each human was also recorded. The distance to the target of interest (sugar, host, resting site, breeding site) and time of day could dictate vector flight paths. Capping the number of eggs per larval habitat required the current number of eggs to be recorded.

There were also considerations relating the method of agent simulation and size of the time step. Where compartmental models inspired an ABM, daily time steps were most common [15]. To simulate host infections using distributions for variables, a temporal resolution of 1–5 days was sufficient to represent changing disease characteristics [32]. However, if agent actions were explicitly modelled, hourly status updates were most common, with some models tracking vectors as frequently as each second [21, 31].

Agency and elements modelled

The complexity of malaria transmission prevents any one model simulating all transmission factors in depth. In practice, each framework is centred on a few core concepts; these are generally sources of heterogeneity that motivate ABM use. These concepts can be categorized as host, vector, parasite, environment, and intervention factors. The core components of a model appeared to dictate the characteristics of each model framework and this relationship is discussed in the following sections. Figure 3 outlines some of the factors that vary between transmission scenarios, and therefore may be suited to modelling with an ABM.

Fig. 3

Diagram outlining factors influencing malaria transmission that have been modelled by ABMs. Factors pertaining to humans and mosquitoes are in red and blue circles, respectively. Factors about the disease process are within the arrows linking these circles. Factors at the top and bottom of the diagram are environmental factors and interventions, respectively

Model agents were either humans only (49/90 models), mosquitoes only (9/90) or both humans and mosquitoes (32/90) (Additional file 2, column 4). In a number of instances, the agency of vectors was unclear, although had been used in precursor studies. Pizzitutti et al. considered a Plasmodium agent as the infection of a host or vector [30, 37], whereas other models considered these elements as part of a vector or human agent. The AGiLESim framework specified larval habitats as agents, having properties such as location, larval capacity and current egg population [22]. The HYDREMATS simulations of the Bomblies et al. hydrology model [20, 25, 38,39,40,41,42] describe individual larval habitats with many variables and time-dependent formulae, but did not refer to habitats as agents.


Models with only humans as agents commonly focused on the impact of infection characteristics or medical interventions on transmission. For example, the OpenMalaria models from the Swiss Tropical and Public Health Institute had at their core individual humans with varying parasite densities [18]. These models were extended to incorporate individual effects of innate immunity [34], pyrogenic thresholds [35] and the impact of vaccination [43], particularly pre-erythrocytic vaccines [18, 33, 44]. The models from this group accounted for 20 of the 49 human-only studies.

Key simulated human factors included age and host immunity, with immunity either being maternal, [13, 14] acquired [18, 39, 45] or both [34]. Details about immunity varied from general descriptions of the impact of immunity on transmission, to individualized antibody concentrations for sexual-stage parasites [46]. Human behaviours regarding intervention use [30] and treatment decisions have been simulated for both patients [47] and carers [28]. One model simulated the impact of anti-malarial use on HIV-positive pregnant women [48], including disease severity and improvements in birth weight.

The most commonly simulated disease aspects were gametocyte density (incorporated in all OpenMalaria models [18]) and the infectiousness of hosts to mosquitoes [49], with other common elements including fever [35], disease morbidity and mortality [24] and disease severity [24, 50].

Of 81 papers with human agents, 50 stratified individuals by age (Additional file 2, column 13), either to vary disease profiles or report results. Commonly, age was used in the calculation of human biting rates, because of its strong correlation with body surface area [18]. Other age-varying factors were adaptive and maternal immunity and duration of infection [17]. If interventions were age-specific, for example vaccination programmes [33], stratification was used to assess the impact at different levels of the community [44].


Models with mosquitoes as agents most commonly assessed interventions impacting vector mortality, such as habitat removal [30], IRS [16] or ITN [51]. Detailed models of vector life cycle and ecology were common, with mosquito behaviours, such as feeding habits [31, 52], movement patterns [30, 53] and biting frequency [54], modelled in detail. Vector movement was often modelled as random motion, but CO2 gradients were also considered [20, 55,56,57]. Simulations such as the EMOD framework modelled the egg population as a cohort for their progression to adulthood [19], and two models simulated eggs as individual agents [58, 59].

In most publications, only female mosquitoes were simulated [60, 61], as the sex responsible for transmission. Male mosquitoes were generally modelled when a complete population was required, for example to assess vector control interventions (e.g. [21]). Males were specifically required to examine the effects of introducing gene drive mosquitoes into an environment, such as “driving-Y” modified male vectors [62]. Although dominated by An. gambiae representation, models also varied by Anopheles species, including An. vagus [63], An. stephensi [64], An. arabiensis [62], and An. darlingi [30, 37]. One An. gambiae model [65] was adapted to An. vagus vectors [63] to simulate transmission in Bangladesh, and several models incorporated multiple Anopheles populations [14, 56, 66,67,68], were adaptable across species [51, 58, 59, 69], or were nonspecific in Anopheles species [60, 70, 71].

When both hosts and vectors were simulated, models focused on mosquito life cycles [19], population dynamics [58], physical environments [55] and interventions for vector control [21]. Human agents were often included to relate Anopheles populations and malaria transmission, in models with vector dynamics as a core component. When models also included a spatial component, malaria transmission generally required vector and host to be co-located. These simulations used a ‘decision tree’ to represent the timing of movements and other necessary actions [21, 30, 31] (e.g. Fig. 2). Simulation of both hosts and vectors at the individual level was regularly used to assess interventions directed at mosquitoes and their effects on malaria transmission [19, 72].


Of 54 models that specified a malaria parasite (see Additional file 2, column five), 51 modelled the dynamics of Plasmodium falciparum, two simulated both P. falciparum and Plasmodium vivax [30, 37], and one Plasmodium berghei [64]. The in-host modelling of P. falciparum dynamics by Molineaux et al. [12, 18, 73] was used in 21 papers. In some studies, Plasmodium infections within an individual were modelled as agents [30, 37]; when a Plasmodium agent reached specific points in its life-cycle, the vector or host agent would change infectious status.

One study characterized the infectious reservoir in humans according to parasite characteristics [68], with parasite densities known to be linked to host infectiousness. Studies investigated the impacts of infection not only on hosts, but also vector behaviours, such as altered biting rates [56, 57]. Parasite biology was also considered, including parasite strains [17, 74], PfHRP2 status [75, 76] and recrudescence of P. vivax malaria [30]; the latter model considered the disease to be a submodule of the overall simulation framework. Three models simulated multiple parasite clones [15, 17, 77], and four allowed antigenic parasite variation [15, 78,79,80,81,82], particularly to capture antimalarial resistance.

The choice of parasite, much like the Anopheles species, was often based on the dominant species in the target location. The dominance of P. falciparum malaria simulation reflects the attention paid to it, which is largely due to the historic relative burden of disease. Additionally, a number of interventions simulated were specific to P. falciparum, in particular the pre-erythrocytic RTS,S vaccine [14, 44, 83].

Modelling different Plasmodium species may require changes to model structure, for example to account for the recrudescence seen in P. vivax but not P. falciparum malaria. Pizzitutti et al. [30, 37] incorporated these differences by adding parameters governing recurrence time and risk of P. falciparum-triggered P. vivax recurrence. The parameters representing infection vary between Plasmodium species, based on prior evidence. Beyond these changes, the decision-tree model structures dictating the actions of agents were not changed.

Environment and spatial modelling

Core environmental aspects related to vector activity were included, such as water sources for oviposition (egg-laying) [20], houses for blood meal locations [72], and meteorological data to account for seasonal patterns in transmission [58, 84]. Rainfall and temperature data were regularly used; model questions included the impact of wet season lengths [42] and hysteresis [40] on vector populations.

Thirty of the 90 models explicitly incorporated a spatial framework into their model (Fig. 1 and Additional file 2, column 11; explored further in Additional file 4). Seventeen models investigated malaria transmission in Africa, by simulating a specific location (e.g. [20, 85]) or using a hypothetical landscape representing a typical African village (e.g. [31, 53]). Two models represented transmission in the Amazon region [30, 37] and one each in Papua New Guinea [17], France [71], Haiti [86] and Bangladesh [63]. Spatial models also simulated local vector activity in a generic physical environment [51, 57, 58].

While many spatial models had detailed representations of a known landscape [25], this was not always the case, particularly with early spatial models [22, 86]. Three methods of spatial construction predominated: grid-based systems, patches and continuous landscapes. Additionally, Depinay et al. [58] explicitly linked vector locations to houses, but did not assign spatial coordinates, whilst Karl et al. [17] generated a probability of transmission based on the distance between host and vector locations to incorporate spatial biting factors.

In a grid-based landscape, the modelled area was most commonly divided into squares of a constant size, with elements of the landscape residing within a square (e.g. [22]). The HYDREMATS hydrology/entomology models were the exception, with grid size varying depending on proximity to human habitats, and complexity of water sources in the area [20]. There was interplay between the physical properties of the study area and the construction of the model. For example, the size of grid squares were related to parameters such as typical distance moved by mosquitoes in one time step, or larval habitat size [20]. The overall simulation area, typically a square between 600 m and 3000 m in length, was commonly justified based on Anopheles dispersal behaviour and typical village sizes being modelled [31]. Six models did simulate much larger areas [16, 17, 37, 62, 67, 85]; the last of these models used 1 km grids to describe the environment, yet characterized movement and interventions on much finer spatial scales [67]. Distances from which mosquitoes can target hosts or oviposition sites were guided by previous field data where possible [30, 31].

The additional spatial methods addressed opposing limitations of grid-based landscapes, namely total area and spatial precision. A patch framework was used by Silal et al. [16] and Eckhoff et al. [62] to represent areas of differing transmission characteristics (e.g. changing EIRs within districts), without explicitly modelling the environment. This method allowed for simulation of broader spatial areas. Rateb et al. [86] achieved this by representing Haiti as a set of “microenvironments” that host agents were impacted by. Conversely, models with a continuous landscape were used to allow variation in the location of houses, larval habitats, sugar sources and resting sites, and to model the impact of different physical environments [21, 31]. Each object had a specified location, and mosquitoes could sense objects within a certain circular distance, as opposed to a particular grid square. However, even when using a village-sized continuous landscape, models from Zhu et al. [21, 31, 36] still mapped random mosquito flight using a grid-based system, albeit on a 1 m × 1 m scale.

Some models also impute spatial results from non-spatial models by alternate means. The Malaria Atlas Project [87] and the Markham Seasonality Index [88] have generated methods to predict EIRs by location, and researchers have combined these tools with model outputs to generate maps of estimated disease burden across sub-Saharan Africa (SSA) [27, 54, 69, 78, 89, 90]. Models also used SSA rainfall data to replicate seasonality patterns [69, 91] and population estimates for host locations over large areas [91, 92]. In one instance, a spatial model of hypothetical locations was simulated in conjunction with geographic information systems (GIS) technology, to represent transmission in specific locations [61].


Fifty-eight studies assessed the impact of at least one intervention on malaria transmission, mosquito prevalence or EIR (Additional file 2, column 6). The majority of papers assessed multiple interventions, 21 assessed interventions in combination, and one investigated the removal of current malaria strategies [52]. Interventions could broadly be divided into those targeted at the human host (e.g. pharmacological) or the vector. Model structure tended to support simulation of the intervention being tested. For example, larval habitat elimination [51] or toxic sugar baits [21] were modelled spatially, requiring vectors as agents; whereas vaccination [18, 33, 44] and drug administration [93] studies involved human agents, in order to capture individual immune responses.

Interventions were included in models for two reasons: for assessment of effectiveness [18], or to replicate the ‘baseline’ characteristics of regional transmission [26], to allow additional interventions to be explored. In general, baseline scenarios included ITN [54] and “case management” [26], as defined by the current interventions in the location being modelled.

Simulation of interventions usually took one of three approaches. First, a specific intervention was defined to have a certain effect when in use, for example a medication having a fixed efficacy, and the impact on transmission is measured. Second, where environments were explicitly modelled, interventions were usually defined by their physical impact. For example, from a baseline scenario, an intervention might alter individual larval habitats [51] or remove sugar sources to reduce vector feeding [31], and simulations highlighted the change in the output of interest, such as disease burden or vector populations. Third, hypothetical interventions were described by their impacts, to target a specific aspect of malaria transmission [65]. The hypothetical impact may approximate a pre-existing intervention, such as halving human biting rates for a fixed period [94], or evaluate novel targets that may guide future control techniques.


The ease with which ABMs allow for heterogeneity and model complexity enables some less often considered transmission factors to be simulated. Examples include anthropophily (i.e. vector preference for human hosts) [19], sugar sources (e.g. fruit trees) [21, 31], travellers [16, 95], the impacts of larval densities in oviposition sites [59, 65], and the pharmacokinetics and pharmacodynamics of anti-malarials [75, 80, 91, 96].

Data use, model outputs and analysis

Parameter estimation and robustness of results

Most models used past literature or simulations to determine baseline parameter values (Additional file 2, column 8). Due to the nature of certain inputs, many parameters cannot or have not been estimated in field studies, and consequently authors used expert knowledge to select these values. Five papers provided no evidence that prior knowledge was used to estimate parameters. These models generally introduced a new modelling technique, or assessed hypothetical scenarios, as opposed to claiming accurate simulation of real-world situations.

Fifty-two papers used models that were either previously calibrated to data or presented calibration as a component of their work. Calibration techniques included the use of calibration vectors, least squares, maximum likelihood functions, and visual estimations. Twelve papers mentioned using Bayesian techniques for model fitting, with five providing credible intervals alongside point estimates of parameters [14, 23, 32, 50, 89, 91, 92, 97]. Of note, no papers calibrated all model parameters to data, which necessarily excludes certain parameter combinations that could produce accurate calibrations.

Validation was performed on models used in 31 papers. All major model frameworks were reported as validated. Validation techniques were rarely explained in detail. When described, validation was most commonly performed by running a calibrated simulation, and comparing model outputs to a dataset not used for calibration [98]. Methods of comparing models to data were rarely explained; use of a square distance function [37], log likelihoods [16] and docking techniques [99] were outlined. Successful model validation was often used to justify extending a model framework to include interventions or to assess their potential impact in the location of interest. A number of studies concluded that their model did not accurately fit the data used for validation [98], suggesting incomplete data may explain any discrepancies.

Forty-five papers either explicitly mentioned sensitivity analysis or described parameter variation and comparison of outputs. The techniques used were generally informal, with methods used rarely explained in detail, and reporting of results was uncommon. Where sensitivity analysis was explained, it involved altering calibrated baseline parameter values by a fixed percentage and assessing changes to outputs. Formal techniques employed included Latin Hypercube sampling [70, 72], regression tree analysis [72], and one-way and probabilistic sensitivity analysis [28].

Optimization and cost-effectiveness analysis

A number of papers referenced optimization use, including in titles [21, 45, 85] and keywords [100], and seven papers performed analysis for the clear purpose of optimization (Additional file 2, column 16). Methods included changing the timing [94, 95], location [21] or implementation of interventions [45, 85]. Five studies assessed between nine and 144 scenarios each. In contrast, two studies described formal optimization methods [69, 97], whereby one or more objective functions are optimized using a defined search algorithm. By changing the combination and coverage of interventions in a variety of baseline settings, these studies compared 98,784 and 306,000 scenarios, respectively.

CEA was performed in 11 studies, with six of these using the OpenMalaria platform to assess various medical interventions [101,102,103] or pre-erythrocytic vaccine RTS,S/SA02A [26, 27, 104] (Additional file 2, column 14). CEA studies used known local costs of the physical intervention and the costs of implementation [27, 104] to determine the financial impact of interventions. These values were compared to benchmarks, such as the World Health Organization (WHO) standard of $150 per life-year gained [28], to determine the potential value of interventions. No studies investigated the maximum impact on outputs for a fixed cost, instead assessing cost-effectiveness of interventions for a specific level of impact.

When determining how to scale interventions for CEA, one method employed was incremental cost-effectiveness ratios (ICERs) [27, 103], in which a baseline scenario was compared to the introduction of a number of interventions, to determine the net cost per case averted of each method up to a certain threshold. Interventions are added iteratively in a manner that does not perfectly optimize their delivery but has the lowest additional cost for each new intervention block.

Ensemble modelling and model comparison

Where multiple models analyse the same transmission scenario, ensemble modelling can be used to highlight robustness of results, or to assess the appropriateness of model structure. Similar to the role of sensitivity analysis in evaluating the sensitivity of model outputs to parameter values, ensemble modelling can assess aspects of both parametric and structural uncertainty. A precursor to this was performed by Ross et al., extending a baseline model to five variants to assess hypotheses regarding treatment and illness patterns [105]. In addition to the individual studies in this review, 12 publications included ensemble modelling of OpenMalaria model variants [106,107,108,109,110,111,112,113,114,115,116,117], and three studies performed “consensus modelling” across different frameworks [118,119,120] (see Additional file 5).

The OpenMalaria ensembles used from six to 14 model variants, selected from 30 tested alternatives. The 16 rejected models were either very similar to an accepted choice or could not be parameterized to the dataset used. For the three articles that compared outputs across the ensemble, the selected models were in agreement for most outputs [106, 107, 113]. The most common divergence was differing intervention efficacies predicted by models with in-built transmission heterogeneity, compared to those with uniform transmission levels between individuals. Interestingly, even though the fourteen models often reached the same conclusion, the estimated parameter values from calibration varied significantly across the ensemble [107].

The nine additional OpenMalaria ensembles aggregated the outputs from model variants, to investigate interventions such as vaccination [108, 116], seasonal malaria chemoprevention [115], mass test-and-treat strategies [109], and long-lasting insecticidal nets (LLINs) and IRS [110,111,112, 114]. One study predicted changes to disease outcomes and transmission upon the removal of vector control interventions from elimination programmes [117]. The aggregation of model outputs for each parameter set was performed to reduce the uncertainty arising from structural differences in models. In this way, the ensemble is considered as one model and its outputs interpreted as such.

Unlike the ensemble modelling above, when comparing inter-model variation, simulation outputs could not be aggregated. Despite the varied approaches, the consensus modelling largely drew consistent findings on the relationship between available malaria prevalence data and clinical incidence [118], and on the impacts of vaccination [119] and mass drug administration (MDA) [120].


Mathematical modelling plays an important role in malaria elimination, and agent-based approaches make a major contribution to these efforts. The extension of compartmental models to their early ABM equivalents arose from the need to understand malaria transmission at the individual level. The result is a rich array of model families and simulation techniques, adapted to a range of key issues in transmission and control.

In general, three core themes emerged regarding justification of ABM use. First, the greater importance of stochasticity in low-transmission settings, particularly settings approaching elimination, requires an alternative approach to traditional compartmental methods. Second, attempts to eliminate local transmission require discrete population simulations to incorporate spatially explicit environments at increasingly fine resolutions. Third, heterogeneities in disease progression and severity on the individual patient level result in varying efficacy of drug and vaccine interventions, which may be difficult to capture within a compartmental framework. These three arguments stem from a common point: a compartmental structure, based on averaging over a population, has limitations when that average does not adequately represent the individuals.

In addressing these issues, the benefits of agent-based techniques in this space are evident. Many papers in this review explicitly aimed to fill the knowledge gap regarding intervention use in low transmission environments. Most projects provided outputs robust at multiple transmission intensities, highlighting the flexibility of ABMs in low-prevalence settings. The HYDREMATS framework was used in multiple locations, incorporating environmental factors such as temperature and rainfall at different times, [25, 38, 39]. The OpenMalaria models progressed from assessing the force of infection of malaria transmission [32], to estimating cost-effectiveness of a vaccination programme [27]. Given the similarities of compartmental models of malaria to the original Ross–McDonald framework [6], it is clear that the depth and flexibility of agent-based methods are allowing new insights into malaria transmission and prevention.

The variation in the models described above highlights the difficulties in developing a standardized style of ABM for use in malaria epidemiology. However, this is arguably a major advantage, with the abundance of techniques allowing for the flexibility desired when transitioning from solely using compartmental models. Instead of suggesting a “gold standard” approach, it may be preferable to ensure the model style used is appropriate for the question at hand. For example, OpenMalaria’s early modelling of gametocyte densities did not use vector agents [18], but successfully provided insights into risks of fever, morbidity and mortality of patients [24]. The EMOD models initially described host-vector interactions without spatial consideration [19], but added this capability when required to assess interventions [62, 85]. Therefore, while not every model incorporated every aspect of malaria epidemiology, each was tailored to the research question at hand.

Conversely, if modelling groups are considering extending their model frameworks, particularly to influence policy, there is potential to draw from the features of one another. For example, HYDREMATS currently includes human and mosquito agents, while the characteristics of human infection are more detailed in the OpenMalaria simulations. Therefore, the time variability of individual gametocyte density, probabilities of fever, morbidity and mortality, and the infectivity of hosts to vectors used in the OpenMalaria framework could be adapted into HYDREMATS to more realistically replicate disease transmission. However, in neither of these simulations do humans move, whereas this process is explicitly simulated in Zhu et al. [31] and Pizzitutti et al. [37] to better represent vector-host feeding patterns. Pizzitutti et al. [30] and EMOD [19] include human behavioural reactions to biting rates (i.e. time-dependent intervention use) and the probabilities of successful blood meals, respectively. These five models have components that simulate vector, egg and human populations, effects of climate on larval habitats, anthropophily, ITN, IRS, larval habitat removal, vaccination, anti-malarial use, attractive toxic sugar baits, and rates of human disease. As each framework provides insights into key components of malaria transmission, all of which are important in guiding elimination strategies.

To some extent, combining model structures across research teams can be considered an extension of the use of “submodules” already undertaken by larger modelling frameworks. The HYDREMATS team have successfully integrated detailed larval habitat and entomologic models [55], and OpenMalaria now includes upwards of seven modules of human disease states and interventions [100]. Modular projects such as HYDREMATS, EMOD [85] and OpenMalaria have provided insights into transmission dynamics, vector populations, disease severity, and the contributors to these factors. Given the importance of comprehensive modelling to guide policy decisions, the potential for combining the strengths of validated models to enhance decision-making capabilities of ABMs could be explored.

A key target for modelling low-transmission settings is a focus on spatial representation and heterogeneity. ABMs can shift spatial heterogeneity from a typically “patch-based” compartmental framework into a continuous space, by having explicit locations for environmental objects, dwellings and agents. These detailed descriptions of the landscape are coupled with local knowledge of physical characteristics (such as host/vector movement patterns) to simulate malaria transmission, ecology, and the impact of interventions based on their location. These insights include the distances between larval habitats and houses to effectively reduce malaria transmission [53], and the impact on systematic versus random location of attractive toxic sugar baits (ATSB) on mosquito abundance [21]. This style of intervention inherently requires spatial modelling, although interventions such as ITNs and IRS have been modelled in both spatial and non-spatial simulations. Examining these interventions in physical environments allows the impact of factors such as vector movement and the proximity of unprotected individuals to be measured.

Potential extensions to spatial models include varying elimination strategies across a landscape, and increasing the size of the geographical area modelled. For example, consider a small community near a local water source, and a nearby larger population with better access to healthcare. A model could implement larval source management at the local water sources, whilst increasing access to vaccines and anti-malarials in the healthcare centre. Human movement dynamics [121] could be incorporated to assess the relative effectiveness of each interventions across both populations. This style of modelling may more accurately represent the manner in which interventions are implemented at the local level. Regarding model areas, most simulations of real-world environments only covered the area of a specific village, with sizes ranging from 600 m × 600 m to 3000 m × 3000 m. These spatial ranges have been limited by computational power, but this limitation will continue to decrease over time. There also may be a lack of access to consistent geospatial data over larger areas; however, many models only included spatial data on physical habitats, which would be collected in a similar manner over larger areas. If it is deemed useful to model malaria over a wider area, techniques from other fields may be used, such as probability modelling of invasive species, which has been performed for an area of over 35,000 km2 [122].

Regarding the locations of malaria modelling, there is an understandable focus on SSA, which was responsible for 88% of the global malaria burden in 2015 [2]. However, there has also been a recent increase in attention on South-East Asia (SEA), which is responsible for 10% of global cases and has emerging issues with drug resistance [2]. Moreover, challenges such as artemisinin resistance and insecticide resistance are more prominent in these areas. Despite this, approximately half of all malaria cases outside Africa in 2015 were due to P. vivax [2], while Plasmodium knowlesi malaria transmission is increasing in locations, such as Malaysia [123]. The methods of ABM construction used in SSA and SEA, and for P. falciparum malaria, suggest transferability to other regions and Plasmodium species, which will be important as data availability from these areas improves and attention turns to global elimination.

Many ABMs reviewed here had an emphasis on informing policy and explicitly aimed to understand specific programmatic questions (e.g. [16, 28, 78, 83]). To reliably inform public health decisions, there must be confidence in the assumptions guiding model creation, in particular regarding choice of parameter values. The methods for estimating key parameters varied greatly across the literature. Parameter justification was not always clear [13, 22, 51, 65]; when explained, models generally calibrated a range of parameters to existing data, or provided references for their choice of fixed values. Further, no previous models have calibrated all parameters to data. Importantly, ensemble modelling of the OpenMalaria variants calibrated 14 variants to the same dataset, but parameter values varied significantly across between models [107]. Given that uncertainty remains even after calibration, it is important to apply a systematic and comprehensive approach to parameter estimation before using models for predicting parameter impact or forecasting.

Whilst field data exists for a range of parameters, researchers must be pragmatic about the possibility of adequately calibrating complex ABMs, particularly when data is required in resource-poor settings. The fixing of well-established values can reduce the parameter space to be searched using calibration methods. Alternatively, techniques such as Markov chain Monte Carlo (MCMC) can search the entire parameter space, or more precisely that part of parameter space that has non-negligible posterior probability. MCMC has already been used in malaria ensemble modelling [108, 116, 118] and parameter estimation by Griffin et al. [14], as well as in modelling of other infectious diseases [124, 125]. Approaches such as MCMC and approximate Bayesian computation are increasing in popularity as including uncertainty in model parameters becomes more common [126].

An arguably more pressing area of need for development is optimization, with methods for agent-based models still in their infancy. In this review, studies that reported optimization usually simulated a suite of different interventions, or the same interventions at different levels of coverage or timing. Cost-effectiveness analysis was typically approached in the same manner. Whilst conclusions were provided as to the most effective simulation approach, true optimization was rare [69, 97], using formal techniques to identify parameter values that optimize one or more objective functions. Given the role of ABMs in modelling interventions in low-transmission settings, formal optimization techniques are important for enhancing the ability of models to guide policy.

Deterministic models have already been used as the basis of optimization of interventions for various infection diseases [125, 127, 128]. Strategies for the optimization of interventions within ABMs appear less common, possibly due to the high computational burden of finding consistent minima in the presence of stochasticity. It is difficult to define how to best approach optimization from an agent-based standpoint. A systematic review of ABMs for optimization problems [129] highlighted techniques used for disciplines such as scheduling, supply chain management, energy systems planning and transportation and logistics. As is likely the case regarding spatial methods, optimization of malaria transmission modelling (and infectious disease simulation more broadly) may benefit from adapting approaches outside the field to a new context.

Increased clarity in model reporting would be of great benefit to both the creators of ABMs and their audience. While many papers included detailed supplementary materials for additional results, project descriptions and calibration, validation, sensitivity analysis and optimization techniques, the intricacies of these techniques often unclear. A protocol exists for the description of ABMs [130], and models that used it [21, 31, 76] were simple to understand and appeared easily replicable by external groups. Further transparency includes sharing of the mathematics and code [76] of models. These small steps in documentation would allow for increased verification and validation of models, as well as increasing opportunities for collaboration between modelling groups.

Beyond individual models, ensemble modelling is an important tool for generating robust conclusions about malaria transmission. A review of ebola models advocated for an ensemble modelling approach that adequately compares state-of-the-art models, but also allows for model diversity [131]. For malaria models built for similar purposes (for example, to estimate certain parameters of interest, or to predict the value of an intervention), both inter- and intra-model comparison has been conducted. In some cases, two or three interventions have been simultaneously assessed using this approach [110, 114, 115]. The next steps for ensemble and consensus modelling may include more interventions, spatial modelling, and developing techniques to determine which models are most appropriate when ensemble members differ in parameter estimates or outputs.


As malaria transmission continues to decline and interventions become more nuanced, agent-based modelling will have an increasingly important role to play in elimination programmes. A variety of techniques have already been developed, and models are increasingly tailored to the question of interest. The flexibility of ABMs is a key feature, with progressive model extensions and fine-grain spatiotemporal simulation two clear examples.

The breadth of model frameworks makes it difficult to develop guidelines for ABM construction, but this should be considered a strength. The correct approach to agent-based modelling is likely a “horses for courses” approach, whereby the question at hand guides the development of the model; current approaches may be used to guide such choices. As models transition from tools for analysing an epidemic to guiding policy directions, it is important to be aware of the current literature and techniques available. Existing model frameworks cover many transmission factors, so cross-collaboration may bring about larger models that can simulate many interventions and provide outputs regarding vector populations, host disease progression and the success of elimination strategies. However, using larger models may create difficulties with simulation. Improved techniques in parameter estimation and optimization could enhance the role ensemble modelling currently plays in evaluating interventions for specific geographical and transmission contexts.

With the ever-increasing computing power available to researchers, detailed ABMs that accurately reflect the biology of malaria transmission are increasingly feasible on a fine spatial resolution over large geographical regions. As such, agent-based modelling will be an important tool for helping to inform malaria elimination strategies over the coming years.



agent-based model


sub-Saharan Africa


insecticide-treated bed nets


artemisinin-based combination therapy


indoor residual spraying


Global Malaria Eradication Programme


National Malaria Control Programme


cost-effectiveness analysis


entomological inoculation rate


World Health Organization


mass drug administration


incremental cost-effectiveness ratio


long-lasting insecticidal net


South-East Asia


Markov chain Monte Carlo


  1. 1.

    Feachem RG, Phillips AA, Targett GA, Snow RW. Call to action: priorities for malaria elimination. Lancet. 2010;376:1517–21.

  2. 2.

    WHO. World Malaria Report 2015. Geneva: World Health Organization; 2015. Accessed 4 Oct 2017.

  3. 3.

    Ross R. Report on the prevention of malaria in Mauritius. London: Waterlow; 1908.

  4. 4.

    Macdonald G. Epidemiological basis of malaria control. Bull World Health Organ. 1956;15:613–26.

  5. 5.

    Macdonald G, Cuellar CB, Foll CV. The dynamics of malaria. Bull World Health Organ. 1968;38:743–55.

  6. 6.

    Reiner RC, Perkins AT, Barker CM, Niu T, Chaves LF, Ellis AM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J R Soc Interface. 2013;10:20120921.

  7. 7.

    Chitnis N, Schapira A, Smith DL, Smith T, Hay SI, Steketee R. Mathematical modelling to support malaria control and elimination. Progress & impact series, vol 5. Geneva: Roll Back Malaria; 2010.

  8. 8.

    Maude RJ, Pontavornpinyo W, Saralamba S, Aguas R, Yeung S, Dondorp AM, et al. The last man standing is the most resistant: eliminating artemisinin-resistant malaria in Cambodia. Malar J. 2009;8:31.

  9. 9.

    Acevedo MA, Prosper O, Lopiano K, Ruktanonchai N, Caughlin TT, Martcheva M, et al. Spatial heterogeneity, host movement and mosquito-borne disease transmission. PLoS ONE. 2015;10:e0127552.

  10. 10.

    Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9:e01001165.

  11. 11.

    Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche PC, Ioannidis JPA, et al. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate health care interventions: explanation and elaboration. PLoS Med. 2009;6:e1000100.

  12. 12.

    Dietz K, Raddatz G, Molineaux L. Mathematical model of the first wave of Plasmodium falciparum asexual parasitemia in non-immune and vaccinated individuals. Am J Trop Med Hyg. 2006;75(2 suppl):46–55.

  13. 13.

    Gurarie D, McKenzie FE. A stochastic model of immune-modulated malaria infection and disease in children. Math Biosci. 2007;210:576–97.

  14. 14.

    Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, et al. Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010;7:e1000324.

  15. 15.

    Gatton ML, Cheng Q. Interrupting malaria transmission: quantifying the impact of interventions in regions of low to moderate transmission. PLoS ONE. 2010;5:e15149.

  16. 16.

    Silal SP, Little F, Barnes KI, White LJ. Predicting the impact of border control on malaria transmission: a simulated focal screen and treat campaign. Malar J. 2015;12(14):268.

  17. 17.

    Karl S, White MT, Milne GJ, Gurarie D, Hay SI, Barry AE, et al. Spatial effects on the multiplicity of Plasmodium falciparum infections. PLoS ONE. 2016;11:e0164054.

  18. 18.

    Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, et al. Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: overview. Am J Trop Med Hyg. 2006;75(2 suppl):1–10.

  19. 19.

    Eckhoff PA. A malaria transmission-directed model of mosquito life cycle and ecology. Malar J. 2011;10:303.

  20. 20.

    Bomblies A, Duchemin JB, Eltahir EAB. Hydrology of malaria: model development and application to a Sahelian village. Water Resour Res. 2008;44:1–26.

  21. 21.

    Zhu L, Marshall JM, Qualls WA, Schlein Y, McManus JW, Arheart KL, et al. Modelling optimum use of attractive toxic sugar bait stations for effective malaria vector control in Africa. Malar J. 2015;14:492.

  22. 22.

    Arifin SMN, Davis GJ, Zhou Y. A spatial agent-based model of malaria: model verification and effects of spatial heterogeneity. Int J Agent Technol Syst. 2011;3:17–34.

  23. 23.

    Ross A, Killeen G, Smith T. Relationships between host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg. 2006;75(2 suppl):32–7.

  24. 24.

    Ross A, Maire N, Molineaux L, Smith T. An epidemiologic model of severe morbidity and mortality caused by Plasmodium falciparum. Am J Trop Med Hyg. 2006;75(2 suppl):63–73.

  25. 25.

    Bomblies A, Eltahir EAB. Assessment of the impact of climate shifts on malaria transmission in the Sahel. EcoHealth. 2009;6:426–37.

  26. 26.

    Tediosi F, Hutton G, Maire N, Smith TA, Ross A, Tanner M. Predicting the cost-effectiveness of introducing a pre-erythrocytic malaria vaccine into the expanded program on immunization in tanzania. Am J Trop Med Hyg. 2006;75(2 suppl):119–30.

  27. 27.

    Maire N, Shillcutt SD, Walker DG, Tediosi F, Smith TA. Cost-effectiveness of the introduction of a pre-erythrocytic malaria vaccine into the expanded program on immunization in sub-Saharan Africa: analysis of uncertainties using a stochastic individual-based simulation model of Plasmodium falciparum malaria. Value Health. 2011;14:1028–38.

  28. 28.

    Phillips V, Njau J, Li S, Kachur P. Simulations show diagnostic testing for malaria in young African children can be cost-saving or cost-effective. Health Aff. 2015;34:1196–203.

  29. 29.

    McKenzie FE, Wong RC, Bossert WH. Discrete-event simulation models of Plasmodium falciparum malaria. Simulation. 1998;71:250–61.

  30. 30.

    Pizzitutti F, Pan W, Barbieri A, Miranda JJ, Feingold B, Guedes GR, et al. A validated agent-based model to study the spatial and temporal heterogeneities of malaria incidence in the rainforest environment. Malar J. 2015;14:514.

  31. 31.

    Zhu L, Qualls WA, Marshall JM, Arheart KL, DeAngelis DL, McManus JW, et al. A spatial individual-based model predicting a great impact of copious sugar sources and resting sites on survival of Anopheles gambiae and malaria parasite transmission. Malar J. 2015;14:59.

  32. 32.

    Smith T, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, et al. Relationship between the entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006;75(2 suppl):11–8.

  33. 33.

    Maire N, Aponte JJ, Ross A, Thompson R, Alonso P, Utzinger J, et al. Modeling a field trial of the RTS, S/AS02A malaria vaccine. Am J Trop Med Hyg. 2006;75(2 suppl):104–10.

  34. 34.

    Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L. A model for natural immunity to asexual blood stages of Plasmodium falciparum malaria in endemic areas. Am J Trop Med Hyg. 2006;75(2 Suppl):19–31.

  35. 35.

    Smith T, Ross A, Maire N, Rogier C, Trape J-F, Molineaux L. An epidemiologic model of the incidence of acute illness in Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006;75(2 suppl):56–62.

  36. 36.

    Zhu L, Müller GC, Marshall JM, Arheart KL, Qualls WA, Hlaing WM, et al. Is outdoor vector control needed for malaria elimination? An individual-based modelling study. Malar J. 2017;16:266.

  37. 37.

    Pizzitutti F, Pan W, Feingold B, Zaitchik B, Álvarez CA, Mena CF. Out of the net: an agent-based model to study human movements influence on local-scale malaria transmission. PLoS ONE. 2018;13:e0193493.

  38. 38.

    Bomblies A, Duchemin J-B, Eltahir EAB. A mechanistic approach for accurate simulation of village scale malaria transmission. Malar J. 2009;8:223–34.

  39. 39.

    Yamana TK, Bomblies A, Laminou IM, Duchemin J-B, Eltahir EAB. Linking environmental variability to village-scale malaria transmission using a simple immunity model. Parasit Vectors. 2013;6:226.

  40. 40.

    Yamana TK, Qiu X, Eltahir EAB. Hysteresis in simulations of malaria transmission. Adv Water Resour. 2017;108:416–22.

  41. 41.

    Endo N, Eltahir EAB. Environmental determinants of malaria transmission around the Koka Reservoir in Ethiopia. GeoHealth. 2018;2:104–15.

  42. 42.

    Endo N, Eltahir EAB. Environmental determinants of malaria transmission in African villages. Malar J. 2016;15:578.

  43. 43.

    Penny MA, Maire N, Studer A, Schapira A, Smith TA. What should vaccine developers ask? Simulation of the effectiveness of malaria vaccines. PLoS ONE. 2008;3:e3193.

  44. 44.

    Maire N, Tediosi F, Ross A, Smith T. Predictions of the epidemiologic impact of introducing a pre-erythrocytic vaccine into the expanded program on immunization in sub-Saharan Africa. Am J Trop Med Hyg. 2006;75(2 suppl):111–8.

  45. 45.

    Nguyen TD, Olliaro P, Dondorp AM, Baird JK, Lam HM, Farrar J, et al. Optimum population-level use of artemisinin combination therapies: a modelling study. Lancet Glob Health. 2015;3:e758.

  46. 46.

    Ouédraogo AL, Eckhoff PA, Luty AJF, Roeffen W, Sauerwein RW, Bousema T, et al. Modeling the impact of Plasmodium falciparum sexual stage immunity on the composition and dynamics of the human infectious reservoir for malaria in natural settings. PLoS Pathog. 2018;14:e1007034.

  47. 47.

    Ross A, Smith T. Interpreting malaria age-prevalence and incidence curves: a simulation study of the effects of different types of heterogeneity. Malar J. 2010;9:132–40.

  48. 48.

    Choi SE, Brandeau ML, Bendavid E. Cost-effectiveness of malaria preventive treatment for HIV-infected pregnant women in sub-Saharan Africa. Malar J. 2017;16:1–10.

  49. 49.

    Killeen GF, Ross A, Smith T. Infectiousness of malaria-endemic human populations to vectors. Am J Trop Med Hyg. 2006;75(2 Suppl):38–45.

  50. 50.

    Griffin JT, Hollingsworth TD, Reyburn H, Drakeley CJ, Riley EM, Ghani AC. Gradual acquisition of immunity to severe malaria with increasing exposure. Proc R Soc B Biol Sci. 2015;282:20142657.

  51. 51.

    Arifin SMN, Madey GR, Collins FH. Examining the impact of larval source management and insecticide-treated nets using a spatial agent-based model of Anopheles gambiae and a landscape generator tool. Malar J. 2013;12:290–313.

  52. 52.

    Chitnis N, Hardy D, Smith T. A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012;74:1098–124.

  53. 53.

    Gu W, Novak RJ. Agent-based modelling of mosquito foraging behaviour for malaria control. Trans R Soc Trop Med Hyg. 2009;103:1105–18.

  54. 54.

    Cairns ME, Walker PGT, Okell LC, Griffin JT, Garske T, Asante KP, et al. Seasonality in malaria transmission: implications for case-management with long-acting artemisinin combination therapy in sub-Saharan Africa. Malar J. 2015;14:321.

  55. 55.

    Bomblies A. Agent-based modeling of malaria vectors: the importance of spatial simulation. Parasit Vectors. 2014;7:308.

  56. 56.

    Shcherbacheva A, Haario H, Killeen GF. Modeling host-seeking behavior of African malaria vector mosquitoes in the presence of long-lasting insecticidal nets. Math Biosci. 2018;295:36–47.

  57. 57.

    Shcherbacheva A, Haario H. The impact of household size on malaria reduction in relation with alterations in mosquito behavior by malaria parasite. J Multi-Valued Log Soft Comput. 2017;29:455–68.

  58. 58.

    Depinay JMO, Mbogo CM, Killeen G, Knols B, Beier J, Carlson J, et al. A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004;3:29.

  59. 59.

    Gentile JE, Rund SSC, Madey GR. Modelling sterile insect technique to control the population of Anopheles gambiae. Malar J. 2015;14:92.

  60. 60.

    Mckenzie FE, Killeen GF, Beier JC, Bossert WH. Seasonality, parasite diversity, and local extinctions in plasmodium falciparum malaria. Ecology. 2001;82:2673–81.

  61. 61.

    Arifin S, Arifin R, Pitts D, Rahman M, Nowreen S, Madey G, et al. Landscape epidemiology modeling using an agent-based model and a geographic information system. Land. 2015;4:378–412.

  62. 62.

    Eckhoff PA, Wenger EA, Godfray HCJ, Burt A. Impact of mosquito gene drive on malaria elimination in a computational model with explicit spatial and temporal dynamics. Proc Natl Acad Sci USA. 2016;114:e255–64.

  63. 63.

    Alam MSMZ, Niaz Arifin SM, Al-Amin HM, Alam MSMZ, Rahman MS. A spatial agent-based model of Anopheles vagus for malaria epidemiology: examining the impact of vector control interventions. Malar J. 2017;16:432.

  64. 64.

    Churcher TS, Dawes EJ, Sinden RE, Christophides GK, Koella JC, Basáñez M-G. Population biology of malaria within the mosquito: density-dependent processes and potential implications for transmission-blocking interventions. Malar J. 2010;9:311.

  65. 65.

    Arifin SMN, Zhou Y, Davis GJ, Gentile JE, Madey GR, Collins FH. An agent-based model of the population dynamics of Anopheles gambiae. Malar J. 2014;13:424.

  66. 66.

    Eckhoff P. Mathematical models of within-host and transmission dynamics to determine effects of malaria interventions in a variety of transmission settings. Am J Trop Med Hyg. 2013;88:817–27.

  67. 67.

    Gerardin J, Bever CA, Bridenbecker D, Hamainza B, Silumbe K, Miller JM, et al. Effectiveness of reactive case detection for malaria elimination in three archetypical transmission settings: a modelling study. Malar J. 2017;16:248.

  68. 68.

    Gerardin J, Ouédraogo AL, McCarthy KA, Eckhoff PA, Wenger EA. Characterization of the infectious reservoir of malaria with an agent-based model calibrated to age-stratified parasite densities and infectiousness. Malar J. 2015;14:231.

  69. 69.

    Walker PGT, Griffin JT, Ferguson NM, Ghani AC. Estimating the most efficient allocation of interventions to achieve reductions in Plasmodium falciparum malaria burden and transmission in Africa: a modelling study. Lancet Glob Health. 2016;4:e474.

  70. 70.

    Ferreira CP, Lyra SP, Azevedo F, Greenhalgh D, Massad E. Modelling the impact of the long-term use of insecticide-treated bed nets on Anopheles mosquito biting time. Malar J. 2017;16:373.

  71. 71.

    Linard C, Ponçon N, Fontenille D, Lambin EF. A multi-agent simulation to assess the risk of malaria re-emergence in southern France. Ecol Modell. 2009;220:160–74.

  72. 72.

    Gu W, Novak RJ. Predicting the impact of insecticide-treated bed nets on malaria transmission: the devil is in the detail. Malar J. 2009;8:256.

  73. 73.

    Molineaux L, Diebner HH, Eichner M, Collins WE, Jeffery GM, Dietz K. Plasmodium falciparum parasitaemia described by a new mathematical model. Parasitology. 2001;122(Pt 4):379–91.

  74. 74.

    McKenzie FE, Bossert WH. An integrated model of Plasmodium falciparum dynamics. J Theor Biol. 2005;232:411–26.

  75. 75.

    Gatton ML, Dunn J, Chaudhry A, Ciketic S, Cunningham J, Cheng Q. Implications of parasites lacking Plasmodium falciparum histidine-rich protein 2 on malaria morbidity and control when rapid diagnostic tests are used for diagnosis. J Infect Dis. 2017;215:1156–66.

  76. 76.

    Watson OJ, Slater HC, Verity R, Parr JB, Mwandagalirwa MK, Tshefu A, et al. Modelling the drivers of the spread of Plasmodium falciparum hrp2 gene deletions in sub-Saharan Africa. Elife. 2017;6:e25008.

  77. 77.

    Klein EY. The impact of heterogeneous transmission on the establishment and spread of antimalarial drug resistance. J Theor Biol. 2014;340:177–85.

  78. 78.

    McCarthy KA, Wenger EA, Huynh GH, Eckhoff PA. Calibration of an intrahost malaria model and parameter ensemble evaluation of a pre-erythrocytic vaccine. Malar J. 2015;14:6.

  79. 79.

    Gurarie D, Karl S, Zimmerman PA, King CH, St Pierre TG, Davis TME. Mathematical modeling of malaria infection with innate and adaptive immunity in individuals and agent-based communities. PLoS ONE. 2012;7:e34040.

  80. 80.

    Gerardin J, Eckhoff P, Wenger EA. Mass campaigns with antimalarial drugs: a modelling comparison of artemether-lumefantrine and DHA-piperaquine with and without primaquine as tools for malaria control and elimination. BMC Infect Dis. 2015;15:144.

  81. 81.

    Eckhoff PA. Malaria parasite diversity and transmission intensity affect development of parasitological immunity in a mathematical model. Malar J. 2012;11:419.

  82. 82.

    Wenger EA, Eckhoff PA. A mathematical model of the impact of present and future malaria vaccines. Malar J. 2013;12:126.

  83. 83.

    Sauboin CJ, Van Bellinghen L-A, Van De Velde N, Van Vlaenderen I. Potential public health impact of RTS, S malaria candidate vaccine in sub-Saharan Africa: a modelling study. Malar J. 2015;14:524.

  84. 84.

    Stryker JJ, Bomblies A. The impacts of land use change on malaria vector abundance in a water-limited, highland region of Ethiopia. EcoHealth. 2012;9:455–70.

  85. 85.

    Gerardin J, Bever CA, Hamainza B, Miller JM, Eckhoff PA, Wenger EA. Optimal population-level infection detection strategies for malaria control and elimination in a spatial model of malaria transmission. PLoS Comput Biol. 2016;12:e1004707.

  86. 86.

    Rateb F, Pavard B, Bellamine-BenSaoud N, Merelo JJ, Arenas MG. Modeling malaria with multi-agent systems. Int J Intell Inf Technol. 2005;1:17–27.

  87. 87.

    Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, et al. A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009;6:e1000048.

  88. 88.

    Markham CG. Seasonaility of precipitation in the United States. Ann Assoc Am Geogr. 1970;60:593–7.

  89. 89.

    Griffin JT, Ferguson NM, Ghani AC. Estimates of the changing age-burden of Plasmodium falciparum malaria disease in sub-Saharan Africa. Nat Commun. 2014;5:3136.

  90. 90.

    Slater HC, Griffin JT, Ghani AC, Okell LC. Assessing the potential impact of artemisinin and partner drug resistance in sub-Saharan Africa. Malar J. 2016;15:10.

  91. 91.

    Okell LC, Cairns M, Griffin JT, Ferguson NM, Tarning J, Jagoe G, et al. Contrasting benefits of different artemisinin combination therapies as first-line malaria treatments using model-based cost-effectiveness analysis. Nat Commun. 2014;5:5606.

  92. 92.

    Griffin JT, Bhatt S, Sinka ME, Gething PW, Lynch M, Patouillard E, et al. Potential for reduction of burden and local elimination of malaria by reducing Plasmodium falciparum malaria transmission: a mathematical modelling study. Lancet Infect Dis. 2016;16:465–72.

  93. 93.

    Bretscher MT, Griffin JT, Ghani AC, Okell LC. Modelling the benefits of long-acting or transmission-blocking drugs for reducing Plasmodium falciparum transmission by case management or by mass treatment. Malar J. 2017;16:341.

  94. 94.

    Gu W, Killeen GF, Mbogo CM, Regens JL, Githure JI, Beier JC. An individual-based model of Plasmodium falciparum malaria transmission on the coast of Kenya. Trans R Soc Trop Med Hyg. 2003;97:43–50.

  95. 95.

    Gerardin J, Bertozzi-Villa A, Eckhoff PA, Wenger EA. Impact of mass drug administration campaigns depends on interaction with seasonal human movement. Int Health. 2018;10:252–7.

  96. 96.

    Slater HC, Walker PGT, Bousema T, Okell LC, Ghani AC. The potential impact of adding ivermectin to a mass treatment intervention to reduce malaria transmission: a modelling study. J Infect Dis. 2014;210:1972–80.

  97. 97.

    Winskill P, Walker PG, Griffin JT, Ghani AC. Modelling the cost-effectiveness of introducing the RTS, S malaria vaccine relative to scaling up other malaria interventions in sub-Saharan Africa. BMJ Glob Health. 2017;2:e000090.

  98. 98.

    Okell LC, Griffin JT, Kleinschmidt I, Hollingsworth TD, Churcher TS, White MJ, et al. The potential contribution of mass treatment to the control of Plasmodium falciparum malaria. PLoS ONE. 2011;6:e20179.

  99. 99.

    Arifin SMN, Davis GJ, Zhou Y, Madey GR. Verification and validation by docking: a case study of agent-based models of Anopheles gambiae. In: Proceedings of the 2010 Summer Computer Simulation Conference; 2010. p. 1–8.

  100. 100.

    Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, et al. Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008;135:1507–16.

  101. 101.

    Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, et al. An approach to model the costs and effects of case management of Plasmodium falciparum malaria in sub-Saharan Africa. Am J Trop Med Hyg. 2006;75(2 Suppl):90–103.

  102. 102.

    Ross A, Maire N, Sicuri E, Smith T, Conteh L. Determinants of the cost-effectiveness of intermittent preventive treatment for malaria in infants and children. PLoS ONE. 2011;6:e18391.

  103. 103.

    Crowell V, Briët OJ, Hardy D, Chitnis N, Maire N, Di Pasquale A, et al. Modelling the cost-effectiveness of mass screening and treatment for reducing Plasmodium falciparum malaria burden. Malar J. 2013;12:4.

  104. 104.

    Tediosi F, Maire N, Penny M, Studer A, Smith TA. Simulation of the cost-effectiveness of malaria vaccines. Malar J. 2009;8:127.

  105. 105.

    Ross A, Penny M, Maire N, Studer A, Carneiro I, Schellenberg D, et al. Modelling the epidemiological impact of intermittent preventive treatment against malaria in infants. PLoS ONE. 2008;3:e2661.

  106. 106.

    Briët OJT, Hardy D, Smith TA. Importance of factors determining the effective lifetime of a mass, long-lasting, insecticidal net distribution: a sensitivity analysis. Malar J. 2012;11:20.

  107. 107.

    Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, et al. Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012;9:e1001157.

  108. 108.

    Penny MA, Pemberton-Ross P, Smith TA. The time-course of protection of the RTS, S vaccine against malaria infections and clinical disease. Malar J. 2015;14:437.

  109. 109.

    Stuckey EM, Miller JM, Littrell M, Chitnis N, Steketee R. Operational strategies of anti-malarial drug campaigns for malaria elimination in Zambia’s southern province: a simulation study. Malar J. 2016;15:148.

  110. 110.

    Stuckey EM, Stevenson JC, Cooke MK, Owaga C, Marube E, Oando G, et al. Simulation of malaria epidemiology and control in the highlands of western Kenya. Malar J. 2012;11:357.

  111. 111.

    Briët OJT, Penny MA, Hardy D, Awolola TS, Van Bortel W, Corbel V, et al. Effects of pyrethroid resistance on the cost effectiveness of a mass distribution of long- lasting insecticidal nets: a modelling study. Malar J. 2013;12:77.

  112. 112.

    Briët OJT, Chitnis N. Effects of changing mosquito host searching behaviour on the cost effectiveness of a mass distribution of long-lasting, insecticidal nets: a modelling study. Malar J. 2013;12:215.

  113. 113.

    Briët OJT, Penny MA. Repeated mass distributions and continuous distribution of long-lasting insecticidal nets: modelling sustainability of health benefits from mosquito nets, depending on case management. Malar J. 2013;12:401.

  114. 114.

    Stuckey EM, Stevenson J, Galactionova K, Baidjoe AY, Bousema T, Odongo W, et al. Modeling the cost effectiveness of malaria control interventions in the highlands of western Kenya. PLoS ONE. 2014;9:e107700.

  115. 115.

    Pemberton-Ross P, Smith TA, Hodel EM, Kay K, Penny MA. Age-shifting in malaria incidence as a result of induced immunological deficit: a simulation study. Malar J. 2015;14:287.

  116. 116.

    Penny MA, Galactionova K, Tarantino M, Tanner M, Smith TA. The public health impact of malaria vaccine RTS, S in malaria endemic Africa: country-specific predictions using 18 month follow-up Phase III data and simulation models. BMC Med. 2015;13:170.

  117. 117.

    Yukich J, Chitnis N. When can malaria control and elimination programs safely reduce vector control efforts?. A simulation study. Geneva: World Health Organization; 2015.

  118. 118.

    Cameron E, Battle KE, Bhatt S, Weiss DJ, Bisanzio D, Mappin B, et al. Defining the relationship between infection prevalence and clinical incidence of Plasmodium falciparum malaria. Nat Commun. 2015;6:8170.

  119. 119.

    Penny MA, Verity R, Bever CA, Sauboin C, Galactionova K, Flasche S, et al. Public health impact and cost-effectiveness of the RTS, S/AS01 malaria vaccine: a systematic comparison of predictions from four mathematical models. Lancet. 2016;387:367–75.

  120. 120.

    Brady OJ, Slater HC, Pemberton-Ross P, Wenger E, Maude RJ, Ghani AC, et al. Role of mass drug administration in elimination of Plasmodium falciparum malaria: a consensus modelling study. Lancet Glob Health. 2017;5:e680–7.

  121. 121.

    Deville P, Linard C, Martin S, Gilbert M, Stevens FR, Gaughan AE, et al. Dynamic population mapping using mobile phone data. Proc Natl Acad Sci USA. 2014;111:15888–93.

  122. 122.

    Keith JM, Spring D. Agent-based Bayesian approach to monitoring the progress of invasive species eradication programs. Proc Natl Acad Sci US. 2013;110:13428–33.

  123. 123.

    Barber BE, William T, Dhararaj P, Anderios F, Grigg MJ, Yeo TW, et al. Epidemiology of Plasmodium knowlesi malaria in north-east Sabah, Malaysia: family clusters and wide age distribution. Malar J. 2012;11:401.

  124. 124.

    Lekone PE, Finkenstädt BF. Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study. Biometrics. 2006;62:1170–7.

  125. 125.

    Kerr CC, Stuart RM, Gray RT, Shattock AJ, Fraser-hurt N, Benedikt C, et al. Optima: a model for HIV epidemic analysis, program prioritization, and resource optimization. J Acquir Immune Defic Syndr. 2015;69:365–76.

  126. 126.

    Heesterbeek H, Anderson RM, Andreasen V, Bansal S, De Angelis D, Dye C, et al. Modeling infectious disease dynamics in the complex landscape of global health. Science. 2015;347:4339.

  127. 127.

    Kelly SL, Shattock AJ, Kerr CC, Stuart RM, Papoyan A, Grigoryan T, et al. Optimizing HIV/AIDS resources in Armenia: increasing ART investment and examining HIV programmes for seasonal migrant labourers. J Int AIDS Soc. 2016;19:1–7.

  128. 128.

    Scott N, Hussain SA, Martin-Hughes R, Fowkes FJI, Kerr CC, Pearson R, et al. Maximizing the impact of malaria funding through allocative efficiency: using the right interventions in the right locations. Malar J. 2017;16:368.

  129. 129.

    Barbati M, Bruno G, Genovese A. Applications of agent-based models for optimization problems: a literature review. Expert Syst Appl. 2012;39:6020–8.

  130. 130.

    Grimm V, Berger U, DeAngelis DL, Polhill JG, Giske J, Railsback SF. The ODD protocol: a review and first update. Ecol Modell. 2010;221:2760–8.

  131. 131.

    Chretien JP, Riley S, George DB. Mathematical modeling of the West Africa ebola epidemic. Elife. 2015;4:1–15.

Download references

Authors’ contributions

JK and JF conceived the project and jointly supervised the work undertaken for this manuscript. JSR and RM contributed to the study design and edited the manuscript. NRS designed and performed the literature review, analysed the results and drafted the manuscript. JF drafted the introduction of the manuscript. JF, JMT, MG and JK cross-checked papers, edited the manuscript and provided the structure of the paper. All authors read and approved the final manuscript.


Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Not applicable.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


N.R.S. is a recipient of an Australian Government Research Training Program (RTP) Scholarship. Funding was provided by the National Health and Medical Research Council (NHMRC) of Australia including an Early Career Research Fellowship (APP1037722) and Project Grant (APP1125788) to J.S.R. Burnet Institute received funding from the NHMRC Independent Research Institutes Infrastructure Support Scheme, and the Victorian State Government Operational Infrastructure Support Scheme. J.M.T. is a recipient of an Early Career Fellowship (APP1142638) from the NHMRC. Mahidol-Oxford Tropical Medicine Research Unit is funded by the Wellcome Trust of Great Britain.

Funding bodies had no role in the design of the study, collection, analysis, and interpretation of data or in writing the manuscript.

Publisher’s Note


Author information

Correspondence to Neal R. Smith.

Additional files


Additional file 1. Systematic review search strategies.


Additional file 2. Key characteristics of individual-based models of malaria.


Additional file 3. Adapted PRISMA search flow diagram of study selection.


Additional file 4. Overview of elements of spatial models.


Additional file 5. Overview of ensemble and comparative models.

Rights and permissions


Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark


  • Malaria
  • Infectious disease transmission
  • Agent-based model
  • Individual-based model
  • Review

